title: Survey of Inter-satellite Communication for Small Satellite Systems: Physical Layer to Network Layer View

abstract: Small satellite systems enable whole new class of missions for navigation, communications, remote sensing and scientific research for both civilian and military purposes. As individual spacecraft are limited by the size, mass and power constraints, mass-produced small satellites in large constellations or clusters could be useful in many science missions such as gravity mapping, tracking of forest fires, finding water resources, etc. Constellation of satellites provide improved spatial and temporal resolution of the target. Small satellite constellations contribute innovative applications by replacing a single asset with several very capable spacecraft which opens the door to new applications. With increasing levels of autonomy, there will be a need for remote communication networks to enable communication between spacecraft. These space based networks will need to configure and maintain dynamic routes, manage intermediate nodes, and reconfigure themselves to achieve mission objectives. Hence, inter-satellite communication is a key aspect when satellites fly in formation. In this paper, we present the various researches being conducted in the small satellite community for implementing inter-satellite communications based on the Open System Interconnection (OSI) model. This paper also reviews the various design parameters applicable to the first three layers of the OSI model, i.e., physical, data link and network layer. Based on the survey, we also present a comprehensive list of design parameters useful for achieving inter-satellite communications for multiple small satellite missions. Specific topics include proposed solutions for some of the challenges faced by small satellite systems, enabling operations using a network of small satellites, and some examples of small satellite missions involving formation flying aspects.


title: Asymmetric de Finetti Theorem for Infinite-dimensional Quantum Systems

abstract: The de Finetti representation theorem for continuous variable quantum system is first developed to approximate an N-partite continuous variable quantum state with a convex combination of independent and identical subsystems, which requires the original state to obey permutation symmetry conditioned on successful experimental verification on k of N subsystems. We generalize the de Finetti theorem to include asymmetric bounds on the variance of canonical observables and biased basis selection during the verification step. Our result thereby enables application of infinite-dimensional de Finetti theorem to situations where two conjugate measurements obey different statistics, such as the security analysis of quantum key distribution protocols based on squeezed state against coherent attack.


title: Finite-Dimensional Representations constructed from Random Walks

abstract: Given a $1$-cocycle $b$ with coefficients in an orthogonal representation, we show that any finite dimensional summand of $b$ is cohomologically trivial if and only if $\| b(X_n) \|^2/n$ tends to a constant in probability, where $X_n$ is the trajectory of the random walk $(G,\mu)$. As a corollary, we obtain sufficient conditions for $G$ to satisfy Shalom's property $H_{\mathrm{FD}}$. Another application is a convergence to a constant in probability of $\mu^{*n}(e) -\mu^{*n}(g)$, $n\gg m$, normalized by its average with respect to $\mu^{*m}$, for any finitely generated amenable group without infinite virtually Abelian quotients. Finally, we show that the harmonic equivariant mapping of $G$ to a Hilbert space obtained as an $U$-ultralimit of normalized $\mu^{*n}- g \mu^{*n}$ can depend on the ultrafilter $U$ for some groups.


title: Determination of the proton-to-helium ratio in cosmic rays at ultra-high energies from the tail of the $X_{\rm max}$ distribution

abstract: We present a method to determine the proton-to-helium ratio in cosmic rays at ultra-high energies. It makes use of the exponential slope, $\Lambda$, of the tail of the $X_{\rm max}$ distribution measured by an air shower experiment. The method is quite robust with respect to uncertainties from modeling hadronic interactions and to systematic errors on $X_{\rm max}$ and energy, and to the possible presence of primary nuclei heavier than helium. Obtaining the proton-to-helium ratio with air shower experiments would be a remarkable achievement. To quantify the applicability of a particular mass-sensitive variable for mass composition analysis despite hadronic uncertainties we introduce as a metric the `analysis indicator' and find an improved performance of the $\Lambda$ method compared to other variables currently used in the literature. The fraction of events in the tail of the $X_{\rm max}$ distribution can provide additional information on the presence of nuclei heavier than helium in the primary beam.


title: Lowest Landau level on a cone and zeta determinants

abstract: We consider the integer QH state on Riemann surfaces with conical singularities, with the main objective of detecting the effect of the gravitational anomaly directly from the form of the wave function on a singular geometry. We suggest the formula expressing the normalisation factor of the holomorphic state in terms of the regularized zeta determinant on conical surfaces and check this relation for some model geometries. We also comment on possible extensions of this result to the fractional QH states.


title: Efficient Algorithms for Scheduling Moldable Tasks

abstract: We study the problem of scheduling $n$ independent moldable tasks on $m$ processors that arises in large-scale parallel computations. When tasks are monotonic, the best known result is a $(\frac{3}{2}+\epsilon)$-approximation algorithm for makespan minimization with a complexity linear in $n$ and polynomial in $\log{m}$ and $\frac{1}{\epsilon}$ where $\epsilon$ is arbitrarily small. We propose a new perspective of the existing speedup models: the speedup of a task $T_{j}$ is linear when the number $p$ of assigned processors is small (up to a threshold $\delta_{j}$) while it presents monotonicity when $p$ ranges in $[\delta_{j}, k_{j}]$; the bound $k_{j}$ indicates an unacceptable overhead when parallelizing on too many processors. The generality of this model is proved to be between the classic monotonic and linear-speedup models. For any given integer $\delta\geq 5$, let $u=\left\lceil \sqrt[2]{\delta} \right\rceil-1\geq 2$. In this paper, we propose a $\frac{1}{\theta(\delta)} (1+\epsilon)$-approximation algorithm for makespan minimization where $\theta(\delta) = \frac{u+1}{u+2}\left( 1- \frac{k}{m} \right)$ ($m\gg k$). As a by-product, we also propose a $\theta(\delta)$-approximation algorithm for throughput maximization with a common deadline.


title: Sewing Greenberger-Horne-Zeilinger states with a quantum zipper

abstract: A mechanism of a chiral spin wave rotation is introduced to systematically generate mesoscopic Greenberger-Horne-Zeilinger states.


title: A Lower Bound on the Renormalized Nelson Model

abstract: We provide explicit lower bounds for the ground-state energy of the renormalized Nelson model in terms of the coupling constant $\alpha$ and the number of particles $N$, uniform in the meson mass and valid even in the massless case. In particular, for any number of particles $N$ and large enough $\alpha$ we provide a bound of the form $-C\alpha^2 N^3\log^2(\alpha N)$, where $C$ is an explicit positive numerical constant; and if $\alpha$ is sufficiently small, we give one of the form $-C\alpha^2 N^3\log^2 N$ for $N \geq 2$, and $-C\alpha^2$ for $N = 1$. Whereas it is known that the renormalized Hamiltonian of the Nelson model is bounded below (as realized by E. Nelson) and implicit lower bounds have been given elsewhere (as in a recent work by Gubinelli, Hiroshima, and L\"{o}rinczi), ours seem to be the first fully explicit lower bounds with a reasonable dependence on $\alpha$ and $N$. We emphasize that the logarithmic term in the bounds above is probably an artifact in our calculations, since one would expect that the ground-state energy should behave as $-C\alpha^2 N^3$ for large $N$ or $\alpha$, as in the polaron model of H. Fr\"{o}hlich.


title: The $C \: ^1 \Sigma ^+$, $A \: ^1 \Sigma ^+$, and $b \: ^3 \Pi_{0^+}$ states of LiRb

abstract: We present the first spectroscopic studies of the $C \ ^1\Sigma^+$ electronic state and the $A \ ^1\Sigma^+$ - $b \ ^3\Pi_{0^+}$ complex in $^7$Li - $^{85}$Rb. Using resonantly-enhanced, two-photon ionization, we observed $v = 7$, 9, 12, 13 and $26-44$ of the $C \ ^1\Sigma^+$ state. We augment the REMPI data with a form of depletion spectra in regions of dense spectral lines. The $A \ ^1\Sigma^+$ - $b \ ^3\Pi_{0^+}$ complex was observed with depletion spectroscopy, depleting to vibrational levels $v=0 \rightarrow 29$ of the $A \ ^1\Sigma^+$ state and $v=8 \rightarrow 18$ of the $b \ ^3\Pi_{0^+}$ state. For all three series, we determine the term energy and vibrational constants. Finally, we outline several possible future projects based on the data presented here.


title: Additive Classical Capacity of Quantum Channels Assisted by Noisy Entanglement

abstract: We give a capacity formula for the classical information transmission over a noisy quantum channel, with separable encoding by the sender and limited resources provided by the receiver's pre-shared ancilla. Instead of a pure state, we consider the signal-ancilla pair in a mixed state, purified by a "witness". Thus, the signal-witness correlation limits the resource available from the signal-ancilla correlation. Our formula characterizes the utility of different forms of resources, including noisy or limited entanglement assistance, for classical communication. With separable encoding, the sender's signals across multiple channel uses are still allowed to be entangled, yet our capacity formula is additive. In particular, for generalized covariant channels our capacity formula has a simple closed-form. Moreover, our additive capacity formula upper bounds the general coherent attack's information gain in various two-way quantum key distribution protocols. For Gaussian protocols, the additivity of the formula indicates that the collective Gaussian attack is the most powerful.


title: Multi-graded nilpotent tuples

abstract: We discuss multi-graded nilpotent tuples of multi-graded vector spaces which are a generalization of graded nilpotent pairs. The multi-grading yields a natural notion of a shape of such tuple and our main interest is to answer the question "Is the number of multi-graded nilpotent tuples of a fixed shape, up to base change in the homogeneous components, finite?" Our methods make use of a translation to the class of so-called "Multi-staircase algebras" and we classify their representation types.


title: Capacity Estimates via comparison with TRO channels

abstract: A ternary ring of operators (TRO) in finite dimensions is a diagonal sum of spaces of rectangular matrices. TRO as operator space corresponds to quantum channels that are diagonal sums of partial traces, which we call TRO channels. TRO channels admits simple, single-letter capacity formula. Using operator space and complex interpolation techniques, we give perturbative capacities estimates for a wider class of quantum channels by comparison to TRO channels. Our estimates applies mainly for quantum and private capacity and also strong converse rates. The examples includes random unitary from group representations which in general are non-degradable channels.


title: Distinguishing quantum states using Clifford orbits

abstract: It is a fundamental property of quantum mechanics that information is lost as a result of performing measurements. Indeed, with every quantum measurement one can associate a number -- its POVM norm constant -- that quantifies how much the distinguishability of quantum states degrades in the worst case as a result of the measurement. This raises the obvious question which measurements preserve the most information in these sense of having the largest norm constant. While a number of near-optimal schemes have been found (e.g. the uniform POVM, or complex projective 4-designs), they all seem to be difficult to implement in practice. Here, we analyze the distinguishability of quantum states under measurements that are orbits of the Clifford group. The Clifford group plays an important role e.g. in quantum error correction, and its elements are considered simple to implement. We find that the POVM norm constants of Clifford orbits depend on the effective rank of the states that should be distinguished, as well as on a quantitative measure of the "degree of localization in phase space" of the vectors in the orbit. The most important Clifford orbit is formed by the set of stabilizer states. Our main result implies that stabilizer measurements are essentially optimal for distinguishing pure quantum states. As an auxiliary result, we use the methods developed here to prove new entropic uncertainty relations for stabilizer measurements. This paper is based on a very recent analysis of the representation theory of tensor powers of the Clifford group.


title: $h^\ast$-polynomials of zonotopes

abstract: The Ehrhart polynomial of a lattice polytope $P$ encodes information about the number of integer lattice points in positive integral dilates of $P$. The $h^\ast$-polynomial of $P$ is the numerator polynomial of the generating function of its Ehrhart polynomial. A zonotope is any projection of a higher dimensional cube. We give a combinatorial description of the $h^\ast$-polynomial of a lattice zonotope in terms of refined descent statistics of permutations and prove that the $h^\ast$-polynomial of every lattice zonotope has only real roots and therefore unimodal coefficients. Furthermore, we present a closed formula for the $h^\ast$-polynomial of a zonotope in matroidal terms which is analogous to a result by Stanley (1991) on the Ehrhart polynomial. Our results hold not only for $h^\ast$-polynomials but carry over to general combinatorial positive valuations. Moreover, we give a complete description of the convex hull of all $h^\ast$-polynomials of zonotopes in a given dimension: it is a simplicial cone spanned by refined Eulerian polynomials.


title: Quantitative estimate on singularities in isoperimetric clusters

abstract: We prove a quantitative estimate on the number of certain singularities in almost minimizing clusters. In particular, we consider the singular points belonging to the lowest stratum of the Federer-Almgren stratification (namely, where each tangent cone does not split a $\R$) with maximal density. As a consequence we obtain an estimate on the number of triple junctions in $2$-dimensional clusters and on the number of tetrahedral points in $3$ dimensions, that in turn implies that the boundaries of volume-constrained minimizing clusters form at most a finite number of equivalence classes modulo homeomorphism of the boundary, provided that the prescribed volumes vary in a compact set. The method is quite general and applies also to other problems: for instance, to count the number of singularities in a codimension 1 area-minimizing surface in $\R^8$.


title: Properties of the First-order Fermi acceleration in fast magnetic reconnection driven by turbulence in collisional MHD flows

abstract: Fast magnetic reconnection may occur in different astrophysical sources, producing flare-like emission and particle acceleration. Currently, this process is being studied as an efficient mechanism to accelerate particles via a first-order Fermi process. In this work we analyse the acceleration rate and the energy distribution of test particles injected in three-dimensional magnetohydrodynamical (MHD) domains with large-scale current sheets where reconnection is made fast by the presence of turbulence. We study the dependence of the particle acceleration time with the relevant parameters of the embedded turbulence, i.e., the Alfv\'en speed $V_{\rm A}$, the injection power $P_{\rm inj}$ and scale $k_{\rm inj}$ ($k_{\rm inj} = 1/l_{\rm inj}$). We find that the acceleration time follows a power-law dependence with the particle kinetic energy: $t_{acc} \propto E^{\alpha}$, with $0.2 < \alpha < 0.6$ for a vast range of values of $c / V_{\rm A} \sim 20 - 1000$. The acceleration time decreases with the Alfv\'en speed (and therefore with the reconnection velocity) as expected, having an approximate dependence $t_{acc} \propto (V_{\rm A} / c)^{-\kappa}$, with $\kappa \sim 2.1- 2.4$ for particles reaching kinetic energies between $1 - 100 \, m_p c^2$, respectively. Furthermore, we find that the acceleration time is only weakly dependent on the $P_{\rm inj}$ and $l_{\rm inj}$ parameters of the turbulence. The particle spectrum develops a high-energy tail which can be fitted by a hard power-law already in the early times of the acceleration, in consistency with the results of kinetic studies of particle acceleration by magnetic reconnection in collisionless plasmas.


title: f-Biharmonic and bi-f-harmonic submanifolds of generalized space forms

abstract: We study f-biharmonic and bi-f-harmonic submanifolds in both generalized complex and Sasakian space forms. We prove necessary and sufficient condition for f-biharmonicity and bi-f-harmonicity in the general case and many particular cases. Some non-existence results are also obtained.


title: The range and valence of a real Smirnov function

abstract: We give a complete description of the possible ranges of real Smirnov functions (quotients of two bounded analytic functions on the open unit disk where the denominator is outer and such that the radial boundary values are real almost everywhere on the unit circle). Our techniques use the theory of unbounded symmetric Toeplitz operators, some general theory of unbounded symmetric operators, classical Hardy spaces, and an application of the uniformization theorem. In addition, we completely characterize the possible valences for these real Smirnov functions when the valence is finite. To do so we construct Riemann surfaces we call disk trees by welding together copies of the unit disk and its complement in the Riemann sphere. We also make use of certain trees we call valence trees that mirror the structure of disk trees.


title: On normality of f.pk-structures on g-manifolds

abstract: We consider higher dimensional generalisations of normal almost contact structures, the so called f.pk-structures where parallelism spans a Lie algebra g (f.pk-g-structures). Two types of these structures are discussed. In the first case, we construct an almost complex structure on a product manifold mirroring K-structures. We show that the natural normality condition can be satisfied only when g is abelian. The second case we consider is when the Lie algebra in question is 3-dimensional, but the almost complex structure on a product is constructed in a different manner. In both cases the normality conditions are expressed in terms of the structure tensors.


title: On the number of factorizations of an integer

abstract: Let $f(n)$ denote the number of unordered factorizations of a positive integer $n$ into factors larger than $1$. We show that the number of distinct values of $f(n)$, less than or equal to $x$, is at most $\exp \left( C \sqrt{\frac{\log x}{\log \log x}} \left( 1 + o(1) \right) \right)$, where $C=2\pi\sqrt{2/3}$ and $x$ is sufficiently large. This improves upon a previous result of the first author and F. Luca.


title: A short proof of Minkowski's second theorem and its application to some open questions

abstract: In this paper we present a surprisingly short proof of Minkowski's second theorem. The author hopes there is no mistake in it, though the argument seems to be too plain to contain one. Also, we apply the main construction of the proof to some problems concerning the anomaly of a convex body, the density of the densest lattice packing, lattice point enumerator, and Ehrhart polynomial.


title: Gyromagnetic factor of rotating disks of electrically charged dust in general relativity

abstract: We calculated the dimensionless gyromagnetic ratio ("$g$-factor") of self-gravitating, uniformly rotating disks of dust with a constant specific charge $\epsilon$. These disk solutions to the Einstein-Maxwell equations depend on $\epsilon$ and a "relativity parameter" $\gamma$ ($0<\gamma\le 1$) up to a scaling parameter. Accordingly, the $g$-factor is a function $g=g(\gamma,\epsilon)$. The Newtonian limit is characterized by $\gamma \ll 1$, whereas $\gamma\to 1$ leads to a black-hole limit. The $g$-factor, for all $\epsilon$, approaches the values $g=1$ as $\gamma\to 0$ and $g=2$ as $\gamma\to 1$.


title: Effects of binary stellar populations on direct collapse black hole formation

abstract: The critical Lyman--Werner flux required for direct collapse blackholes (DCBH) formation, or J$_{crit}$, depends on the shape of the irradiating spectral energy distribution (SED). The SEDs employed thus far have been representative of {{realistic}} single stellar populations. We study the effect of binary stellar populations on the formation of DCBH, as a result of their contribution to the Lyman--Werner radiation field. Although binary populations with ages $>$ 10 Myr yield a larger LW photon output, we find that the corresponding values of J$_{crit}$ can be up to 100 times higher than single stellar populations. We attribute this to the shape of the binary SEDs as they produce a sub--critical rate of H$^-$ photodetaching 0.76 eV photons as compared to single stellar populations, reaffirming the role that H$^-$ plays in DCBH formation. This further corroborates the idea that DCBH formation is better understood in terms of a critical region in the H$_2$--H$^-$ photo--destruction rate parameter space, rather than a single value of LW flux.


title: The blazar sequence 2.0

abstract: I discuss the spectral energy distribution (SED) of all blazars with redshift detected by the {\it Fermi} satellite and listed in the 3LAC catalog. I will update the so called "blazar sequence" from the phenomenological point of view, with no theory or modelling. I will show that: i) pure data show that jet and accretion power are related; ii) the updated blazar sequence maintains the properties of the old version, albeit with a less pronounced dominance of the $\gamma$--ray emission; iii) at low bolometric luminosities, two different type of objects have the same high energy power: low black hole mass flat spectrum radio quasars and high mass BL Lacs. Therefore, at low luminosities, there is a very large dispersion of SED shapes; iv) in low power BL Lacs, the contribution of the host galaxy is important. Remarkably, the luminosity distribution of the host galaxies of BL Lacs are spread in a very narrow range; v) a simple sum of two smoothly joining power laws can describe the blazar SEDs very well.


title: Inequalities for Quadratic Operator Perspective of Convex Functions and Bounded Linear Operators in Hilbert Spaces

abstract: In this paper we introduce the concept of quadratic operator perspective for a continuous function {\Phi} defined on the positive semi-axis of real numbers. This generalize the quadratic weighted operator geometric mean and the quadratic relative operator entropy. Some inequalities for this perspective of convex functions are established. Applications for quadratic weighted operator geometric mean and quadratic relative operator entropy are also provided.


title: CP-violating phase on magnetized toroidal orbifolds

abstract: We study the CP-violating phase of the quark sector in the $U(8)$ flavor model on $T^2/Z_N \, (N=2,3,4,6)$ with non-vanishing magnetic fluxes, where properties of possible origins of the CP violation are investigated minutely. In this system, a non-vanishing value is mandatory in the real part of the complex modulus parameter $\tau$ of the two-dimensional torus. On $T^2$ without orbifolding, underlying discrete flavor symmetries severely restrict the form of Yukawa couplings and it is very difficult to reproduce the observed pattern in the quark sector including the CP-violating phase $\delta_{\rm CP}$. In cases of multiple Higgs doublets emerging on $T^2/Z_2$, the mass matrices of the zero-mode fermions can be written in the Gaussian textures by choosing appropriate configurations of vacuum expectation values of the Higgs fields. When such Gaussian textures of mass matrices are realized, we show that all of the quark profiles, which are mass hierarchies among the quarks, quark mixing angles, and $\delta_{\rm CP}$ can be simultaneously realized.


title: Comprehending particle production at RHIC and LHC energies using global measurements

abstract: The centrality dependence of the charged-particle multiplicity densities ($dN_{ch}/d\eta$) and transverse energy densities ($dE_{T}/d\eta$) are investigated using the two-component Glauber approach for broad range of energies in heavy ion collisions at RHIC and LHC. A comprehensive study shows that the data is well described within the framework of two component model which includes the contribution of "soft processes" and "hard processes" for different centrality classes and energies. The data at two different energies are compared by means of the ratio of $dN_{ch}/d\eta$ (and $dE_{T}/d\eta$) to see the interplay of energy and $x$-scaling.


title: Spectral line intensity and polarization in gas-dusty medium

abstract: We study the standard spectral line radiative transfer equation for media consisting of resonant atoms and non-resonant components (the dust grains and the atoms without considered spectral transition). Our goal is to study the intensity and polarization of the resonance radiation emerging from semi-infinite atmosphere. Using the known technique of resolvent matrices, we obtain the exact solution of vectorial radiative transfer equation for various sources of non-polarized radiation in semi-infinite atmosphere. Homogeneous, linear increasing and exponentially decreasing sources are considered. Recall, that homogeneous sources correspond to the isothermal atmosphere and the exponential ones correspond to the incident radiation with different angles of inclination.


title: Expansion under climate change: the genetic consequences

abstract: Range expansion and range shifts are crucial population responses to climate change. Genetic consequences are not well understood but are clearly coupled to ecological dynamics that, in turn, are driven by shifting climate conditions. We model a population with a deterministic reaction-- diffusion model coupled to a heterogeneous environment that develops in time due to climate change. We decompose the resulting travelling wave solution into neutral genetic components to analyse the spatio-temporal dynamics of its genetic structure. Our analysis shows that range expansions and range shifts under slow climate change preserve genetic diversity. This is because slow climate change creates range boundaries that promote spatial mixing of genetic components. Mathematically , the mixing leads to so-called pushed travelling wave solutions. This mixing phenomenon is not seen in spatially homogeneous environments, where range expansion reduces genetic diversity through gene surfing arising from pulled travelling wave solutions. However, the preservation of diversity is diminished when climate change occurs too quickly. Using diversity indices, we show that fast expansions and range shifts erode genetic diversity more than slow range expansions and range shifts. Our study provides analytical insight into the dynamics of travelling wave solutions in heterogeneous environments.


title: Representations of $p$-convolution algebras on $L^q$-spaces

abstract: For a nontrivial locally compact group $G$, and $p\in [1,\infty)$, consider the Banach algebras of $p$-pseudofunctions, $p$-pseudomeasures, $p$-convolvers, and the full group $L^p$-operator algebra. We show that these Banach algebras are operator algebras if and only if $p=2$. More generally, we show that for $q\in [1,\infty)$, these Banach algebras can be represented on an $L^q$-space if and only if one of the following holds: (a) $p=2$ and $G$ is abelian; or (b) $|\frac 1p - \frac 12|=|\frac 1q - \frac 12|$. This result can be interpreted as follows: for $p,q\in [1,\infty)$, the $L^p$- and $L^q$-representation theories of a group are incomparable, except in the trivial cases when they are equivalent. As an application, we show that, for distinct $p,q\in [1,\infty)$, if the $L^p$ and $L^q$ crossed products of a topological dynamical system are isomorphic, then $\frac 1p + \frac 1q=1$. In order to prove this, we study the following relevant aspects of $L^p$-crossed products: existence of approximate identities, duality with respect to $p$, and existence of canonical isometric maps from group algebras into their multiplier algebras.


title: The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: higher-order correlations revealed by germ-grain Minkowski Functionals

abstract: We probe the higher-order clustering of the galaxies in the final data release (DR12) of the Sloan Digital Sky Survey Baryon Oscillation Spectroscopic Survey (BOSS) using the method of germ-grain Minkowski Functionals (MFs). Our sample consists of 410,615 BOSS galaxies from the northern Galactic cap in the redshift range 0.450--0.595. We show the MFs to be sensitive to contributions up to the six-point correlation function for this data set. We ensure with a custom angular mask that the results are more independent of boundary effects than in previous analyses of this type. We extract the higher-order part of the MFs and quantify the difference to the case without higher-order correlations. The resulting $\chi^{2}$ value of over 10,000 for a modest number of degrees of freedom, O(200), indicates a 100-sigma deviation and demonstrates that we have a highly significant signal of the non-Gaussian contributions to the galaxy distribution. This statistical power can be useful in testing models with differing higher-order correlations. Comparing the galaxy data to the QPM and MultiDark-Patchy mocks, we find that the latter better describes the observed structure. From an order-by-order decomposition we expect that, for example, already a reduction of the amplitude of the MD-Patchy mock power spectrum by 5% would remove the remaining tension.


title: Shaping of the inner Oort cloud by Planet Nine

abstract: We present a numerical simulation of the dynamical interaction between the proposed Planet Nine and a debris disk around the Sun for $4{\rm Gyr}$, accounting for the secular perturbation of the four giant planets in two scenarios: (a) an initially thin circular disk around the Sun (b) inclined and eccentric disk. We show, in both scenarios, that Planet Nine governs the dynamics in between $1000{\rm -5000{\rm AU}}$ and forms spherical structure in the inner part ($\sim1000{\rm AU}$) and inclined disk. This structure is the outcome of mean motion resonances and secular interaction with Planet Nine. We compare the morphology of this structure with the outcome from a fly-by encounter of a star with the debris disk and show distinct differences between the two cases. We predict that this structure serves as a source of comets and calculate the resulting comet production rate to be detectable.


title: The Infrared-Gamma-Ray Connection: A WISE View of the Extragalactic Gamma-Ray Sky

abstract: Using data from the WISE all-sky survey we discovered that the non-thermal infrared (IR) emission of blazars, the largest known population of extragalactic gamma-ray sources, has peculiar spectral properties. In this work, we confirm and strengthen our previous analyses using the latest available releases of both the WISE and the Fermi source catalogs. We also show that there is a tight correlation between the mid-IR colors and the gamma-ray spectral index of Fermi blazars. We name this correlation "the infrared--gamma-ray connection". We discuss how this connection links both the emitted powers and the spectral shapes of particles accelerated in jets arising from blazars over ten decades in energy. Based on this evidence, we argue that the infrared--gamma-ray connection is stronger than the well known radio--gamma-ray connection.


title: Composite fermion duality for half-filled multicomponent Landau Levels

abstract: We study the interplay of particle-hole symmetry and fermion-vortex duality in multicomponent half-filled Landau levels, such as quantum Hall gallium arsenide bilayers and graphene. For the $\nu{=}1/2{+}1/2$ bilayer, we show that particle-hole-symmetric interlayer Cooper pairing of composite fermions leads to precisely the same phase as the electron exciton condensate realized in experiments. This equivalence is easily understood by applying the recent Dirac fermion formulation of $\nu{=}1/2$ to two components. It can also be described by Halperin-Lee-Read composite fermions undergoing interlayer $p_x{+}ip_y$ pairing. An RG analysis showing strong instability to interlayer pairing at large separation $d\rightarrow \infty$ demonstrates that two initially-decoupled composite Fermi liquids can be smoothly tuned into the conventional bilayer exciton condensate without encountering a phase transition. We also discuss multicomponent systems relevant to graphene, derive related phases including a $Z_2$ gauge theory with spin-half visons, and argue for symmetry-enforced gaplessness under full SU$(N_f)$ flavor symmetry when the number of components $N_f$ is even.


title: Theory of quantum kagome ice and vison zero modes

abstract: We derive an effective $Z_2$ gauge theory to describe the quantum kagome ice (QKI) state that has been observed by Carrasquilla $\textit{et. al.}$ in Monte Carlo studies of the $S = 1/2$ kagome XYZ model in a Zeeman field. The numerical results on QKI are consistent with, but do not confirm or rule out, the hypothesis that it is a $Z_2$ spin liquid. Our effective theory allows us to explore this hypothesis and make a striking prediction for future numerical studies, namely that symmetry-protected vison zero modes arise at lattice disclination defects, leading to a Curie defect term in the spin susceptibility, and a characteristic $(N_{dis} - 1) \ln 2$ contribution to the entropy, where $N_{dis}$ is the number of disclinations. Only the $Z_2$ Ising symmetry is required to protect the vison zero modes. This is remarkable because a unitary $Z_2$ symmetry cannot be responsible for symmetry-protected degeneracies of local degrees of freedom. We also discuss other signatures of symmetry fractionalization in the $Z_2$ spin liquid, and phase transitions out of the $Z_2$ spin liquid to nearby ordered phases.


title: Realizing topological surface states in a lower-dimensional flat band

abstract: The anomalous surface states of symmetry protected topological (SPT) phases are usually thought to be only possible in conjunction with the higher dimensional topological bulk. However, it has recently been realized that a class of anomalous SPT surface states can be realized in the same dimension if symmetries are allowed to act in a nonlocal fashion. An example is the particle-hole symmetric half filled Landau level, which effectively realizes the anomalous surface state of a 3D chiral Topological Insulator (class AIII). A dual description in terms of Dirac composite fermions has also been discussed. Here we explore generalizations of these constructions to multicomponent quantum Hall states. Our results include a duality mapping of the bilayer case to composite bosons with Kramers degeneracy and the possibility of a particle hole symmetric integer quantum Hall state when the number of components is a multiple of eight. Next, we make a further extension by half filling other classes of topological bands and imposing particle hole symmetry. When applied to time-reversal invariant topological insulators we realize a different chiral class (CII) topological surface state. Notably, half-filling a 3D TI band allows for the realization of the surface of the otherwise inaccessible 4D topological insulator, which supports an anomalous 3D Dirac cone. Surface topological orders equivalent to the 3D Dirac cone (from the global anomaly standpoint) are constructed and connections to Witten's SU(2) anomaly are made. These observations may also be useful for numerical simulations of topological surface states and of Dirac fermions without fermion doubling.


title: Superfluid helium II as the QCD vacuum

abstract: We study the winding number susceptibility < I^2> in a superfluid system and the topological susceptibility < Q^2> in QCD. We argue that both correlation functions exhibit similar structures, including the generation of the contact terms. We discuss the nature of the contact term in superfluid system and argue that it has exactly the same origin as in QCD, and it is related to the long distance physics which cannot be associated with conventional microscopical degrees of freedom such as phonons and rotons. We emphasize that the conceptual similarities between superfluid system and QCD may lead, hopefully, to a deeper understanding of the topological features of a superfluid system as well as the QCD vacuum.

categories: cond-mat.stat-mech hep-lat hep-ph hep-th


title: Quantum critical properties of a metallic spin density wave transition

abstract: We report on numerically exact determinantal quantum Monte Carlo simulations of the onset of spin-density wave (SDW) order in itinerant electron systems captured by a sign-problem-free two-dimensional lattice model. Extensive measurements of the SDW correlations in the vicinity of the phase transition reveal that the critical dynamics of the bosonic order parameter are well described by a dynamical critical exponent z = 2, consistent with Hertz-Millis theory, but are found to follow a finite-temperature dependence that does not fit the predicted behavior of the same theory. The presence of critical SDW fluctuations is found to have a strong impact on the fermionic quasiparticles, giving rise to a dome-shaped superconducting phase near the quantum critical point. In the superconducting state we find a gap function that has an opposite sign between the two bands of the model and is nearly constant along the Fermi surface of each band. Above the superconducting $T_c$ our numerical simulations reveal a nearly temperature and frequency independent self energy causing a strong suppression of the low-energy quasiparticle spectral weight in the vicinity of the hot spots on the Fermi surface. This indicates a clear breakdown of Fermi liquid theory around these points.


title: PKDGRAV3: Beyond Trillion Particle Cosmological Simulations for the Next Era of Galaxy Surveys

abstract: We report on the successful completion of a 2 trillion particle cosmological simulation to z=0 run on the Piz Daint supercomputer (CSCS, Switzerland), using 4000+ GPU nodes for a little less than 80h of wall-clock time or 350,000 node hours. Using multiple benchmarks and performance measurements on the US Oak Ridge National Laboratory Titan supercomputer, we demonstrate that our code PKDGRAV3, delivers, to our knowledge, the fastest time-to-solution for large-scale cosmological N-body simulations. This was made possible by using the Fast Multipole Method in conjunction with individual and adaptive particle time steps, both deployed efficiently (and for the first time) on supercomputers with GPU-accelerated nodes. The very low memory footprint of PKDGRAV3 allowed us to run the first ever benchmark with 8 trillion particles on Titan, and to achieve perfect scaling up to 18000 nodes and a peak performance of 10 Pflops.


title: The dust content of galaxies from z = 0 to z = 9

abstract: We study the dust content of galaxies from $z=0$ to $z=9$ in semi-analytic models of galaxy formation that include new recipes to track the production and destruction of dust. We include condensation of dust in stellar ejecta, the growth of dust in the interstellar medium (ISM), the destruction of dust by supernovae and in the hot halo, and dusty winds and inflows. The rate of dust growth in the ISM depends on the metallicity and density of molecular clouds. Our fiducial model reproduces the relation between dust mass and stellar mass from $z=0$ to $z=7$, the number density of galaxies with dust masses less than $10^{8.3}\,\rm{M}_\odot$, and the cosmic density of dust at $z=0$. The model accounts for the double power-law trend between dust-to-gas (DTG) ratio and gas-phase metallicity of local galaxies and the relation between DTG ratio and stellar mass. The dominant mode of dust formation is dust growth in the ISM, except for galaxies with $M_*<10^7\,\rm{M}_\odot$, where condensation of dust in supernova ejecta dominates. The dust-to-metal ratio of galaxies depends on the gas-phase metallicity, unlike what is typically assumed in cosmological simulations. Model variants including higher condensation efficiencies, a fixed timescale for dust growth in the ISM, or no growth at all reproduce some of the observed constraints, but fail to simultaneously reproduce the shape of dust scaling relations and the dust mass of high-redshift galaxies.


title: A framework for testing leptonic unitarity by neutrino oscillation experiments

abstract: If leptonic unitarity is violated by new physics at an energy scale much lower than the electroweak scale, which we call low-scale unitarity violation, it has different characteristic features from those expected in unitarity violation at high-energy scales. They include maintaining flavor universality and absence of zero-distance flavor transition. We present a framework for testing such unitarity violation at low energies by neutrino oscillation experiments. Starting from the unitary 3 active plus $N$ (arbitrary integer) sterile neutrino model we show that by restricting the active-sterile and sterile-sterile neutrino mass squared differences to $\gtrsim$ 0.1 eV$^2$ the oscillation probability in the $(3+N)$ model becomes insensitive to details of the sterile sector, providing a nearly model-independent framework for testing low-scale unitarity violation. Yet, the presence of the sterile sector leaves trace as a constant probability leaking term, which distinguishes low-scale unitarity violation from the high-scale one. The non-unitary mixing matrix in the active neutrino subspace is common for the both cases. We analyze how severely the unitarity violation can be constrained in $\nu_{e}$-row by taking a JUNO-like setting to simulate medium baseline reactor experiments. Possible modification of the features of the $(3+N)$ model due to matter effect is discussed to first order in the matter potential.


title: SED-dependent Galactic Extinction Prescription for Euclid and Future Cosmological Surveys

abstract: The outcome of upcoming cosmological surveys will depend on the accurate estimates of photometric redshifts. In the framework of the implementation of the photo-z algorithm for Euclid, we are exploring new avenues to improve template-fitting methods. The paper focusses on the prescription of the extinction of source light by dust in the Milky Way. Since Galactic extinction strongly correlates with wavelength and photometry is commonly obtained in broad-band filters, the amount of absorption depends on the source SED, a point often neglected as the SED is not known a-priori. A consequence of this is that the observed E(B-V) (=A_B-A_V) will be different from the E(B-V) used to normalise the absorption law k_lambda (=A_lambda/E(B-V)). Band-pass corrections are required to renormalise the law for a given SED. We assess the band-pass corrections of a range of SEDs and find they vary by up to 20%. We investigate how dust-to-reddening scaling factors depend of the sources used for their calibration. We derive scaling factors from the color excesses of z<0.4 SDSS red galaxies and show that band-pass corrections predict the observed differences. Extinction is then estimated for a range of SEDs and filters relevant to Euclid and other cosmological ground-based surveys. For high extinction line-of-sights (E(B-V)>0.1, ~8% of the Euclid survey), the variations in corrections can be ~0.1mag in the `bluer' optical filters and ~0.04mag in the NIR filters. An inaccurate correction of extinction critically affects photo-z. In particular, for high extinctions and z<0.5, the bias (mean D_z=z_phot-z_real) exceeds 0.2%(1+z), the precision required by weak-lensing analyses. Additional uncertainty on the MW extinction law further reduces the photo-z precision. We propose a new prescription of Galactic absorption for template-fitting algorithms that takes into consideration the dependence of extinction with SED.


title: BAT AGN Spectroscopic Survey-III. An observed link between AGN Eddington ratio and narrow emission line ratios

abstract: We investigate the observed relationship between black hole mass ($M_{\rm BH}$), bolometric luminosity ($L_{\rm bol}$), and Eddington ratio (${\lambda}_{\rm Edd}$) with optical emission line ratios ([NII] {\lambda}6583/H{\alpha}, [SII] {\lambda}{\lambda}6716,6731/H{\alpha}, [OI] {\lambda}6300/H{\alpha}, [OIII] {\lambda}5007/H{\beta}, [NeIII] {\lambda}3869/H{\beta}, and HeII {\lambda}4686/H{\beta}) of hard X-ray-selected AGN from the BAT AGN Spectroscopic Survey (BASS). We show that the [NII] {\lambda}6583/H{\alpha} ratio exhibits a significant correlation with ${\lambda}_{\rm Edd}$ ($R_{\rm Pear}$ = -0.44, $p$-value=$3\times10^{-13}$, {\sigma} = 0.28 dex), and the correlation is not solely driven by $M_{\rm BH}$ or $L_{\rm bol}$. The observed correlation between [NII] {\lambda}6583/H{\alpha} ratio and $M_{\rm BH}$ is stronger than the correlation with $L_{\rm bol}$, but both are weaker than the ${\lambda}_{\rm Edd}$ correlation. This implies that the large-scale narrow lines of AGN host galaxies carry information about the accretion state of the AGN central engine. We propose that the [NII] {\lambda}6583/H{\alpha} is a useful indicator of Eddington ratio with 0.6 dex of rms scatter, and that it can be used to measure ${\lambda}_{\rm Edd}$ and thus $M_{\rm BH}$ from the measured $L_{\rm bol}$, even for high redshift obscured AGN. We briefly discuss possible physical mechanisms behind this correlation, such as the mass-metallicity relation, X-ray heating, and radiatively driven outflows.


title: A Testable Conspiracy: Simulating Baryonic Effects on Self-Interacting Dark Matter Halos

abstract: We investigate the response of self-interacting dark matter (SIDM) halos to the growth of galaxy potentials using idealized simulations, each run in tandem with standard collisionless Cold Dark Matter (CDM). We find a greater diversity in the SIDM halo profiles compared to the CDM halo profiles. If the stellar gravitational potential strongly dominates in the central parts of a galaxy, then SIDM halos can be as dense as CDM halos on observable scales. For extreme cases with highly compact disks core collapse can occur, leading to SIDM halos that are denser and cuspier than their CDM counterparts. If the stellar potential is not dominant, then SIDM halos retain constant density cores with densities far below CDM predictions. When a disk potential is present, the inner SIDM halo becomes \em{more flattened} in the disk plane than the CDM halo. These results are in excellent quantitative agreement with the predictions of Kaplinghat et al. (2014). We also simulated a galaxy cluster halo with a central stellar distribution similar to the brightest central galaxy of the cluster A2667. A SIDM halo simulated with cross section over mass $\sigma/m = 0.1\ \mathrm{cm^2 g^{-1}}$ provides a good match to the measured dark matter density profile of A2667, while an adiabatically-contracted CDM halo is denser and cuspier. The cored profile of the same halo simulated with $\sigma/m = 0.5\ \mathrm{cm^2 g^{-1}}$ is not dense enough to match A2667. Our findings are in agreement with previous results that $\sigma/m \gtrsim 0.1\ \mathrm{cm^2 g^{-1}}$ is disfavored for dark matter collision velocities in excess of about 1500 km/s. More generally, the predictive cross-talk between baryonic potentials and SIDM density distributions offers new directions for constraining SIDM cross sections in massive galaxies where baryons are dynamically important.


title: Positivity, Grassmannian Geometry and Simplex-like Structures of Scattering Amplitudes

abstract: This article revisits and elaborates the significant role of positive geometry of momentum twistor Grassmannian for planar N=4 SYM scattering amplitudes. First we establish the fundamentals of positive Grassmannian geometry for tree amplitudes. Then we formulate this subject around these four major facets: 1. Deriving the tree and 1-loop BCFW recursion relations solely from positivity, after introducing the simple building blocks called positive components. 2. Applying Grassmannian geometry and Pluecker coordinates to determine the signs in N$^2$MHV homological identities. Most of them in fact reflect the secret incarnation of the simple 6-term NMHV identity. 3. Deriving the stacking positivity relation, which is powerful for parameterizing the matrix representative in terms of positive variables, of a given Grassmannian geometric configuration. 4. Introducing a highly refined formalism for tree BCFW recursion relation, which reveals its two-fold simplex-like structures. First, the BCFW contour in terms of (reduced) Grassmannian geometry representatives is dissected into a triangle-shape sum, as only a small fraction of the sum needs to be explicitly identified. Second, this fraction can be further dissected, according to different growing modes and parameters. The growing modes possess the shapes of solid simplices of various dimensions, with which infinite number of BCFW cells can be entirely captured by the characteristic objects called fully-spanning cells. We find that for a given k, beyond n=4k+1 there is no more new fully-spanning cell, which signifies the essential termination of the recursive growth of BCFW cells. As n increases beyond the termination point, the BCFW contour simply replicates itself according to the simplex-like patterns, which enables us to master all BCFW cells once for all without actually identifying most of them.


title: Coarse-grained hidden entropy production in partially inaccessible quantum jump trajectories

abstract: We consider an open quantum system for which only a subset of all possible transitions are accessible, while the remaining ones are hidden from direct observation. Using a modification of the notion of quantum jump trajectories we introduce the coarse-grained hidden entropy, which quantifies the entropy production in the hidden subsystem conditioned on our observations of the visible part. The entropy production consisting of the sum of visible and coarse-grained hidden entropy is shown to satisfy an integral fluctuations theorem. Depending on the nature of the continuously occurring measurement and feedback processes between the interacting subsystems, the visible entropy can take negative values in which case the hidden systems acts as a Maxwell's demon. Our results allow us to study quantum-mechanical effects in autonomous systems, such as the autonomous quantum Maxwell's demon we introduce as an illustrative example.


title: The impact of protocluster environments at z = 1.6

abstract: We investigate the effects of dense environments on galaxy evolution by examining how the properties of galaxies in the z = 1.6 protocluster Cl 0218.3-0510 depend on their location. We determine galaxy properties using spectral energy distribution fitting to 14-band photometry, including data at three wavelengths that tightly bracket the Balmer and 4000A breaks of the protocluster galaxies. We find that two-thirds of the protocluster galaxies, which lie between several compact groups, are indistinguishable from field galaxies. The other third, which reside within the groups, differ significantly from the intergroup galaxies in both colour and specific star formation rate. We find that the fraction of red galaxies within the massive protocluster groups is twice that of the intergroup region. These excess red galaxies are due to enhanced fractions of both passive galaxies (1.7 times that of the intergroup region) and dusty star-forming galaxies (3 times that of the intergroup region). We infer that some protocluster galaxies are processed in the groups before the cluster collapses. These processes act to suppress star formation and change the mode of star formation from unobscured to obscured.


title: Reconstructing the three-dimensional local dark matter velocity distribution

abstract: Directionally sensitive dark matter (DM) direct detection experiments present the only way to observe the full three-dimensional velocity distribution of the Milky Way halo local to Earth. In this work we compare methods for extracting information about the local DM velocity distribution from a set of recoil directions and energies in a range of hypothetical directional and non-directional experiments. We compare a model independent empirical parameterisation of the velocity distribution based on an angular discretisation with a model dependent approach which assumes knowledge of the functional form of the distribution. The methods are tested under three distinct halo models which cover a range of possible phase space structures for the local velocity distribution: a smooth Maxwellian halo, a tidal stream and a debris flow. In each case we use simulated directional data to attempt to reconstruct the shape and parameters describing each model as well as the DM particle properties. We find that the empirical parametrisation is able to make accurate unbiased reconstructions of the DM mass and cross section as well as capture features in the underlying velocity distribution in certain directions without any assumptions about its true functional form. We also find that by extracting directionally averaged velocity parameters with this method one can discriminate between halo models with different classes of substructure.


title: Bilocal theory and gravity I

abstract: With the aim of investigating the relation between gravity and non-locality at the classical level, we study a bilocal scalar field model. Bilocality introduces new (internal) degrees of freedom that seem to reproduce gravity. We show that the equations of motion of the massless branch of the free bilocal model match those of linearized gravity. We also discuss higher orders of perturbation theory, where there is self-interaction in both gravity and the bilocal field sectors.


title: A marked correlation function for constraining modified gravity models

abstract: Future large scale structure surveys will provide increasingly tight constraints on our cosmological model. These surveys will report results on the distance scale and growth rate of perturbations through measurements of Baryon Acoustic Oscillations and Redshift-Space Distortions. It is interesting to ask: what further analyses should become routine, so as to test as-yet-unknown models of cosmic acceleration? Models which aim to explain the accelerated expansion rate of the Universe by modifications to General Relativity often invoke screening mechanisms which can imprint a non-standard density dependence on their predictions. This suggests density-dependent clustering as a `generic' constraint. This paper argues that a density-marked correlation function provides a density-dependent statistic which is easy to compute and report and requires minimal additional infrastructure beyond what is routinely available to such survey analyses. We give one realization of this idea and study it using low order perturbation theory. We encourage groups developing modified gravity theories to see whether such statistics provide discriminatory power for their models.


title: Exoplanet Orbital Eccentricities Derived From LAMOST-Kepler Analysis

abstract: The nearly circular (mean eccentricity <e>~0.06) and coplanar (mean mutual inclination <i>~3 deg) orbits of the Solar System planets motivated Kant and Laplace to put forth the hypothesis that planets are formed in disks, which has developed into the widely accepted theory of planet formation. Surprisingly, the first several hundred extrasolar planets (mostly Jovian) discovered using the Radial Velocity (RV) technique are commonly on eccentric orbits (<e> ~ 0.3). This raises a fundamental question: Are the Solar System and its formation special? The Kepler mission has found thousands of transiting planets dominated by sub-Neptunes, but most of their orbital eccentricities remain unknown. By using the precise spectroscopic host star parameters from the LAMOST observations, we measure the eccentricity distributions for a large (698) and homogeneous Kepler planet sample with transit duration statistics. Nearly half of the planets are in systems with single transiting planets (singles), while the other half are multiple-transiting planets (multiples). We find an eccentricity dichotomy: on average, Kepler singles are on eccentric orbits with <e>~0.3, while the multiples are on nearly circular (<e> = 0.04^{+0.03}_{-0.04}) and coplanar (<i> = 1.4^{+0.8}_{-1.1} deg) orbits similar to the Solar System planets. Our results are consistent with previous studies of smaller samples and individual systems. We also show that Kepler multiples and solar system objects follow a common relation <e>~(1-2)x<i> between mean eccentricities and mutual inclinations. The prevalence of circular orbits and the common relation may imply that the solar system is not so atypical in the galaxy after all.


title: Global Instability of Exo-Moon System Triggered by Photo-Evaporation

abstract: Many exoplanets have been found in orbits close to their host stars and thus they are subject to the effects of photo-evaporation. Previous studies have shown that a large portion of exoplanets detected by the Kepler mission have been significantly eroded by photo-evaporation. In this paper, we numerically study the effects of photo-evaporation on the orbital evolution of a hypothesized moon system around a planet. We find that photo-evaporation is crucial to the stability of the moon system. Photo-evaporation can erode the atmosphere of the planet thus leading to significant mass loss. As the planet loses mass, its Hill radius shrinks and its moons increase their orbital semi-major axes and eccentricities. When some moons approach their critical semi-major axes, global instability of the moon system would be triggered, which usually ends up with two, one or even zero surviving moons. Some lost moons could escape from the moon system to become a new planet orbiting the star or run away further to become a free-floating object in the Galaxy. Given the destructive role of photo-evaporation, we speculate that exo-moons are less common for close-in planets ($<0.1$ AU), especially those around M-type stars, because they are more X-ray luminous and thus enhancing photo-evaporation. The lessons we learn in this study may be helpful for the target selection of on-going/future exomoon searching programs.


title: The Fate of Debris in the Pluto-Charon System

abstract: The Pluto-Charon system has come into sharper focus following the fly by of New Horizons. We use N-body simulations to probe the unique dynamical history of this binary dwarf planet system. We follow the evolution of the debris disc that might have formed during the Charon-forming giant impact. First, we note that in-situ formation of the four circumbinary moons is extremely difficult if Charon undergoes eccentric tidal evolution. We track collisions of disc debris with Charon, estimating that hundreds to hundreds of thousands of visible craters might arise from 0.3--5 km radius bodies. New Horizons data suggesting a dearth of these small craters may place constraints on the disc properties. While tidal heating will erase some of the cratering history, both tidal and radiogenic heating may also make it possible to differentiate disc debris craters from Kuiper belt object craters. We also track the debris ejected from the Pluto-Charon system into the Solar System; while most of this debris is ultimately lost from the Solar System, a few tens of 10--30 km radius bodies could survive as a Pluto-Charon collisional family. Most are plutinos in the 3:2 resonance with Neptune, while a small number populate nearby resonances. We show that migration of the giant planets early in the Solar System's history would not destroy this collisional family. Finally, we suggest that identification of such a family would likely need to be based on composition as they show minimal clustering in relevant orbital parameters.


title: Asymmetric Orbital Distribution near Mean Motion Resonance: Application to Planets Observed by Kepler and Radial Velocities

abstract: Many multiple-planet systems have been found by the Kepler transit survey and various radial velocity (RV) surveys. Kepler planets show an asymmetric feature, namely, there are small but significant deficits/excesses of planet pairs with orbital period spacing slightly narrow/wide of the exact resonance, particularly near the first order mean motion resonance (MMR), such as 2:1 and 3:2 MMR. Similarly, if not exactly the same, an asymmetric feature (pileup wide of 2:1 MMR) is also seen in RV planets, but only for massive ones. We analytically and numerically study planets' orbital evolutions near and in the MMR. We find that their orbital period ratios could be asymmetrically distributed around the MMR center regardless of dissipation. In the case of no dissipation, Kepler planets' asymmetric orbital distribution could be partly reproduced for 3:2 MMR but not for 2:1 MMR, implying that dissipation might be more important to the latter. The pileup of massive RV planets just wide of 2:1 MMR is found to be consistent with the scenario that planets formed separately then migrated toward the MMR. The location of the pileup infers a K value of 1-100 on the order of magnitude for massive planets, where K is the damping rate ratio between orbital eccentricity and semimajor axis during planet migration.


title: Non-Unitarity, sterile neutrinos, and Non-Standard neutrino Interactions

abstract: The simplest Standard Model extension to explain neutrino masses involves the addition of right-handed neutrinos. At some level, this extension will impact neutrino oscillation searches. In this work we explore the differences and similarities between the case in which these neutrinos are kinematically accessible (sterile neutrinos) or not (mixing matrix non-unitarity). We clarify apparent inconsistencies in the present literature when using different parametrizations to describe these effects and recast both limits in the popular neutrino non-standard interaction (NSI) formalism. We find that, in the limit in which sterile oscillations are averaged out at the near detector, their effects at the far detector coincide with non-unitarity at leading order, even in presence of a matter potential. We also summarize the present bounds existing in both limits and compare them with the expected sensitivities of near future facilities taking the DUNE proposal as a benchmark. We conclude that non-unitarity effects are too constrained to impact present or near future neutrino oscillation facilities but that sterile neutrinos can play an important role at long baseline experiments. The role of the near detector is also discussed in detail.


title: Ultraviolet and X-ray variability of active galactic nuclei with Swift

abstract: We analyse a sample of 21 active galactic nuclei (AGN) using data from the Swift satellite to study the variability properties of the population in the X-ray, UV and optical band. We find that the variable part of the UV-optical emission has a spectrum consistent with a powerlaw, with an average index of $-2.21\pm0.13$, as would be expected from central illumination of a thin disc (index of -7/3). We also calculate the slope of a powerlaw from UV to X-ray variable emission, $\alpha_{\rm OX,Var}$; the average for this sample is $\alpha_{\rm OX,Var} = -1.06 \pm 0.04$. The anticorrelation of $\alpha_{\rm OX}$ with the UV luminosity, $L_{\rm UV}$, previously found in the average emission is also present in the variable part: $\alpha_{\rm OX,Var} = (-0.177 \pm 0.083) {\rm log} (L_{\rm \nu,Var} (2500\,\AA)) + (3.88 \pm 2.33)$. Correlated variability between the emission in X-rays and UV is detected significantly for 9 of the 21 sources. All these cases are consistent with the UV lagging the X-rays, as would be seen if the correlated UV variations were produced by the reprocessing of X-ray emission. The observed UV lags are tentatively longer than expected for a standard thin disc.


title: Phase diagram and correlation functions of the anisotropic imperfect Bose gas in $d$ dimensions

abstract: We study an anisotropic variant of the $d$-dimensional imperfect Bose gas, where the asymptotic behaviour of the dispersion $\epsilon_{\bf k}$ at vanishing momentum $\bf{k}$ may differ from the standard quadratic form. The analysis reveals the key role of the shift exponent $\psi$ governing the asymptotic behaviour of the critical temperature $T_c(\mu)$ as a function of the chemical potential $\mu$ at $T_c\to 0$. We argue that the universality classes of Bose-Einstein condensation admitted by the model may be classified according to the allowed values of $\psi$ so that spatial dimensionality has only an indirect impact on the transition properties. We analyse the correlation function of the model and discuss its asymptotics depending on the direction. Both for the perfect and imperfect anisotropic Bose gases, the correlation function $\chi({\bf x})$ at $T>T_c$ turns out to show either exponential decay or exponentially damped oscillatory behaviour depending on the orientation of ${\bf x}$ with respect to the dispersion anisotropies.

categories: cond-mat.stat-mech cond-mat.quant-gas math-ph math.MP


title: Relation between proton decay and PMNS phase in the minimal SUSY $SO(10)$ GUT

abstract: Proton decay is one of the most important predictions of the grand unified theory (GUT). In the supersymmetric (SUSY) GUT, proton decays via the dimension-five operators need to be suppressed. In the $SO(10)$ model where ${\bf 10}+\overline{\bf 126}$ Higgs fields couple to fermions, neutrino oscillation parameters including the CP-violating Pontecorvo-Maki-Nakagawa-Sakata (PMNS) phase can be related to the Yukawa couplings to generate the dimension-five operators in the unified framework. We show how the suppressed proton decay depends on the PMNS phase, and stress the importance of the precise measurements of the PMNS phase as well as the neutrino 23-mixing angle. These become especially important if the SUSY particles are found around less than a few TeV at LHC and proton decays are observed at Hyper-Kamiokande and DUNE experiments in the near future.


title: Graphical modelling of multivariate spatial point processes with continuous marks

abstract: This paper is the second in a series of papers which combine graphical modelling and marked spatial point patterns. Extending the previous results of \cite Eckardt (2016a), we introduce a marked spatial dependence graph model which depicts the global dependence structure of quantitatively marked multi-type points that occur in space based on the marked conditional partial spectral coherence. Most beneficial, no structural assumption with respect to the characteristics in the data are to be made prior to analysis. This approach presents a computationally efficient method of pattern recognition in highly structured and high dimensional multi-type spatial point processes where also quantitative marks are available. Unlike all previous methods, our new model permits the simultaneous analysis of all multivariate conditional interrelations. The new technique is illustrated analysing the diameter at breast hight of 37 different tree species recorded at 10053 locations in Duke Forest.


title: Benchmarking the Graphulo Processing Framework

abstract: Graph algorithms have wide applicablity to a variety of domains and are often used on massive datasets. Recent standardization efforts such as the GraphBLAS specify a set of key computational kernels that hardware and software developers can adhere to. Graphulo is a processing framework that enables GraphBLAS kernels in the Apache Accumulo database. In our previous work, we have demonstrated a core Graphulo operation called \textit{TableMult} that performs large-scale multiplication operations of database tables. In this article, we present the results of scaling the Graphulo engine to larger problems and scalablity when a greater number of resources is used. Specifically, we present two experiments that demonstrate Graphulo scaling performance is linear with the number of available resources. The first experiment demonstrates cluster processing rates through Graphulo's TableMult operator on two large graphs, scaled between $2^{17}$ and $2^{19}$ vertices. The second experiment uses TableMult to extract a random set of rows from a large graph ($2^{19}$ nodes) to simulate a cued graph analytic. These benchmarking results are of relevance to Graphulo users who wish to apply Graphulo to their graph problems.


title: The Herschel Exploitation of Local Galaxy Andromeda (HELGA) VII: A SKIRT radiative transfer model and insights on dust heating

abstract: The radiation of stars heats dust grains in the diffuse interstellar medium and in star-forming regions in galaxies. Modelling this interaction provides information on dust in galaxies, a vital ingredient for their evolution. It is not straightforward to identify the stellar populations heating the dust, and to link attenuation to emission on a sub-galactic scale. Radiative transfer models are able to simulate this dust-starlight interaction in a realistic, three-dimensional setting. We investigate the dust heating mechanisms on a local and global galactic scale, using the Andromeda galaxy (M31) as our laboratory. We perform a series of panchromatic radiative transfer simulations of Andromeda with our code SKIRT. The high inclination angle of M31 complicates the 3D modelling and causes projection effects. However, the observed morphology and flux density are reproduced fairly well from UV to sub-millimeter wavelengths. Our model reveals a realistic attenuation curve, compatible with previous, observational estimates. We find that the dust in M31 is mainly (91 % of the absorbed luminosity) heated by the evolved stellar populations. The bright bulge produces a strong radiation field and induces non-local heating up to the main star-forming ring at 10 kpc. The relative contribution of unevolved stellar populations to the dust heating varies strongly with wavelength and with galactocentric distance.The dust heating fraction of unevolved stellar populations correlates strongly with NUV-r colour and specific star formation rate. These two related parameters are promising probes for the dust heating sources at a local scale.


title: Scalar field dark matter in clusters of galaxies

abstract: One alternative to the cold dark matter (CDM) paradigm is the scalar field dark matter (SFDM) model, which assumes dark matter is a spin-0 ultra-light scalar field (SF) with a typical mass $m\sim10^{-22}\mathrm{eV}/c^2$ and positive self-interactions. Due to the ultra-light boson mass, the SFDM could form Bose-Einstein condensates (BEC) in the very early Universe, which are interpreted as the dark matter haloes. Although cosmologically the model behaves as CDM, they differ at small scales: SFDM naturally predicts fewer satellite haloes, cores in dwarf galaxies and the formation of massive galaxies at high redshifts. The ground state (or BEC) solution at zero temperature suffices to describe low-mass galaxies but fails for larger systems. A possible solution is adding finite-temperature corrections to the SF potential which allows combinations of excited states. In this work, we test the finite-temperature multistate SFDM solution at galaxy cluster scales and compare our results with the Navarro-Frenk-White (NFW) and BEC profiles. We achieve this by fitting the mass distribution of 13 \textit{Chandra} X-ray clusters of galaxies, excluding the region of the brightest cluster galaxy. We show that the SFDM model accurately describes the clusters' DM mass distributions offering an equivalent or better agreement than the NFW profile. The complete disagreement of the BEC model with the data is also shown. We conclude that the theoretically motivated multistate SFDM profile is an interesting alternative to empirical profiles and ad hoc fitting-functions that attempt to couple the asymptotic NFW decline with the inner core in SFDM.


title: Colouring squares of claw-free graphs

abstract: Is there some absolute $\varepsilon > 0$ such that for any claw-free graph $G$, the chromatic number of the square of $G$ satisfies $\chi(G^2) \le (2-\varepsilon) \omega(G)^2$, where $\omega(G)$ is the clique number of $G$? Erd\H{o}s and Ne\v{s}et\v{r}il asked this question for the specific case of $G$ the line graph of a simple graph and this was answered in the affirmative by Molloy and Reed. We show that the answer to the more general question is also yes, and moreover that it essentially reduces to the original question of Erd\H{o}s and Ne\v{s}et\v{r}il.


title: Squared chromatic number without claws or large cliques

abstract: Let $G$ be a claw-free graph on $n$ vertices with clique number $\omega$, and consider the chromatic number $\chi(G^2)$ of the square $G^2$ of $G$. Writing $\chi'_s(d)$ for the supremum of $\chi(L^2)$ over the line graphs $L$ of simple graphs of maximum degree at most $d$, we prove that $\chi(G^2)\le \chi'_s(\omega)$ for $\omega \in \{3,4\}$. For $\omega=3$, this implies the sharp bound $\chi(G^2) \leq 10$. For $\omega=4$, this implies $\chi(G^2)\leq 22$, which is within $2$ of the conjectured best bound. This work is motivated by a strengthened form of a conjecture of Erd\H{o}s and Ne\v{s}et\v{r}il.


title: Voigt transformations in retrospect: missed opportunities?

abstract: The teaching of modern physics often uses the history of physics as a didactic tool. However, as in this process the history of physics is not something studied but used, there is a danger that the history itself will be distorted in, as Butterfield calls it, a "Whiggish" way, when the present becomes the measure of the past. It is not surprising that reading today a paper written more than a hundred years ago, we can extract much more of it than was actually thought or dreamed by the author himself. We demonstrate this Whiggish approach on the example of Woldemar Voigt's 1887 paper. From the modern perspective, it may appear that this paper opens a way to both the special relativity and to its anisotropic Finslerian generalization which came into the focus only recently, in relation with the Cohen and Glashow's very special relativity proposal. With a little imagination, one can connect Voigt's paper to the notorious Einstein-Poincare priority dispute, which we believe is a Whiggish late time artifact. We use the related historical circumstances to give a broader view on special relativity, than it is usually anticipated.


title: Negative curves on symmetric blowups of the projective plane, resurgences and Waldschmidt constants

abstract: The Klein and Wiman configurations are highly symmetric configurations of lines in the projective plane arising from complex reflection groups. One noteworthy property of these configurations is that all the singularities of the configuration have multiplicity at least three. In this paper we study the surface X obtained by blowing up the projective plane in the singular points of one of these line configurations. We study invariant curves on X in detail, with a particular emphasis on curves of negative self-intersection. We use the representation theory of the stabilizers of the singular points to discover several invariant curves of negative self-intersection on X, and use these curves to study Nagata-type questions for linear series on X. The homogeneous ideal I of the collection of points in the configuration is an example of an ideal where the symbolic cube of the ideal is not contained in the square of the ideal; ideals with this property are seemingly quite rare. The resurgence and asymptotic resurgence are invariants which were introduced to measure such failures of containment. We use our knowledge of negative curves on X to compute the resurgence of I exactly. We also compute the asymptotic resurgence and Waldschmidt constant exactly in the case of the Wiman configuration of lines, and provide estimates on both for the Klein configuration.


title: Projective Curvature Tensors of Second Type Almost Geodesic Mappings

abstract: We consider equitorsion second type almost geodesic mappings of a non-symmetric affine connection space in this article. Using different computational methods, we obtained some invariants of these mappings. Last generalized Thomas projective parameter and Weyl projective tensor as invariants of a second type almost geodesic mapping of a non-symmetric affine connection space are further generalized here.


title: Detection of Faults in Power System Using Wavelet Transform and Independent Component Analysis

abstract: Uninterruptible power supply is the main motive of power utility companies that motivate them for identifying and locating the different types of faults as quickly as possible to protect the power system prevent complete power black outs using intelligent techniques. Thus, the present research work presents a novel method for detection of fault disturbances based on Wavelet Transform (WT) and Independent Component Analysis (ICA). The voltage signal is taken offline under fault conditions and is being processed through wavelet and ICA for detection. The time-frequency resolution from WT transform detects the fault initiation instant in the signal. Again, a performance index is calculated from independent component analysis under fault condition which is used to detect the fault disturbance in the voltage signal. The proposed approach is tested to be robust enough under various operating scenarios like without noise, with 20-dB noise and variation in frequency. Further, the detection study is carried out using a performance index, energy content, by applying the existing Fourier transform (FT), short time Fourier transform (STFT) and the proposed wavelet transform. Fault disturbances are detected if the energy calculated in each scenario is greater than the corresponding threshold value. The fault detection study is simulated in MATLAB/Simulink for a typical power system.


title: Strain-dependent internal parameters in hyperelastic biological materials

abstract: The behavior of hyperelastic energies depending on an internal parameter, which is a function of the deformation gradient, is discussed. As an example, the analysis of two models where the parameter describes the activation of a tetanized skeletal muscle tissue is presented. In those models, the activation parameter depends on the strain and it is shown the importance of considering the derivative of the parameter with respect to the strain in order to capture the proper stress-strain relations.


title: Thermodynamics of Giant Planet Formation: Shocking Hot Surfaces on Circumplanetary Disks

abstract: The luminosity of young giant planets can inform about their formation and accretion history. The directly imaged planets detected so far are consistent with the "hot-start" scenario of high entropy and luminosity. If nebular gas passes through a shock front before being accreted into a protoplanet, the entropy can be substantially altered. To investigate this, we present high resolution, 3D radiative hydrodynamic simulations of accreting giant planets. The accreted gas is found to fall with supersonic speed in the gap from the circumstellar disk's upper layers onto the surface of the circumplanetary disk and polar region of the protoplanet. There it shocks, creating an extended hot supercritical shock surface. This shock front is optically thick, therefore, it can conceal the planet's intrinsic luminosity beneath. The gas in the vertical influx has high entropy which when passing through the shock front decreases significantly while the gas becomes part of the disk and protoplanet. This shows that circumplanetary disks play a key role in regulating a planet's thermodynamic state. Our simulations furthermore indicate that around the shock surface extended regions of atomic - sometimes ionized - hydrogen develop. Therefore circumplanetary disk shock surfaces could influence significantly the observational appearance of forming gas-giants.


title: Predictions for the detection of Tidal Streams with Gaia using Great Circle Methods

abstract: The Gaia astrometric mission may offer an unprecedented opportunity to discover new tidal streams in the Galactic halo. To test this, we apply nGC3, a great-circle-cell count method that combines position and proper motion data to identify streams, to ten mock Gaia catalogues of K giants and RR Lyrae stars constructed from cosmological simulations of Milky Way analogues. We analyse two sets of simulations, one using a combination of $N$-body and semi-analytical methods which has extremely high resolution, the other using hydro-dynamical methods, which captures the dynamics of baryons, including the formation of an in situ halo. These ten realisations of plausible Galactic merger histories allow us to assess the potential for the recovery of tidal streams in different Milky Way formation scenarios. We include the Gaia~selection function and observational errors in these mock catalogues. We find that the nGC3 method has a well-defined detection boundary in the space of stream width and projected overdensity, that can be predicted based on direct observables alone. We predict that about 4-13 dwarf galaxy streams can be detected in a typical Milky Way-mass halo with Gaia+nGC3, with an estimated efficiency of $>$80\% inside the detection boundary. The progenitors of these streams are in the mass range of the classical dwarf galaxies and may have been accreted as early as redshift $\sim3$. Finally, we analyse how different possible extensions of the Gaia mission will improve the detection of tidal streams.


title: Spin equilibrium in strongly-magnetized accreting stars

abstract: Strongly magnetized accreting stars are often hypothesized to be in `spin equilibrium' with their surrounding accretion flows, which requires that the accretion rate changes more slowly than it takes the star to reach spin equilibrium. This is not true for most magnetically accreting stars, which have strongly variable accretion outbursts on time-scales much shorter than the time it would take to reach spin equilibrium. This paper examines how accretion outbursts affect the time a star takes to reach spin equilibrium and its final equilibrium spin period. I consider several different models for angular momentum loss -- either carried away in an outflow, lost to a stellar wind, or transferred back to the accretion disc (the `trapped disc'). For transient sources, the outflow scenario leads to significantly longer times to reach spin equilibrium ($\sim$10x), and shorter equilibrium spin periods than would be expected from spin equilibrium arguments, while the `trapped disc' does not. The results suggest that disc trapping plays a significant role in the spin evolution of strongly magnetic stars, with some caveats for young stellar objects.


title: An empirical analysis of the Ebola outbreak in West Africa

abstract: The data for the Ebola outbreak that occurred in 2014-2016 in three countries of West Africa are analysed within a common framework. The analysis is made using the results of an agent based Susceptible-Infected-Removed (SIR) model on a Euclidean network, where nodes at a distance $l$ are connected with probability $P(l) \propto l^{-\delta }$, $\delta$ determining the range of the interaction, in addition to nearest neighbors. The cumulative (total) density of infected population here has the form $R(t) = \frac{a\exp(t/T)}{1+c\exp(t/T)}$, where the parameters depend on $\delta$ and the infection probability $q$. This form is seen to fit well with the data. Using the best fitting parameters, the time at which the peak is reached is estimated and is shown to be consistent with the data. We also show that in the Euclidean model, one can choose $\delta$ and $q$ values which reproduce the data for the three countries qualitatively. These choices are correlated with population density, control schemes and other factors. Comparing the real data and the results from the model one can also estimate the size of the actual population susceptible to the disease. Rescaling the real data a reasonably good quantitative agreement with the simulation results is obtained.


title: Efficient Offline Waveform Design Using Quincunx/Hexagonal Time-Frequency Lattices For 5G Systems

abstract: Conventional OFDM, adopted in LTE-A systems, cannot provide the quality of service requirements sought in 5G systems because of extreme natural channel impairments caused by higher Doppler spreads and unexpected artificial impairments caused by multi-source transmission, to be brought by 5G, and by synchronization relaxation for closed-loop signaling overhead reduction in some 5G applications. These severe impairments induce a strong loss of orthogonality between subcarriers and OFDM symbols and, therefore, lead to a dramatic increase in ICI and ISI. To be well armed against these dramatic impairments, we, in the present paper, optimize the transmit/receive waveforms for FBMC systems, with hexagonal time-frequency lattices, operating over severe doubly dispersive channels, accounting for both natural and artificial impairments. For this, we exploit the POPS paradigm, recently proposed for rectangular time-frequency lattices, to design offline waveforms maximizing the SINR for hexagonal time-frequency lattices. We show that FBMC, with hexagonal lattices, offers a strong improvement in SIR with respect to conventional OFDM and an improvement of 1dB with respect to POPS-FBMC, with classical rectangular lattices. Furthermore, we show that the hexagonal POPS-FBMC brings more robustness to frequency synchronization errors and offers a 10dB reduction in OOB emissions with respect to rectangular POPS-FBMC.


title: Enhanced Dynamical Stability with Harmonic Slip-stacking

abstract: We develop a configuration of radio-frequency (rf) cavities to dramatically improve the performance of slip-stacking. Slip-stacking is an accumulation technique used at Fermilab to nearly double proton intensity by maintaining two beams of different momenta in the same storage ring. The two particle beams are longitudinally focused in the Recycler by two 53 MHz 100 kV rf cavities with a small frequency difference between them. We propose an additional 106 MHz 20 kV rf cavity with a frequency at the double the average of the upper and lower main rf frequencies. We show the harmonic rf cavity cancels out the resonances generated between the two main rf cavities and we derive the relationship between the harmonic rf voltage and the main rf voltage. We find the area factors that can be used to calculate the available phase space area for any set of beam parameters without individual simulation. We establish Booster beam quality requirements to achieve 99\% slip-stacking efficiency. We measure the longitudinal distribution of the Booster beam and use it to generate a realistic beam model for slip-stacking simulation. We demonstrate that the harmonic rf cavity can not only reduce particle loss during slip-stacking, but also reduce the final longitudinal emittance.


title: J-Frame Sequences in Krein Space

abstract: Let $\{f_n:n\in\mathbb{N}\}$ be a $J$-frame for a Krein space ${\textbf{\textit{K}}}$ and $P_M$ be a $J$-orthogonal projection from ${\textbf{\textit{K}}}$ onto a subspace $M$. In this article we find sufficient conditions under which $\{P_M(f_n):n\in\mathbb{N}\}$ is a $J$-frame for $P_M\textbf{\textit{K}}$ and $\{(I-P_M)f_n\}_{n\in{\mathbb{N}}}$ is a $J$-frame for $(I-P_M)\textbf{\textit{K}}$. We also introduce $J$-frame sequence for a Krein space ${\textbf{\textit{K}}}$ and study some properties of $J$-frame sequence analogues to Hilbert space frame theory.


title: Tight J-frames in Krein space and the associated J-frame potential

abstract: Motivated by the idea of $J$-frame for a Krein space $\textbf{\textit{K}}$, introduced by Giribet \textit{et al.} (J. I. Giribet, A. Maestripieri, F. Mart\'inez Per\'{i}a, P. G. Massey, \textit{On frames for Krein spaces}, J. Math. Anal. Appl. (1), {\bf 393} (2012), 122--137.), we introduce the notion of $\zeta-J$-tight frame for a Krein space $\textbf{\textit{K}}$. In this paper we characterize $J$-orthonormal basis for $\textbf{\textit{K}}$ in terms of $\zeta-J$-Parseval frame. We show that a Krein space is richly supplied with $\zeta-J$-Parseval frames. We also provide a necessary and sufficient condition when the linear sum of two $\zeta-J$-Parseval frames is again a $\zeta-J$-Parseval frame. We then generalize the notion of $J$-frame potential in Krein space from Hilbert space frame theory. Finally we provided a necessary and sufficient condition for a $J$-frame potential of the corresponding $\zeta-J$-tight frame to be minimum.


title: The Role of Star-Formation and AGN in Dust Heating of z=0.3-2.8 Galaxies - II. Informing IR AGN fraction estimates through simulations

abstract: A key question in extragalactic studies is the determination of the relative roles of stars and AGN in powering dusty galaxies at $z\sim$1-3 where the bulk of star-formation and AGN activity took place. In Paper I, we present a sample of $336$ 24$\mu$m-selected (Ultra)Luminous Infrared Galaxies, (U)LIRGs, at $z \sim 0.3$-$2.8$, where we focus on determining the AGN contribution to the IR luminosity. Here, we use hydrodynamic simulations with dust radiative transfer of isolated and merging galaxies, to investigate how well the simulations reproduce our empirical IR AGN fraction estimates and determine how IR AGN fractions relate to the UV-mm AGN fraction. We find that: 1) IR AGN fraction estimates based on simulations are in qualitative agreement with the empirical values when host reprocessing of the AGN light is considered; 2) for star-forming galaxy-AGN composites our empirical methods may be underestimating the role of AGN, as our simulations imply $>$50% AGN fractions, $\sim$3$\times$ higher than previous estimates; 3) 6% of our empirically classified "SFG" have AGN fractions $\gtrsim$50%. While this is a small percentage of SFGs, if confirmed, would imply the true number density of AGN may be underestimated; 4) this comparison depends on the adopted AGN template -- those that neglect the contribution of warm dust lower the empirical fractions by up to 2$\times$; and 5) the IR AGN fraction is only a good proxy for the intrinsic UV-mm AGN fraction when the extinction is high ($A_V\gtrsim 1$ or up to and including coalescence in a merger).


title: Task Specific Adversarial Cost Function

abstract: The cost function used to train a generative model should fit the purpose of the model. If the model is intended for tasks such as generating perceptually correct samples, it is beneficial to maximise the likelihood of a sample drawn from the model, Q, coming from the same distribution as the training data, P. This is equivalent to minimising the Kullback-Leibler (KL) distance, KL[Q||P]. However, if the model is intended for tasks such as retrieval or classification it is beneficial to maximise the likelihood that a sample drawn from the training data is captured by the model, equivalent to minimising KL[P||Q]. The cost function used in adversarial training optimises the Jensen-Shannon entropy which can be seen as an even interpolation between KL[Q||P] and KL[P||Q]. Here, we propose an alternative adversarial cost function which allows easy tuning of the model for either task. Our task specific cost function is evaluated on a dataset of hand-written characters in the following tasks: Generation, retrieval and one-shot learning.


title: Inverted anisotropy in the in-plane dissipation of BSCCO tapes: preliminary results

abstract: The electrical properties of superconducting tapes and coatings in the direction transverse to the long dimension of the composite has been rarely studied. However, transverse dissipation can eventually determine the behavior of a transmission line in the case of failure due to the presence of transversal cracks, and is also fundamental in the AC regime. In this paper we present a preliminary experimental study of the electrical transport properties along the transverse direction of BSCCO-metal tapes, and compare them with those measured along the long axis of the material. In spite of the fact that the tapes under study are not multi-filamentary, our experiments suggest that there is a measurable anisotropy of the transport properties between the longitudinal and transverse directions.


title: Learning Genomic Representations to Predict Clinical Outcomes in Cancer

abstract: Genomics are rapidly transforming medical practice and basic biomedical research, providing insights into disease mechanisms and improving therapeutic strategies, particularly in cancer. The ability to predict the future course of a patient's disease from high-dimensional genomic profiling will be essential in realizing the promise of genomic medicine, but presents significant challenges for state-of-the-art survival analysis methods. In this abstract we present an investigation in learning genomic representations with neural networks to predict patient survival in cancer. We demonstrate the advantages of this approach over existing survival analysis methods using brain tumor data.


title: Universal hierarchical structure of quasiperiodic eigenfunctions

abstract: We determine exact exponential asymptotics of eigenfunctions and of corresponding transfer matrices of the almost Mathieu operators for all frequencies in the localization regime. This uncovers a universal structure in their behavior, governed by the continued fraction expansion of the frequency, explaining some predictions in physics literature. In addition it proves the arithmetic version of the frequency transition conjecture. Finally, it leads to an explicit description of several non-regularity phenomena in the corresponding non-uniformly hyperbolic cocycles, which is also of interest as both the first natural example of some of those phenomena and, more generally, the first non-artificial model where non-regularity can be explicitly studied.


title: A Bayesian Risk Approach to Data-driven Stochastic Optimization: Formulations and Asymptotics

abstract: A large class of stochastic programs involve optimizing an expectation taken with respect to an underlying distribution that is unknown in practice. One popular approach to addressing the distributional uncertainty, known as the distributionally robust optimization (DRO), is to hedge against the worst case over an uncertainty set of candidate distributions. However, it has been observed that inappropriate construction of the uncertainty set can sometimes result in over-conservative solutions. To explore the middle ground between optimistically ignoring the distributional uncertainty and pessimistically fixating on the worst-case scenario, we propose a Bayesian risk optimization (BRO) framework for parametric underlying distributions, which is to optimize a risk functional applied to the posterior distribution of an unknown distribution parameter. Of our particular interest are four risk functionals: mean, mean-variance, value-at-risk, and conditional value-at-risk. To unravel the implication of BRO, we establish the consistency of objective functions and optimal solutions, as well as the asymptotic normality of objective functions and optimal values. More importantly, our analysis reveals a hidden interpretation: the objectives of BRO can be approximately viewed as a weighted sum of posterior mean objective and the (squared) half-width of the true objective's confidence interval.


title: Low and High Energy Asymptotic Behavior of Electroweak Corrections in Polarized $e^+ e^- \to \mu^+ \mu^-$ Process

abstract: Electroweak radiative corrections will play a major role in the analysis of several upcoming ultra-precision experiments such as Belle-II, so is crucial to make sure that they are fully under control. The article outlines the recent developments in the theoretical and computational approaches to one-loop (NLO) electroweak radiative corrections to the parity-violating asymmetry in $e^- e^+ \to \mu^- \mu^+ (\gamma)$ process with longitudinally polarized electrons. We derive asymptotic expressions for low and high energy regions (well below or above $Z$-resonance, correspondingly) and analyze the leading contributions. For most of energy regions, our results are in good agreement with precise computer-algebra based calculation and can used as a quicker alternative.


title: Deep Reinforcement Learning for Mention-Ranking Coreference Models

abstract: Coreference resolution systems are typically trained with heuristic loss functions that require careful tuning. In this paper we instead apply reinforcement learning to directly optimize a neural mention-ranking model for coreference evaluation metrics. We experiment with two approaches: the REINFORCE policy gradient algorithm and a reward-rescaled max-margin objective. We find the latter to be more effective, resulting in significant improvements over the current state-of-the-art on the English and Chinese portions of the CoNLL 2012 Shared Task.


title: Singular elliptic equation involving the GJMS operator on the standard unit sphere

abstract: Given a Riemannian compact manifold (M,g) of dimension n>4, we have proven in [1] under some conditions that the equation : Pg(u) = Bu +Au2+Cu (1) where Pg is the GJMS-operator, n = dim(M) > 2k, A, B and C are smooth positive functions on M, p > 1 and 2] denotes the critical Sobolev admits twodistinct positive solutions. The proof of this result is essentially based on the given smooth function ' > 0 with norm k'kPg = 1 fulfilling some conditions ( see Theorem 3 in [1]). In this note we construct an example of such function on the unit standard sphere (Sn; h). Con- sequently the conditions of the Theorem are improved in the case of (Sn; h)