title: Graph Product Multilayer Networks: Spectral Properties and Applications

abstract: This paper aims to establish theoretical foundations of graph product multilayer networks (GPMNs), a family of multilayer networks that can be obtained as a graph product of two or more factor networks. Cartesian, direct (tensor), and strong product operators are considered, and then generalized. We first describe mathematical relationships between GPMNs and their factor networks regarding their degree/strength, adjacency, and Laplacian spectra, and then show that those relationships can still hold for nonsimple and generalized GPMNs. Applications of GPMNs are discussed in three areas: predicting epidemic thresholds, modeling propagation in nontrivial space and time, and analyzing higher-order properties of self-similar networks. Directions of future research are also discussed.


title: Strict s-numbers of the Volterra operator

abstract: For Volterra operator $V\colon L^1(0,1) \to C[0,1]$ and summation operator $\sigma\colon \ell^1 \to c$, we obtain exact values of Approximation, Gelfand, Kolmogorov, Mityagin and Isomorphism numbers.


title: Review of QCD, QGP, Heavy Quark Meson Production Enhancement and Suppression

abstract: This review of the Quantum Chromodynamics (QCD),the early universe Cosmoloical Phase Transition from the Quark-Gluon Plasma (QGP) to our present universe (QCDPT), Relativistic Heavy Ion Collisions (RHIC) which can produce the QGP, the possible detection of the QGP produced by the production of mixed hybrid heavy quark mesons. We also review recent studies of the production of mixed heavy quark hybrids via RHIC and heavy quark meson supression in p-Pb and Pb-Pb collisions.


title: HD 209458b in New Light: Evidence of Nitrogen Chemistry, Patchy Clouds and Sub-Solar Water

abstract: Interpretations of exoplanetary transmission spectra have been undermined by apparent obscuration due to clouds/hazes. Debate rages on whether weak H$_2$O features seen in exoplanet spectra are due to clouds or inherently depleted oxygen. Assertions of solar H$_2$O abundances have relied on making a priori model assumptions, e.g. chemical/radiative equilibrium. In this work, we attempt to address this problem with a new retrieval paradigm for transmission spectra. We introduce POSEIDON, a two-dimensional atmospheric retrieval algorithm including generalised inhomogeneous clouds. We demonstrate that this prescription allows one to break vital degeneracies between clouds and prominent molecular abundances. We apply POSEIDON to the best transmission spectrum presently available, for the hot Jupiter HD 209458b, uncovering new insights into its atmosphere at the day-night terminator. We extensively explore the parameter space with an unprecedented 10$^8$ models, spanning the continuum from fully cloudy to cloud-free atmospheres, in a fully Bayesian retrieval framework. We report the first detection of nitrogen chemistry (NH$_3$ and/or HCN) in an exoplanet atmosphere at 3.7-7.7$\sigma$ confidence, non-uniform cloud coverage at 4.5-5.4$\sigma$, high-altitude hazes at $>$3$\sigma$, and sub-solar H$_2$O at $\gtrsim$3-5$\sigma$, depending on the assumed cloud distribution. We detect NH$_3$ at 3.3$\sigma$ and 4.9$\sigma$ for fully cloudy and cloud-free scenarios, respectively. For the model with the highest Bayesian evidence, we constrain H$_2$O at 5-15 ppm (0.01-0.03$\times$ solar) and NH$_3$ at 0.01-2.7 ppm, strongly suggesting disequilibrium chemistry and cautioning against equilibrium assumptions. Our results herald new promise for retrieving cloudy atmospheres using high-precision HST and JWST spectra.


title: Strong Clustering of Lyman Break Galaxies around Luminous Quasars at z~4

abstract: In the standard picture of structure formation, the first massive galaxies are expected to form at the highest peaks of the density field, which constitute the cores of massive proto-clusters. Luminous quasars (QSOs) at z~4 are the most strongly clustered population known, and should thus reside in massive dark matter halos surrounded by large overdensities of galaxies, implying a strong QSO-galaxy cross-correlation function. We observed six z~4 QSO fields with VLT/FORS exploiting a novel set of narrow band filters custom designed to select Lyman Break Galaxies (LBGs) in a thin redshift slice of Delta_z~0.3, mitigating the projection effects that have limited the sensitivity of previous searches for galaxies around z>~4 QSOs. We find that LBGs are strongly clustered around QSOs, and present the first measurement of the QSO-LBG cross-correlation function at z~4, on scales of 0.1<~R<~9 Mpc/h (comoving). Assuming a power law form for the cross-correlation function xi=(r/r0_QG)^gamma, we measure r0_QG=8.83^{+1.39}_{-1.51} Mpc/h for a fixed slope of gamma=2.0. This result is in agreement with the expected cross-correlation length deduced from measurements of the QSO and LBG auto-correlation function, and assuming a linear bias model. We also measure a strong auto-correlation of LBGs in our QSO fields finding r0_GG=21.59^{+1.72}_{-1.69} Mpc/h for a fixed slope of gamma=1.5, which is ~4 times larger than the LBG auto-correlation length in random fields, providing further evidence that QSOs reside in overdensities of LBGs. Our results qualitatively support a picture where luminous QSOs inhabit exceptionally massive (M_halo>10^12 M_sun) dark matter halos at z~4.


title: Self-induced dust traps: overcoming planet formation barriers

abstract: Planet formation is thought to occur in discs around young stars by the aggregation of small dust grains into much larger objects. The growth from grains to pebbles and from planetesimals to planets is now fairly well understood. The intermediate stage has however been found to be hindered by the radial-drift and fragmentation barriers. We identify a powerful mechanism in which dust overcomes both barriers. Its key ingredients are i) backreaction from the dust onto the gas, ii) grain growth and fragmentation, and iii) large-scale gradients. The pile-up of growing and fragmenting grains modifies the gas structure on large scales and triggers the formation of pressure maxima, in which particles are trapped. We show that these self-induced dust traps are robust: they develop for a wide range of disc structures, fragmentation thresholds and initial dust-to-gas ratios. They are favored locations for pebbles to grow into planetesimals, thus opening new paths towards the formation of planets.


title: Testing strong-field gravity with tidal Love numbers

abstract: The tidal Love numbers (TLNs) encode the deformability of a self-gravitating object immersed in a tidal environment and depend significantly both on the object's internal structure and on the dynamics of the gravitational field. An intriguing result in classical general relativity is the vanishing of the TLNs of black holes. We extend this result in three ways, aiming at testing the nature of compact objects: (i) we compute the TLNs of exotic compact objects, including different families of boson stars, gravastars, wormholes, and other toy models for quantum corrections at the horizon scale. In the black-hole limit, we find a universal logarithmic dependence of the TLNs on the location of the surface; (ii) we compute the TLNs of black holes beyond vacuum general relativity, including Einstein-Maxwell, Brans-Dicke and Chern-Simons gravity; (iii) We assess the ability of present and future gravitational-wave detectors to measure the TLNs of these objects, including the first analysis of TLNs with LISA. Both LIGO, ET and LISA can impose interesting constraints on boson stars, while LISA is able to probe even extremely compact objects. We argue that the TLNs provide a smoking gun of new physics at the horizon scale, and that future gravitational-wave measurements of the TLNs in a binary inspiral provide a novel way to test black holes and general relativity in the strong-field regime.


title: A Massive Shell of Supernova-formed Dust in SNR G54.1+0.3

abstract: While theoretical dust condensation models predict that most refractory elements produced in core-collapse supernovae (SNe) efficiently condense into dust, a large quantity of dust has so far only been observed in SN 1987A. We present the analysis of Spitzer Space Telescope, Herschel Space Observatory, Stratospheric Observatory for Infrared Astronomy (SOFIA), and AKARI observations of the infrared (IR) shell surrounding the pulsar wind nebula in the supernova remnant G54.1+0.3. We attribute a distinctive spectral feature at 21 $\mu$m to a magnesium silicate grain species that has been invoked in modeling the ejecta-condensed dust in Cas A, which exhibits the same spectral signature. If this species is responsible for producing the observed spectral feature and accounts for a significant fraction of the observed IR continuum, we find that it would be the dominant constituent of the dust in G54.1+0.3, with possible secondary contributions from other compositions, such as carbon, silicate, or alumina grains. The smallest mass of SN-formed dust required by our models is 1.1 $\pm$ 0.8 $\rm M_{\odot}$. We discuss how these results may be affected by varying dust grain properties and self-consistent grain heating models. The spatial distribution of the dust mass and temperature in G54.1+0.3 confirms the scenario in which the SN-formed dust has not yet been processed by the SN reverse shock and is being heated by stars belonging to a cluster in which the SN progenitor exploded. The dust mass and composition suggest a progenitor mass of 16$-$27 $\rm M_{\odot}$ and imply a high dust condensation efficiency, similar to that found for Cas A and SN 1987A. The study provides another example of significant dust formation in a Type IIP SN and sheds light on the properties of pristine SN-condensed dust.


title: First-principles study of the dynamic Jahn-Teller distortion of the neutral vacancy in diamond

abstract: First-principles density functional theory methods are used to investigate the structure, energetics, and vibrational motions of the neutral vacancy defect in diamond. The measured optical absorption spectrum demonstrates that the tetrahedral $T_d$ point group symmetry of pristine diamond is maintained when a vacancy defect is present. This is shown to arise from the presence of a dynamic Jahn-Teller distortion that is stabilised by large vibrational anharmonicity. Our calculations further demonstrate that the dynamic Jahn-Teller-distorted structure of $T_d$ symmetry is lower in energy than the static Jahn-Teller distorted tetragonal $D_{2d}$ vacancy defect, in agreement with experimental observations. The tetrahedral vacancy structure becomes more stable with respect to the tetragonal structure by increasing temperature. The large anharmonicity arises mainly from quartic vibrations, and is associated with a saddle point of the Born-Oppenheimer surface and a minimum in the free energy. This study demonstrates that the behaviour of Jahn-Teller distortions of point defects can be calculated accurately using anharmonic vibrational methods. Our work will open the way for first-principles treatments of dynamic Jahn-Teller systems in condensed matter.


title: A Hamiltonian approach for the Thermodynamics of AdS black holes

abstract: In this work we study the Thermodynamics of D-dimensional Schwarzschild-anti de Sitter (SAdS) black holes. The minimal Thermodynamics of the SAdS spacetime is briefly discussed, highlighting some of its strong points and shortcomings. The minimal SAdS Thermodynamics is extended within a Hamiltonian approach, by means of the introduction of an additional degree of freedom. We demonstrate that the cosmological constant can be introduced in the thermodynamic description of the SAdS black hole with a canonical transformation of the Schwarzschild problem, closely related to the introduction of an anti-de Sitter thermodynamic volume. The treatment presented is consistent, in the sense that it is compatible with the introduction of new thermodynamic potentials, and respects the laws of black hole Thermodynamics. By demanding homogeneity of the thermodynamic variables, we are able to construct a new equation of state that completely characterizes the Thermodynamics of SAdS black holes. The treatment naturally generates phenomenological constants that can be associated with different boundary conditions in underlying microscopic theories. A whole new set of phenomena can be expected from the proposed generalization of SAdS Thermodynamics.

categories: hep-th cond-mat.stat-mech gr-qc math-ph math.MP


title: Sensitivity estimate of the MACE gamma ray telescope

abstract: The MACE (Major Atmospheric Cherenkov Experiment) is a 21 m diameter gamma-ray telescope which is presently being installed at Hanle in Ladakh, India (32^0 46^' 46^" N, 78^0 58^' 35^" E) at an altitude of 4270 m a.s.l. Once operational, it will become the highest altitude very high energy (VHE) gamma-ray telescope in the world based on Imaging Atmospheric Cherenkov Technique (IACT). In the present work, we discuss the sensitivity estimate of the MACE telescope by using a substantially large Monte Carlo simulation database at 5^0 zenith angle. The sensitivity of MACE telescope is estimated by carrying out the gamma-hadron segregation using the Random Forest method. It is estimated that the MACE telescope will have an analysis energy threshold of 38 GeV for image intensities above 50 photoelectrons. The integral sensitivity for point like sources with Crab Nebula-like spectrum above 38 GeV is ~2.7% of Crab Nebula flux at 5 sigma statistical significance level in 50 hrs of observation.


title: Herschel and Hubble study of a lensed massive dusty starbursting galaxy at $z\sim3$

abstract: We present the results of combined deep Keck/NIRC2, HST/WFC3 near-infrared and Herschel far infrared observations of an extremely star forming dusty lensed galaxy identified from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS J133542.9+300401). The galaxy is gravitationally lensed by a massive WISE identified galaxy cluster at $z\sim1$. The lensed galaxy is spectroscopically confirmed at $z=2.685$ from detection of $\rm {CO (1 \rightarrow 0)}$ by GBT and from detection of $\rm {CO (3 \rightarrow 2)}$ obtained with CARMA. We use the combined spectroscopic and imaging observations to construct a detailed lens model of the background dusty star-forming galaxy (DSFG) which allows us to study the source plane properties of the target. The best-fit lens model provide magnification of $\mu_{\rm star}=2.10\pm0.11$ and $\mu_{\rm dust}=2.02\pm0.06$ for the stellar and dust components respectively. Multi-band data yields a magnification corrected star formation rate of $1900(\pm200)\,M_{\odot}{\rm yr^{-1}}$ and stellar mass of $6.8_{-2.7}^{+0.9}\times10^{11}\,M_{\odot}$ consistent with a main sequence of star formation at $z\sim2.6$. The CO observations yield a molecular gas mass of $8.3(\pm1.0)\times10^{10}\,M_{\odot}$, similar to the most massive star-forming galaxies, which together with the high star-formation efficiency are responsible for the intense observed star formation rates. The lensed DSFG has a very short gas depletion time scale of $\sim40$ Myr. The high stellar mass and small gas fractions observed indicate that the lensed DSFG likely has already formed most of its stellar mass and could be a progenitor of the most massive elliptical galaxies found in the local Universe.


title: The link between galaxy and black hole growth in the EAGLE simulation

abstract: We investigate the connection between the star formation rate (SFR) of galaxies and their central black hole accretion rate (BHAR) using the EAGLE cosmological hydrodynamical simulation. We find, in striking concurrence with recent observational studies, that the <SFR>--BHAR relation for an AGN selected sample produces a relatively flat trend, whilst the <BHAR>--SFR relation for a SFR selected sample yields an approximately linear trend. These trends remain consistent with their instantaneous equivalents even when both SFR and BHAR are time-averaged over a period of 100~Myr. There is no universal relationship between the two growth rates. Instead, SFR and BHAR evolve through distinct paths that depend strongly on the mass of the host dark matter halo. The galaxies hosted by haloes of mass M200 $\lesssim 10^{11.5}$Msol grow steadily, yet black holes (BHs) in these systems hardly grow, yielding a lack of correlation between SFR and BHAR. As haloes grow through the mass range $10^{11.5} \lesssim$ M200 $\lesssim 10^{12.5 }$Msol BHs undergo a rapid phase of non-linear growth. These systems yield a highly non-linear correlation between the SFR and BHAR, which are non-causally connected via the mass of the host halo. In massive haloes (M200 $\gtrsim 10^{12.5}$Msol) both SFR and BHAR decline on average with a roughly constant scaling of SFR/BHAR $\sim 10^{3}$. Given the complexity of the full SFR--BHAR plane built from multiple behaviours, and from the large dynamic range of BHARs, we find the primary driver of the different observed trends in the <SFR>--BHAR and <BHAR>--SFR relationships are due to sampling considerably different regions of this plane.


title: Argyres-Douglas matter and N=2 dualities

abstract: We study S duality of four dimensional N=2 Argyres-Douglas (AD) theory engineered from 6d A_{N-1} (2,0) theory. We find a (p,q) sequence of SCFTs, here (p,q) is co-prime and class S theory defined on sphere corresponds to class (0,1) theory. We represent these theories by a sphere with marked points, and S duality is interpreted as different pants decompositions of the same punctured sphere. The weakly coupled gauge theory description involves gauging AD matter which is represented by three punctured sphere.


title: Numerical simulations and infrared spectro-interferometry reveal the wind collision region in gamma2 Velorum

abstract: Colliding stellar winds in massive binary systems have been studied through their radio, optical lines and strong X-ray emission for decades. More recently, near-infrared spectrointerferometric observations have become available in a few systems, but isolating the contribution from the individual stars and the wind collision region still remains a challenge. In this paper, we study the colliding wind binary gamma2 Velorum and aim at identifying the wind collision zone from infrared interferometric data, which provide unique spatial information to determine the wind properties. Our analysis is based on multi-epoch VLTI/AMBER data that allows us to separate the spectral components of both stars. First, we determine the astrometric solution of the binary and confirm previous distance measurements. We then analyse the spectra of the individual stars, showing that the O star spectrum is peculiar within its class. Then, we perform three-dimensional hydrodynamic simulations of the system from which we extract model images, visibility curves and closure phases which can be directly compared with the observed data. The hydrodynamic simulations reveal the 3D spiral structure of the wind collision region, which results in phase-dependent emission maps. Our model visibility curves and closure phases provide a good match when the wind collision region accounts for 3 to 10 per cent of the total flux in the near infrared. The dialogue between hydrodynamic simulations, radiative transfer models and observations allows us to fully exploit the observations. Similar efforts will be crucial to study circumstellar environments with the new generation of VLTI instruments like GRAVITY and MATISSE.


title: Emergence and complexity in theoretical models of self-organized criticality

abstract: In this thesis we present few theoretical studies of the models of self-organized criticality. Following a brief introduction of self-organized criticality, we discuss three main problems. The first problem is about growing patterns formed in the abelian sandpile model (ASM). The patterns exhibit proportionate growth where different parts of the pattern grow in same rate, keeping the overall shape unchanged. This non-trivial property, often found in biological growth, has received increasing attention in recent years. In this thesis, we present a mathematical characterization of a large class of such patterns in terms of discrete holomorphic functions. In the second problem, we discuss a well known model of self-organized criticality introduced by Zhang in 1989. We present an exact analysis of the model and quantitatively explain an intriguing property known as the emergence of quasi-units. In the third problem, we introduce an operator algebra to determine the steady state of a class of stochastic sandpile models.


title: Textual Entailment with Structured Attentions and Composition

abstract: Deep learning techniques are increasingly popular in the textual entailment task, overcoming the fragility of traditional discrete models with hard alignments and logics. In particular, the recently proposed attention models (Rockt\"aschel et al., 2015; Wang and Jiang, 2015) achieves state-of-the-art accuracy by computing soft word alignments between the premise and hypothesis sentences. However, there remains a major limitation: this line of work completely ignores syntax and recursion, which is helpful in many traditional efforts. We show that it is beneficial to extend the attention model to tree nodes between premise and hypothesis. More importantly, this subtree-level attention reveals information about entailment relation. We study the recursive composition of this subtree-level entailment relation, which can be viewed as a soft version of the Natural Logic framework (MacCartney and Manning, 2009). Experiments show that our structured attention and entailment composition model can correctly identify and infer entailment relations from the bottom up, and bring significant improvements in accuracy.


title: Bosonic integer quantum Hall effect as topological pumping

abstract: Based on a quasi-one-dimensional limit of quantum Hall states on a thin torus, we construct a model of interaction-induced topological pumping which mimics the Hall response of the bosonic integer quantum Hall (BIQH) state. The quasi-one-dimensional counterpart of the BIQH state is identified as the Haldane phase composed of two-component bosons which form effective spin-$1$ degrees of freedom. An adiabatic change between the Haldane phase and trivial Mott insulators constitute {\it off-diagonal} topological pumping in which the translation of the lattice potential for one component induces a current in the other. The mechanism of this pumping is interpreted in terms of changes in polarizations between symmetry-protected quantized values.


title: Ferroelectric or non-ferroelectric: why so many materials exhibit ferroelectricity on the nanoscale

abstract: Ferroelectric materials have remained one of the foci of condensed matter physics and materials science for over 50 years. In the last 20 years, the development of voltage-modulated scanning probe microscopy techniques, exemplified by Piezoresponse force microscopy (PFM) and associated time and voltage spectroscopies, opened a pathway to explore these materials on a single-digit nanometer level. Consequently, domain structures, walls and polarization dynamics can now be imaged in real space. More generally, PFM has allowed studying electromechanical coupling in a broad variety of materials ranging from ionics to biological systems. It can also be anticipated that the recent Nobel prize in molecular electromechanical machines will result in rapid growth in interest in PFM as a method to probe their behavior on single device and device assembly levels. However, the broad introduction of PFM also resulted in a growing number of reports on nearly ubiquitous presence of ferroelectric-like phenomena including remnant polar states and electromechanical hysteresis loops in materials which are non-ferroelectric in the bulk, or in cases where size effects are expected to suppress ferroelectricity. While in certain cases plausible physical mechanisms can be suggested, there is remarkable similarity in observed behaviors, irrespective of the materials system. In this review, we summarize the basic principles of PFM, briefly discuss the features of ferroelectric surfaces salient to PFM imaging and spectroscopy, and summarize existing reports on ferroelectric like responses in non-classical ferroelectric materials. We further discuss possible mechanisms behind observed behaviors, and possible experimental strategies for their identification.


title: Diophantine approximation in prescribed degree

abstract: We investigate approximation to a given real number by algebraic numbers and algebraic integers of prescribed degree. We deal with both best and uniform approximation, and highlight the similarities and differences compared with the intensely studied problem of approximation by algebraic numbers (and integers) of bounded degree. We establish the answer to a question of Bugeaud concerning approximation to transcendental real numbers by quadratic irrational numbers, and thereby we refine a result of Davenport and Schmidt from 1967. We also obtain several new characterizations of Liouville numbers, and certain new insights on inhomogeneous Diophantine approximation. As an auxiliary side result, we provide an upper bound for the number of certain linear combinations of two given relatively prime integer polynomials with a linear factor. We conclude with several open problems.


title: Quenching vs. Quiescence: forming realistic massive ellipticals with a simple starvation model

abstract: The decrease in star formation (SF) and the morphological change necessary to produce the $z=0$ elliptical galaxy population are commonly ascribed to a sudden quenching event, which is able to rid the central galaxy of its cold gas reservoir in a short time. Following this event, the galaxy is able to prevent further SF and stay quiescent via a maintenance mode. We test whether such a quenching event is truly necessary using a simple model of quiescence. In this model, hot gas (all gas above a temperature threshold) in a $\sim10^{12} M_{\odot}$ halo mass galaxy at redshift $z\sim3$ is prevented from cooling. The cool gas continues to form stars at a decreasing rate and the galaxy stellar mass, morphology, velocity dispersion and position on the color magnitude diagram (CMD) proceed to evolve. By $z=0$, the halo mass has grown to $10^{13} M_{\odot}$ and the galaxy has attained characteristics typical of an observed $z=0$ elliptical galaxy. Our model is run in the framework of a cosmological, smooth particle hydrodynamic code which includes SF, early stellar feedback, supernova feedback, metal cooling and metal diffusion. Additionally, we post-process our simulations with a radiative transfer code to create a mock CMD. In contrast to previous assumptions that a pure "fade away" model evolves too slowly to account for the sparsity of galaxies in the "green valley", we demonstrate crossing times of $<$1 Gyr. We conclude that no sudden quenching event is necessary to produce such rapid colour transitions.


title: A note on spaces of continuous functions on compact scattered spaces

abstract: In 1959, Pelczynski and Semadeni proved a theorem in which they gave some equivalent conditions for a compact Hausdorff space to be scattered. The purpose of the current note is that to clarify the meaning of the subtle term "conditionally weakly sequentially compact" they used as the basis for the proof of their theorem. Unfortunately, the term now is taken over by a similar but subtle concept that may cause a serious problem.


title: Performance Evaluation of Spectrum Mobility in Multi-homed Mobile IPv6 Cognitive Radio Cellular Networks

abstract: Technological developments alongside VLSI achievements enable mobile devices to be equipped with multiple radio interfaces which is known as multihoming. On the other hand, the combination of various wireless access technologies, known as Next Generation Wireless Networks (NGWNs) has been introduced to provide continuous connection to mobile devices in any time and location. Cognitive radio networks as a part of NGWNs aroused to overcome spectrum inefficiency and spectrum scarcity issues. In order to provide seamless and ubiquitous connection across heterogeneous wireless access networks in the context of cognitive radio networks, utilizing Mobile IPv6 is beneficial. In this paper, a mobile device equipped with two radio interfaces is considered in order to evaluate performance of spectrum handover in terms of handover latency. The analytical results show that the proposed model can achieve better performance compared to other related mobility management protocols mainly in terms of handover latency.


title: Topological phases of parafermionic chains with symmetries

abstract: We study the topological classification of parafermionic chains in the presence of a modified time reversal symmetry that satisfies ${\cal T}^2=1 $. Such chains can be realized in one dimensional structures embedded in fractionalized two dimensional states of matter, e.g. at the edges of a fractional quantum spin Hall system, where counter propagating modes may be gapped either by back-scattering or by coupling to a superconductor. In the absence of any additional symmetries, a chain of $\mathbb{Z}_m$ parafermions can belong to one of several distinct phases. We find that when the modified time reversal symmetry is imposed, the classification becomes richer. If $m $ is odd, each of the phases splits into two subclasses. We identify the symmetry protected phase as a Haldane phase that carries a Kramers doublet at each end. When $m $ is even, each phase splits into four subclasses. The origin of this split is in the emergent Majorana fermions associated with even values of $m$. We demonstrate the appearance of such emergent Majorana zero modes in a system where the constituents particles are either fermions or bosons.


title: Pruned cellular free resolutions of monomial ideals

abstract: Using discrete Morse theory, we give an algorithm that prunes the excess of information in the Taylor resolution and constructs a new cellular free resolution for an arbitrary monomial ideal. The pruned resolution is not simplicial in general, but we can slightly modify our algorithm in order to obtain a simplicial resolution. We also show that the Lyubeznik resolution fits into our pruning strategy. We finally use our methods to give a different approach to the theory of splitting of monomial ideals.


title: On utmost multiplicity of hierarchical stellar systems

abstract: According to theoretical considerations, multiplicity of hierarchical stellar systems can reach, depending on masses and orbital parameters, several hundred, while observational data confirm existence of at most septuple (seven-component) systems. In this study, we cross-match very high multiplicity (six and more components) stellar systems in modern catalogues of visual double and multiple stars, to find candidates to hierarchical systems among them. After cross-matching with catalogues of closer binaries (eclipsing, spectroscopic, etc.), some of their components were found to be binary/multiple themselves, which increases the system's degree of multiplicity. Optical pairs, known from literature or filtered by the authors, are flagged and excluded from the statistics. We have compiled a list of potentially very high multiplicity hierarchical systems that contains 10~objects. Their multiplicity does not exceed 12, and we discuss a number of ways to explain the lack of extremely high multiplicity systems.


title: Fuzzy Dark Matter from Infrared Confining Dynamics

abstract: A very light boson of mass $\mathcal{O}(10^{-22})$ eV may potentially be a viable dark matter (DM) candidate which can avoid phenomenological problems associated with cold DM. Such "fuzzy DM (FDM)" may naturally be an axion with a decay constant $f_a \sim 10^{16} \div 10^{18}$ GeV, and a mass $m_a \sim \mu^2/f_a$ with $\mu\sim 10^2$ eV. Here we propose a concrete model where $\mu$ arises as a dynamical scale from infrared confining dynamics, analogous to QCD. Our model is an alternative to the usual approach of generating $\mu$ through string theoretic instanton effects. We outline the features of this scenario that result from various cosmological constraints. We find that those constraints are suggestive of a period of mild of inflation, perhaps from a strong first order phase transition, that reheats the Standard Model (SM) sector only. A typical prediction of our scenario, broadly speaking, is a larger effective number of neutrinos compared to the SM value $N_{\text{eff}} \approx 3$, as inferred from precision measurements of the cosmic microwave background. Some of the new degrees of freedom may be identified as "sterile neutrinos," which may be required to explain certain neutrino oscillation anomalies. Hence, aspects of our scenario could be testable in terrestrial experiments, which is a novelty of our FDM model.


title: An architecture for efficient gravitational wave parameter estimation with multimodal linear surrogate models

abstract: The recent direct observation of gravitational waves has further emphasized the desire for fast, low-cost, and accurate methods to infer the parameters of gravitational wave sources. Due to expense in waveform generation and data handling, the cost of evaluating the likelihood function limits the computational performance of these calculations. Building on recently developed surrogate models and a novel parameter estimation pipeline, we show how to quickly generate the likelihood function as an analytic, closed-form expression. Using a straightforward variant of a production-scale parameter estimation code, we demonstrate our method using surrogate models of effective-one-body and numerical relativity waveforms. Our study is the first time these models have been used for parameter estimation and one of the first ever parameter estimation calculations with multi-modal numerical relativity waveforms, which include all l <= 4 modes. Our grid-free method enables rapid parameter estimation for any waveform with a suitable reduced-order model. The methods described in this paper may also find use in other data analysis studies, such as vetting coincident events or the computation of the coalescing-compact-binary detection statistic.


title: How many elements of a Coxeter group have a unique reduced expression?

abstract: Let $(W,R)$ be an arbitrary Coxeter system. We determine the number of elements of $W$ that have a unique reduced expression.


title: Pure states of maximum uncertainty with respect to a given POVM

abstract: One of the differences between classical and quantum world is that in the former we can always perform a measurement that gives certain outcomes for all pure states, while such a situation is not possible in the latter. The degree of randomness of the distribution of the measurement outcomes can be quantified by the Shannon entropy. While it is well known that this entropy, as a function of quantum states, needs to be minimized by some pure states, we would like to address the question how 'badly' can we end by choosing initially any pure state, i.e., which pure states produce the maximal amount of uncertainty under given measurement. We find these maximizers for all highly symmetric POVMs in dimension 2, and for all SIC-POVMs in any dimension.


title: Learning causal effects from many randomized experiments using regularized instrumental variables

abstract: Scientific and business practices are increasingly resulting in large collections of randomized experiments. Analyzed together, these collections can tell us things that individual experiments in the collection cannot. We study how to learn causal relationships between variables from the kinds of collections faced by modern data scientists: the number of experiments is large, many experiments have very small effects, and the analyst lacks metadata (e.g., descriptions of the interventions). Here we use experimental groups as instrumental variables (IV) and show that a standard method (two-stage least squares) is biased even when the number of experiments is infinite. We show how a sparsity-inducing l0 regularization can --- in a reversal of the standard bias--variance tradeoff in regularization --- reduce bias (and thus error) of interventional predictions. Because we are interested in interventional loss minimization we also propose a modified cross-validation procedure (IVCV) to feasibly select the regularization parameter. We show, using a trick from Monte Carlo sampling, that IVCV can be done using summary statistics instead of raw data. This makes our full procedure simple to use in many real-world applications.


title: Parton saturation scaling function for exclusive production of vector mesons and Deeply Virtual Compton Scattering

abstract: We provide a universal expression of cross sections for the exclusive vector meson production and Deeply Virtual Compton Scattering (DVCS) in photon-proton and photon-nucleus interactions based on the geometric scaling phenomenon. The theoretical parameterization based on the scaling property depends only on the single variable $\tau_A = Q^2/Q_{\mathrm{sat}}^2$, where the saturation scale, $Q_{\mathrm{sat}}$, drives the energy dependence and the corresponding nuclear effects. This phenomenological result describes all available data from DESY-HERA for $\rho,\,\phi,\,J/\psi$ production and DVCS measurements. A discussion on the size of nuclear shadowing corrections on photon-nucleus interaction is also done.


title: The Dem@Care Experiments and Datasets: a Technical Report

abstract: The objective of Dem@Care is the development of a complete system providing personal health services to people with dementia, as well as medical professionals and caregivers, by using a multitude of sensors, for context-aware, multi-parametric monitoring of lifestyle, ambient environment, and health parameters. Multi-sensor data analysis, combined with intelligent decision making mechanisms, will allow an accurate representation of the person's current status and will provide the appropriate feedback, both to the person and the associated caregivers, enhancing the standard clinical workflow. Within the project framework, several data collection activities have taken place to assist technical development and evaluation tasks. In all these activities, particular attention has been paid to adhere to ethical guidelines and preserve the participants' privacy. This technical report describes shorty the (a) the main objectives of the project, (b) the main ethical principles and (c) the datasets that have been already created.


title: More lessons from the six box toy experiment

abstract: Following a paper in which the fundamental aspects of probabilistic inference were introduced by means of a toy experiment, details of the analysis of simulated long sequences of extractions are shown here. In fact, the striking performance of probability-based inference and forecasting, compared to those obtained by simple `rules', might impress those practitioners who are usually underwhelmed by the philosophical foundation of the different methods. The analysis of the sequences also shows how the smallness of the probability of what has been actually observed, given the hypotheses of interest, is irrelevant for the purpose of inference.


title: Tropical limit and the micro-macro correspondence

abstract: Tropical mathematics is used to establish a correspondence between certain microscopic and macroscopic objects in statistical models. Tropical algebra gives a common framework for macrosystems (subsets) and their elementary constituents (elements) that is well-behaved with respect to composition. This kind of connection is studied with maps that preserve a monoid structure. The approach highlights an underlying order relation that is explored through the concepts of filter and ideal. Particular attention is paid to asymmetry and duality between max- and min-criteria. Physical implementations are presented through simple examples in thermodynamics and non-equilibrium physics. The phenomenon of ultrametricity, the notion of tropical equilibrium and the role of ground energy in non-equilibrium models are discussed. Tropical symmetry, i.e. idempotence, is investigated.


title: The analytical solution of the problem on plasma oscillations in half-space with diffusion boundary conditions

abstract: The boundary problem about behaviour (oscillations) of the electronic plasmas with arbitrary degree of degeneration of electronic gas in half-space with diffusion boundary conditions is analytically solved. The kinetic equation of Vlasov - Boltzmann with integral of collisions of type BGK (Bhatnagar, Gross, Krook) and Maxwell equation for electric field are applied. Distribution function for electrons and electric field in plasma in the form of expansion under eigen solutions of the initial system of equations are received. Coefficients of these expansions are found by means of the boundary conditions.


title: A Comparison Between Physics-based and Polytropic MHD Models for Stellar Coronae and Stellar Winds of Solar Analogs

abstract: The development of the Zeeman-Doppler Imagine (ZDI) technique has provided synoptic observations of surface magnetic fields of low-mass stars. This led the Stellar Astrophysics community to adopt modeling techniques that have been used in solar physics using solar magnetograms. However, many of these techniques have been neglected by the solar community due to their failure to reproduce solar observations. Nevertheless, some of these techniques are still used to simulate the coronae and winds of solar analogs. Here we present a comparative study between two MHD models for the solar corona and solar wind. The first type of model is a polytropic wind model, and the second is the physics-based AWSOM model. We show that while the AWSOM model consistently reproduces many solar observations, the polytropic model fails to reproduce many of these observations and in the cases it does, its solutions are unphysical. Our recommendation is that polytropic models, which are used to estimate mass-loss rates and other parameters of solar analogs, must be first calibrated with solar observations. Alternatively, these models can be calibrated with models that capture more detailed physics of the solar corona (such as the AWSOM model), and that can reproduce solar observations in a consistent manner. Without such a calibration, the results of the polytropic models cannot be validated, but they can be wrongly used by others.


title: Alternate Product Adjacencies in Digital Topology

abstract: We study properties of Cartesian products of digital images, using a variety of adjacencies that have appeared in the literature.


title: Transverse single-spin asymmetries in $\ell \,p^\uparrow \to h \,X$ within a TMD approach: Role of quasireal photon exchange

abstract: We present an updated study of transverse single-spin asymmetries for the inclusive large-$P_T$ processes $\ell \, p^\uparrow \to h\, X$ and $\ell\, p^\uparrow \to {\rm jet}\,X$, within a transverse momentum-dependent approach, including the contribution of quasireal (Weizs\"acker-Williams) photons. In the spirit of a unified transverse momentum-dependent scheme, predictions are obtained adopting the Sivers and transversity distributions and the Collins fragmentation functions as extracted from fits to the azimuthal asymmetries measured in semi-inclusive deep inelastic scattering and $e^+e^-$ annihilation processes. The description of the available data is extremely good, showing a clear general improvement with respect to the previous leading-order analysis. Predictions for unpolarized cross sections and single-spin asymmetries for ongoing and future experiments are also given.


title: Thick Subcategories of the stable category of modules over the exterior algebra

abstract: We study thick subcategories defined by modules of complexity one in $\underline{\md}R$, where $R$ is the exterior algebra in $n+1$ indeterminates.


title: The Diverse Club: The Integrative Core of Complex Networks

abstract: A complex system can be represented and analyzed as a network, where nodes represent the units of the network and edges represent connections between those units. For example, a brain network represents neurons as nodes and axons between neurons as edges. In many networks, some nodes have a disproportionately high number of edges. These nodes also have many edges between each other, and are referred to as the rich club. In many different networks, the nodes of this club are assumed to support global network integration. However, another set of nodes potentially exhibits a connectivity structure that is more advantageous to global network integration. Here, in a myriad of different biological and man-made networks, we discover the diverse club--a set of nodes that have edges diversely distributed across the network. The diverse club exhibits, to a greater extent than the rich club, properties consistent with an integrative network function--these nodes are more highly interconnected and their edges are more critical for efficient global integration. Moreover, we present a generative evolutionary network model that produces networks with a diverse club but not a rich club, thus demonstrating that these two clubs potentially evolved via distinct selection pressures. Given the variety of different networks that we analyzed--the c. elegans, the macaque brain, the human brain, the United States power grid, and global air traffic--the diverse club appears to be ubiquitous in complex networks. These results warrant the distinction and analysis of two critical clubs of nodes in all complex systems.


title: SPIRITS: Uncovering Unusual Infrared Transients With Spitzer

abstract: We present an ongoing, systematic search for extragalactic infrared transients, dubbed SPIRITS --- SPitzer InfraRed Intensive Transients Survey. In the first year, using Spitzer/IRAC, we searched 190 nearby galaxies with cadence baselines of one month and six months. We discovered over 1958 variables and 43 transients. Here, we describe the survey design and highlight 14 unusual infrared transients with no optical counterparts to deep limits, which we refer to as SPRITEs (eSPecially Red Intermediate Luminosity Transient Events). SPRITEs are in the infrared luminosity gap between novae and supernovae, with [4.5] absolute magnitudes between -11 and -14 (Vega-mag) and [3.6]-[4.5] colors between 0.3 mag and 1.6 mag. The photometric evolution of SPRITEs is diverse, ranging from < 0.1 mag/yr to > 7 mag/yr. SPRITEs occur in star-forming galaxies. We present an in-depth study of one of them, SPIRITS 14ajc in Messier 83, which shows shock-excited molecular hydrogen emission. This shock may have been triggered by the dynamic decay of a non-hierarchical system of massive stars that led to either the formation of a binary or a proto-stellar merger.


title: A Rough Path Perspective on Renormalization

abstract: We develop the algebraic theory of rough path translation. Particular attention is given to the case of branched rough paths, whose underlying algebraic structure (Connes-Kreimer, Grossman-Larson) makes it a useful model case of a regularity structure in the sense of Hairer. Pre-Lie structures are seen to play a fundamental rule which allow a direct understanding of the translated (i.e. renormalized) equation under consideration. This construction is also novel with regard to the algebraic renormalization theory for regularity structures due to Bruned--Hairer--Zambotti (2016), the links with which are discussed in detail.


title: Measurements of minimum entropy at K-shell photoelectric effect provides the radius and shape of the proton

abstract: Constructive resonances with minimum entropy are derived from measurements of photoelectric effect cross sections at atomic K-shell thresholds. A resonance region with optimal constructive interference is given by a principal wave length {\lambda} of the order of Bohr atom radius. Our study shows that the proton shape is not a sphere but it has an elliptical volumetric shape with two equal axes and the other different. Resonance waves allowed us the possibility to measure both proton radius and shape through an interference term. This last, was a necessary condition in order to have an effective cross section maximum at threshold. The minimum entropy means minimum shape deformation and it was found to be (0.830 +/- 0.015) fm and the average proton radius was found to be (0.851 +/- 0.201) fm.


title: Perfect Sequences and Arrays over the Unit Quaternions

abstract: We introduce several new constructions for perfect periodic autocorrelation sequences and arrays over the unit quaternions. This paper uses both mathematical proofs and com- puter experiments to prove the (bounded) array constructions have perfect periodic auto- correlation. Furthermore, the first sequence construction generates odd-perfect sequences of unbounded lengths, with good ZCZ.


title: Regularized Reduced Order Models for a Stochastic Burgers Equation

abstract: In this paper, we study the numerical stability of reduced order models for convection-dominated stochastic systems in a relatively simple setting: a stochastic Burgers equation with linear multiplicative noise. Our preliminary results suggest that, in a convection-dominated regime, standard reduced order models yield inaccurate results in the form of spurious numerical oscillations. To alleviate these oscillations, we use the Leray reduced order model, which increases the numerical stability of the standard model by smoothing (regularizing) the convective term with an explicit spatial filter. The Leray reduced order model yields significantly better results than the standard reduced order model and is more robust with respect to changes in the strength of the noise.


title: Adaptive Real-Time Software Defined MIMO Visible Light Communications using Spatial Multiplexing and Spatial Diversity

abstract: In this paper, we experimentally demonstrate a real-time software defined multiple input multiple output (MIMO) visible light communication (VLC) system employing link adaptation of spatial multiplexing and spatial diversity. Real-time MIMO signal processing is implemented by using the Field Programmable Gate Array (FPGA) based Universal Software Radio Peripheral (USRP) devices. Software defined implantation of MIMO VLC can assist in enabling an adaptive and reconfigurable communication system without hardware changes. We measured the error vector magnitude (EVM), bit error rate (BER) and spectral efficiency performance for single carrier M-QAM MIMO VLC using spatial diversity and spatial multiplexing. Results show that spatial diversity MIMO VLC improves error performance at the cost of spectral efficiency that spatial multiplexing should enhance. We propose the adaptive MIMO solution that both modulation schema and MIMO schema are dynamically adapted to the changing channel conditions for enhancing the error performance and spectral efficiency. The average error-free spectral efficiency of adaptive 2x2 MIMO VLC achieved 12 b/s/Hz over 2 meters indoor dynamic transmission.


title: Poisson distribution for gaps between sums of two squares and level spacings for toral point scatterers

abstract: We investigate the level spacing distribution for the quantum spectrum of the square billiard. Extending work of Connors--Keating, and Smilansky, we formulate an analog of the Hardy--Littlewood prime $k$-tuple conjecture for sums of two squares, and show that it implies that the spectral gaps, after removing degeneracies and rescaling, are Poisson distributed. Consequently, by work of Rudnick and Uebersch\"ar, the level spacings of arithmetic toral point scatterers, in the weak coupling limit, are also Poisson distributed. We also give numerical evidence for the conjecture and its implications.

categories: math-ph math.AP math.MP math.NT math.SP nlin.CD


title: Examples of renormalized SDEs

abstract: We demonstrate two examples of stochastic processes whose lifts to geometric rough paths require a renormalisation procedure to obtain convergence in rough path topologies. Our first example involves a physical Brownian motion subject to a magnetic force which dominates over the friction forces in the small mass limit. Our second example involves a lead-lag process of discretised fractional Brownian motion with Hurst parameter $H \in (1/4,1/2)$, in which the stochastic area captures the quadratic variation of the process. In both examples, a renormalisation of the second iterated integral is needed to ensure convergence of the processes, and we comment on how this procedure mimics negative renormalisation arising in the study of singular SPDEs and regularity structures.


title: On the Decay of Strong Magnetization in Global Disc Simulations with Toroidal Fields

abstract: Strong magnetization in accretion discs could resolve a number of outstanding issues related to stability and state transitions in low-mass X-ray binaries. However, it is unclear how real discs become strongly magnetized and, even if they do, whether they can remain in such a state. In this paper, we address the latter issue through a pair of global disc simulations. Here, we only consider cases of initially purely toroidal magnetic fields contained entirely within a compact torus. We find that, over only a few tens of orbital periods, the magnetization of an initially strongly magnetized disc, $P_\mathrm{mag}/P_\mathrm{gas} \ge 10$, drops to $\lesssim 0.1$, similar to the steady-state value reached in initially weakly magnetized discs. This is consistent with recent shearing box simulations with initially strong toroidal fields, the robust conclusion being that strongly magnetized toroidal fields can not be locally self-sustaining. These results appear to leave net poloidal flux or extended radial fields as the only avenues for establishing strongly magnetized discs, ruling out the thermal collapse scenario.


title: An Approach to Non-Abelian Cyclotomic Fields

abstract: We mainly study a polynomial $f_{1,n}(x)=x^{n-1} + 2x^{n-2} + 3x^{n-3} + \cdots + kx^{n-k} + \cdots + (n-1)x + n$ over $\mathbb{Z}$ and the Galois group of the minimal splitting field. First, we show that an arbitrary root $\alpha_{n}$ of $f_{1,n}(x)$ satisfies $|\alpha_{n}|\to 1$ ($n\to \infty$), and discuss the irreducibility of $f_{1,n}(x)$ over $\mathbb{Z}$ for several type $n$. After that, we show that the Galois group of $f_{1,n}(x)$ is Symmetric group $S_{n-1}$ for several type $n$. Although those roots of $f_{1,n}(x)=0$ don't draw an exact circle, it looks like a circle on complex plane. Moreover by considering that Galois groups of $f_{1,n}(x)$ are not abelian in many cases, we call such extension fields over $\mathbb{Q}$ "Non-Abelian Cycrotomic Fields" here.


title: The World's First Real-Time Testbed for Massive MIMO: Design, Implementation, and Validation

abstract: This paper sets up a framework for designing a massive multiple-input multiple-output (MIMO) testbed by investigating hardware (HW) and system-level requirements such as processing complexity, duplexing mode and frame structure. Taking these into account, a generic system and processing partitioning is proposed which allows flexible scaling and processing distribution onto a multitude of physically separated devices. Based on the given HW constraints such as maximum number of links and maximum throughput for peer-to-peer interconnections combined with processing capabilities, the framework allows to evaluate modular HW components. To verify our design approach, we present the LuMaMi (Lund University Massive MIMO) testbed which constitutes the first reconfigurable real-time HW platform for prototyping massive MIMO. Utilizing up to 100 base station antennas and more than 50 Field Programmable Gate Arrays, up to 12 user equipments are served on the same time/frequency resource using an LTE-like Orthogonal Frequency Division Multiplexing time-division duplex-based transmission scheme. Proof-of-concept tests with this system show that massive MIMO can simultaneously serve a multitude of users in a static indoor and static outdoor environment utilizing the same time/frequency resource.


title: Noncommutative cosmological model in the presence of a phantom fluid

abstract: We study noncommutative classical Friedmann-Robertson-Walker cosmological models. The constant curvature of the spatial sections can be positive ($k=1$), negative ($k=-1$) or zero ($k=0$). The matter is represented by a perfect fluid with negative pressure, phantom fluid, which satisfies the equation of state $p =\alpha \rho$ , with $\alpha < - 1$, where $p$ is the pressure and $\rho$ is the energy density. We use Schutz's formalism in order to write the perfect fluid Hamiltonian. The noncommutativity is introduced by nontrivial Poisson brackets between few variables of the models. In order to recover a description in terms of commutative variables, we introduce variables transformations that depend on a noncommutative parameter ($\gamma$). The main motivation for the introduction of the noncommutativity is trying to explain the present accelerated expansion of the universe. We obtain the dynamical equations for these models and solve them. The solutions have four constants: $\gamma$, a parameter associated with the fluid energy $C$, $k$, $\alpha$ and the initial conditions of the models variables. For each value of $\alpha$, we obtain different equations of motion. Then, we compare the evolution of the universe between the present noncommutative models and the corresponding commutative ones ($\gamma \to 0$). The results show that $\gamma$ is very useful for describing an accelerating universe. We estimate the value of $\gamma$, for the present conditions of the Universe. Then, using that value of $\gamma$, in one of the noncommutative cosmological models, we compute the amount of time this universe would take to reach the {\it big rip}.


title: K\"ahler groups and subdirect products of surface groups

abstract: We present a construction that produces infinite classes of K\"ahler groups that arise as fundamental groups of fibres of maps to higher dimensional tori. Following the work of Delzant and Gromov, there is great interest in knowing which subgroups of direct products of surface groups are K\"ahler. We apply our construction to obtain new classes of irreducible, coabelian K\"ahler subgroups of direct products of $r$ surface groups. These cover the full range of possible finiteness properties of irreducible subgroups of direct products of $r$ surface groups: For any $r\geq 3$ and $2\leq k \leq r-1$, our classes of subgroups contain K\"ahler groups that have a classifying space with finite $k$-skeleton while not having a classifying space with finitely many $(k+1)$-cells. We also address the converse question of finding constraints on K\"ahler subdirect products of surface groups and, more generally, on homomorphisms from K\"ahler groups to direct products of surface groups. We show that if a K\"ahler subdirect product of $r$ surface groups admits a classifying space with finite $k$-skeleton for $k>\frac{r}{2}$, then it is virtually the kernel of an epimorphism from a direct product of surface groups onto a free abelian group of even rank.


title: Effect of Cell-Selection on the Effective Fading Distribution in a Downlink K-tier HetNet

abstract: This paper characterizes the statistics of effective fading gain in multi-tier cellular networks with strongest base station (BS) cell association policy. First, we derive the probability of association with the $n$-th nearest BS in the $k$-th tier. Next, we use this result to derive the probability density function (PDF) of the channel fading gain (effective fading) experienced by the user when associating with the strongest BS. Interestingly, our results show that the effective channel gain distribution solely depends upon the original channel fading and the path-loss exponent. Moreover, we show that in the case of Nakagami-$m$ fading channels (Gamma distribution), the distribution of the effective fading is also Gamma but with a gain of $\frac{\alpha}{2}$ in the shape parameter, where $\alpha$ is the path-loss exponent.


title: Singular Limit of BSDEs and Optimal Control of two Scale Stochastic Systems in Infinite Dimensional Spaces

abstract: In this paper we study, by probabilistic techniques, the convergence of the value function for a two-scale, infinite-dimensional, stochastic controlled system as the ratio between the two evolution speeds diverges. The value function is represented as the solution of a backward stochastic differential equation (BSDE) that it is shown to converge towards a reduced BSDE. The noise is assumed to be additive both in the slow and the fast equations for the state. Some non degeneracy condition on the slow equation are required. The limit BSDE involves the solution of an ergodic BSDE and is itself interpreted as the value function of an auxiliary stochastic control problem on a reduced state space.


title: Quaternions in collective dynamics

abstract: We introduce a model of multi-agent dynamics for self-organised motion; individuals travel at a constant speed while trying to adopt the averaged body attitude of their neighbours. The body attitudes are represented through unitary quaternions. We prove the correspondance with the model presented in a previous work by the three first authors where the body attitudes are represented by rotation matrices. Differently from this previous work, the individual based model (IBM) introduced here is based on nematic (rather than polar) alignment. From the IBM, the kinetic and macroscopic equations are derived. The benefit of this approach is twofold: firstly, it allows for a better understanding of the macroscopic equations obtained and, secondly, these equations are prone to numerical studies, which is key for applications.


title: Efficient Multi-Dimensional Mapping Using QAM Constellations for BICM-ID

abstract: Bit-interleaved coded modulation with iterative decoding (BICM-ID) offers very good error performance over additive white Gaussian noise (AWGN) and fading channels if it uses a wisely designed signal mapping. Further, error performance of BICM-ID can significantly be improved by employing multi-dimensional (MD) modulation. However, suitable MD mappings are obtained by computer search techniques except for MD modulations that use smaller constellation e.g., binary phase shift keying (BPSK), quadrature phase shift keying (QPSK) and $8$-ary phase shift keying ($8$-PSK) as basic modulation. The alphabet size of MD modulations increases exponentially as the order of the basic modulation increases and computer search becomes intractable. In this paper, we propose a systematic mapping method for MD modulations. The innovativeness of our proposed method is that it generates MD mappings using $16$- and $64$-quadrature amplitude modulation (QAM) very efficiently. The presented numerical results show that the proposed method improves bit error rate (BER) of BICM-ID.


title: The Dynamics of Wave-Particle Duality

abstract: Both classical and wave-mechanical monochromatic waves may be treated in terms of exact ray-trajectories (encoded in the structure itself of Helmholtz-like equations) whose mutual coupling is the one and only cause of any diffraction and interference process. In the case of Wave Mechanics, de Broglie's merging of Maupertuis's and Fermat's principles (see Section 3) provides, without resorting to the probability-based guidance-laws and flow-lines of the Bohmian theory, the simple law addressing particles along the Helmholtz rays of the relevant matter waves. The purpose of the present research was to derive the exact Hamiltonian ray-trajectory systems concerning, respectively, classical electromagnetic waves, non-relativistic matter waves and relativistic matter waves. We faced then, as a typical example, the numerical solution of non-relativistic wave-mechanical equation systems in a number of numerical applications, showing that each particle turns out to "dances a wave-mechanical dance" around its classical trajectory, to which it reduces when the ray-coupling is neglected. Our approach reaches the double goal of a clear insight into the mechanism of wave-particle duality and of a reasonably simple computability. We finally compared our exact dynamical approach, running as close as possible to Classical Mechanics, with the hydrodynamic Bohmian theory, based on fluid-like "guidance laws".


title: The $N$-Body Problem in Spaces with Uniformly Varying Curvature

abstract: We generalize the curved $N$-body problem to spheres and hyperbolic spheres whose curvature $\kappa$ varies in time. Unlike in the particular case when the curvature is constant, the equations of motion are non-autonomous. We first briefly consider the analogue of the Kepler problem and then investigate the homographic orbits for any number of bodies, proving the existence of several such classes of solutions on spheres. Allowing the curvature to vary in time offers some insight into the effect of an expanding universe, in the context the curved $N$-body problem, when $\kappa$ satisfies Hubble's law. The study of these equations also opens the possibility of finding new connections between classical mechanics and general relativity.


title: Gunrock: GPU Graph Analytics

abstract: For large-scale graph analytics on the GPU, the irregularity of data access and control flow, and the complexity of programming GPUs, have presented two significant challenges to developing a programmable high-performance graph library. "Gunrock", our graph-processing system designed specifically for the GPU, uses a high-level, bulk-synchronous, data-centric abstraction focused on operations on a vertex or edge frontier. Gunrock achieves a balance between performance and expressiveness by coupling high performance GPU computing primitives and optimization strategies with a high-level programming model that allows programmers to quickly develop new graph primitives with small code size and minimal GPU programming knowledge. We characterize the performance of various optimization strategies and evaluate Gunrock's overall performance on different GPU architectures on a wide range of graph primitives that span from traversal-based algorithms and ranking algorithms, to triangle counting and bipartite-graph-based algorithms. The results show that on a single GPU, Gunrock has on average at least an order of magnitude speedup over Boost and PowerGraph, comparable performance to the fastest GPU hardwired primitives and CPU shared-memory graph libraries such as Ligra and Galois, and better performance than any other GPU high-level graph library.


title: The Large D Limit of Planar Diagrams

abstract: We show that in $\text{O}(D)$ invariant matrix theories containing a large number $D$ of complex or Hermitian matrices, one can define a $D\rightarrow\infty$ limit for which the sum over planar diagrams truncates to a tractable, yet non-trivial, sum over melon diagrams. In particular, results obtained recently in SYK and tensor models can be generalized to traditional, string-inspired matrix quantum mechanical models of black holes.


title: The Parsec-Scale Morphology of Southern GPS sources

abstract: Multi-frequency, multi-epoch ATCA observations of a sample of AGN resulted in the identification of 9 new candidate Giga-hertz Peaked Spectrum (GPS) sources. Here we present Long Baseline Array observations at 4.8 GHz of the four candidates with no previously published VLBI image, and consider these together with previously published VLBI images of the other five sources. We find core-jet or compact double morphologies dominate, with further observations required to distinguish between these two possibilities for some sources. One of the nine candidates, PKS 1831-711, displays appreciable variability, suggesting its GPS spectrum is more ephemeral in nature. We focus in particular on the apparent relationship between a narrow spectral width and "compact double" parsec-scale morphology, finding further examples, but also exceptions to this trend. An examination of the VLBI morphologies high-redshift (z>3) sub-class of GPS sources suggests that core-jet morphologies predominate in this class.


title: Follower, Predecessor, and Extender Set Sequences of $\beta$-Shifts

abstract: Given a one-dimensional shift $X$, let $|F_X(n)|$ be the number of follower sets of words of length $n$ in $X$, and $|P_X(n)|$ be the number of predecessor sets of words of length $n$ in $X$. We call the sequence $\{|F_X(n)|\}_{n \in \mathbb{N}}$ the follower set sequence of the shift $X$, and $\{|P_X(n)|\}_{n \in \mathbb{N}}$ the predecessor set sequence of the shift $X$. Extender sets are a generalization of follower sets, and we define the extender set sequence similarly. In this paper, we examine achievable differences in limiting behavior of follower, predecessor, and extender set sequences. This is done through the classical $\beta$-shifts. We show that the follower set sequences of $\beta$-shifts must grow at most linearly in $n$, while the predecessor and extender set sequences may demonstrate exponential growth rate in $n$, depending on choice of $\beta$.


title: The Non-Split Bessel Model on GSp(2n) as an Iwahori-Hecke Algebra Module

abstract: We realize the non-split Bessel model of Novodvorsky and Piatetski-Shapiro as a generalized Gelfand-Graev representation of GSp(4), as suggested by Kawanaka. With uniqueness of the model already established by Novodvorsky and Piatetski-Shapiro, we establish existence of a Bessel model for unramified principal series representations. We then connect the Iwahori-fixed vectors in the Bessel model to a linear character of the Hecke algebra of GSp(4) following the method outlined more generally by Brubaker, Bump, and Friedberg. We use this connection to calculate the image of Iwahori-fixed vectors of unramified principal series in the model, and ultimately provide an explicit alternator expression for the spherical vector in the model. We show that the resulting alternator expression matches previous results of Bump, Friedberg, and Furusawa. We offer the conjecture that a generalized Bessel model on GSp(2n) retains the uniqueness property estabilished in the case when $n=2$; assuming that this conjecture holds, we extend all of the previously mentioned results to the case where $n>2$, including existence of the model for unramified principal series representations.


title: Fast quantum computation at arbitrarily low energy

abstract: One version of the energy-time uncertainty principle states that the minimum time $T_{\perp}$ for a quantum system to evolve from a given state to any orthogonal state is $h/(4 \Delta E)$ where $\Delta E$ is the energy uncertainty. A related bound called the Margolus-Levitin theorem states that $T_{\perp} \geq h/(2 E)$ where E is the expectation value of energy and the ground energy is taken to be zero. Many subsequent works have interpreted $T_{\perp}$ as defining a minimal time for an elementary computational operation and correspondingly a fundamental limit on clock speed determined by a system's energy. Here we present local time-independent Hamiltonians in which computational clock speed becomes arbitrarily large relative to E and $\Delta E$ as the number of computational steps goes to infinity. We argue that energy considerations alone are not sufficient to obtain an upper bound on computational speed, and that additional physical assumptions such as limits to information density and information transmission speed are necessary to obtain such a bound.


title: Global Energetics of Solar Flares: V. Energy Closure in Flares and Coronal Mass Ejections

abstract: In this study we synthesize the results of four previous studies on the global energetics of solar flares and associated coronal mass ejections (CMEs), which include magnetic, thermal, nonthermal, and CME energies in 399 solar M and X-class flare events observed during the first 3.5 years of the Solar Dynamics Observatory (SDO) mission. Our findings are: (1) The sum of the mean nonthermal energy of flare-accelerated particles ($E_{\mathrm{nt}}$), the energy of direct heating ($E_{\mathrm{dir}}$), and the energy in coronal mass ejections ($E_{\mathrm{CME}}$), which are the primary energy dissipation processes in a flare, is found to have a ratio of $(E_{\mathrm{nt}}+E_{\mathrm{dir}}+ E_{\mathrm{CME}})/E_{\mathrm{mag}} = 0.87 \pm 0.18$, compared with the dissipated magnetic free energy $E_{\mathrm{mag}}$, which confirms energy closure within the measurement uncertainties and corroborates the magnetic origin of flares and CMEs; (2) The energy partition of the dissipated magnetic free energy is: $0.51\pm0.17$ in nonthermal energy of $\ge 6$ keV electrons, $0.17\pm0.17$ in nonthermal $\ge 1$ MeV ions, $0.07\pm0.14$ in CMEs, and $0.07\pm0.17$ in direct heating; (3) The thermal energy is almost always less than the nonthermal energy, which is consistent with the thick-target model; (4) The bolometric luminosity in white-light flares is comparable with the thermal energy in soft X-rays (SXR); (5) Solar Energetic Particle (SEP) events carry a fraction $\approx 0.03$ of the CME energy, which is consistent with CME-driven shock acceleration; and (6) The warm-target model predicts a lower limit of the low-energy cutoff at $e_c \approx 6$ keV, based on the mean differential emission measure (DEM) peak temperature of $T_e=8.6$ MK during flares. This work represents the first statistical study that establishes energy closure in solar flare/CME events.


title: The Width Distribution of Loops and Strands in the Solar Corona -- Are we Hitting Rock Bottom ?

abstract: In this study we analyze {\sl Atmospheric Imaging Assembly (AIA)} and Hi-C images in order to investigate absolute limits for the finest loop strands. We develop a model of the occurrence-size distribution function of coronal loop widths, characterized by a lower limit of widths $w_{min}$, a peak width $w_p$, a peak occurrence number $n_p$, and a power law slope $a$. Our data analysis includes automated tracing of curvi-linear features with the OCCULT-2 code, automated sampling of the cross-sectional widths of coronal loops, and fitting of the theoretical size distribution to the observed distribution. With Monte-Carlo simulations and variable pixel sizes $\Delta x$ we derive a first diagnostic criterion to discriminate whether the loop widths are unresolved $(w_p/\Delta x \approx 2.5\pm0.2)$, or fully resolved (if $w_p/\Delta x > 2.7$). For images with resolved loop widths we can apply a second diagnostic criterion that predicts the lower limit of loop widths, $w_{min} \approx 3 (\Delta x_{crit}-0.37")$ as a function of the critical resolution $\Delta x_{crit}$. We find that the loop widths are marginally resolved in AIA images, but are fully resolved in Hi-C images, where our model predicts a lower limit of loop widths at $w_{min} \approx 100$ km and a most frequent (peak) value at $w_p \approx 300$ km, in agreement with recent results of Brooks et al. This result agrees with the statistics of photospheric granulation sizes and thus supports coronal heating mechanisms operating on the macroscopic scale of photospheric magneto-convection, rather than nanoflare heating models with unresolved microscopic scales.


title: A local to global principle for densities over function fields

abstract: Let $d$ be a positive integer and $\mathbb H$ be an integrally closed subring of a global function field $F$. The purpose of this paper is to provide a general sieve method to compute densities of subsets of $\mathbb H^d$ defined by local conditions. The main advantage of the method relies on the fact that one can use results from measure theory to extract density results over $\mathbb H^{d}$. Using this method we are able to compute the density of the set of polynomials with coefficients in $\mathbb H$ which give rise to "good" totally ramified extensions of the global function field $F$. As another application, we give a closed expression for the density of rectangular unimodular matrices with coefficients in $\mathbb H$ in terms of the $L$-polynomial of the function field.


title: Early Science with the Large Millimetre Telescope: Molecules in the Extreme Outflow of a proto-Planetary Nebula

abstract: Extremely high velocity emission likely related to jets is known to occur in some proto-Planetary Nebulae. However, the molecular complexity of this kinematic component is largely unknown. We observed the known extreme outflow from the proto-Planetary Nebula IRAS 16342-3814, a prototype water fountain, in the full frequency range from 73 to 111 GHz with the RSR receiver on the Large Millimetre Telescope. We detected the molecules SiO, HCN, SO, and $^{13}$CO. All molecular transitions, with the exception of the latter are detected for the first time in this source, and all present emission with velocities up to a few hundred km s$^{-1}$. IRAS 16342-3814 is therefore the only source of this kind presenting extreme outflow activity simultaneously in all these molecules, with SO and SiO emission showing the highest velocities found of these species in proto-Planetary Nebulae. To be confirmed is a tentative weak SO component with a FWHM $\sim$ 700 km s$^{-1}$. The extreme outflow gas consists of dense gas (n$_{\rm H_2} >$ 10$^{4.8}$--10$^{5.7}$ cm$^{-3}$), with a mass larger than $\sim$ 0.02--0.15 M$_{\odot}$. The relatively high abundances of SiO and SO may be an indication of an oxygen-rich extreme high velocity gas.


title: Tetravalent $2$-arc-transitive Cayley graphs on non-abelian simple groups

abstract: A graph Gamma is said to be 2-arc-transitive if its full automorphism group Aut(\Gamma) has a single orbit on ordered paths of length 2, and for G\leq Aut(\Gamma), \Gamma is G-regular if G is regular on the vertex set of \Gamma. Let G be a finite non-abelian simple group and let \Gamma be a connected tetravalent 2-arc-transitive G-regular graph. In 2004, Fang, Li and Xu proved that either G\unlhd \Aut(\Gamma) or G is one of 22 possible candidates. In this paper, the number of candidates is reduced to 7, and for each candidate G, it is shown that \Aut(\Gamma) has a normal arc-transitive non-abelian simple subgroup T such that G\leq T and the pair (G,T) is explicitly given


title: Vietoris-type Topologies on Hyperspaces

abstract: We introduce a new Vietoris-type hypertopology by means of the upper-Vietoris-type hypertopology defined by G. Dimov and D. Vakarelov [On Scott consequence systems, Fundamenta Informaticae, 33 (1998), 43-70] (it was called there {\em Tychonoff-type hypertopology}) and the lower-Vietoris-type hypertopology introduced by E. Ivanova-Dimova [Lower-Vietoris-type topologies on hyperspaces, Topology Appl. (to appear)]. We study this new Vietoris-type hypertopology and show that it is, in general, different from the Vietoris topology. Also, some of the results of E. Michael [Topologies on spaces of subsets, Trans. Amer. Math. Soc. 71 (1951), 152-182] about hyperspaces with Vietoris topology are extended to analogous results for hyperspaces with Vietoris-type topology. We obtain as well some results about hyperspaces with Vietoris-type topology which concern some problems analogous to those regarded by H.-J. Schmidt [Hyperspaces of quotient and subspaces. I. Hausdorff topological spaces, Math. Nachr. 104 (1981), 271-280].


title: Row bounds needed to justifiably express flagged Schur functions with Gessel-Viennot determinants

abstract: Let $\lambda$ be a partition with no more than $n$ parts. Let $\beta$ be a weakly increasing $n$-tuple with entries from $\{ 1, ... , n \}$. The flagged Schur function in the variables $x_1, ... , x_n$ that is indexed by $\lambda$ and $\beta$ has been defined to be the sum of the content weight monomials for the semistandard Young tableaux of shape $\lambda$ whose values are row-wise bounded by the entries of $\beta$. Gessel and Viennot gave a determinant expression for the flagged Schur function indexed by $\lambda$ and $\beta$; this could be done since the pair $(\lambda, \beta)$ satisfied their "nonpermutable" condition for the sequence of terminals of an $n$-tuple of lattice paths that they used to model the tableaux. We generalize flagged Schur functions by dropping the requirement that $\beta$ be weakly increasing. Then for each $\lambda$ we give a condition on the entries of $\beta$ for the pair $(\lambda, \beta)$ to be nonpermutable that is both necessary and sufficient. When the parts of $\lambda$ are not distinct there will be multiple row bound $n$-tuples $\beta$ that will produce the same set of tableaux. We accordingly group the bounding $\beta$ into equivalence classes and identify the most efficient $\beta$ in each class for the determinant computation. We recently showed that many other sets of objects that are indexed by $n$ and $\lambda$ are enumerated by the number of these efficient $n$-tuples. We called these counts "parabolic Catalan numbers". It is noted that the $GL(n)$ Demazure characters (key polynomials) indexed by 312-avoiding permutations can also be expressed with these determinants.


title: Localization and Stationary Phase Approximation on Supermanifolds

abstract: Given an odd vector field $Q$ on a supermanifold $M$ and a $Q$-invariant density $\mu$ on $M$, under certain compactness conditions on $Q$, the value of the integral $\int_{M}\mu$ is determined by the value of $\mu$ on any neighborhood of the vanishing locus $N$ of $Q$. We present a formula for the integral in the case where $N$ is a subsupermanifold which is appropriately non-degenerate with respect to $Q$. In the process, we discuss the linear algebra necessary to express our result in a coordinate independent way. We also extend stationary phase approximation and the Morse-Bott Lemma to supermanifolds.


title: High-energy emission as signature of magnetic field amplification in Neutron Star Mergers

abstract: The merger of a binary neutron star system is suggested as the central engine of short gamma-ray bursts (sGRBs). For the merger process, simulations predict that magnetic field is amplified beyond magnetar field strength by Kelvin-Helmholtz instabilities. With the Large Area Telescope (LAT), bursts have been detected that show a temporally extended component in coincidence with a short-lasting peak at the end of the prompt phase. We show that the presence of these LAT components in a sGRB could provide evidence of magnetic field amplification in the neutron star merger.


title: Efficient asymptotic variance reduction when estimating volatility in high frequency data

abstract: This paper shows how to carry out efficient asymptotic variance reduction when estimating volatility in the presence of stochastic volatility and microstructure noise with the realized kernels (RK) from [Barndorff-Nielsen et al., 2008] and the quasi-maximum likelihood estimator (QMLE) studied in [Xiu, 2010]. To obtain such a reduction, we chop the data into B blocks, compute the RK (or QMLE) on each block, and aggregate the block estimates. The ratio of asymptotic variance over the bound of asymptotic efficiency converges as B increases to the ratio in the parametric version of the problem, i.e. 1.0025 in the case of the fastest RK Tukey-Hanning 16 and 1 for the QMLE. The impact of stochastic sampling times and jump in the price process is examined carefully. The finite sample performance of both estimators is investigated in simulations, while empirical work illustrates the gain in practice.


title: Towards entanglement entropy with UV cutoff in conformal nets

abstract: We consider the entanglement entropy for a spacetime region and its spacelike complement in the framework of algebraic quantum field theory. For a M\"obius covariant local net satisfying a certain nuclearity property, we consider the von Neumann entropy for type I factors between local algebras and introduce an entropic quantity. Then we implement a cutoff on this quantity with respect to the conformal Hamiltonian and show that it remains finite as the distance of two intervals tends to zero. We compare our definition to others in the literature.


title: Pentavalent symmetric graphs admitting vertex-transitive non-abelian simple groups

abstract: A graph \Gamma is said to be {\em symmetric} if its automorphism group \Aut(\Gamma) is transitive on the arc set of \Gamma. Let $G$ be a finite non-abelian simple group and let \Gamma be a connected pentavalent symmetric graph such that G\leq \Aut(\Gamma). In this paper, we show that if $G$ is transitive on the vertex set of \Gamma, then either G\unlhd \Aut(\Gamma) or \Aut(\Gamma) contains a non-abelian simple normal subgroup $T$ such that $G\leq T$ and $(G,T)$ is one of $58$ possible pairs of non-abelian simple groups. In particular, if $G$ is arc-transitive, then $(G,T)$ is one of $17$ possible pairs, and if $G$ is regular on the vertex set of \Gamma, then $(G,T)$ is one of $13$ possible pairs, which improves the result on pentavalent symmetric Cayley graph given by Fang, Ma and Wang in 2011.


title: Current-phase relations of few-mode InAs nanowire Josephson junctions

abstract: Gate-tunable semiconductor nanowires with superconducting leads have great potential for quantum computation and as model systems for mesoscopic Josephson junctions. The supercurrent, $I$, versus the phase, $\phi$, across the junction is called the current-phase relation (CPR). It can reveal not only the amplitude of the critical current, but also the number of modes and their transmission. We measured the CPR of many individual InAs nanowire Josephson junctions, one junction at a time. Both the amplitude and shape of the CPR varied between junctions, with small critical currents and skewed CPRs indicating few-mode junctions with high transmissions. In a gate-tunable junction, we found that the CPR varied with gate voltage: Near the onset of supercurrent, we observed behavior consistent with resonant tunneling through a single, highly transmitting mode. The gate dependence is consistent with modeled subband structure that includes an effective tunneling barrier due to an abrupt change in the Fermi level at the boundary of the gate-tuned region. These measurements of skewed, tunable, few-mode CPRs are promising both for applications that require anharmonic junctions and for Majorana readout proposals.


title: GPU Multisplit: an extended study of a parallel algorithm

abstract: Multisplit is a broadly useful parallel primitive that permutes its input data into contiguous buckets or bins, where the function that categorizes an element into a bucket is provided by the programmer. Due to the lack of an efficient multisplit on GPUs, programmers often choose to implement multisplit with a sort. One way is to first generate an auxiliary array of bucket IDs and then sort input data based on it. In case smaller indexed buckets possess smaller valued keys, another way for multisplit is to directly sort input data. Both methods are inefficient and require more work than necessary: the former requires more expensive data movements while the latter spends unnecessary effort in sorting elements within each bucket. In this work, we provide a parallel model and multiple implementations for the multisplit problem. Our principal focus is multisplit for a small (up to 256) number of buckets. We use warp-synchronous programming models and emphasize warp-wide communications to avoid branch divergence and reduce memory usage. We also hierarchically reorder input elements to achieve better coalescing of global memory accesses. On a GeForce GTX 1080 GPU, we can reach a peak throughput of 18.93 Gkeys/s (or 11.68 Gpairs/s) for a key-only (or key-value) multisplit. Finally, we demonstrate how multisplit can be used as a building block for radix sort. In our multisplit-based sort implementation, we achieve comparable performance to the fastest GPU sort routines, sorting 32-bit keys (and key-value pairs) with a throughput of 3.0 G keys/s (and 2.1 Gpair/s).


title: On the number of gapped repeats with arbitrary gap

abstract: For any functions $f(x)$, $g(x)$ from $\mathbb {N}$ to $\mathbb {R}$ we call repeats $uvu$ such that $g(|u|)\le |v|\le f(|u|)$ as {\it $f,g$-gapped repeats}. We study the possible number of $f,g$-gapped repeats in words of fixed length~$n$. For quite weak conditions on $f(x)$, $g(x)$ we obtain an upper bound on this number which is linear to~$n$.


title: Heavy quarks within the electroweak multiplet

abstract: Standard-model fields and their associated electroweak Lagrangian are equivalently expressed in a shared spin basis. The scalar-vector terms are written with scalar-operator components acting on quark-doublet elements, and shown to be parametrization-invariant. Such terms, and the t- and b-quark Yukawa terms are linked by the identification of the common mass-generating Higgs operating upon the other fields, after acquiring a vacuum expectation value $v$. Thus, the customary vector masses are related to the fermions', fixing the t-quark mass $m_t$ with the relation $m^2_t+m^2_b=v^2/2$ either for maximal hierarchy, or given the b-quark mass $m_b$, implying $m_t \simeq 173.9$ GeV, for $v=246$ GeV. A sum rule is derived for all quark masses that generalizes this restriction. An interpretation follows that electroweak bosons and heavy quarks belong in a multiplet.


title: Liquid-like thermal conduction in a crystalline solid

abstract: A solid conducts heat through both transverse and longitudinal acoustic phonons, but a liquid employs only longitudinal vibrations. Here, we report that the crystalline solid AgCrSe2 has liquid-like thermal conduction. In this compound, Ag atoms exhibit a dynamic duality that they are exclusively involved in intense low-lying transverse acoustic phonons while they also undergo local fluctuations inherent in an order-to-disorder transition occurring at 450 K. As a consequence of this extreme disorder-phonon coupling, transverse acoustic phonons become damped as approaching the transition temperature, above which they are not defined anymore because their lifetime is shorter than the relaxation time of local fluctuations. Nevertheless, the damped longitudinal acoustic phonon survives for thermal transport. This microscopic insight might reshape the fundamental idea on thermal transport properties of matter and facilitates the optimization of thermoelectrics.


title: Graphene and its elemental analogue: A molecular dynamics view of fracture phenomenon

abstract: Graphene and some graphene like two dimensional materials; hexagonal boron nitride (hBN) and silicene have unique mechanical properties which severely limit the suitability of conventional theories used for common brittle and ductile materials to predict the fracture response of these materials. This study revealed the fracture response of graphene, hBN and silicene nanosheets under different tiny crack lengths by molecular dynamics (MD) simulations using LAMMPS. The useful strength of these large area two dimensional materials are determined by their fracture toughness. Our study shows a comparative analysis of mechanical properties among the elemental analogues of graphene and suggested that hBN can be a good substitute for graphene in terms of mechanical properties. We have also found that the pre-cracked sheets fail in brittle manner and their failure is governed by the strength of the atomic bonds at the crack tip. The MD prediction of fracture toughness shows significant difference with the fracture toughness determined by Griffth's theory of brittle failure which restricts the applicability of Griffith's criterion for these materials in case of nano-cracks. Moreover, the strengths measured in armchair and zigzag directions of nanosheets of these materials implied that the bonds in armchair direction has the stronger capability to resist crack propagation compared to zigzag direction.


title: Bessel bridge representation for heat kernel in hyperbolic space

abstract: This article shows a Bessel bridge representation for the transition density of Brownian motion on the Poincare space. This transition density is also referred to as the heat kernel on the hyperbolic space in differential geometry literature. The representation recovers the well-known closed form expression for the heat kernel on hyperbolic space in dimension three. However, the newly derived bridge representation is different from the McKean kernel in dimension two and from the Gruet's formula in higher dimensions. The methodology is also applicable to the derivation of an analogous Bessel bridge representations for heat kernel on a Cartan-Hadamard radially symmetric space and for the transition density of hyperbolic Bessel process.


title: Low Complexity Receiver for Uplink SCMA System via Expectation Propagation

abstract: Sparse code multiple access (SCMA) scheme is considered to be one promising non-orthogonal multiple access technology for the future fifth generation (5G) communications. Due to the sparse nature, message passing algorithm (MPA) has been used as the receiver to achieve close to maximum likelihood (ML) detection performance with much lower complexity. However, the complexity order of MPA is still exponential with the size of codebook and the degree of signal superposition on a given resource element. In this paper, we propose a novel low complexity iterative receiver based on expectation propagation algorithm (EPA), which reduces the complexity order from exponential to linear. Simulation results demonstrate that the proposed EPA receiver achieves nearly the same block error rate (BLER) performance as the conventional message passing algorithm (MPA) receiver with orders less complexity.


title: Seifert manifolds admitting partially hyperbolic diffeomorphisms

abstract: We characterize which 3-dimensional Seifert manifolds admit transitive partially hyperbolic diffeomorphisms. In particular, a circle bundle over a higher-genus surface admits a transitive partially hyperbolic diffeomorphism if and only if it admits an Anosov flow.


title: Measurement of higher-order multipole amplitudes in $\psi(3686)\rightarrow\gamma\chi_{c1,2}$ with $\chi_{c1,2}\to\gamma J/\psi$ and search for the transition $\eta_{c}(2S)\to\gamma J/\psi$

abstract: Using 106~million $\psi(3686)$ events collected with the BESIII detector, we measure multipole amplitudes for the decay $\psi(3686)\rightarrow\gamma\chi_{c1,2}\to\gamma\gamma J/\psi$ beyond the dominant electric-dipole amplitudes. The normalized magnetic-quadrupole amplitude for $\psi(3686)\rightarrow\gamma\chi_{c1,2}\rightarrow\gamma\gamma J/\psi$ and the normalized electric-octupole amplitudes for $\psi(3686)rightarrow\gamma\chi_{c2}$,~$\chi_{c2}\rightarrow\gamma J/\psi$ are determined. The M2 amplitudes for $\psi(3686)\rightarrow\gamma\chi_{c1}$ and $\chi_{c1,2}\rightarrow\gamma J/\psi$ are found to differ significantly from zero and are consistent with theoretical predictions. We also obtain the ratios of M2 contributions of $\psi(3686)$ and $J/\psi$ decays to $\chi_{c1,2}$, $b_{2}^{1}/b_{2}^{2} = 1.35\pm0.72$ and $a_{2}^{1}/a_{2}^{2} = 0.617\pm0.083$, which agree well with theoretical expectations. By considering the multipole contributions of $\chi_{c1,2}$, we measure the product branching fractions for the cascade decays $\psi(3686)\rightarrow\gamma\chi_{c0,1,2}\to\gamma\gamma J/\psi$ and search for the process $\eta_{c}(2S)\to\gamma J/\psi$ through $\psi(3686)\rightarrow\gamma\eta_{c}(2S)$. The product branching fraction for $\psi(3686)\rightarrow\gamma\chi_{c0}\to\gamma\gamma J/\psi$ is 3$\sigma$ larger than published measurements, while those of $\psi(3686)rightarrow\gamma\chi_{c1,2}\to\gamma\gamma J/\psi$ are consistent. No significant signal for the decay $\psi(3686)\to\gamma \eta_c(2S)\to\gamma \gamma J/\psi$ is observed, and the upper limit of the product branching fraction at the 90\% confidence level is determined.


title: Enhancing quantum control by bootstrapping a quantum processor of 12 qubits

abstract: Accurate and efficient control of quantum systems is one of the central challenges for quantum information processing. Current state-of-the-art experiments rarely go beyond 10 qubits and in most cases demonstrate only limited control. Here we demonstrate control of a 12-qubit system, and show that the system can be employed as a quantum processor to optimize its own control sequence by using measurement-based feedback control (MQFC). The final product is a control sequence for a complex 12-qubit task: preparation of a 12-coherent state. The control sequence is about 10% more accurate than the one generated by the standard (classical) technique, showing that MQFC can correct for unknown imperfections. Apart from demonstrating a high level of control over a relatively large system, our results show that even at the 12-qubit level, a quantum processor can be a useful lab instrument. As an extension of our work, we propose a method for combining the MQFC technique with a twirling protocol, to optimize the control sequence that produces a desired Clifford gate.


title: Generalized Minimum Distance Estimators in Linear Regression with Dependent Errors

abstract: This paper discusses minimum distance estimation method in the linear regression model with dependent errors which are strongly mixing. The regression parameters are estimated through the minimum distance estimation method, and asymptotic distributional properties of the estimators are discussed. A simulation study compares the performance of the minimum distance estimator with other well celebrated estimator. This simulation study shows the superiority of the minimum distance estimator over another estimator. KoulMde (R package) which was used for the simulation study is available online. See section 4 for the detail.


title: A renormalization method for three-boson system with a triboson field

abstract: A renormalization method that introduces an auxiliary field to represent a three-body bound state is studied in a three-boson system with a triboson field. A cutoff dependence in the three-boson system emerges as a limit cycle, and the cyclic singularity is renormalized by employing two methods: a standard method and the auxiliary field method. For each method, different sets of diagrams are involved for renormalization, and thus we numerically study three quantities: counter term for renormalization, scattering length of $s$-wave boson-diboson scattering, and normalized wavefunction. We confirm that the two methods would lead to the same result.