title: The case for 100 GeV bino dark matter: A dedicated LHC tri-lepton search

abstract: Global fit studies performed in the pMSSM and the photon excess signal originating from the Galactic Center seem to suggest compressed electroweak supersymmetric spectra with a $\sim$100 GeV bino-like dark matter particle. We find that these scenarios are not probed by traditional electroweak supersymmetry searches at the LHC. We propose to extend the ATLAS and CMS electroweak supersymmetry searches with an improved strategy for bino-like dark matter, focusing on chargino plus next-to-lightest neutralino production, with a subsequent decay into a tri-lepton final state. We explore the sensitivity for pMSSM scenarios with $\Delta m = m_{\rm NLSP} - m_{\rm LSP} \sim (5 - 50)$ GeV in the $\sqrt{s} = 14$ TeV run of the LHC. Counterintuitively, we find that the requirement of low missing transverse energy increases the sensitivity compared to the current ATLAS and CMS searches. With 300 fb$^{-1}$ of data we expect the LHC experiments to be able to discover these supersymmetric spectra with mass gaps down to $\Delta m \sim 9$ GeV for DM masses between 40 and 140 GeV. We stress the importance of a dedicated search strategy that targets precisely these favored pMSSM spectra.


title: NEXT: In-Network Nonconvex Optimization

abstract: We study nonconvex distributed optimization in multi-agent networks with time-varying (nonsymmetric) connectivity. We introduce the first algorithmic framework for the distributed minimization of the sum of a smooth (possibly nonconvex and nonseparable) function - the agents' sum-utility - plus a convex (possibly nonsmooth and nonseparable) regularizer. The latter is usually employed to enforce some structure in the solution, typically sparsity. The proposed method hinges on successive convex approximation techniques while leveraging dynamic consensus as a mechanism to distribute the computation among the agents: each agent first solves (possibly inexactly) a local convex approximation of the nonconvex original problem, and then performs local averaging operations. Asymptotic convergence to (stationary) solutions of the nonconvex problem is established. Our algorithmic framework is then customized to a variety of convex and nonconvex problems in several fields, including signal processing, communications, networking, and machine learning. Numerical results show that the new method compares favorably to existing distributed algorithms on both convex and nonconvex problems.


title: Mean field limit of interacting filaments and vector valued non linear PDEs

abstract: Families of $N$ interacting curves are considered, with long range, mean field type, interaction. A family of curves defines a 1-current, concentrated on the curves, analog of the empirical measure of interacting point particles. This current is proved to converge, as $N$ goes to infinity, to a mean field current, solution of a nonlinear, vector valued, partial differential equation. In the limit, each curve interacts with the mean field current and two different curves have an independence property if they are independent at time zero. This set-up is inspired from vortex filaments in turbulent fluids, although for technical reasons we have to restrict to smooth interaction, instead of the singular Biot-Savart kernel. All these results are based on a careful analysis of a nonlinear flow equation for 1-currents, its relation with the vector valued PDE and the continuous dependence on the initial conditions.


title: Investigation of pionic contribution in the lepton and anti-lepton production cross section in p-Cu and p-Pt collision

abstract: For detailed explanation of the experimental results of lepton production cross section in hadronic collisions such as nucleon-nucleon or nucleon-nuclei, it is of great importance to use quarks and sea quarks distribution function inside free and bound nucleons. In this paper the role of pion cloud inside the nucleus in the structure function of Cu and Pt nuclei and the EMC ratio of these nuclei were investigated by using harmonic oscillator model.


title: Randomization and sparsity in huge-scale optimization on the Mirror Descent example

abstract: We investigate different randomizations for mirror descent method. We try to propose such a randomization that allows us to use sparsity of the problem as much as it possible. In the paper one can also find a generalization of randomizaed mirror descent for the convex optimization problems with functional restrictions.


title: Dynamic Heterogeneity In The Monoclinic Phase Of CCl$_4$

abstract: Carbon tetrachloride (CCl$_4$) is one of the simplest compounds having a translationally stable monoclinic phase while exhibiting a rich rotational dynamics below 226 K. Recent nuclear quadrupolar resonance (NQR) experiments revealed that the dynamics of CCl$_4$ is similar to that of the other members of the isostructural series CBr$_{n}$Cl$_{4-n}$, suggesting that the universal relaxation features of canonical glasses such as $\alpha$- and $\beta$-relaxation are also present in non-glass formers. Using molecular dynamics (MD) simulations we studied the rotational dynamics in the monoclinic phase of CCl$_4$. The molecules undergo $C3$ type jump-like rotations around each one of the four C-Cl bonds. The rotational dynamics is very well described with a master equation using as the only input the rotational rates measured from the simulated trajectories. It is found that the heterogeneous dynamics emerges from faster and slower modes associated with different rotational axes, which have fixed orientations relative to the crystal and are distributed among the four non-equivalent molecules of the unit cell.


title: A Decentralized Second-Order Method with Exact Linear Convergence Rate for Consensus Optimization

abstract: This paper considers decentralized consensus optimization problems where different summands of a global objective function are available at nodes of a network that can communicate with neighbors only. The proximal method of multipliers is considered as a powerful tool that relies on proximal primal descent and dual ascent updates on a suitably defined augmented Lagrangian. The structure of the augmented Lagrangian makes this problem non-decomposable, which precludes distributed implementations. This problem is regularly addressed by the use of the alternating direction method of multipliers. The exact second order method (ESOM) is introduced here as an alternative that relies on: (i) The use of a separable quadratic approximation of the augmented Lagrangian. (ii) A truncated Taylor's series to estimate the solution of the first order condition imposed on the minimization of the quadratic approximation of the augmented Lagrangian. The sequences of primal and dual variables generated by ESOM are shown to converge linearly to their optimal arguments when the aggregate cost function is strongly convex and its gradients are Lipschitz continuous. Numerical results demonstrate advantages of ESOM relative to decentralized alternatives in solving least squares and logistic regression problems.


title: Revisiting Zariski Main Theorem from a constructive point of view

abstract: This paper deals with the Peskine version of Zariski Main Theorem published in 1965 and discusses some applications. It is written in the style of Bishop's constructive mathematics. Being constructive, each proof in this paper can be interpreted as an algorithm for constructing explicitly the conclusion from the hypothesis. The main non-constructive argument in the proof of Peskine is the use of minimal prime ideals. Essentially we substitute this point by two dynamical arguments; one about gcd's, using subresultants, and another using our notion of strong transcendence. In particular we obtain algorithmic versions for the Multivariate Hensel Lemma and the structure theorem of quasi-finite algebras.


title: Stepwise Positional-Orientational Order and the Multicritical-Multistructural Global Phase Diagram of the s=3/2 Ising Model from Renormalization-Group Theory

abstract: The spin-3/2 Ising model, with nearest-neighbor interactions only, is the prototypical system with two different ordering species, with concentrations regulated by a chemical potential. Its global phase diagram, obtained in d=3 by renormalization-group theory in the Migdal-Kadanoff approximation or equivalently as an exact solution of a d=3 hierarchical lattice, with flows subtended by 40 different fixed points, presents a very rich structure containing eight different ordered and disordered phases, with more than fourteen different types of phase diagrams in temperature and chemical potential. It exhibits phases with orientational and/or positional order. It also exhibits quintuple phase transition reentrances. Universality of critical exponents is conserved across different renormalization-group flow basins, via redundant fixed points. One of the phase diagrams contains a plastic crystal sequence, with positional and orientational ordering encountered consecutively as temperature is lowered. The global phase diagram also contains double critical points, first-order and critical lines between two ordered phases, critical endpoints, usual and unusual (inverted) bicritical points, tricritical points, multiple tetracritical points, and zero-temperature criticality and bicriticality. The 4-state Potts permutation-symmetric subspace is contained in this model.


title: Effects of neutron-star dynamic tides on gravitational waveforms within the effective-one-body approach

abstract: Extracting the unique information on ultradense nuclear matter from the gravitational waves emitted by merging, neutron-star binaries requires robust theoretical models of the signal. We develop a novel effective-one-body waveform model that includes, for the first time, dynamic (instead of only adiabatic) tides of the neutron star, as well as the merger signal for neutron-star--black-hole binaries. We demonstrate the importance of the dynamic tides by comparing our model against new numerical-relativity simulations of nonspinning neutron-star--black-hole binaries spanning more than 24 gravitational-wave cycles, and to other existing numerical simulations for double neutron-star systems. Furthermore, we derive an effective description that makes explicit the dependence of matter effects on two key parameters: tidal deformability and fundamental oscillation frequency.


title: Separating Putative Pathogens from Background Contamination with Principal Orthogonal Decomposition: Evidence for Leptospira in the Ugandan Neonatal Septisome

abstract: Neonatal sepsis (NS) is responsible for over a 1 million yearly deaths worldwide. In the developing world NS is often treated without an identified microbial pathogen. Amplicon sequencing of the bacterial 16S rRNA gene can be used to identify organisms that are difficult to detect by routine microbiological methods. However, contaminating bacteria are ubiquitous in both hospital settings and research reagents, and must be accounted for to make effective use of these data. In the present study, we sequenced the bacterial 16S rRNA gene obtained from blood and cerebrospinal fluid (CSF) of 80 neonates presenting with NS to the Mbarara Regional Hospital in Uganda. Assuming that patterns of background contamination would be independent of pathogenic microorganism DNA, we applied a novel quantitative approach using principal orthogonal decomposition to separate background contamination from potential pathogens in sequencing data. We designed our quantitative approach contrasting blood, CSF, and control specimens, and employed a variety of statistical random matrix bootstrap hypotheses to estimate statistical significance. These analyses demonstrate that Leptospira appears present in some infants presenting within 48 hr of birth, indicative of infection in utero, and up to 28 days of age, suggesting environmental exposure. This organism cannot be cultured in routine bacteriological settings, and is enzootic in the cattle that the rural peoples of western Uganda often live in close proximity. Our findings demonstrate that statistical approaches to remove background organisms common in 16S sequence data can reveal putative pathogens in small volume biological samples from newborns. This computational analysis thus reveals an important medical finding that has the potential to alter therapy and prevention efforts in a critically ill population.


title: Dynamical Selection of Critical Exponents

abstract: In renormalized field theories there are in general one or few fixed points which are accessible by the renormalization-group flow. They can be identified from the fixed-point equations. Exceptionally, an infinite family of fixed points exists, parameterized by a scaling exponent $\zeta$, itself function of a non-renormalizing parameter. Here we report a different scenario with an infinite family of fixed points of which seemingly only one is chosen by the renormalization-group flow. This dynamical selection takes place in systems with an attractive interaction ${\cal V}(\phi)$, as in standard $\phi^4$ theory, but where the potential $\cal V$ at large $\phi$ goes to zero, as e.g. the attraction by a defect.


title: Virtual Machine Warmup Blows Hot and Cold

abstract: Virtual Machines (VMs) with Just-In-Time (JIT) compilers are traditionally thought to execute programs in two phases: the initial warmup phase determines which parts of a program would most benefit from dynamic compilation, before JIT compiling those parts into machine code; subsequently the program is said to be at a steady state of peak performance. Measurement methodologies almost always discard data collected during the warmup phase such that reported measurements focus entirely on peak performance. We introduce a fully automated statistical approach, based on changepoint analysis, which allows us to determine if a program has reached a steady state and, if so, whether that represents peak performance or not. Using this, we show that even when run in the most controlled of circumstances, small, deterministic, widely studied microbenchmarks often fail to reach a steady state of peak performance on a variety of common VMs. Repeating our experiment on 3 different machines, we found that at most 43.5% of <VM, benchmark> pairs consistently reach a steady state of peak performance.


title: Robust flux error estimation of an unfitted Nitsche method for high-contrast interface problems

abstract: We prove an optimal error estimate for the flux variable for a stabilized unfitted Nitsche finite element method applied to an elliptic interface problem with discontinuous constant coefficients. Our result shows explicitly that this error estimate is totally independent of the diffusion coefficients


title: Emergent smectic order in simple active particle models

abstract: Novel "smectic-P" behavior, in which self-propelled particles form rows and move on average along them, occurs generically within the orientationally-ordered phase of simple models that we simulate. Both apolar (head-tail symmetric) and polar (head-tail asymmetric) models with aligning and repulsive interactions exhibit slow algebraic decay of smectic order with system size up to some finite length scale, after which faster decay occurs. In the apolar case, this scale is that of an undulation instability of the rows. In the polar case, this instability is absent, but traveling fluctuations disrupt the rows in large systems and motion and smectic order may spontaneously globally rotate. These observations agree with a new hydrodynamic theory which we present here. Variants of our models also exhibit active smectic "A" and "C" order, with motion orthogonal and oblique to the layers respectively.


title: A Borel-reducibility Counterpart of Shelah's Main Gap Theorem

abstract: We study the Borel-reducibility of isomorphism relations of complete first order theories and show the consistency of the following: For all such theories T and T', if T is classifiable and T' is not, then the isomorphism of models of T' is strictly above the isomorphism of models of T with respect to Borel-reducibility. In fact, we can also ensure that a range of equivalence relations modulo various non-stationary ideals are strictly between those isomorphism relations. The isomorphism relations are considered on models of some fixed uncountable cardinality obeying certain restrictions.


title: Young and middle age pulsar light-curve morphology: Comparison of Fermi observations with gamma-ray and radio emission geometries

abstract: Thanks to the huge amount of gamma-ray pulsar photons collected by the Fermi Large Area Telescope since June 2008, it is now possible to constrain gamma-ray geometrical models by comparing simulated and observed light-curve morphological characteristics. We assumed vacuum-retarded dipole pulsar magnetic field and tested simulated and observed morphological light-curve characteristics in the framework of two pole emission geometries, Polar Cap (PC), radio, and Slot Gap (SG), and Outer Gap (OG)/One Pole Caustic (OPC) emission geometries. We compared simulated and observed/estimated light-curve morphological parameters as a function of observable and non-observable pulsar parameters. The PC model gives the poorest description of the LAT pulsar light-curve morphology. The OPC best explains both the observed gamma-ray peak multiplicity and shape classes. The OPC and SG models describe the observed gamma-ray peak-separation distribution for low- and high-peak separations, respectively. This suggests that the OPC geometry best explains the single-peak structure but does not manage to describe the widely separated peaks predicted in the framework of the SG model as the emission from the two magnetic hemispheres. The OPC radio-lag distribution shows higher agreement with observations suggesting that assuming polar radio emission, the gamma-ray emission regions are likely to be located in the outer magnetosphere. The larger agreement between simulated and LAT estimations in the framework of the OPC suggests that the OPC model best predicts the observed variety of profile shapes. The larger agreement between observations and the OPC model jointly with the need to explain the abundant 0.5 separated peaks with two-pole emission geometries, calls for thin OPC gaps to explain the single-peak geometry but highlights the need of two-pole caustic emission geometry to explain widely separated peaks.


title: Standing Waves Braneworlds

abstract: The class of non-stationary braneworld models generated by the coupled gravitational and scalar fields is reviewed. The model represents a brane in a space-time with single time and one large (infinite) and several small (compact) space-like extra dimensions. In some particular cases the model has the solutions corresponding to the bulk gravi-scalar standing waves bounded by the brane. Pure gravitational localization mechanism of matter particles on the node of standing waves, where the brane is placed, is discussed. Cosmological applications of the model also considered.


title: Quantum theory of the Generalised Uncertainty Principle

abstract: We extend significantly previous works on the Hilbert space representations of the Generalized Uncertainty Principle (GUP) in 3+1 dimensions of the form $[X_i,P_j] = i F_{ij}$ where $ F_{ij} = f(P^2) \delta_{ij} + g(P^2) P_i P_j $ for any functions $f$. However, we restrict our study to the case of commuting $X$'s. We focus in particular on the symmetries of the theory, and the minimal length that emerge in some cases. We first show that, at the algebraic level, there exists an unambiguous mapping between the GUP with a deformed quantum algebra and a quadratic Hamiltonian into a standard, Heisenberg algebra of operators and an aquadratic Hamiltonian, provided the boost sector of the symmetries is modified accordingly. The theory can also be mapped to a completely standard Quantum Mechanics with standard symmetries, but with momentum dependent position operators. Next, we investigate the Hilbert space representations of these algebraically equivalent models, and focus, specifically on whether they exhibit a minimal length. We carry the functional analysis of the various operators involved, and show that the appearance of a minimal length critically depends on the relationship between the generators of translations and the physical momenta. In particular, because this relationship is preserved by the algebraic mapping presented in this paper, when a minimal length is present in the standard GUP, it is also present in the corresponding Aquadratic Hamiltonian formulation, despite the perfectly standard algebra of this model. In general, a minimal length requires bounded generators of translations, i.e. a specific kind of quantization of space, and this depends on the precise shape of the function $f$ defined previously. This result provides an elegant and unambiguous classification of which universal quantum gravity corrections lead to the emergence of a minimal length.


title: Emergence of transverse spin in optical modes of semiconductor nanowires

abstract: The transverse spin angular momentum of light has recently received tremendous attention as it adds a new degree of freedom for controlling light-matter interactions. In this work we demonstrate the generation of transverse spin angular momentum by the weakly-guided mode of semiconductor nanowires. The evanescent field of these modes in combination with the transversality condition rigorously accounts for the occurrence of transverse spin angular momentum. The intriguing and nontrivial spin properties of optical modes in semiconductor nanowires are of high interest for a broad range of new applications including chiral optical trapping, quantum information processing, and nanophotonic circuitry.


title: Integral formulae for codimension-one foliated Finsler manifolds

abstract: We study extrinsic geometry of a codimension-one foliation ${\cal F}$ of a closed Finsler space $(M,F)$, in particular, of a Randers space $(M,\alpha+\beta)$. Using a unit vector field $\nu$ orthogonal (in the Finsler sense) to the leaves of ${\cal F}$ we define a new Riemannian metric $g$ on $M$, which for Randers case depends nicely on $(\alpha,\beta)$. For that $g$ we derive several geometric invariants of ${\cal F}$ (e.g. the Riemann curvature and the shape operator) in terms of $F$, then under natural assumptions on $\beta$ which simplify derivations, we express them in terms of corresponding invariants arising from $\alpha$ and $\beta$. Using our approach (2012), we produce the integral formulae for ${\cal F}$ on $(M, F)$ and $(M, \alpha+\beta)$, which relate integrals of mean curvatures with those involving algebraic invariants obtained from the shape operator of a foliation, and the Riemann curvature in the direction $\nu$. They generalize the formulae by Brito, Langevin and Rosenberg, which show that total mean curvatures (of arbitrary order $k$) for codimension-one foliations on a closed $(m+1)$-dimensional manifold of constant curvature $K$ don't depend on a choice of ${\cal F}$.


title: Cosmology with velocity dispersion counts: an alternative to measuring cluster halo masses

abstract: The evolution of galaxy cluster counts is a powerful probe of several fundamental cosmological parameters. A number of recent studies using this probe have claimed tension with the cosmology preferred by the analysis of the Planck primary CMB data, in the sense that there are fewer clusters observed than predicted based on the primary CMB cosmology. One possible resolution to this problem is systematic errors in the absolute halo mass calibration in cluster studies, which is required to convert the standard theoretical prediction (the halo mass function) into counts as a function of the observable (e.g., X-ray luminosity, Sunyaev-Zel'dovich flux, optical richness). Here we propose an alternative strategy, which is to directly compare predicted and observed cluster counts as a function of the one-dimensional velocity dispersion of the cluster galaxies. We argue that the velocity dispersion of groups/clusters can be theoretically predicted as robustly as mass but, unlike mass, it can also be directly observed, thus circumventing the main systematic bias in traditional cluster counts studies. With the aid of the BAHAMAS suite of cosmological hydrodynamical simulations, we demonstrate the potential of the velocity dispersion counts for discriminating even similar $\Lambda$CDM models. These predictions can be compared with the results from existing redshift surveys such as the highly-complete Galaxy And Mass Assembly (GAMA) survey, and upcoming wide-field spectroscopic surveys such as the Wide Area Vista Extragalactic Survey (WAVES) and the Dark Energy Survey Instrument (DESI).


title: Coulomb and quenching effects in small nanoparticle-based spasers

abstract: We study numerically the effect of mode mixing and direct dipole-dipole interactions between gain molecules on spasing in a small composite nanoparticles with a metallic core and a dye-doped dielectric shell. By combining Maxwell-Bloch equations with Green's function formalism, we calculate lasing frequency and threshold population inversion for various gain densities in the shell. We find that gain coupling to nonresonant plasmon modes has a negligible effect on spasing threshold. In contrast, the direct dipole-dipole coupling, by causing random shifts of gain molecules' excitation frequencies, hinders reaching the spasing threshold in small systems. We identify a region of parameter space in which spasing can occur considering these effects.


title: Static electric multipole susceptibilities of the relativistic hydrogen-like atom in the ground state: Application of the Sturmian expansion of the generalized Dirac-Coulomb Green function

abstract: The ground state of the Dirac one-electron atom, placed in a weak, static electric field of definite $2^{L}$-polarity, is studied within the framework of the first-order perturbation theory. The Sturmian expansion of the generalized Dirac-Coulomb Green function [R. Szmytkowski, J. Phys. B 30 (1997) 825, erratum: 30 (1997) 2747] is used to derive closed-form analytical expressions for various far-field and near-nucleus static electric multipole susceptibilities of the atom. The far-field multipole susceptibilities --- the polarizabilities $\alpha_{L}$, electric-to-magnetic cross-susceptibilities $\alpha_{\mathrm{E}L\to\mathrm{M}(L\mp1)}$ and electric-to-toroidal-magnetic cross-susceptibilities $\alpha_{\mathrm{E}L\to\mathrm{T}L}$ --- are found to be expressible in terms of one or two non-terminating generalized hypergeometric functions ${}_{3}F_{2}$ with the unit argument. Counterpart formulas for the near-nucleus multipole susceptibilities --- the electric nuclear shielding constants $\sigma_{\mathrm{E}L\to\mathrm{E}L}$, near-nucleus electric-to-magnetic cross-susceptibilities $\sigma_{\mathrm{E}L\to\mathrm{M}(L\mp1)}$ and near-nucleus electric-to-toroidal-magnetic cross-susceptibilities $\sigma_{\mathrm{E}L\to\mathrm{T}L}$ --- involve terminating ${}_{3}F_{2}(1)$ series and for each $L$ may be rewritten in terms of elementary functions. Exact numerical values of the far-field dipole, quadrupole, octupole and hexadecapole susceptibilities are provided for selected hydrogenic ions. Analytical quasi-relativistic approximations, valid to the second order in $\alpha Z$, where $\alpha$ is the fine-structure constant and $Z$ is the nuclear charge number, are derived for all types of the far-field and near-nucleus susceptibilities considered in the paper.


title: Consensus, Cohesion and Connectivity

abstract: Social life clusters into groups held together by ties that also transmit information. When collective problems occur, group members use their ties to discuss what to do and to establish an agreement, to be reached quick enough to prevent discounting the value of the group decision. The speed at which a group reaches consensus can be predicted by the algebraic connectivity of the network, which also imposes a lower bound on the group's cohesion. This specific measure of connectivity is put to the test by re-using experimental data, which confirm the prediction.


title: EMFS: Repurposing SMTP and IMAP for Data Storage and Synchronization

abstract: Cloud storage has become a massive and lucrative business, with companies like Apple, Microsoft, Google, and Dropbox providing hundreds of millions of clients with synchronized and redundant storage. These services often command price-to-storage ratios significantly higher than the market rate for physical storage, as well as increase the surface area for data leakage. In place of this consumer-unfriendly status quo, I propose using widely available, well standardized email protocols like SMTP and IMAP in conjunction with free email service providers to store, synchronize, and share files across discrete systems.


title: Magneto-optical polarization rotation in a ladder-type atomic system for tunable offset locking

abstract: We demonstrate an easily tunable locking scheme for stabilizing frequency-sum of two lasers on a two-photon ladder transition based on polarization rotation in warm rubidium vapors induced by magnetic field and circularly polarized drive field. Unprecedented tunability of the two-photon offset frequency is due to strong splitting and shifting of magnetic states in external field. In our experimental setup we achieve two-photon detuning of up to 700 MHz.


title: On the number of limit cycles for perturbed pendulum equations

abstract: We consider perturbed pendulum-like equations on the cylinder of the form $ \ddot x+\sin(x)= \varepsilon \sum_{s=0}^{m}{Q_{n,s} (x)\, \dot x^{s}}$ where $Q_{n,s}$ are trigonometric polynomials of degree $n$, and study the number of limit cycles that bifurcate from the periodic orbits of the unperturbed case $\varepsilon=0$ in terms of $m$ and $n$. Our first result gives upper bounds on the number of zeros of its associated first order Melnikov function, in both the oscillatory and the rotary regions. These upper bounds are obtained expressing the corresponding Abelian integrals in terms of polynomials and the complete elliptic functions of first and second kind. Some further results give sharp bounds on the number of zeros of these integrals by identifying subfamilies which are shown to be Chebyshev systems.


title: Chiral magnetic conductivity and surface states of Weyl semimetals in topological insulator ultra-thin film multilayer

abstract: We investigate an ultra-thin film of topological insulator (TI) multilayer as a model for a three-dimensional (3D) Weyl semimetal. We introduce tunneling parameters $t_S$, $t_\perp$, and $t_D$, where the former two parameters couple layers of the same thin film at small and large momenta, and the latter parameter couples neighbouring thin film layers along the $z$-direction. The Chern number is computed in each topological phase of the system and we find that for $t_S,t_D>0$, the tunneling parameter $t_\perp$ changes from positive to negative as the system transits from Weyl semi-metallic phase to insulating phases. We further study the chiral magnetic effect (CME) of the system in the presence of a time dependent magnetic field. We compute the low-temperature dependence of the chiral magnetic conductivity and show that it captures three distinct phases of the system separated by plateaus. Furthermore, we propose and study a 3D lattice model of Porphyrin thin film, an organic material known to support topological Frenkel exciton edge states. We show that this model exhibits a 3D Weyl semi-metallic phase and also supports a 2D Weyl semi-metallic phase. We further show that this model recovers that of 3D Weyl semimetal in topological insulator thin film multilayer. Thus, paving the way for simulating a 3D Weyl semimetal in topological insulator thin film multilayer. We obtain the surface states (Fermi arcs) in the 3D model and the chiral edge states in the 2D model and analyze their topological properties.


title: Stock loans with liquidation

abstract: We derive a "semi-analytic" solution for a stock loan in which the lender forces liquidation when the loan-to-collateral ratio drops beneath a certain threshold. We use this to study the sensitivity of the contract to model parameters.


title: Augmented kludge waveforms and Gaussian process regression for EMRI data analysis

abstract: Extreme-mass-ratio inspirals (EMRIs) will be an important type of astrophysical source for future space-based gravitational-wave detectors. There is a trade-off between accuracy and computational speed for the EMRI waveform templates required in the analysis of data from these detectors. We discuss how the systematic error incurred by using faster templates may be reduced with improved models such as augmented kludge waveforms, and marginalised over with statistical techniques such as Gaussian process regression.


title: On pattern matching with k mismatches and few don't cares

abstract: We consider the problem of pattern matching with $k$ mismatches, where there can be don't care or wild card characters in the pattern. Specifically, given a pattern $P$ of length $m$ and a text $T$ of length $n$, we want to find all occurrences of $P$ in $T$ that have no more than $k$ mismatches. The pattern can have don't care characters, which match any character. Without don't cares, the best known algorithm for pattern matching with $k$ mismatches has a runtime of $O(n\sqrt{k \log k})$. With don't cares in the pattern, the best deterministic algorithm has a runtime of $O(nk polylog m)$. Therefore, there is an important gap between the versions with and without don't cares. In this paper we give an algorithm whose runtime increases with the number of don't cares. We define an {\em island} to be a maximal length substring of $P$ that does not contain don't cares. Let $q$ be the number of islands in $P$. We present an algorithm that runs in $O(n\sqrt{k\log m}+n\min\{\sqrt[3]{qk\log^2 m},\sqrt{q\log m}\})$ time. If the number of islands $q$ is $O(k)$ this runtime becomes $O(n\sqrt{k\log m})$, which essentially matches the best known runtime for pattern matching with $k$ mismatches without don't cares. If the number of islands $q$ is $O(k^2)$, this algorithm is asymptotically faster than the previous best algorithm for pattern matching with $k$ mismatches with don't cares in the pattern.


title: Asteroid 4 Vesta: dynamical and collisional evolution during the Late Heavy Bombardment

abstract: Vesta is the only currently identified asteroid for which we possess samples, which revealed us that the asteroid is differentiated and possesses a relatively thin basaltic crust that survived to the evolution of the asteroid belt and the Solar System. However, little is know about the effects of past events like the Late Heavy Bombardment on this crust. We address this gap in our knowledge by simulating the LHB in the different dynamical scenarios proposed for the migration of the giant planets in the broad framework of the Nice Model. The results of simulations generate information about produced crater population, surface saturation, mass loss and mass gain of Vesta and number of energetic or catastrophic impacts during LHB. Our results reveal that planet-planet scattering is a dynamically favourable migration mechanism for the survival of Vesta and its crust. The number of impacts on Vesta estimated as due to the LHB is $31\pm5$, i.e. about 5 times larger than the number of impacts that would have occurred in an unperturbed main belt in the same time interval. The contribution of a possible extended belt, instead, is quite limited and can be quantified in $2\pm1$ impacts. The chance of energetic and catastrophic impacts is less than 10\% and is compatible with the absence of giant craters dated back to 4 Ga ago and with the survival of the asteroid during the LHB. The mass loss translates in the erosion of $3-5$ meters of the crust, consistently with the global survival of the basaltic crust of Vesta confirmed by the Dawn mission. Our analysis revealed that the contribution of the LHB to the cratering of Vesta' surface is not significant and is actually erased by the crater population produced by the following 4 Ga of collisional evolution of the asteroid, in agreement with the data provided by the Dawn mission.


title: Ordering multiple soft gluon emissions

abstract: We present an expression for the QCD amplitude for a general hard scattering process with any number of soft gluon emissions, to one-loop accuracy. The amplitude is written in two different but equivalent ways: as a product of operators ordered in dipole transverse momentum and as a product of loop-expanded currents. We hope that these results will help in the development of an all-orders algorithm for multiple emissions that includes the full colour structure and both the real and imaginary contributions to the amplitude.


title: Influence of cooling on dynamics of buoyant jet

abstract: The Rayleigh--Taylor instability which is responsible for the occurrence of narrow upward jets are studied in the scope of the nonhydrostatic model with horizontally--nonuniform density and the Newtonian cooling. As analysis shows, the total hierarchy of instabilities in this model consists of three regimes -- collapse, algebraic instability, and inertial motion. Realization of these stages, mutual transitions and interference depend on a ratio between two characteristic time scales -- collapse time and cooling time.


title: Ab initio phonon coupling and optical response of hot electrons in plasmonic metals

abstract: Ultrafast laser measurements probe the non-equilibrium dynamics of excited electrons in metals with increasing temporal resolution. Electronic structure calculations can provide a detailed microscopic understanding of hot electron dynamics, but a parameter-free description of pump-probe measurements has not yet been possible, despite intensive research, because of the phenomenological treatment of electron-phonon interactions. We present ab initio predictions of the electron-temperature dependent heat capacities and electron-phonon coupling coefficients of plasmonic metals. We find substantial differences from free-electron and semi-empirical estimates, especially in noble metals above transient electron temperatures of 2000 K, because of the previously-neglected strong dependence of electron-phonon matrix elements on electron energy. We also present first-principles calculations of the electron-temperature dependent dielectric response of hot electrons in plasmonic metals, including direct interband and phonon-assisted intraband transitions, facilitating complete theoretical predictions of the time-resolved optical probe signatures in ultrafast laser experiments.


title: Generalized Kato Decomposition For Operator Matrices and SVEP

abstract: In this paper, we show that for a bounded linear operator $T$, the corresponding generalized Kato decomposition spectrum $\sigma_{gK}(T)$ satisfies the equality $\sigma_{gD}(T)=\sigma_{gK}(T)\cup (S(T)\cup S(T^*))$ where $\sigma_{gD} (T ) $ is the generalized Drazin spectrum of $T$ and $S(T )$ (resp., $S(T^*)$ is the set where T (resp., $T^*$) fails to have SVEP. As application, we give sufficient conditions which assure that the generalized Kato decomposition spectrum of an upper triangular operator matrices is the union of its diagonal entries spectra. Moreover, some applications are given.


title: Pucci eigenvalues on geodesic balls

abstract: We study the eigenvalue problem for the Riemannian Pucci operator on geodesic balls. We establish upper and lower bounds for the principal Pucci eigenvalues depending on the curvature, extending Cheng's eigenvalue comparison theorem for the Laplace-Beltrami operator. For manifolds with bounded sectional curvature, we prove Cheng's bounds hold for Pucci eigenvalues on geodesic balls of radius less than the injectivity radius. For manifolds with Ricci curvature bounded below, we prove Cheng's upper bound holds for Pucci eigenvalues on certain small geodesic balls. We also prove that the principal Pucci eigenvalues of an $O(n)$-invariant hypersurface immersed in $\mathbb{R}^{n+1}$ with one smooth boundary component are smaller than the eigenvalues of an $n$-dimensional Euclidean ball with the same boundary.


title: Ultrafast Depopulation of a Quantum Dot by LA-phonon-assisted Stimulated Emission

abstract: We demonstrate ultrafast incoherent depopulation of a quantum dot from above to below the transparency point using LA-phonon-assisted emission stimulated by a redshifted laser pulse. The QD is turned from a weakly vibronic system into a strongly vibronic one by laser driving which enables the phonon-assisted relaxation between the excitonic components of two dressed states. The depopulation is achieved within a laser pulse-width-limited time of 20 ps and exhibits a broad tuning range of a few meV. Our experimental results are well reproduced by path-integral calculations.


title: How to improve accuracy for DFA technique

abstract: This paper extends the existing literature on empirical estimation of the confidence intervals associated to the Detrended Fluctuation Analysis (DFA). We used Montecarlo simulation to evaluate the confidence intervals. Varying the parameters in DFA technique, we point out the relationship between those and the standard deviation of H. The parameters considered are the finite time length L, the number of divisors d used and the values of those. We found that all these parameters play a crucial role, determining the accuracy of the estimation of H.


title: Analytical solution for multisingular vortex Gaussian beams: The mathematical theory of scattering modes

abstract: We present a novel procedure to solve the Schr\"odinger equation, which in optics is the paraxial wave equation, with an initial multisingular vortex Gaussian beam. This initial condition has a number of singularities in a plane transversal to propagation embedded in a Gaussian beam. We use the scattering modes, which are solutions of the paraxial wave equation that can be combined straightforwardly to express the initial condition and therefore permit to solve the problem. To construct the scattering modes one needs to obtain a particular set of polynomials, which play an analogous role than Laguerre polynomials for Laguerre-Gaussian modes. We demonstrate here the recurrence relations needed to determine these polynomials. To stress the utility and strength of the method we solve first the problem of an initial Gaussian beam with two positive singularities and a negative one embedded in. We show that the solution permits one to obtain analytical expressions. These can used to obtain closed expressions for meaningful quantities, like the distance at which the positive and negative singularities merge, closing the loop of a vortex line. Furthermore, we present an example of calculation of an specific discrete-Gauss state, which is the solution of the diffraction of a Laguerre-Gauss state showing definite angular momentum (that is, a highly charged vortex) by a thin diffractive element showing certain discrete symmetry. We show that thereby this problem is solved in a much simpler way than using the previous procedure based in the integral Fresnel diffraction method.

categories: physics.optics cond-mat.quant-gas math-ph math.MP quant-ph


title: The ecological and evolutionary energetics of hunter-gatherer residential mobility

abstract: Residential mobility is deeply entangled with all aspects of hunter-gatherer life ways, and is therefore an issue of central importance in hunter-gatherer studies. Hunter-gatherers vary widely in annual rates of residential mobility, and understanding the sources of this variation has long been of interest to anthropologists and archaeologists. Since mobility is, to a large extent, driven by the need for a continuous supply of food, a natural framework for addressing this question is provided by the metabolic theory of ecology. This provides a powerful framework for formulating formal testable hypotheses concerning evolutionary and ecological constraints on the scale and variation of hunter-gatherer residential mobility. We evaluate these predictions using extant data and show strong support for the hypotheses. We show that the overall scale of hunter-gatherer residential mobility is predicted by average human body size, and the limited capacity of mobile hunter-gatherers to store energy internally. We then show that the majority of variation in residential mobility observed across cultures is predicted by energy availability in local ecosystems. Our results demonstrate that large-scale evolutionary and ecological processes, common to all plants and animals, constrain hunter-gatherers in predictable ways as they move through territories to effectively exploit resources over the course of a year. Moreover, our results extend the scope of the metabolic theory of ecology by showing how it successfully predicts variation in the behavioral ecology of populations within a species.


title: Spinor Lie derivatives and Fermion stress-energies

abstract: Stress-energies for Fermi fields are derived from the principle of general covariance. This is done by developing a notion of Lie derivatives of spinors along arbitrary vector fields. A substantial theory of such derivatives was first introduced by Kosmann; here I show how an apparent conflict in the literature on this is due to a difference in the definitions of spinors, and that tracking the Lie derivative of the Infeld-van der Waerden symbol, as well as the spinor fields under consideration, gives a fuller picture of the geometry and leads to the Fermion stress-energy. The differences in the definitions of spinors do not affect the results here, but could matter in certain quantum-gravity programs and for spinor transformations under discrete symmetries.


title: A Subluminal Relativistic Magnetohydrodynamics Scheme with ADER-WENO Predictor and Multidimensional Riemann Solver-Based Corrector

abstract: The relativistic magnetohydrodynamics (RMHD) set of equations has recently seen increased use in astrophysical computations. Even so, RMHD codes remain fragile. The reconstruction can sometimes yield superluminal velocities in certain parts of the mesh. In this paper we present a reconstruction strategy that overcomes this problem by making a single conservative to primitive transformation per cell followed by higher order WENO reconstruction on a carefully chosen set of primitives that guarantee subluminal reconstruction of the flow variables. For temporal evolution via a predictor step we also present second, third and fourth order accurate ADER methods that keep the velocity subluminal during the predictor step. The RMHD system also requires the magnetic field to be evolved in a divergence-free fashion. In the treatment of classical numerical MHD the analogous issue has seen much recent progress with the advent of multidimensional Riemann solvers. By developing multidimensional Riemann solvers for RMHD, we show that similar advances extend to RMHD. As a result, the face-centered magnetic fields can be evolved much more accurately using the edge-centered electric fields in the corrector step. Those edge-centered electric fields come from a multidimensional Riemann solver for RMHD which we present in this paper. In this paper we also develop several new test problems for RMHD. We show that RMHD vortices can be designed that propagate on the computational mesh as self-preserving structures. These RMHD vortex test problems provide a means to do truly multidimensional accuracy testing for RMHD codes. Several other stringent test problems are presented. We show the importance of resolution in certain test problems. Our tests include a demonstration that RMHD vortices are stable when they interact with shocks.

categories: astro-ph.HE astro-ph.IM gr-qc physics.comp-ph


title: Singular values for products of two coupled random matrices: hard edge phase transition

abstract: Consider the product $GX$ of two rectangular complex random matrices coupled by a constant matrix $\Omega$, where $G$ can be thought to be a Gaussian matrix and $X$ is a bi-invariant polynomial ensemble. We prove that the squared singular values form a biorthogonal ensemble in Borodin's sense, and further that for $X$ being Gaussian the correlation kernel can be expressed as a double contour integral. When all but finitely many eigenvalues of $\Omega^{} \Omega^{*}$ are equal, the corresponding correlation kernel is shown to admit a phase transition phenomenon at the hard edge in four different regimes as the coupling matrix changes. Specifically, the four limiting kernels in turn are the Meijer G-kernel for products of two independent Gaussian matrices, a new critical and interpolating kernel, the perturbed Bessel kernel and the finite coupled product kernel associated with $GX$. In the special case that $X$ is also a Gaussian matrix and $\Omega$ is scalar, such a product has been recently investigated by Akemann and Strahov. We also propose a Jacobi-type product and prove the same transition.


title: The Vertical Current Approximation Nonlinear Force-Free Field Code - Description, Performance Tests, and Measurements of Magnetic Energies Dissipated in Solar Flares

abstract: In this work we provide an updated description of the Vertical Current Approximation Nonlinear Force-Free Field (VCA-NLFFF) code, which is designed to measure the evolution of the potential, nonpotential, free energies, and the dissipated magnetic energies during solar flares. This code provides a complementary and alternative method to existing traditional NLFFF codes. The chief advantages of the VCA-NLFFF code over traditional NLFFF codes are the circumvention of the unrealistic assumption of a force-free photosphere in the magnetic field extrapolation method, the capability to minimize the misalignment angles between observed coronal loops (or chromospheric fibril structures) and theoretical model field lines, as well as computational speed. In performance tests of the VCA-NLFFF code, by comparing with the NLFFF code of Wiegelmann (2004), we find agreement in the potential, nonpotential, and free energy within a factor of about 1.3, but the Wiegelmann code yields in the average a factor of 2 lower flare energies. The VCA-NLFFF code is found to detect decreases in flare energies in most X, M, and C-class flares. The successful detection of energy decreases during a variety of flares with the VCA-NLFFF code indicates that current-driven twisting and untwisting of the magnetic field is an adequate model to quantify the storage of magnetic energies in active regions and their dissipation during flares. - The VCA-NLFFF code is also publicly available in the Solar SoftWare (SSW).


title: Pseudo semi B-Fredholm and Generalized Drazin invertible operators Through Localized SVEP

abstract: In this paper, we define and study the pseudo upper and lower semi B-Fredholm of bounded operators in a Banach space. In particular, we prove equality up to $S(T)$ between the left generalized Drazin spectrum and the pseudo upper semi B-Fredholm spectrum, $S(T )$ is the set where $T$ fails to have the SVEP. Also, we prove that both of the left and the right generalized Drazin operators are invariant under additive commuting power finite rank perturbations and some perturbations for the pseudo upper and lower semi B-Fredholm operators are given. As applications, we investigate some classes of operators as the supercyclic and multiplier operators. Furthermore, we investigate the left and the right generalized Drazin invertibility of upper triangular operator matrices by giving sufficient conditions which assure that the left and the right generalized Drazin spectrum or the pseudo upper and lower semi B-Fredholm of an upper triangular operator matrices is the union of its diagonal entries spectra.


title: Radiative Breaking of the Minimal Supersymmetric Left-Right Model

abstract: We study a variation to the SUSY Left-Right symmetric model based on the gauge group $SU(3)_c\times SU(2)_L\times SU(2)_R\times U(1)_{BL}$. Beyond the quark and lepton superfields we only introduce a second Higgs bidoublet to produce realistic fermion mass matrices. This model does not include any $SU(2)_R$ triplets. We calculate renormalization group evolutions of soft SUSY parameters at the one-loop level down to low energy. We find that an $SU(2)_R$ slepton doublet acquires a negative mass squared at low energies, so that the breaking of $SU(2)_R\times U(1)_{BL}\rightarrow U(1)_Y$ is realized by a non-zero vacuum expectation value of a right-handed sneutrino. Small neutrino masses are produced through neutrino mixings with gauginos. Mass limits on the $SU(2)_R\times U(1)_{BL}$ sector are obtained by direct search results at the LHC as well as lepton-gaugino mixing bounds from the LEP precision data.


title: The K2-ESPRINT Project IV: A Hot Jupiter in a Prograde Orbit with a Possible Stellar Companion

abstract: We report on the detection and early characterization of a hot Jupiter in a 3-day orbit around K2-34 (EPIC~212110888), a metal-rich F-type star located in the K2 Cycle 5 field. Our follow-up campaign involves precise radial velocity (RV) measurements and high-contrast imaging using multiple facilities. The absence of a bright nearby source in our high-contrast data suggests that the transit-like signals are not due to light variations from such a companion star. Our intensive RV measurements show that K2-34b has a mass of $1.773\pm0.086M_J$, confirming its status as a planet. We also detect the Rossiter-McLaughlin effect for K2-34b and show that the system has a good spin-orbit alignment ($\lambda=-1_{-9}^{+10}$ degrees). High-contrast images obtained by the HiCIAO camera on the Subaru 8.2-m telescope reveal a faint companion candidate ($\Delta m_H=6.19\pm 0.11$ mag) at a separation of $0\farcs36$. Follow-up observations are needed to confirm that the companion candidate is physically associated with K2-34. K2-34b appears to be an example of a typical "hot Jupiter," albeit one which can be precisely characterized using a combination of K2 photometry and ground-based follow-up.


title: Online Ski Rental for ON/OFF Scheduling of Energy Harvesting Base Stations

abstract: The co-existence of small cell base stations (SBSs) with conventional macrocell base station is a promising approach to boost the capacity and coverage of cellular networks. However, densifying the network with a viral deployment of SBSs can significantly increase energy consumption. To reduce the reliance on unsustainable energy sources, one can adopt self-powered SBSs that rely solely on energy harvesting. Due to the uncertainty of energy arrival and the finite capacity of energy storage systems, self-powered SBSs must smartly optimize their ON and OFF schedule. In this paper, the problem of ON/OFF scheduling of self-powered SBSs is studied, in the presence of energy harvesting uncertainty with the goal of minimizing the operational costs consisted of energy consumption and transmission delay of a network. For the original problem, we show an algorithm can solve the problem in the illustrative case. To reduce the complexity of the original problem, an approximation is proposed. To solve the approximated problem, a novel approach based on the ski rental framework, a powerful online optimization tool, is proposed. Using this approach, each SBS can effectively decide on its ON/OFF schedule autonomously, without any prior information on future energy arrivals. By using competitive analysis, a deterministic online algorithm (DOA) and a randomized online algorithm (ROA) are developed. ROA is shown to achieve the optimal competitive ratio in the approximation problem. Simulation results show that, compared to a baseline approach, the ROA can yield performance gains reaching up to 15.6% in terms of reduced total energy consumption of SBSs and up to 20.6% in terms of per-SBS network delay reduction. The results shed light on the fundamental aspects that impact the ON time of SBSs while demonstrating that the proposed ROA can reduce up to 69.9% the total cost compared to a baseline approach.


title: Optical polarization and intervalley scattering in single layers of MoS2 and MoSe2

abstract: Single layers of MoS2 and MoSe2 were optically pumped with circularly polarized light and an appreciable polarization was initialized as the pump energy was varied. The circular polarization of the emitted photoluminescence was monitored as function of the difference between the excitation energy and the A-exciton emission at the K-point of the Brillouin zone. Our results show a threshold of twice the LA phonon energy, specific to the material, above which phonon-assisted intervalley scattering causes depolarization. In both materials this lead to almost complete depolarization within ~100 meV above the threshold energy. We identify the extra kinetic energy of the exciton (independent of whether it is neutral or charged) as the key parameter for presenting a unifying picture of the depolarization process.


title: On the convergence of monotone Hurwitz generating functions

abstract: Monotone Hurwitz numbers were introduced by the authors as a combinatorially natural desymmetrization of the Hurwitz numbers studied in enumerative algebraic geometry. Over the course of several papers, we developed the structural theory of monotone Hurwitz numbers and demonstrated that it is in many ways parallel to that of their classical counterparts. In this note, we identify an important difference between the monotone and classical worlds: fixed-genus generating functions for monotone double Hurwitz numbers are absolutely summable, whereas those for classical double Hurwitz numbers are not. This property is crucial for applications of monotone Hurwitz theory in analysis. We quantify the growth rate of monotone Hurwitz numbers in fixed genus by giving universal upper and lower bounds on the radii of convergence of their generating functions.


title: Directed percolation process in the presence of velocity fluctuations: Effect of compressibility and finite correlation time

abstract: The direct bond percolation process (Gribov process) is studied in the presence of random velocity fluctuations generated by the Gaussian self-similar ensemble with finite correlation time. We employ the renormalization group in order to analyze a combined effect of the compressibility and finite correlation time on the long-time behavior of the phase transition between an active and an absorbing state. The renormalization procedure is performed to the one-loop order. Stable fixed points of the renormalization group and their regions of stability are calculated in the one-loop approximation within the three-parameter $(\eps,y,\eta)$-expansion. Different regimes corresponding to the rapid-change limit and frozen velocity field are discussed, and their fixed points' structure is determined in numerical fashion.


title: Burstiness and fractional diffusion on complex networks

abstract: Many dynamical processes on real world networks display complex temporal patterns as, for instance, a fat-tailed distribution of inter-events times, leading to heterogeneous waiting times between events. In this work, we focus on distributions whose average inter-event time diverges, and study its impact on the dynamics of random walkers on networks. The process can naturally be described, in the long time limit, in terms of Riemann-Liouville fractional derivatives. We show that all the dynamical modes possess, in the asymptotic regime, the same power law relaxation, which implies that the dynamics does not exhibit time-scale separation between modes, and that no mode can be neglected versus another one, even for long times. Our results are then confirmed by numerical simulations.


title: Universal features of cluster numbers in percolation

abstract: The number of clusters per site $n(p)$ in percolation at the critical point $p = p_c$ is not itself a universal quantity---it depends upon the lattice and percolation type (site or bond). However, many of its properties, including finite-size corrections, scaling behavior with $p$, and amplitude ratios, show various degrees of universal behavior. Some of these are universal in the sense that the behavior depends upon the shape of the system, but not lattice type. Here, we elucidate the various levels of universality for elements of $n(p)$ both theoretically and by carrying out extensive studies on several two- and three-dimensional systems, by high-order series analysis, Monte-Carlo simulation, and exact enumeration. We find many new results, including precise values for $n(p_c)$ for several systems, a clear demonstration of the singularity in $n''(p)$, and metric scale factors. We make use of the matching polynomial of Sykes and Essam to find exact relations between properties for lattices and matching lattices. We propose a criterion for an absolute metric factor $b$ based upon the singular behavior of the scaling function, rather than a relative definition of the metric that has previously been used.


title: Frobenius semisimplicity for convolution morphisms

abstract: This article concerns properties of mixed $\ell$-adic complexes on varieties over finite fields, related to the action of the Frobenius automorphism. We establish a fiberwise criterion for the semisimplicity and Frobenius semisimplicity of the direct image complex under a proper morphism of varieties over a finite field. We conjecture that the direct image of the intersection complex on the domain is always semisimple and Frobenius semisimple; this conjecture would imply that a strong form of the decomposition theorem of Beilinson-Bernstein-Deligne-Gabber is valid over finite fields. We prove our conjecture for (generalized) convolution morphisms associated with partial affine flag varieties for split connected reductive groups over finite fields, and we prove allied Frobenius semisimplicity results for the intersection cohomology groups of twisted products of Schubert varieties. We offer two proofs for these results: one is based on the paving by affine spaces of the fibers of certain convolution morphisms, the other involves a new schematic theory of "big cells" adapted to partial affine flag varieties, and combines Deligne's theory of weights with a suitable contracting $\mathbb G_m$-action on those big cells. Both proofs rely on our general result that the intersection complex of the image of a proper map of varieties over a finite field is a direct summand of the direct image of the intersection complex of the domain. With suitable reformulations, the main results are valid over any algebraically closed ground field.


title: Autonomous Agent Behaviour Modelled in PRISM -- A Case Study

abstract: Formal verification of agents representing robot behaviour is a growing area due to the demand that autonomous systems have to be proven safe. In this paper we present an abstract definition of autonomy which can be used to model autonomous scenarios and propose the use of small-scale simulation models representing abstract actions to infer quantitative data. To demonstrate the applicability of the approach we build and verify a model of an unmanned aerial vehicle (UAV) in an exemplary autonomous scenario, utilising this approach.


title: Cosmic ray air showers from sphalerons

abstract: The discovery of the Higgs boson marks a key ingredient to establish the electroweak structure of the Standard Model. Its non-abelian gauge structure gives rise to, yet unobserved, non-perturbative baryon and lepton number violating processes. We propose to use cosmic ray air showers, as measured at the Pierre Auger Observatory, to set a limit on the hadronic production cross section of sphalerons. We identify several observables to discriminate between sphaleron and QCD induced air showers.


title: Hybrid Beamforming for Massive MIMO Backhaul (Working Title)

abstract: The uplink where both the transmitter and receiver can use a large antenna array is considered. This is proposed as a method of antenna offloading and connecting small cell access points (SCAP) in a Two-Tier cellular network. Due to having a limited number of RF-chains, hybrid beamformers are designed where phase-only processing is done at the RF-band, followed by digital processing at the baseband. The proposed receiver is a row combiner that clusters sufficiently correlated antenna elements, and its' performance is compared against random projection via a Discrete Fourier Transform (DFT) matrix. The analogue to the row combiner is a column spreader, which is dependent on the transmit correlation, and repeats the transmitted signal over antenna elements that are correlated. A key benefit of this approach is to reduce the number of phase shifters used, while outperforming the DFT scheme. When only the transmitter has correlation and is RF-chain limited, the baseband precoding vectors are shown to be the eigenvectors of the effective transmit correlation matrix. Depending on the channel correlation, this matrix can be approximated to have a tridiagonal Toeplitz structure with the proposed column spreader (CS). The resulting eigenvalues have a closed form solution which allows us to characterize the sum rate of the system. Most interestingly, the associated eigenvectors do not require knowledge of the effective transmit correlation matrix to be calculated using an Eigenvalue Decomposition (EVD) method.


title: Analysis of the rational Krylov subspace projection method for large-scale algebraic Riccati equations

abstract: In the numerical solution of the algebraic Riccati equation $A^* X + X A - X BB^* X + C^* C =0$, where $A$ is large, sparse and stable, and $B$, $C$ have low rank, projection methods have recently emerged as a possible alternative to the more established Newton-Kleinman iteration. In spite of convincing numerical experiments, a systematic matrix analysis of this class of methods is still lacking. We derive new relations for the approximate solution, the residual and the error matrices, giving new insights into the role of the matrix $A-BB^*X$ and of its approximations in the numerical procedure. The new results provide theoretical ground for recently proposed modifications of projection methods onto rational Krylov subspaces.


title: A Quantitative Analysis of Localized Robustness of MYCN in Neuroblastoma

abstract: The amplification of the gene MYCN (V-myc myelocytomatosis viral-valeted oncogene, neuroblastoma derived) has been a well-documented indicator for poor prognosis in neuroblastoma, a childhood cancer. Unfortunately, there has been limited success in understanding MYCN functionality in the landscape of neuroblastoma and more importantly, given that MYCN has been deemed undruggable, the need to potentially illuminate key opportunities that indirectly target MYCN is of great interest. To this end, this work employs an emerging quantitative technique from network science, namely network curvature, to quantify the biological robustness of MYCN and its surrounding neighborhood. In particular, when amplified in Stage IV cancer, MYCN exhibits higher curvature (more robust) than those samples with under expressed MYCN levels. When examining the surrounding neighborhood, the above argument still holds for network curvature, but is lost when only analyzing differential expression - a common technique amongst oncologists and computational/molecular biologists. This finding points to the problem (and possible solution) of drug targeting in the context of complexity and indirect cell signaling affects that have often been obfuscated through traditional techniques.


title: Fast Computation of Minimal Interpolation Bases in Popov Form for Arbitrary Shifts

abstract: We compute minimal bases of solutions for a general interpolation problem, which encompasses Hermite-Pad\'e approximation and constrained multivariate interpolation, and has applications in coding theory and security. This problem asks to find univariate polynomial relations between $m$ vectors of size $\sigma$; these relations should have small degree with respect to an input degree shift. For an arbitrary shift, we propose an algorithm for the computation of an interpolation basis in shifted Popov normal form with a cost of $\mathcal{O}\tilde{~}(m^{\omega-1} \sigma)$ field operations, where $\omega$ is the exponent of matrix multiplication and the notation $\mathcal{O}\tilde{~}(\cdot)$ indicates that logarithmic terms are omitted. Earlier works, in the case of Hermite-Pad\'e approximation and in the general interpolation case, compute non-normalized bases. Since for arbitrary shifts such bases may have size $\Theta(m^2 \sigma)$, the cost bound $\mathcal{O}\tilde{~}(m^{\omega-1} \sigma)$ was feasible only with restrictive assumptions on the shift that ensure small output sizes. The question of handling arbitrary shifts with the same complexity bound was left open. To obtain the target cost for any shift, we strengthen the properties of the output bases, and of those obtained during the course of the algorithm: all the bases are computed in shifted Popov form, whose size is always $\mathcal{O}(m \sigma)$. Then, we design a divide-and-conquer scheme. We recursively reduce the initial interpolation problem to sub-problems with more convenient shifts by first computing information on the degrees of the intermediate bases.


title: A survey of the cold molecular gas in gravitationally lensed star-forming galaxies at z>2

abstract: Using the Australia Telescope Compact Array (ATCA), we conducted a survey of CO J=1-0 and J=2-1 line emission towards strongly lensed high-redshift dusty star forming galaxies (DSFGs) previously discovered with the South Pole Telescope (SPT). Our sample comprises 17 sources that had CO-based spectroscopic redshifts obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) and the Atacama Pathfinder Experiment (APEX). We detect all sources with known redshifts in either CO J=1-0 or J=2-1. Twelve sources are detected in the 7-mm continuum. The derived CO luminosities imply gas masses in the range (0.5-11)x10^{10} M_sun and gas depletion timescales <200 Myr, using a CO to gas mass conversion factor alpha_CO=0.8 M_sun (K km/s pc^2)^{-1}. Combining the CO luminosities and dust masses, along with a fixed gas-to-dust ratio, we derive alpha_CO factors in the range 0.4-1.8, similar to what is found in other starbursting systems. We find small scatter in alpha_CO values within the sample, even though inherent variations in the spatial distribution of dust and gas in individual cases could bias the dust-based alpha_CO estimates. We find that lensing magnification factors based on the CO linewidth to luminosity relation (mu_CO) are highly unreliable, but particularly when mu<5. Finally, comparison of the gas and dynamical masses suggest that the average molecular gas fraction stays relatively constant at z=2-5 in the SPT DSFG sample.


title: Distributed Subgradient Projection Algorithm over Directed Graphs

abstract: We propose a distributed algorithm, termed the Directed-Distributed Projected Subgradient (D-DPS), to solve a constrained optimization problem over a multi-agent network, where the goal of agents is to collectively minimize the sum of locally known convex functions. Each agent in the network owns only its local objective function, constrained to a commonly known convex set. We focus on the circumstance when communications between agents are described by a directed network. The D-DPS augments an additional variable for each agent, to overcome the asymmetry caused by the directed communication network. The convergence analysis shows that D-DPS converges at a rate of $O(\frac{\ln k}{\sqrt{k}})$, where k is the number of iterations.


title: A representation stability theorem for VI-modules

abstract: A VI-module gives rise to a sequence of representations of the finite general linear groups. We prove that the sequence obtained from any finitely generated VI-module over an algebraically closed field of characteristic zero is representation stable.


title: Ordinary X-rays from Three Extraordinary Millisecond Pulsars: XMM-Newton Observations of PSRs J0337+1715, J0636+5129, and J0645+5158

abstract: We present the first X-ray observations of three recently discovered millisecond pulsars (MSPs) with interesting characteristics: PSR J0337+1715, PSR J0636+5129, and PSR J0645+5158. PSR J0337+1715 is a fast-spinning, bright, and so-far unique MSP in a hierarchical triple system with two white dwarf (WD) companions. PSR J0636+5129 is a MSP in a very tight 96-min orbit with a low-mass, 8 $M_J$ companion. PSR J0645+5158 is a nearby, isolated MSP with a very small duty cycle (1-2%), which has led to its inclusion in high-precision pulsar timing programs. Using data from XMM-Newton, we have analyzed X-ray spectroscopy for these three objects, as well as optical/ultraviolet photometry for PSR J0337+1715. The X-ray data for each are largely consistent with expectations for most MSPs with regards to the ratios of thermal and non-thermal emission. We discuss the implications of these data on the pulsar population, and prospects for future observations of these pulsars.


title: Aspects of electron-phonon interactions with strong forward scattering in FeSe Thin Films on SrTiO$_3$ substrates

abstract: Mono- and multilayer FeSe thin films grown on SrTiO$_\mathrm{3}$ and BiTiO$_\mathrm{3}$ substrates exhibit a greatly enhanced superconductivity over that found in bulk FeSe. A number of proposals have been advanced for the mechanism of this enhancement. One possibility is the introduction of a cross-interface electron-phonon ($e$-$ph$) interaction between the FeSe electrons and oxygen phonons in the substrates that is peaked in the forward scattering (small ${\bf q}$) direction due to the two-dimensional nature of the interface system. Motivated by this, we explore the consequences of such an interaction on the superconducting state and electronic structure of a two-dimensional system using Migdal-Eliashberg theory. This interaction produces not only deviations from the expectations of conventional phonon-mediated pairing but also replica structures in the spectral function and density of states, as probed by angle-resolved photoemission spectroscopy, scanning tunneling microscopy/spectroscopy, and quasi-particle interference imaging. We also discuss the applicability of Migdal-Eliashberg theory for a situation where the \ep interaction is peaked at small momentum transfer and in the FeSe/STO system.


title: Low Temperature Asymptotics in Spherical Mean Field Spin Glasses

abstract: In this paper, we study the low temperature limit of the spherical Crisanti-Sommers variational problem. We identify the $\Gamma$-limit of the Crisanti-Sommers functionals, thereby establishing a rigorous variational problem for the ground state energy of spherical mixed $p$-spin glasses. As an application, we compute moderate deviations of the corresponding minimizers in the low temperature limit. In particular, for a large class of models this yields moderate deviations for the overlap distribution. We then analyze the ground state energy problem. We show that this variational problem is dual to an obstacle-type problem. This duality is at the heart of our analysis. We present the regularity theory of the optimizers of the primal and dual problems. This culminates in a simple method for constructing a finite dimensional space in which these optimizers live for any model. As a consequence of these results, we unify independent predictions of Crisanti-Leuzzi and Auffinger-Ben Arous regarding the 1RSB phase in this limit. We find that the "positive replicon eigenvalue" and "pure-like" conditions are together necessary for optimality, but that neither are themselves sufficient, answering a question of Auffinger and Ben Arous in the negative. We end by proving that these conditions completely characterize the 1RSB phase in $2+p$-spin models.


title: Solubility of Rock in Steam Atmospheres of Planets

abstract: Extensive experimental studies show all major rocky elements (Si, Mg, Fe, Ni, Ca, Al, Na, K) dissolve in steam to a greater or lesser extent. We use these results to compute chemical equilibrium abundances of rocky element-bearing gases in steam atmospheres equilibrated with silicate magma oceans. Rocky elements partition into steam atmospheres as volatile hydroxide gases and via reaction with HF or HCl as volatile halide gases in much larger amounts than expected from their vapor pressures over volatile-free solid or molten rock at the same temperature. We compute the extent of fractional vaporization by defining gas to magma partition coefficients and show Earth's sub-solar Si to Mg bulk elemental ratio may be due to loss of a primordial steam atmosphere. We conclude hot rocky exoplanets that are undergoing or have undergone escape of steam atmospheres may experience fractional vaporization and loss of Si, Mg, Fe, Ni, Ca, Al, Na, and K. This loss may modify their bulk composition, density, heat balance, and internal structure.


title: Low Starting Electron Beam Current in Degenerate Band Edge Oscillators

abstract: We propose a new principle of operation in vacuum electron-beam-based oscillators that leads to a low beam current for starting oscillations. The principle is based on super synchronous operation of an electron beam interacting with four degenerate electromagnetic modes in a slow-wave structure (SWS). The four mode super synchronous regime is associated with a very special degeneracy condition in the dispersion diagram of a cold periodic SWS called degenerate band edge (DBE). This regime features a giant group delay in the finitelength SWS and low starting-oscillation beam current. The starting beam current is at least an order of magnitude smaller compared to a conventional backward wave oscillator (BWO) of the same length. As a representative example we consider a SWS conceived by a periodically-loaded metallic waveguide supporting a DBE, and investigate starting-oscillation conditions using Pierce theory generalized to coupled transmission lines (CTL). The proposed super synchronism regime can be straightforwardly adapted to waveguide geometries others than the periodically-loaded waveguide considered here since DBE is a general property that can be realized in a variety of structures.


title: On boundary behavior of mappings in terms of prime ends

abstract: A boundary behavior of mappings, which are closely related with Sobolev and Orlicz--Sobolev classes in the plane and in the space, is investigated. There are obtained theorems on boundary behavior of classes mentioned above.


title: Detection and localization of change points in temporal networks with the aid of stochastic block models

abstract: A framework based on generalized hierarchical random graphs (GHRGs) for the detection of change points in the structure of temporal networks has recently been developed by Peel and Clauset [1]. We build on this methodology and extend it to also include the versatile stochastic block models (SBMs) as a parametric family for reconstructing the empirical networks. We use five different techniques for change point detection on prototypical temporal networks, including empirical and synthetic ones. We find that none of the considered methods can consistently outperform the others when it comes to detecting and locating the expected change points in empirical temporal networks. With respect to the precision and the recall of the results of the change points, we find that the method based on a degree-corrected SBM has better recall properties than other dedicated methods, especially for sparse networks and smaller sliding time window widths.

categories: cs.SI physics.data-an physics.soc-ph stat.ME


title: Haag's theorem in renormalised quantum field theories

abstract: We review a package of no-go results in axiomatic quantum field theory with Haag's theorem at its centre. Since the concept of operator-valued distributions in this framework comes very close to what we believe canonical quantum fields are about, these results are of consequence to quantum field theory: they suggest the seeming absurdity that this highly victorious theory is incapable of describing interactions. We single out unitarity of the interaction picture's intertwiner as the most salient provision of Haag's theorem and critique canonical perturbation theory to argue that renormalisation bypasses Haag's theorem by violating this very assumption.


title: A Coherent Polariton Laser

abstract: The semiconductor polariton laser promises a new source of coherent light, which, compared to conventional semiconductor photon lasers, has input-energy threshold orders of magnitude lower. However, intensity stability, a defining feature of a coherent state, has remained poor. Intensity noise at many times of the shot-noise of a coherent state has persisted, which has been attributed to multiple mechanisms that are difficult to separate in conventional polariton systems. The large intensity noise in turn limited the phase coherence. These limit the capability of the polariton laser as a source of coherence light. Here, we demonstrate a polariton laser with shot-noise limited intensity stability, as expected of a fully coherent state. This is achieved by using an optical cavity with high mode selectivity to enforce single-mode lasing, suppress condensate depletion, and establish gain saturation. The absence of spurious intensity fluctuations moreover enabled measurement of a transition from exponential to Gaussian decay of the phase coherence of the polariton laser. It suggests large self-interaction energies in the polariton condensate, exceeding the laser bandwidth. Such strong interactions are unique to matter-wave laser and important for nonlinear polariton devices. The results will guide future development of polariton lasers and nonlinear polariton devices.


title: Analytic torsion, dynamical zeta functions, and the Fried conjecture

abstract: We prove the equality of the analytic torsion and the value at zero of a Ruelle dynamical zeta function associated with an acyclic unitarily flat vector bundle on a closed locally symmetric reductive manifold. This solves a conjecture of Fried. This article should be read in conjunction with an earlier paper by Moscovici and Stanton.


title: Long-term behaviour in a chemotaxis-fluid system with logistic source

abstract: We consider the coupled chemotaxis Navier-Stokes model with logistic source terms \[ n_t + u\cdot \nabla n = \Delta n - \chi \nabla \cdot (n \nabla c) + \kappa n - \mu n^2\] \[ c_t + u\cdot \nabla c = \Delta c - nc\] \[ u_t + (u\cdot \nabla)u = \Delta u +\nabla P + n\nabla \Phi + f, \quad\qquad \nabla \cdot u=0 \] in a bounded, smooth domain $\Omega\subset \mathbb{R}^3$ under homogeneous Neumann boundary conditions for $n$ and $c$ and homogeneous Dirichlet boundary conditions for $u$ and with given functions $f\in L^\infty(\Omega\times(0,\infty))$ satisfying certain decay conditions and $\Phi\in C^{1+\beta}(\bar\Omega)$ for some $\beta\in(0,1)$. We construct weak solutions and prove that after some waiting time they become smooth and finally converge to the semi-trivial steady state $(\frac{\kappa}{\mu},0,0)$. Keywords: chemotaxis, Navier-Stokes, logistic source, boundedness, large-time behaviour


title: Shintani zeta functions and a refinement of Gross's leading term conjecture

abstract: We introduce the notion of Shintani data, which axiomatizes algebraic aspects of Shintani zeta functions. We develop the general theory of Shintani data, and show that the order of vanishing part of Gross's conjecture follows from the existence of a Shintani datum. Recently, Dasgupta and Spiess proved the order of vanishing part of Gross's conjecture under certain conditions. We give an alternative proof of their result by constructing a certain Shintani datum. We also propose a refinement of Gross's leading term conjecture by using the theory of Shintani data. Out conjecture gives a conjectural construction of localized Rubin-Stark elements which can be regarded as a higher rank generalization of the conjectural construction of Gross-Stark units due to Dasgupta and Dasgupta-Spiess.


title: Emergence of core-peripheries in networks

abstract: A number of important transport networks, such as the airline and trade networks of the world, exhibit a characteristic core-periphery structure, wherein a few nodes are highly interconnected and the rest of the network frays into a tree. Mechanisms underlying the emergence of core-peripheries, however, remain elusive. Here, we demonstrate that a simple pruning process based on removal of underutilized links and redistribution of loads can lead to the emergence of core-peripheries. Links are assumed beneficial if they either carry a sufficiently large load or are essential for global connectivity. This incentivized redistribution process is controlled by a single parameter which balances connectivity and profit. The obtained networks exhibit a highly resilient and connected core with a frayed periphery. The balanced network shows a higher resilience than the World Airline Network or the World Trade Network, revealing a pathway towards robust structural features through pruning.


title: On the Structure and Efficient Computation of IsoRank Node Similarities

abstract: The alignment of protein-protein interaction (PPI) networks has many applications, such as the detection of conserved biological network motifs, the prediction of protein interactions, and the reconstruction of phylogenetic trees [1, 2, 3]. IsoRank is one of the first global network alignment algorithms [4, 5, 6], where the goal is to match all (or most) of the nodes of two PPI networks. The IsoRank algorithm first computes a pairwise node similarity metric, and then generates a matching between the two node sets based on this metric. The metric is a convex combination of a structural similarity score (with weight $ \alpha $) and an extraneous amino-acid sequence similarity score for two proteins (with weight $ 1 - \alpha $). In this short paper, we make two contributions. First, we show that when IsoRank similarity depends only on network structure ($\alpha = 1$), the similarity of two nodes is only a function of their degrees. In other words, IsoRank similarity is invariant to any network rewiring that does not affect the node degrees. This result suggests a reason for the poor performance of IsoRank in structure-only ($ \alpha = 1 $) alignment. Second, using ideas from [7, 8], we develop an approximation algorithm that outperforms IsoRank (including recent versions with better scaling, e.g., [9]) by several orders of magnitude in time and memory complexity, despite only a negligible loss in precision.


title: A non-linear subdiffusion model for a cell-cell adhesion in chemotaxis

abstract: The purpose of this work is to propose a non-Markovian and nonlinear model of subdiffusive transport that involves adhesion affects the cells escape rates form position x, with chemotaxis. This leads the escape rates to be dependent on the particles density at the neighbours as well as the chemotactic gradient. We systematically derive subdiffusive fractional master equation, then we consider the diffusive limit of the fractional master equation. We finally solve the resulted fractional subdiffusive master equation stationery and analyse the role of adhesion in the resulted stationary density.


title: Facet-dependent giant spin orbit torque in single crystalline antiferromagnetic Ir-Mn / ferromagnetic permalloy bilayers

abstract: There has been considerable interest in spin-orbit torques for the purpose of manipulating the magnetization of ferromagnetic (FM) films or nano-elements for spintronic technologies. Spin-orbit torques are derived from spin currents created from charge currents in materials with significant spin-orbit coupling that diffuse into an adjacent FM material. There have been intensive efforts to search for candidate materials that exhibit large spin Hall angles, i.e. efficient charge to spin current conversion. Here we report, using spin torque ferromagnetic resonance, the observation of a giant spin Hall angle of up to ~0.35 in (100) oriented single crystalline antiferromagnetic (AF) IrMn3 thin films, coupled to ferromagnetic permalloy layers, and a spin Hall angle that is about three times smaller in (111) oriented films. For the (100) oriented samples we show that the magnitude of the spin Hall angle can be significantly changed by manipulating the populations of the various AF domains through field annealing. Using ab-initio calculations we show that the triangular AF structure of IrMn3 gives rise to a substantial intrinsic spin Hall conductivity that is three times larger for the (100) than for the (111) orientations, consistent with our experimental findings.


title: Sterile Neutrino Fits to Short Baseline Data

abstract: Neutrino oscillation models involving extra mass eigenstates beyond the standard three (3+N) are fit to global short baseline experimental data. We find that 3+1 has a best fit of Delta m^2_41 = 1.75 eV^2 with a Delta chi^2 [null-min] (dof) of 52.34 (3). The 3+2 fit has a Delta chi^2 [null-min] (dof) of 56.99 (7). Bayesian credible intervals are shown for the first time for a 3+1 model. These are found to be in agreement with frequentist intervals. The results of these new fits favor a higher Delta m^2 value than previous studies, which may have an impact on future sterile neutrino searches such as the Fermilab SBN program.


title: Rationality for subclasses of 321-avoiding permutations

abstract: We prove that every proper subclass of the 321-avoiding permutations that is defined either by only finitely many additional restrictions or is well quasi-ordered has a rational generating function. To do so we show that any such class is in bijective correspondence with a regular language. The proof makes significant use of formal languages and of a host of encodings, including a new mapping called the panel encoding that maps languages over the infinite alphabet of positive integers avoiding certain subwords to languages over finite alphabets.


title: Neutronization During Carbon Simmering In Type Ia Supernova Progenitors

abstract: When a Type Ia supernova (SN Ia) progenitor first ignites carbon in its core, it undergoes ${\sim} \,10^{3}-10^{4} \,$yr of convective burning prior to the onset of thermonuclear runaway. This carbon simmering phase is important for setting the thermal profile and composition of the white dwarf. Using the \texttt{MESA} stellar evolution code, we follow this convective burning and examine the production of neutron-rich isotopes. The neutron content of the SN fuel has important consequences for the ensuing nucleosynthesis, and, in particular, for the production of secondary Fe-peak nuclei like Mn and stable Ni. These elements have been observed in the X-ray spectra of SN remnants like Tycho, Kepler, and 3C 397, and their yields can provide valuable insights into the physics of SNe Ia and the properties of their progenitors. We find that weak reactions during simmering can at most generate a neutron excess of ${\approx} \, 3 \times 10^{-4}$. This is ${\approx} \, 8 \times 10^{-4}$ lower than that found in previous studies that do not take the full density and temperature profile of the simmering region into account. Our results imply that the progenitor metallicity is the main contributor to the neutron excess in SN Ia fuel for $Z \gtrsim 1/3 \, Z_{\odot}$. Alternatively, at lower metallicities, this neutron excess provides a floor that should be present in any centrally-ignited SN~Ia scenario.


title: LSS constraints with controlled theoretical uncertainties

abstract: Forecasts and analyses of cosmological observations often rely on the assumption of a perfect theoretical model over a defined range of scales. We explore how model uncertainties and nuisance parameters in perturbative models of the matter and galaxy spectra affect constraints on neutrino mass and primordial non-Gaussianities. We provide a consistent treatment of theoretical errors and argue that their inclusion is a necessary step to obtain realistic cosmological constraints. We find that galaxy surveys up to high redshifts will allow a detection of the minimal neutrino mass and local non-Gaussianity of order unity, but improving the constraints on equilateral non-Gaussianity beyond the CMB limits will be challenging. We argue that similar considerations apply to analyses where theoretical models are based on simulations.


title: Residual-based a posteriori error estimation for the Maxwell's eigenvalue problem

abstract: We present an a posteriori estimator of the error in the L^2-norm for the numerical approximation of the Maxwell's eigenvalue problem by means of N\'ed\'elec finite elements. Our analysis is based on a Helmholtz decomposition of the error and on a superconvergence result between the L^2-orthogonal projection of the exact eigenfunction onto the curl of the N\'ed\'elec finite element space and the eigenfunction approximation. Reliability of the a posteriori error estimator is proved up to higher order terms and local efficiency of the error indicators is shown by using a standard bubble functions technique. The behavior of the a posteriori error estimator is illustrated on a numerical test.


title: Instability of nanometric fluid films on a thermally conductive substrate

abstract: We consider thin fluid films placed on thermally conductive substrates and exposed to time-dependent spatially uniform heat source. The evolution of the films is considered within the long-wave framework in the regime such that both fluid/substrate interaction, modeled via disjoining pressure, and Marangoni forces, are relevant. We analyze the problem by the means of linear stability analysis as well as by time-dependent nonlinear simulations. The main finding is that when self-consistent computation of the temperature field is performed, a complex interplay of different instability mechanisms results. This includes either monotonous or oscillatory dynamics of the free surface. In particular, we find that the oscillatory behavior is absent if the film temperature is assumed to be slaved to the current value of the film thickness. The results are discussed within the context of liquid metal films, but are of relevance to dynamics of any thin film involving variable temperature of the free surface, such that the temperature and the film interface itself evolve on comparable time scales.

categories: cond-mat.soft math-ph math.MP physics.flu-dyn


title: Approximate tangents, harmonic measure, and domains with rectifiable boundaries

abstract: Let $\Omega \subset \mathbb{R}^{n+1}$, $n \geq 1$, be an open and connected set. Set $\mathcal{T}_n$ to be the set of points $\xi \in \partial \Omega$ so that there exists an approximate tangent $n$-plane for $\partial\Omega$ at $\xi$ and $\partial\Omega$ satisfies the weak lower Ahlfors-David $n$-regularity condition at $\xi$. We first show that $\mathcal{T}_n$ can be covered by a countable union of boundaries of bounded Lipschitz domains. Then, letting $\partial^\star \Omega$ be a subset of $\mathcal{T}_n$ where $\Omega$ satisfies an appropriate thickness condition, we prove that $\partial^\star \Omega$ can be covered by a countable union of boundaries of bounded Lipschitz domains contained in $\Omega$. As a corollary we obtain that if $\Omega$ has locally finite perimeter, $\partial\Omega$ is weakly lower Ahlfors-David $n$-regular, and the measure-theoretic boundary coincides with the topological boundary of $\Omega$ up to a set of $\mathcal{H}^n$-measure zero, then $\partial \Omega$ can be covered, up to a set of $\mathcal{H}^n$-measure zero, by a countable union of boundaries of bounded Lipschitz domains that are contained in $\Omega$. This implies that in such domains, $\mathcal{H}^n|_{\partial\Omega}$ is absolutely continuous with respect to harmonic measure.


title: Ensemble Toolkit: Scalable and Flexible Execution of Ensembles of Tasks

abstract: There are many science applications that require scalable task-level parallelism and support for flexible execution and coupling of ensembles of simulations. Most high-performance system software and middleware, however, are designed to support the execution and optimization of single tasks. Motivated by the missing capabilities of these computing systems and the increasing importance of task-level parallelism, we introduce the Ensemble toolkit which has the following application development features: (i) abstractions that enable the expression of ensembles as primary entities, and (ii) support for ensemble-based execution patterns that capture the majority of application scenarios. Ensemble toolkit uses a scalable pilot-based runtime system that decouples workload execution and resource management details from the expression of the application, and enables the efficient and dynamic execution of ensembles on heterogeneous computing resources. We investigate three execution patterns and characterize the scalability and overhead of Ensemble toolkit for these patterns. We investigate scaling properties for up to O(1000) concurrent ensembles and O(1000) cores and find linear weak and strong scaling behaviour.


title: Segmented strings and the McMillan map

abstract: We present new exact solutions describing motions of closed segmented strings in $AdS_3$ in terms of elliptic functions. The existence of analytic expressions is due to the integrability of the classical equations of motion, which in our examples reduce to instances of the McMillan map. We also obtain a discrete evolution rule for the motion in $AdS_3$ of arbitrary bound states of fundamental strings and D1-branes in the test approximation.


title: Lepton Flavor Violating Decays of Neutral Higgses in Extended Mirror Fermion Model

abstract: We perform the one-loop induced charged lepton flavor violating decays of the neutral Higgses in an extended mirror fermion model with non-sterile electroweak-scale right-handed neutrinos and a horizontal $A_4$ symmetry in the lepton sector. We demonstrate that for the 125 GeV scalar $h$ there is tension between the recent LHC result ${\cal B}(h \to \tau \mu) \sim $ 1% and the stringent limits on the rare processes $\mu \to e \gamma$ and $\tau \to (\mu$ or $e) \gamma$ from low energy experiments.