title: Attraction-repulsion transition in the interaction of adatoms and vacancies in graphene

abstract: The interaction of two resonant impurities in graphene has been predicted to have a long-range character with weaker repulsion when the two adatoms reside on the same sublattice and stronger attraction when they are on different sublattices. We reveal that this attraction results from a single energy level. This opens up a possibility of controlling the sign of the impurity interaction via the adjustment of the chemical potential. For many randomly distributed impurities (adatoms or vacancies) this may offer a way to achieve a controlled transition from aggregation to dispersion.


title: Sensitivity analysis for multidimensional and functional outputs

abstract: Let $X:=(X_1, \ldots, X_p)$ be random objects (the inputs), defined on some probability space $(\Omega,{\mathcal{F}}, \mathbb P)$ and valued in some measurable space $E=E_1\times\ldots \times E_p$. Further, let $Y:=Y = f(X_1, \ldots, X_p)$ be the output. Here, $f$ is a measurable function from $E$ to some Hilbert space $\mathbb{H}$ ($\mathbb{H}$ could be either of finite or infinite dimension). In this work, we give a natural generalization of the Sobol indices (that are classically defined when $Y\in\mathbb R$ ), when the output belongs to $\mathbb{H}$. These indices have very nice properties. First, they are invariant. under isometry and scaling. Further they can be, as in dimension $1$, easily estimated by using the so-called Pick and Freeze method. We investigate the asymptotic behaviour of such estimation scheme.


title: Topological lattice field theories from intertwiner dynamics

abstract: We introduce a class of 2D lattice models that describe the dynamics of intertwiners, or, in a condensed matter interpretation, the fusion and splitting of anyons. We identify different families and instances of triangulation invariant, that is, topological, models inside this class. These models give examples for symmetry protected topologically ordered 1D quantum phases with quantum group symmetries. Furthermore the models provide realizations for anyon condensation into a new effective vacuum. We explain the relevance of our findings for the problem of identifying the continuum limit of spin foam and spin net models.


title: Searches for light- and heavy-flavour three-jet resonances in pp collisions at sqrt(s) = 8 TeV

abstract: A search for three-jet hadronic resonance production in pp collisions at a centre-of-mass energy of 8 TeV has been conducted by the CMS Collaboration at the LHC with a data sample corresponding to an integrated luminosity of 19.4 inverse femtobarns. The search method is model independent, and events are selected that have high jet multiplicity and large values of jet transverse momenta. The signal models explored assume R-parity-violating supersymmetric gluino pair production and have final states with either only light-flavour jets or both light- and heavy-flavour jets. No significant deviation is found between the selected events and the expected standard model multijet and t t-bar quark background. For a gluino decaying into light-flavour jets, a lower limit of 650 GeV on the gluino mass is set at a 95% confidence level, and for a gluino decaying into one heavy- and two light-flavour jets, gluino masses between 200 and 835 GeV are, for the first time, likewise excluded.


title: Sur les op\'erations de tores alg\'ebriques de complexit\'e un dans les vari\'et\'es affines

abstract: This thesis is devoted to the study of geometric properties of affine algebraic varieties endowed with an action of an algebraic torus. It comes from three preprints which correspond to the indicated points (1), (2), (3). Let $X$ be an affine variety equipped with an action of the algebraic torus $\mathbb{T}$. The complexity of the $\mathbb{T}$-action on $X$ is the codimension of the general $\mathbb{T}$-orbits. Under the assumption of normality and when the ground field is algebraically closed of characteristic 0, the variety $X$ admits a combinatorial description in terms of convex geometry. This description obtained by Altmann and Hausen in the year 2006 generalizes the classical one for toric varieties. Our purpose is to investigate new problems on the algebraic and geometric properties of the variety $X$ when the $\mathbb{T}$-action on $X$ is of complexity 1. (1) In the first part, a result gives an effective method to determine the integral closure of any affine variety defined over an algebraically field of characteristic 0 with a $\mathbb{T}$-action of complexity 1 in terms of the combinatorial description of Altmann-Hausen. (2) The calculations of the first part suggest a proof of the validity of the presentation of Altmann-Hausen in the case of complexity 1 over an arbitrary ground field. This is done in the second part. (3) In the third part, when the base field is perfect, we classify all the actions of the additive group on $X$ that are normalized by the $\mathbb{T}$-action of complexity 1. This classification generalizes classical works of Flenner and Zaidenberg in the surface case and of Liendo when the base field is algebraically closed of characteristic 0.


title: All-optical sensing of a single-molecule electron spin

abstract: We demonstrate an all-optical method for magnetic sensing of individual molecules in ambient conditions at room temperature. Our approach is based on shallow nitrogen-vacancy (NV) centers near the surface of a diamond crystal, which we use to detect single paramagnetic molecules covalently attached to the diamond surface. The manipulation and readout of the NV centers is all-optical and provides a sensitive probe of the magnetic field fluctuations stemming from the dynamics of the electronic spins of the attached molecules. As a specific example, we demonstrate detection of a single paramagnetic molecule containing a gadolinium (Gd$^{3+}$) ion. We confirm single-molecule resolution using optical fluorescence and atomic force microscopy to co-localize one NV center and one Gd$^{3+}$-containing molecule. Possible applications include nanoscale and in vivo magnetic spectroscopy and imaging of individual molecules.

categories: cond-mat.mes-hall physics.atom-ph physics.bio-ph quant-ph


title: Anisotropic wave propagation in nematic liquid crystals

abstract: Despite the fact that quantitative experimental data have been available for more than forty years now, nematoacoustics still poses intriguing theoretical and experimental problems. In this paper, we prove that the main observed features of acoustic wave propagation through a nematic liquid crystal cell -- namely, the anisotropy of sound velocity and its frequency dependence -- may be plausibly explained by a first-gradient continuum theory characterized by a hyperelastic anisotropic response from an evolving relaxed configuration. We compare and contrast our proposal with a competing theory where the liquid crystal is modeled as an isotropically compressible, anisotropic second-gradient fluid.


title: Strong Stability of Nash Equilibria in Load Balancing Games

abstract: We study strong stability of Nash equilibria in load balancing games of m (m >= 2) identical servers, in which every job chooses one of the m servers and each job wishes to minimize its cost, given by the workload of the server it chooses. A Nash equilibrium (NE) is a strategy profile that is resilient to unilateral deviations. Finding an NE in such a game is simple. However, an NE assignment is not stable against coordinated deviations of several jobs, while a strong Nash equilibrium (SNE) is. We study how well an NE approximates an SNE. Given any job assignment in a load balancing game, the improvement ratio (IR) of a deviation of a job is defined as the ratio between the pre- and post-deviation costs. An NE is said to be a r-approximate SNE (r >= 1) if there is no coalition of jobs such that each job of the coalition will have an IR more than r from coordinated deviations of the coalition. While it is already known that NEs are the same as SNEs in the 2-server load balancing game, we prove that, in the m-server load balancing game for any given m >= 3, any NE is a (5/4)-approximate SNE, which together with the lower bound already established in the literature yields a tight approximation bound. This closes the final gap in the literature on the study of approximation of general NEs to SNEs in load balancing games. To establish our upper bound, we make a novel use of a graph-theoretic tool.


title: The Dilaton Wess-Zumino Action in 6 Dimensions from Weyl Gauging: Local Anomalies and Trace Relations

abstract: We extend a previous analysis on the derivation of the dilaton Wess-Zumino (WZ) action in $d=4$, based on the method of Weyl gauging, to $6$ dimensions. As in the previous case, we discuss the structure of the same action in dimensional regularization using 6-dimensional Weyl invariants, extracting the dilaton interactions in the most general scheme, with the inclusion of the local anomaly terms. As an application, we present the WZ action for the (2,0) tensor multiplet, which has been investigated in the past in the context of the $AdS_7/CFT_6$ holographic anomaly matching. We then extend to $d=6$ the investigation of fully traced correlation functions of EMT's, formerly presented in $d=4$, showing that their hierarchy is functionally related only to the first 6 correlators. We give the explicit expressions of these in the most general scheme, up to rank-4.


title: A Many-Atom Cavity QED System with Homogeneous Atom-Cavity Coupling

abstract: We demonstrate a many-atom-cavity system with a high-finesse dual-wavelength standing wave cavity in which all participating rubidium atoms are nearly identically coupled to a 780-nm cavity mode. This homogeneous coupling is enforced by a one-dimensional optical lattice formed by the field of a 1560-nm cavity mode.


title: Measurement of the dispersion of radiation from a steady cosmological source

abstract: The `missing baryons' of the near universe are believed to be principally in a partially ionized state. Although passing electromagnetic waves are dispersed by the plasma, the effect has hitherto not been utilized as a means of detection because it is generally believed that a successful observation requires the background source to be highly variable, \ie~the class of sources that could potentially deliver a verdict is limited. We argue in two stages that this condition is not necessary. First, by modeling the fluctuations on macroscopic scales as interference between wave packets we show that, in accordance with the ideas advanced by Einstein in 1917, both the behavior of photons as bosons (\ie~the intensity variance has contributions from Poisson and phase noise) and the van-Cittert-Zernike theorem are a consequence of wave-particle duality. Nevertheless, we then point out that in general the variance on some macroscopic timescale $\tau$ consists of (a) a main contributing term $\propto 1/\tau$, plus (b) a small negative term $\propto 1/\tau^2$ due to the finite size of the wave packets. If the radiation passes through a dispersive medium, this size will be enlarged well beyond its vacuum minimum value of $\Delta t \approx 1/\Delta\nu$, leading to a more negative (b) term (while (a) remains unchanged) and hence a suppression of the variance w.r.t. the vacuum scenario. The phenomenon, which is typically at the few parts in 10$^5$ level, enables one to measure cosmological dispersion in principle. Signal-to-noise estimates, along with systematic issues and how to overcome them, will be presented.


title: Gapped symmetric edges of symmetry protected topological phases

abstract: Symmetry protected topological (SPT) phases are gapped quantum phases which host symmetry-protected gapless edge excitations. On the other hand, the edge states can be gapped by spontaneously breaking symmetry. We show that topological defects on the symmetry-broken edge cannot proliferate due to their fractional statistics. A gapped symmetric boundary, however, can be achieved between an SPT phase and certain fractionalized phases by condensing the bound state of a topological defect and an anyon. We demonstrate this by two examples in two dimensions: an exactly solvable model for the boundary between topological Ising paramagnet and double semion model, and a fermionic example about the quantum spin Hall edge. Such a hybrid structure containing both SPT phase and fractionalized phase generally support ground state degeneracy on torus.


title: Gravitational corrections to the scattering of scalar particles

abstract: We evaluated the scattering amplitude of neutral scalar particles at one-loop order in the context of effective field theory of quantum gravity in the presence of a cosmological constant. Our study suggests that quantum gravitational corrections induce an asymptotic freedom behavior to the $\lambda\phi^4$ theory, for a positive cosmological constant. This result hints that a complete theory of quantum gravity can play an important role to avoid the issue of triviality in quantum field theory.


title: How to lift positive Ricci curvature

abstract: We show how to lift positive Ricci and almost non-negative curvatures from an orbit space $M/G$ to the corresponding $G$-manifold, $M$. We apply the results to get new examples of Riemannian manifolds that satisfy both curvature conditions simultaneously.


title: Feasibility of top quark measurements at LHCb and constraints on the large-$x$ gluon PDF

abstract: The forward LHCb acceptance opens interesting possibilities of studying precision Standard Model hard processes in a kinematical region beyond the reach of ATLAS and CMS. In this paper we perform a feasibility study for cross-section measurements of top quark pairs with the LHCb detector, with an analysis of signal and background rates for selected final states, and determine the potential precision achievable at $\sqrt{s} =$ 7 and 14~TeV. We then study the dependence of theoretical uncertainties on the pseudorapidity distribution of top quarks produced in pair production at NLO, and observe that a cross-section measurement at high pseudorapidity has enhanced sensitivity to probe the high-$x$ gluon PDF as compared to measurements in the central-region. Based on simulated pseudodata, the impact of a 14~TeV cross-section measurement on the gluon PDF and charge asymmetry is quantified.


title: The many facets of the (non relativistic) Nuclear Equation of State

abstract: A nucleus is a quantum many body system made of strongly interacting Fermions, protons and neutrons (nucleons). This produces a rich Nuclear Equation of State whose knowledge is crucial to our understanding of the composition and evolution of celestial objects. The nuclear equation of state displays many different features; first neutrons and protons might be treated as identical particles or nucleons, but when the differences between protons and neutrons are spelled out, we can have completely different scenarios, just by changing slightly their interactions. At zero temperature and for neutron rich matter, a quantum liquid gas phase transition at low densities or a quark-gluon plasma at high densities might occur. Furthermore, the large binding energy of the $\alpha$ particle, a Boson, might also open the possibility of studying a system made of a mixture of Bosons and Fermions, which adds to the open problems of the nuclear equation of state.


title: Conditions for Bifurcations in a Non-Autonomous Scalar Differential Equation

abstract: In this paper are provided some sufficient conditions for a non autonomous scalar differential equation to have saddle node, transcritical and pitchfork bifurcations using higher order derivatives.


title: On $r$-cross $t$-intersecting families for weak compositions

abstract: Let $\mathbb N_0$ be the set of non-negative integers, and let $P(n,l)$ denote the set of all weak compositions of $n$ with $l$ parts, i.e., $P(n,l)=\{ (x_1,x_2,\dots, x_l)\in\mathbb N_0^l\ :\ x_1+x_2+\cdots+x_l=n\}$. For any element $\mathbf u=(u_1,u_2,\dots, u_l)\in P(n,l)$, denote its $i$th-coordinate by $\mathbf u(i)$, i.e., $\mathbf u(i)=u_i$. Let $l=\min(l_1,l_2,\dots, l_r)$. Families $\mathcal A_j\subseteq P(n_j,l_j)$ ($j=1,2,\dots, r$) are said to be $r$-cross $t$-intersecting if $\vert \{ i\in [l] \ :\ \mathbf u_1(i)=\mathbf u_2(i)=\cdots=\mathbf u_r(i)\} \vert\geq t$ for all $\mathbf u_j\in \mathcal A_j$. Suppose that $l\geq t+2$. We prove that there exists a constant $n_0=n_0(l_1,l_2,\dots,l_r,t)$ depending only on $l_j$'s and $t$, such that for all $n_j\geq n_0$, if the families $\mathcal A_j\subseteq P(n_j,l_j)$ ($j=1,2,\dots, r$) are $r$-cross $t$-intersecting, then \begin{equation} \prod_{j=1}^r \vert \mathcal{A}_j \vert\leq \prod_{j=1}^r {n_j+l_j-t-1 \choose l_j-t-1}.\notag \end{equation} Moreover, equality holds if and only if there is a $t$-set $T$ of $\{1,2,\dots,l\}$ such that $\mathcal{A}_j=\{\mathbf u\in P(n_j,l_j)\ :\ \mathbf u(i)=0\ {\rm for\ all}\ i\in T\}$ for $j=1,2,\dots, r$.


title: Solar Mini-Dimming Kinematics and Their Positive Correlations with Coronal Mass Ejections and Prominence

abstract: Solar mini-dimmings can be detect in the Extreme Ultra-Violet coro- nal eruptions. Here, sequences of 171 A images taken by Solar Dynamic Observa- tory/Atmospheric Imaging Assembaly on 13 June 2010 are used. In this special day, both of coronal mass ejection and prominence were observed. The average velocities and accelerations of 500 mini-dimmings which were detected using on feature based classifier (Alipour et al 2012) are studied. The large number of mini-dimmings shows positive accelerations in the beginning times as similar as Coronal Mass Ejections. On the the start time of prominence eruptions, the number of mini-dimmings is increased to a maximum values. There is a positive correlation between the kinematics of mini- dimmings and both CME and prominences. This study can be extended to understand the exact relationship of CMEs and mini-dimmings.


title: Ending AIDS in South Africa: How long will it take? How much will it cost?

abstract: South Africa has more people infected with HIV but, by providing access to anti-retroviral therapy (ART), has kept more people alive than any other country. The effectiveness, availability and affordability of potent anti-retroviral therapy (ART) make it possible to contemplate ending the epidemic of HIV/AIDS. We consider what would have happened without ART, the impact of the current roll-out of ART, what might be possible if early treatment becomes available to all, and what could have happened if ART had been provided much earlier in the epidemic. In 2013 the provision of ART has reduced the prevalence of HIV from an estimated 15% to 9% among adults not on ART, the annual incidence from 2% to 0.9%, and the AIDS related deaths from 0.9% to 0.3% p.a. saving 1.5 million lives and USD727M. Regular testing and universal access to ART could reduce the prevalence among adults not on ART in 2023 to 0.06%, annual incidence to 0.05%, and eliminate AIDS deaths. Cumulative costs between 2013 ands 2023 would increase by USD692M only 4% of the total cost of USD17Bn. If a universal testing and early treatment had started in 1998 the prevalence of HIV among adults not on ART in 2013 would have fallen to 0.03%, annual incidence to 0.03%, and saved 2.5 million lives. The cost up to 2013 would have increased by USD18Bn but this would have been cost effective at US$7,200 per life saved. Future surveys of HIV among women attending ante-natal clinics should include testing women for the presence of anti-retroviral drugs, measuring their viral loads, and using appropriate assays for estimating HIV incidence. These data would make it possible to develop better and more reliable estimates of the current state of the epidemic, the success of the current ART programme, levels of viral load suppression for those on ART and the incidence of infection.


title: The dependence of the galaxy mass-metallicity relation on environment and the implied metallicity of the IGM

abstract: We explore the dependence of the galaxy mass-metallicity relation on environment in SDSS, in terms of both over-density and central/satellite dichotomy. We find that at a given stellar mass, there is a strong dependence of metallicity on over-density for star-forming satellites. Instead, for star-forming centrals no correlation is found. Star-forming satellites at different stellar masses form a tight sequence in the average overdensity - metallicity plane, which covers the entire observed range of metallicities and stellar masses. This remarkable result appears to imply that there exists a universal evolutionary path for all star-forming satellites, regardless of their stellar masses. The strong correlation between over-density and metallicity for star-forming satellites indicates that the gas inflow of satellite galaxies is progressively metal-enriched in denser regions. We interpret our results by employing the gas regulator model and find that the metallicity of the enriched inflow of star-forming satellite galaxies strongly increases with increasing over-density, largely independently of stellar mass. If the metallicity of the inflow of star-forming satellites can represent the metallicity of the IGM, then the implied metallicity of the IGM rises from ~ 0.01 solar metallicity in void-like environments to ~ 0.3 solar metallicity in cluster-like environments, in broad agreement with observations. We show that the observed metallicity difference between star-forming centrals and star-forming satellites becoming smaller towards high stellar masses can be simply explained by the mass-independent enriched inflow, without the need to involve any mass-dependent environmental effect on metallicity. Since satellite galaxies account for at least half of the galaxy population, our findings prompt for a revision of many galaxy evolutionary models, which generally assume pristine gas inflows.


title: Quantum Spin Ice: A Search for Gapless Quantum Spin Liquids in Pyrochlore Magnets

abstract: The spin ice materials, including Ho2Ti2O7 and Dy2Ti2O7, are rare earth pyrochlore magnets which, at low temperatures, enter a constrained paramagnetic state with an emergent gauge freedom. Remarkably, the spin ices provide one of very few experimentally realised examples of fractionalization because their elementary excitations can be regarded as magnetic monopoles and, over some temperature range, the spin ice materials are best described as liquids of these emergent charges. In the presence of quantum fluctuations, one can obtain, in principle, a quantum spin liquid descended from the classical spin ice state characterised by emergent photon-like excitations. Whereas in classical spin ices the excitations are akin to electrostatic charges, in the quantum spin liquid these charges interact through a dynamic and emergent electromagnetic field. In this review, we describe the latest developments in the study of such a quantum spin ice, focussing on the spin liquid phenomenology and the kinds of materials where such a phase might be found.


title: Galaxy Assembly Bias: A Significant Source of Systematic Error in the Galaxy-Halo Relationship

abstract: It is common practice for methods that use galaxy clustering to constrain the galaxy-halo relationship, such as the halo occupation distribution (HOD) and/or conditional luminosity function (CLF), to assume that halo mass alone suffices to determine a halo's resident galaxy population. Yet, the clustering strength of cold dark matter halos depends upon halo properties in addition to mass, such as formation time, an effect referred to as assembly bias. If galaxy characteristics are correlated with any of these auxiliary halo properties, the basic assumption of HOD/CLF methods is violated. We estimate the potential for assembly bias to induce systematic errors in inferred halo occupation statistics. We use halo abundance matching and age matching to construct fiducial mock galaxy catalogs that exhibit assembly bias as well as additional mock catalogs with identical HODs, but with assembly bias removed. We fit a parameterized HOD to the projected two-point clustering of mock galaxies in each catalog to assess the systematic errors induced by reasonable levels of assembly bias. In the absence of assembly bias, the inferred HODs generally describe the true underlying HODs well, validating the basic methodology. However, in all of the cases with assembly bias, the inferred HODs have systematic errors that are statistically significant. In most cases, these systematic errors cannot be identified using void statistics as auxiliary observables. We conclude that the galaxy-halo relationship inferred from galaxy clustering should be subject to a non-negligible systematic error induced by assembly bias. Our work suggests that efforts to model and/or constrain assembly bias should be high priorities as it is a threatening source of systematic error in galaxy evolution studies as well as the precision cosmology program.


title: The Dispersal of Protoplanetary Disks

abstract: Protoplanetary disks are the sites of planet formation, and the evolution and eventual dispersal of these disks strongly influences the formation of planetary systems. Disk evolution during the planet-forming epoch is driven by accretion and mass-loss due to winds, and in typical environments photoevaporation by high-energy radiation from the central star is likely to dominate final gas disk dispersal. We present a critical review of current theoretical models, and discuss the observations that are used to test these models and inform our understanding of the underlying physics. We also discuss the role disk dispersal plays in shaping planetary systems, considering its influence on both the process(es) of planet formation and the architectures of planetary systems. We conclude by presenting a schematic picture of protoplanetary disk evolution and dispersal, and discussing prospects for future work.


title: Gravitational and distributed heating effects of a cD galaxy on the hydrodynamical structure of its host cluster

abstract: We investigate the effects of a cD galaxy's gravity and AGN heating of the host galaxy cluster. We consider a standard prescription for the hydrodynamics, with the structures determined by mass continuity, momentum and energy conservation equations in spherical symmetry. The cluster comprises a dark matter halo (DM) and ionised X-ray emitting intracluster gas (ICM), which jointly determine the gravitational potential. The cD galaxy is an additive gravitational potential component. The DM assumes a polytropic equation of state (determined by its microphysics), which could be non-radiative self-interacting particles or more exotically interacting particles. The AGN provides distributed heating, counteracting radiative cooling. Stationary density and velocity dispersion profiles are obtained by numerically integrating the hydrodynamic equations with appropriate boundary conditions. The minimum gas temperature in the cluster core is higher when a cD galaxy is present than when it is absent. The solutions also yield a point-like mass concentration exceeding a minimum mass: presumably the AGN's supermassive black hole (SMBH). Consistency with observed SMBH masses constrains the possible DM equations of state. The constraints are looser when a cD galaxy is present. Distributed (AGN) heating alters cluster global properties, and also reduces the lower limits for the central point-mass, for the preferred DM models in which the dark particles have greater heat capacity than point particles. Eluding these constraints would require dominant non-spherical or anisotropic effects (e.g. bulk rotation, non-radial streaming, asymmetric lumps, or a strong magnetic field).


title: Around 200 new X-ray binary IDs from 13 years of Chandra observations of the M31 center

abstract: We have created 0.3--10 keV, 13 year, unabsorbed luminosity lightcurves for 528 X-ray sources in the central 20' of M31. We have 174 Chandra observations spaced at ~1 month intervals thanks to our transient monitoring program, deeper observations of the M31 nucleus, and some public data from other surveys. We created 0.5--4.5 keV structure functions (SFs) for each source, for comparison with the ensemble structure function of AGN. We find 220 X-ray sources with luminosities > ~1E+35 erg/s that have SFs with significantly more variability than the ensemble AGN SF, and are likely X-ray binaries (XBs). A further 30 X-ray sources were identified as XBs using other methods. We therefore have 250 probable XBs in total, including ~200 new identifications. This result represents great progress over the ~50 XBs and ~40 XB candidates previously identified out of the ~2000 X-ray sources within the D_25 region of M31; it also demonstrates the power of SF analysis for identifying XBs in external galaxies. We also identify a new transient black hole candidate, associated with the M31 globular cluster B128.


title: Quantifying the sensitivity of oscillation experiments to the neutrino mass ordering

abstract: Determining the type of the neutrino mass ordering (normal versus inverted) is one of the most important open questions in neutrino physics. In this paper we clarify the statistical interpretation of sensitivity calculations for this measurement. We employ standard frequentist methods of hypothesis testing in order to precisely define terms like the median sensitivity of an experiment. We consider a test statistic $T$ which in a certain limit will be normal distributed. We show that the median sensitivity in this limit is very close to standard sensitivities based on $\Delta\chi^2$ values from a data set without statistical fluctuations, such as widely used in the literature. Furthermore, we perform an explicit Monte Carlo simulation of the INO, JUNO, LBNE, NOvA, and PINGU experiments in order to verify the validity of the Gaussian limit, and provide a comparison of the expected sensitivities for those experiments.


title: Disentangling a dynamical Higgs

abstract: The pattern of deviations from Standard Model predictions and couplings is different for theories of new physics based on a non-linear realization of the $SU(2)_L\times U(1)_Y$ gauge symmetry breaking and those assuming a linear realization. We clarify this issue in a model-independent way via its effective Lagrangian formulation in the presence of a light Higgs particle, up to first order in the expansions: dimension-six operators for the linear expansion and four derivatives for the non-linear one. Complete sets of pure gauge and gauge-Higgs operators are considered, implementing the renormalization procedure and deriving the Feynman rules for the non-linear expansion. We establish the theoretical relation and the differences in physics impact between the two expansions. Promising discriminating signals include the decorrelation in the non-linear case of signals correlated in the linear one: some pure gauge versus gauge-Higgs couplings and also between couplings with the same number of Higgs legs. Furthermore, anomalous signals expected at first order in the non-linear realization may appear only at higher orders of the linear one, and vice versa. We analyze in detail the impact of both type of discriminating signals on LHC physics.


title: Higgs Production in General 5D Warped Models

abstract: We calculate the production rate of the Higgs boson at the LHC in the context of general 5 dimensional (5D) warped scenarios with spacetime background modified from the usual $AdS_5$, and where all the SM fields, including the Higgs, propagate in the bulk. This modification can alleviate considerably the bounds coming from precision electroweak tests and flavor physics. We evaluate the Higgs production rate and show that it is generically consistent with the current experimental results from the LHC for Kaluza-Klein (KK) masses as low as 2 TeV, unlike in pure $AdS_5$ scenarios, where for the same masses, the Higgs production typically receives corrections too large to be consistent with LHC data. Thus the new pressure on warped models arising from LHC Higgs data is also alleviated in $AdS_5$-modified warped scenarios.


title: Competing regimes of motion in 1D mobile impurities

abstract: We show that a distinguishable mobile impurity inside a one-dimensional many-body state at zero temperature generally does not behave like a quasiparticle (QP). Instead, both the impurities dynamics as well as the ground state of the bath are fundamentally transformed by a diverging number of zero-energy excitations being generated, leading to what we call infrared-dominated (ID) dynamics. Combining analytics and DMRG numerics we provide a general formula for the power law governing ID dynamics at zero momentum, discuss a threshold beyond which quasiparticle dynamics may occur again, and study the competition between the ID and quasiparticle universality classes at larger impurity momenta.


title: Hybrid Bounds on Twisted L-Functions Associated to Modular Forms

abstract: For $f$ a primitive holomorphic cusp form of even weight $k \geq 4$, level $N$, and $\chi$ a Dirichlet character mod $Q$ with $(Q,N)=1$, we establish a new hybrid subconvexity bound for $L(1/2 + it, f_\chi)$, which improves upon all known hybrid bounds. This is done via amplification and taking advantage of a shifted convolution sum of two variables defined and analyzed in a recent paper of Hoffstein and Hulse.


title: Three easy exercises in off-shell string-inspired methods

abstract: Off-shell string-inspired methods (OSSIM) calculate off-shell QCD Green's functions using Schwinger-Feynman proper-time techniques, always in the background field method (BFM) Feynman gauge for technical convenience, and so far only at one loop. We already know that these results are gauge-invariant, because this gauge realizes the prescriptions of the Pinch Technique (PT), a Feynman-graph formulation for any gauge, but the idea of the first exercise is to show this directly in OSSIM. In this exercise we extend proper-time OSSIM beyond the BFM Feynman gauge so that one can apply PT algorithms, and show that the intrinsic PT is equivalent to resolving ambiguities in OSSIM in other gauges. In the second exercise we use forty-year-old rules of the author and Tiktopoulos for expressing loop integrals with numerator momenta directly in terms of Feynman parameters after momentum integration (the goal of OSSIM) and show that these rules elegantly and with economy of effort give rise, at least at one loop, to standard OSSIM algorithms. In the third exercise we apply world-line techniques to the problem of the breaking of adjoint strings, requiring a non-perturbative treatment that in the end reduces to a variant of the Schwinger result for production of electron-positron pairs in an electric field. This generalizes OSSIM to non-perturbative processes.


title: Infrared Properties of z=7 Galaxies from Cosmological Simulations

abstract: Three-dimensional panchromatic dust radiative transfer calculations are performed on a set of 198 galaxies of stellar masses in the range 5x10^8-3x10^10 Msun from a cosmological hydrodynamic simulation (resolved at 29pc/h) at z=7. In a companion paper (Kimm & Cen), the stellar mass and UV luminosity functions, and UV-optical and FUV-NUV colors are shown to be in good agreement with observations, if an SMC-type dust extinction curve is adopted. Here we make useful predictions, self-consistently, of the infrared properties of these z=7 simulated galaxies that can be confronted with upcoming ALMA data. Our findings are as follows. (1) The effective radius in the rest-frame MIPS 70 micron band is in the range of 80-400pc proper for z=7 galaxies with L_FIR=10^{11.3-12}Lsun. (2) The median of the peak wavelength of the far-infrared (FIR) spectral energy distribution is in the range of 45-60 micron, depending on the dust-to-metal ratio. (3) For star formation rate in the range 3-100 Msun/yr the median FIR to bolometric luminosity ratio is 60-90%. (4) The FIR luminosity function displays a power law in the high end with a slope of -3.1 +- 0.4, instead of the usual exponential decline.


title: Higgs characterisation via vector-boson fusion and associated production: NLO and parton-shower effects

abstract: Vector-boson fusion and associated production at the LHC can provide key information on the strength and structure of the Higgs couplings to the Standard Model particles. Using an effective field theory approach, we study the effects of next-to-leading order (NLO) QCD corrections matched to parton shower on selected observables for various spin-0 hypotheses. We find that inclusion of NLO corrections is needed to reduce the theoretical uncertainties on total rates as well as to reliably predict the shapes of the distributions. Our results are obtained in a fully automatic way via FeynRules and MadGraph5_aMC@NLO.


title: QCD analysis of polarized DIS and the SIDIS asymmetry world data and light sea-quark decomposition

abstract: The results of our new QCD analysis of helicity parton distributions of the nucleon at full next-to-leading order (NLO) accuracy in the fixed-flavor number scheme will be presented. Performing a combined QCD fit on the global sets of latest inclusive and the semi-inclusive polarized deep inelastic scattering data, we are able to extract new polarized parton distribution functions (PPDFs) at the input scale $Q_{0}^2=1$ GeV$^2$. Particulary, we have calculated PPDFs considering light sea-quark decomposition and the results are compared with the experimental data and the most precise theoretical models obtained by recent analyses. The latest COMPASS2010 SIDIS data, which were not available for the previous analyses, are employed in the current analysis and the effect of COMPASS SIDIS data is studied in detail. Also the uncertainties of PPDFs are determined using the standard Hessian technique.


title: Linear theory for filtering nonlinear multiscale systems with model error

abstract: We study filtering of multiscale dynamical systems with model error arising from unresolved smaller scale processes. The analysis assumes continuous-time noisy observations of all components of the slow variables alone. For a linear model with Gaussian noise, we prove existence of a unique choice of parameters in a linear reduced model for the slow variables. The linear theory extends to to a non-Gaussian, nonlinear test problem, where we assume we know the optimal stochastic parameterization and the correct observation model. We show that when the parameterization is inappropriate, parameters chosen for good filter performance may give poor equilibrium statistical estimates and vice versa. Given the correct parameterization, it is imperative to estimate the parameters simultaneously and to account for the nonlinear feedback of the stochastic parameters into the reduced filter estimates. In numerical experiments on the two-layer Lorenz-96 model, we find that parameters estimated online, as part of a filtering procedure, produce accurate filtering and equilibrium statistical prediction. In contrast, a linear regression based offline method, which fits the parameters to a given training data set independently from the filter, yields filter estimates which are worse than the observations or even divergent when the slow variables are not fully observed.


title: High-Energy X-rays from J174545.5-285829, the Cannonball: A Candidate Pulsar Wind Nebula Associated with Sgr A East

abstract: We report the unambiguous detection of non-thermal X-ray emission up to 30 keV from the Cannonball, a few-arcsecond long diffuse X-ray feature near the Galactic Center, using the NuSTAR X-ray observatory. The Cannonball is a high-velocity (vproj~500 km/s) pulsar candidate with a cometary pulsar wind nebula (PWN) located ~2' north-east from Sgr A*, just outside the radio shell of the supernova remnant Sagittarius A (Sgr A) East. Its non-thermal X-ray spectrum, measured up to 30 keV, is well characterized by a Gamma~1.6 power-law, typical of a PWN, and has an X-ray luminosity of L(3-30 keV)=1.3e34 erg/s. The spectral and spatial results derived from X-ray and radio data strongly suggest a runaway neutron star born in the Sgr A East supernova event. We do not find any pulsed signal from the Cannonball. The NuSTAR observations allow us to deduce the PWN magnetic field and show that it is consistent with the lower limit obtained from radio observations.


title: Proper holomorphic mappings into $\ell$-concave quadric domains in projective space

abstract: In this paper, we prove a type of partial rigidity result for proper holomorphic mappings of certain $\ell$-concave domains in projective space into model quadratic $\ell$-concave domains. The main technical result is a degree estimate for proper holomorphic mappings into the model domains, provided that the mappings extend to projective space as rational mappings, and the source domain contains algebraic varieties and has a boundary with low CR complexity.


title: Pseudodifferentialoperatoren mit nichtglatten Koeffizienten auf Mannigfaltigkeiten

abstract: A generalization of differential operators are pseudodifferential operators which are used for reasoning about partial differential equations with variable coefficients. A lot of useful properties about classical pseudodifferential operators are already discovered by the works of H. Kumano-Go and L. H\"ormander. In this work, we like to transfer their ideas on those operators which symbols are merely H\"older-continous. We focus on the composition of operators, symbol smoothing and operators on manifolds.


title: Linear Regression without computing pseudo-inverse matrix

abstract: We are presenting a method of linear regression based on Gram-Schmidt orthogonal projection that does not compute a pseudo-inverse matrix. This is useful when we want to make several regressions with random data vectors for simulation purposes.


title: A New Stochastic Model of the Causal Interpretation of Quantum Theory on the Development of the Fundamental Concept of Mass

abstract: In this paper we pose two fundamental ideas on the motion of an elementary particle supporting the internal "spin motion" or $\textit{Zitterbewegung}$ and a particle as concentrated energy. First, the particle moves randomly in a limited area (in a quantum sized-volume) like random vibrating system where the particle will diffuse in a quantum-sized volume when it absorbs or emits the quantized amount of energy. The quantum sized-volume can move too and plays role similar to the carrier amplitude, while the vibrational motion of the fundamental particle with frequency and amplitude represents a modulating signal. The current of diffusion process taking place in the quantum sized-volume will represent the emergence of a spin phenomenon shown by the existence of Clifford algebra. Second, the particle is pure energy concentrated on the surface of 3 dimensional sphere-form (2-manifold without boundary). Afterward we show that by defining the particle mass as an invariant quantity based on the two fundamental ideas, we can derive both the diffusion constant of the particle in the quantum sized-volume as $\beta = \hslash/2m$ and the Schr\"{o}dinger equation. Furthermore, by posing that the vibrational motion of the particle limited on the quantum sized-volume plays fundamental role as time interval unit (proper time or particle clock), we show that the relativistic effects of the particle must represent atomic processes.


title: Split-Family SUSY, U(2)^5 Flavour Symmetry and Neutrino Physics

abstract: In split-family SUSY, one can use a U(2)^3 symmetry to protect flavour observables in the quark sector from SUSY contributions. However, attempts to extend this procedure to the lepton sector by using an analogous U(2)^5 symmetry fail to reproduce the neutrino data without introducing some form of fine-tuning. In this work, we solve this problem by shifting the U(2)^2 symmetry acting on leptons towards the second and third generations. This allows neutrino data to be reproduced without much difficulties, as well as protecting the leptonic flavour observables from SUSY. Key signatures are a $\mu\to e\gamma$ branching ratio possibly observable in the near future, as well as having selectrons as the lightest sleptons.


title: Efficient Regularization of Squared Curvature

abstract: Curvature has received increased attention as an important alternative to length based regularization in computer vision. In contrast to length, it preserves elongated structures and fine details. Existing approaches are either inefficient, or have low angular resolution and yield results with strong block artifacts. We derive a new model for computing squared curvature based on integral geometry. The model counts responses of straight line triple cliques. The corresponding energy decomposes into submodular and supermodular pairwise potentials. We show that this energy can be efficiently minimized even for high angular resolutions using the trust region framework. Our results confirm that we obtain accurate and visually pleasing solutions without strong artifacts at reasonable run times.


title: An Efficiently Solvable Quadratic Program for Stabilizing Dynamic Locomotion

abstract: We describe a whole-body dynamic walking controller implemented as a convex quadratic program. The controller solves an optimal control problem using an approximate value function derived from a simple walking model while respecting the dynamic, input, and contact constraints of the full robot dynamics. By exploiting sparsity and temporal structure in the optimization with a custom active-set algorithm, we surpass the performance of the best available off-the-shelf solvers and achieve 1kHz control rates for a 34-DOF humanoid. We describe applications to balancing and walking tasks using the simulated Atlas robot in the DARPA Virtual Robotics Challenge.


title: First-principles calculation of entropy for liquid metals

abstract: We demonstrate the accurate calculation of entropies and free energies for a variety of liquid metals using an extension of the two phase thermodynamic (2PT) model based on a decomposition of the velocity autocorrelation function into gas-like (hard sphere) and solid-like (harmonic) subsystems. The hard sphere model for the gas-like component is shown to give systematically high entropies for liquid metals as a direct result of the unphysical Lorentzian high-frequency tail. Using a memory function framework we derive a generally applicable velocity autocorrelation and frequency spectrum for the diffusive component which recovers the low frequency (long time) behavior of the hard sphere model while providing for realistic short time coherence and high frequency tails to the spectrum. This approach provides a significant increase in the accuracy of the calculated entropies for liquid metals and is compared to ambient pressure data for liquid sodium, aluminum, gallium, tin, and iron. The use of this method for the determination of melt boundaries is demonstrated with a calculation of the high pressure bcc melt boundary for sodium. With the significantly improved accuracy available with the memory function treatment for softer interatomic potentials, the 2PT model for entropy calculations should find broader application in high energy density science, warm dense matter, planetary science, geophysics, and material science.


title: Towards the event horizon -- the supermassive black hole in the Galactic Center

abstract: The center of our Galaxy hosts the best constrained supermassive black hole in the universe, Sagittarius A* (Sgr A*). Its mass and distance have been accurately determined from stellar orbits and proper motion studies, respectively, and its high-frequency radio, and highly variable near-infrared and X-ray emission originate from within a few Schwarzschild radii of the event horizon. The theory of general relativity (GR) predicts the appearance of a black hole shadow, which is a lensed image of the event horizon. This shadow can be resolved by very long baseline radio interferometry and test basic predictions of GR and alternatives thereof. In this paper we review our current understanding of the physical properties of Sgr A*, with a particular emphasis on the radio properties, the black hole shadow, and models for the emission and appearance of the source. We argue that the Galactic Center holds enormous potential for experimental tests of black hole accretion and theories of gravitation in their strong limits.


title: Extremal Domains for Self-Commutators in the Bergman Space

abstract: In arXiv:1305.5193v1, the authors have shown that Putnam's inequality for the norm of self-commutators can be improved by a factor of $\frac{1}{2}$ for Toeplitz operators with analytic symbol $\varphi$ acting on the Bergman space $A^{2}(\Omega)$. This improved upper bound is sharp when $\varphi(\Omega)$ is a disk. In this paper we show that disks are the only domains for which the upper bound is attained.


title: Nuclear Structure Uncertainties in Parity-Violating Electron Scattering from Carbon 12

abstract: High precision measurements of the parity-violating asymmetry in polarized electron scattering from nuclei can be used to extract information on nuclear and nucleon structure or to determine Standard Model couplings and higher-order radiative corrections. To this end, low uncertainties are also required in the effects that inevitably arise from modeling the underlying nuclear structure. An experimental precision of a few tenths of a percent may be attainable for the asymmetry if the appropriate kinematic range is chosen, as will be discussed here for the case of $^{12}$C. And given this, the dual goal of ascertaining both the sizes of various nuclear structure related effects and of providing estimates of their uncertainties for this particular target will be discussed.


title: Tolman IV solution in the Randall-Sundrum Braneworld

abstract: In the context of the Randall-Sundrum braneworld, the minimal geometric deformation approach (MGD) is used to generate an exact analytic interior solution to four-dimensional effective Einstein's field equations for a spherically symmetric compact distribution. This solution represents the braneworld version of the well known Tolman IV solution in General Relativity. By using this analytic solution, an exhaustive analysis of the braneworld effects on realistic stellar interiors is developed, finding strong evidences in favor of the hypothesis that compactness is reduced due to bulk effects on stellar configurations.


title: Non abelian hydrodynamics and heavy ion collisions

abstract: The goal of the relativistic heavy ion collisions (RHIC) program is to create a state of matter where color degrees of freedom are deconfined. The dynamics of matter in this state, in spite of the complexities of quantum chromodynamics, is largely determined by the conservation laws of energy momentum and color currents. Therefore it is possible to describe its main features in hydrodynamic terms, the very short color neutralization time notwithstanding. In this lecture we shall give a simple derivation of the hydrodynamics of a color charged fluid, by generalizing the usual derivation of hydrodynamics from kinetic theory to the non abelian case.


title: Demonstration of long-lived high power optical waveguides in air

abstract: We demonstrate that femtosecond filaments can set up an extended and robust thermal waveguide structure in air with a lifetime of several milliseconds, making possible the very long range guiding and distant projection of high energy laser pulses and high average power beams. As a proof of principle, we demonstrate guiding of 110 mJ, 7 ns, 532 nm pulses with 90% throughput over ~15 Rayleigh lengths in a 70 cm long air waveguide generated by the long timescale thermal relaxation of an array of femtosecond filaments. The guided pulse was limited only by our available laser energy. In general, these waveguides should be robust against the effects of thermal blooming of extremely high average power laser beams.


title: Search for Natural SUSY with inclusive search strategies at the LHC using the CMS detector

abstract: Natural SUSY suggests the existence of light stop quarks, accessible at the LHC, which are the focus of a dedicated CMS search program. I present two inclusive CMS searches that look for TeV scale colored sparticles in final states with jets, b-tagged jets and missing transverse energy performed using up to 19.4fb-1 of 8TeV LHC proton-proton data. No deviation from the Standard Model was observed in these searches, and the implications for this was shown for several simplified model scenarios and phenomenological MSSM.


title: On the admissibility of certain local system

abstract: A rank one local system on the complement of a hyperplane arrangement is said to be admissible if it satisfies certain non-positivity condition at every resonant edges. It is known that the cohomology of admissible local system can be computed combinatorially. In this paper, we study the structure of the set of all non-admissible local systems in the character torus. We prove that the set of non-admissible local systems forms a union of subtori. The relations with characteristic varieties are also discussed.


title: Hydrodynamics at RHIC and LHC: What have we learned?

abstract: Fluid dynamical description of elementary particle collisions has a long history dating back to the works of Landau and Fermi. Nevertheless, it is during the last 10-15 years when fluid dynamics has become the standard tool to describe the evolution of matter created in ultrarelativistic heavy ion collisions. In this review I briefly describe the hydrodynamical models, what we have learned when analyzing the RHIC and LHC data using these models, and what the latest developments and challenges are.


title: Impacts of demographic catastrophes

abstract: Analysis of demographic catastrophes shows that, with the exception of perhaps only two critical events, they were too weak to influence the growth of human population. These results reinforce the conclusion that the concept of the Epoch of Malthusian Stagnation, the alleged first stage of growth claimed by the Demographic Transition Theory, is not supported by empirical evidence. They show that even if we assume that Malthusian positive checks are capable of suppressing the growth of population their impact was too weak to create the Epoch of Malthusian Stagnation.


title: Efficient Quantum Algorithms for Analyzing Large Sparse Electrical Networks

abstract: Analyzing large sparse electrical networks is a fundamental task in physics, electrical engineering and computer science. We propose two classes of quantum algorithms for this task. The first class is based on solving linear systems, and the second class is based on using quantum walks. These algorithms compute various electrical quantities, including voltages, currents, dissipated powers and effective resistances, in time $\operatorname{poly}(d, c, \operatorname{log}(N), 1/\lambda, 1/\epsilon)$, where $N$ is the number of vertices in the network, $d$ is the maximum unweighted degree of the vertices, $c$ is the ratio of largest to smallest edge resistance, $\lambda$ is the spectral gap of the normalized Laplacian of the network, and $\epsilon$ is the accuracy. Furthermore, we show that the polynomial dependence on $1/\lambda$ is necessary. This implies that our algorithms are optimal up to polynomial factors and cannot be significantly improved.


title: Concrete Categories in Homotopy Type Theory

abstract: We introduce some classes of genuine higher categories in homotopy type theory, defined as well-behaved subcategories of the category of types. We give several examples, and some techniques for showing other things are not examples. While only a small part of what is needed, it is a natural construction, and may be instructive for people seeking to provide a fully general construction.


title: Topologically Massive Yang-Mills Theory and Link Invariants

abstract: Topologically massive Yang-Mills theory is studied in the framework of geometric quantization. Since this theory has a mass gap proportional to the topological mass m, Yang-Mills contribution decays exponentially at very large distances compared to 1/m, leaving a pure Chern-Simons theory with level number k. In this paper, the near Chern-Simons limit is studied where the distance is large enough to give an almost topological theory, with a small contribution from the Yang-Mills term. It is shown that this almost topological theory consists of two copies of Chern-Simons with level number k/2, very similar to the Chern-Simons splitting of topologically massive AdS gravity. Also, gauge invariance of these half-Chern-Simons theories is discussed. As m approaches to infinity, the split parts add up to give the original Chern-Simons term with level k. Reduction of the phase space is discussed in this limit. Finally, a relation between the observables of topologically massive Yang-Mills theory and Chern-Simons theory is shown. One of the two split Chern-Simons pieces is shown to be associated with Wilson loops while the other with 't Hooft loops. This allows one to use skein relations to calculate topologically massive Yang-Mills theory observables in the near Chern-Simons limit.


title: A note on morphisms determined by objects

abstract: We prove that a Hom-finite additive category having determined morphisms on both sides is a dualizing variety. This complements a result by Krause. We prove that in a Hom-finite abelian category having Serre duality, a morphism is right determined by some object if and only if it is an epimorphism. We give a characterization to abelian categories having Serre duality via determined morphisms.


title: A 30 mK, 13.5 T scanning tunneling microscope with two independent tips

abstract: We describe the design, construction, and performance of an ultra-low temperature, high-field scanning tunneling microscope (STM) with two independent tips. The STM is mounted on a dilution refrigerator and operates at a base temperature of 30 mK with magnetic fields of up to 13.5 T. We focus on the design of the two-tip STM head, as well as the sample transfer mechanism, which allows \textit{in situ} transfer from an ultra high vacuum (UHV) preparation chamber while the STM is at 1.5 K. Other design details such as the vibration isolation and rf-filtered wiring are also described. Their effectiveness is demonstrated via spectral current noise characteristics and the root mean square roughness of atomic resolution images. The high-field capability is shown by the magnetic field dependence of the superconducting gap of $\text{Cu}_x\text{Bi}_2\text{Se}_3$. Finally, we present images and spectroscopy taken with superconducting Nb tips with the refrigerator at 35 mK that indicate that the effective temperature of our tips/sample is approximately 184 mK, corresponding to an energy resolution of 16 $\mu$eV.


title: Submodularization for Quadratic Pseudo-Boolean Optimization

abstract: Many computer vision problems require optimization of binary non-submodular energies. We propose a general optimization framework based on local submodular approximations (LSA). Unlike standard LP relaxation methods that linearize the whole energy globally, our approach iteratively approximates the energies locally. On the other hand, unlike standard local optimization methods (e.g. gradient descent or projection techniques) we use non-linear submodular approximations and optimize them without leaving the domain of integer solutions. We discuss two specific LSA algorithms based on "trust region" and "auxiliary function" principles, LSA-TR and LSA-AUX. These methods obtain state-of-the-art results on a wide range of applications outperforming many standard techniques such as LBP, QPBO, and TRWS. While our paper is focused on pairwise energies, our ideas extend to higher-order problems. The code is available online (http://vision.csd.uwo.ca/code/).


title: General relativistic density perturbations

abstract: We investigate a general relativistic mechanism in which spikes generate matter overdensities in the early universe. When the cosmological fluid is tilted, the tilt provides another mechanism in generating matter inhomogeneities. We numerically investigate the effect of a sign change in the tilt, when there is a spike but the tilt does not change sign, and when the spike and the sign change in the tilt coincide. We find that the tilt plays the primary role in generating matter inhomogeneities, and it does so by creating both local overdensities and underdensities. We discuss of the physical implications of the work.


title: Theory of a Directive Optical Leaky Wave Antenna Integrated into a Resonator and Enhancement of Radiation Control

abstract: We provide for the first time the detailed study of the radiation performance of an optical leaky wave antenna (OLWA) integrated into a Fabry-P\'erot resonator. We show that the radiation pattern can be expressed as the one generated by the interference of two leaky waves counter-propagating in the resonator leading to a design procedure for achieving optimized broadside radiation, i.e., normal to the waveguide axis. We thus report a realizable implementation of the OLWA made of semiconductor and dielectric regions. The theoretical modeling is supported by full-wave simulation results, which are found to be in good agreement. We aim to control the radiation intensity in the broadside direction via excess carrier generation in the semiconductor regions. We show that the presence of the resonator can provide an effective way of enhancing the radiation level modulation, which reaches values as high as 13.5 dB, paving the way for novel promising control capabilities that might allow the generation of very fast optical switches, as an example.


title: Simple DFS on the Complement of a Graph and on Partially Complemented Digraphs

abstract: A complementation operation on a vertex of a digraph changes all outgoing arcs into non-arcs, and outgoing non-arcs into arcs. A partially complemented digraph $\widetilde{G}$ is a digraph obtained from a sequence of vertex complement operations on $G$. Dahlhaus et al. showed that, given an adjacency-list representation of $\widetilde{G}$, depth-first search (DFS) on $G$ can be performed in $O(n + \widetilde{m})$ time, where $n$ is the number of vertices and $\widetilde{m}$ is the number of edges in $\widetilde{G}$. To achieve this bound, their algorithm makes use of a somewhat complicated stack-like data structure to simulate the recursion stack, instead of implementing it directly as a recursive algorithm. We give a recursive $O(n+\widetilde{m})$ algorithm that uses no complicated data-structures.


title: Modeling of Atmospheric Flow Around a Coastal Cape: Lee Side Story

abstract: The current research focuses on mesoscale dynamics of the atmospheric circulation around an idealized coastal cape representing typical summertime circulation along the northwest coast of the U.S., studied using a mesoscale coupled ocean-atmosphere modeling system. The orographic wind maximum features a strong NW flow extending a few hundred kilometers downstream and seaward of the cape, which closely follows mesoscale orographic low pressure developed in the lee of the cape. Both wind maximum and the lee trough experience a pronounced diurnal cycle, marked by maximum northwest flow and minimum pressure in the local evening hours (its opposite phase during morning hours), and confirmed by observations from limited buoy and coastal stations. Vertical structure of the atmospheric boundary layer over the coastal ocean on the lee side of the cape indicated the downward propagation of potential temperature and wind features during the course of the day, as opposed to the traditional surface-driven development of the atmospheric boundary layer. Momentum analysis showed the local pressure gradients near the surface vary greatly depending on the relative location about the cape. Strong perturbations in relative vorticity downwind of the cape supported the hypothesis that mountain dynamics have a key role in establishing wind regime around the cape. The sensitivity studies indicated the importance of coastal terrain elevation on the formation and the strength of the lee trough and wind maximum. Higher coastal topography was found to cause greater separation of the upwelling front (jet) from the coast downwind of the cape.


title: A stochastic microscopic model for the dynamics of antigenic variation

abstract: We present a novel model that describes the within-host evolutionary dynamics of parasites undergoing antigenic variation. The approach uses a multi-type branching process with two types of entities defined according to their relationship with the immune system: clans of resistant parasitic cells (i.e. groups of cells sharing the same antigen not yet recognized by the immune system) that may become sensitive, and individual sensitive cells that can acquire a new resistance thus giving rise to the emergence of a new clan. The simplicity of the model allows analytical treatment to determine the subcritical and supercritical regimes in the space of parameters. By incorporating a density-dependent mechanism the model is able to capture additional relevant features observed in experimental data, such as the characteristic parasitemia waves. In summary our approach provides a new general framework to address the dynamics of antigenic variation which can be easily adapted to cope with broader and more complex situations.


title: A Polarization-insensitive and high-efficiency Schottky Photodetector integrated with a silicon ridge waveguide

abstract: We propose a polarization-insensitive and high-efficiency plasmonic silicon Schottky diode for detection of sub-bandgap photons in the optical communication wavelength range through internal photoemission. Our photodiode is based on a hybrid plasmonic silicon waveguide. It has a gold film covering both the top and the sidewalls of a dielectric silicon waveguide with the Schottky contact formed at the gold-silicon interface. An extensive physical model is presented in detail and applied to calculate and analyze the performance of our detector. By comparison with a diode with only top contact of gold, the polarization sensitivity of responsivity is greatly minimized in our photodetector with sidewall coverage of gold. Much higher responsivities for both polarizations are also achieved in a very broad wavelength range of 1.2-1.5 {\mu}m. Moreover, the Schottky contact is only 4 {\mu}m long, leading to a very small dark current. Our design is very promising for practical applications in high-density silicon photonic integration.


title: An invariant of rational homology 3-spheres via vector fields

abstract: We define an invariant of rational homology 3-spheres via vector fields. The construction of our invariant is a generalization of both that of the Kontsevich-Kuperberg-Thurston invariant and that of Watanabe's Morse homotopy invariant, which implies the equivalence of these two invariants.


title: Entropy Production in Simple Special Relativistic Fluids

abstract: It is well known that, in the absence of external forces, simple non-relativistic fluids involve entropy production only through heat conduction and shear viscosity. In this work, it is shown that a number density gradient contributes to the local entropy production of a simple relativistic fluid using special relativistic kinetic theory. Also, the presence of an external field may cause strictly relativistic contributions to the entropy production, a fact not widely recognized. The implications of these effects are thoroughly discussed.


title: The decomposability problem for torsion-free abelian groups is analytic complete

abstract: We discuss the decomposability of torsion-free abelian groups. We show that among computable groups of finite rank this property is $\Sigma^0_3$-complete. However, when we consider groups of infinite rank, it becomes $\Sigma^1_1$-complete, so it cannot be characterized by a first-order formula in the language of arithmetic.


title: Single Spin Asymmetries of Inclusive Hadrons Produced in Electron Scattering from a Transversely Polarized $^3$He Target

abstract: We report the first measurement of target single-spin asymmetries (A$_N$) in the inclusive hadron production reaction, $e~$+$~^3\text{He}^{\uparrow}\rightarrow h+X$, using a transversely polarized $^3$He target. The experiment was conducted at Jefferson Lab in Hall A using a 5.9-GeV electron beam. Three types of hadrons ($\pi^{\pm}$, $\text{K}^{\pm}$ and proton) were detected in the transverse hadron momentum range 0.54 $<p_T<$ 0.74 GeV/c. The range of $x_F$ for pions was -0.29 $<x_F<$ -0.23 and for kaons -0.25 $<x_F<$-0.18. The observed asymmetry strongly depends on the type of hadron. A positive asymmetry is observed for $\pi^+$ and $\text{K}^+$. A negative asymmetry is observed for $\pi^{-}$. The magnitudes of the asymmetries follow $|A^{\pi^-}| < |A^{\pi^+}| < |A^{K^+}|$. The K$^{-}$ and proton asymmetries are consistent with zero within the experimental uncertainties. The $\pi^{+}$ and $\pi^{-}$ asymmetries measured for the $^3$He target and extracted for neutrons are opposite in sign with a small increase observed as a function of $p_T$.


title: A New Discontinuous Galerkin Finite Element Method for Directly Solving the Hamilton-Jacobi Equations

abstract: In this paper, we improve upon the discontinuous Galerkin (DG) method for Hamilton-Jacobi (HJ) equation with convex Hamiltonians in (Y. Cheng and C.-W. Shu, J. Comput. Phys. 223:398-415,2007) and develop a new DG method for directly solving the general HJ equations. The new method avoids the reconstruction of the solution across elements by utilizing the Roe speed at the cell interface. Besides, we propose an entropy fix by adding penalty terms proportional to the jump of the normal derivative of the numerical solution. The particular form of the entropy fix was inspired by the Harten and Hyman's entropy fix (A. Harten and J. M. Hyman. J. Comput. Phys. 50(2):235-269, 1983) for Roe scheme for the conservation laws. The resulting scheme is compact, simple to implement even on unstructured meshes, and is demonstrated to work for nonconvex Hamiltonians. Benchmark numerical experiments in one dimension and two dimensions are provided to validate the performance of the method.


title: Quantum affine $\frak{gl}_n$ via Hecke algebras

abstract: We use the Hecke algebras of affine symmetric groups and their associated Schur algebras to construct a new algebra through a basis, and a set of generators and explicit multiplication formulas of basis elements by generators. We prove that this algebra is isomorphic to the quantum enveloping algebra of the loop algebra of $\mathfrak {gl}_n$. Though this construction is motivated by the work \cite{BLM} by Beilinson--Lusztig--MacPherson for quantum $\frak{gl}_n$, our approach is purely algebraic and combinatorial, independent of the geometric method which seems to work only for quantum $\mathfrak{gl}_n$ and quantum affine $\mathfrak{sl}_n$. As an application, we discover a presentation of the Ringel--Hall algebra of a cyclic quiver by semisimple generators and their multiplications by the defining basis elements.


title: Optimization, Learning, and Games with Predictable Sequences

abstract: We provide several applications of Optimistic Mirror Descent, an online learning algorithm based on the idea of predictable sequences. First, we recover the Mirror Prox algorithm for offline optimization, prove an extension to Holder-smooth functions, and apply the results to saddle-point type problems. Next, we prove that a version of Optimistic Mirror Descent (which has a close relation to the Exponential Weights algorithm) can be used by two strongly-uncoupled players in a finite zero-sum matrix game to converge to the minimax equilibrium at the rate of O((log T)/T). This addresses a question of Daskalakis et al 2011. Further, we consider a partial information version of the problem. We then apply the results to convex programming and exhibit a simple algorithm for the approximate Max Flow problem.


title: Mapping of the Influenza-A Hemagglutinin Serotypes Evolution by the ISSCOR Method

abstract: Analyses and visualizations by the ISSCOR method of influenza virus hemagglutinin genes of different A-subtypes revealed some rather striking temporal relationships between groups of individual gene subsets. Based on these findings we consider application of the ISSCOR-PCA method for analyses of large sets of homologous genes to be a worthwhile addition to a toolbox of genomics - allowing for a rapid diagnostics of trends, and ultimately even aiding an early warning of newly emerging epidemiological threats.


title: Harmonic analysis on homogeneous complex bounded domains and noncommutative geometry

abstract: We define and study a noncommutative Fourier transform on every homogeneous complex bounded domain. We then give an application in noncommutative differential geometry by defining noncommutative Baumslag-Solitar tori.


title: Optimality Conditions in Vector Optimization

abstract: In this paper we obtain second- and first-order optimality conditions of Kuhn-Tucker type and Fritz John one for weak efficiency in the vector problem with inequality constraints. In the necessary conditions we suppose that the objective function and the active constraints are continuously differentiable. We introduce notions of KTSP-invex problem and second-order KTSP-invex one. We obtain that the vector problem is (second-order) KTSP-invex if and only if for every triple $(\bar x,\bar\lambda,\bar\mu)$ with Lagrange multipliers $\bar\lambda$ and $\bar\mu$ for the objective function and constraints respectively, which satisfies the (second-order) necessary optimality conditions, the pair $(\bar x,\bar\mu)$ is a saddle point of the scalar Lagrange function with a fixed multiplier $\bar\lambda$. We introduce notions second-order KT-pseudoinvex-I, second-order KT-pseudo\-invex-II, second-order KT-invex problems. We prove that every second-order Kuhn-Tucker stationary point is a weak global Pareto minimizer (global Pareto minimizer) if and only if the problem is second-order KT-pseudoinvex-I (KT-pseudoinvex-II). It is derived that every second-order Kuhn Tucker stationary point is a global solution of the weighting problem if and only if the vector problem is second-order KT-invex.


title: An Asynchronous Parallel Stochastic Coordinate Descent Algorithm

abstract: We describe an asynchronous parallel stochastic coordinate descent algorithm for minimizing smooth unconstrained or separably constrained functions. The method achieves a linear convergence rate on functions that satisfy an essential strong convexity property and a sublinear rate ($1/K$) on general convex functions. Near-linear speedup on a multicore system can be expected if the number of processors is $O(n^{1/2})$ in unconstrained optimization and $O(n^{1/4})$ in the separable-constrained case, where $n$ is the number of variables. We describe results from implementation on 40-core processors.


title: Estimates on the Growth of Meromorphic Solutions of Linear Differential Equations with density conditions

abstract: We give an alternative and simpler method for getting pointwise estimate of meromorphic solutions of homogeneous linear differential equations with coefficients meromorphic in a finite disk or in the open plane originally obtained by Hayman and the author. In particular, our estimates generally give better upper bounds for higher order derivatives of the meromorphic solutions under consideration, are valid, however, outside an exceptional set of finite logarithmic density. The estimates again show that the growth of meromorphic solutions with a positive deficiency at infinity can be estimated in terms of initial conditions of the solution at or near the origin and the characteristic functions of the coefficients.


title: Large Negative Quantum Renormalization of Excitation Energies in the Spin-1/2 Kagome Lattice Antiferromagnet Cs$_2$Cu$_3$SnF$_{12}$

abstract: Magnetic excitations in the spin-$\frac{1}{2}$ distorted kagome lattice antiferromagnet Cs$_2$Cu$_3$SnF$_{12}$, which has an ordered ground state owing to the strong Dzyaloshinskii-Moriya interaction, were studied using inelastic neutron scattering. Although the spin-wave dispersion can be qualitatively understood in terms of linear spin-wave theory (LSWT), the excitation energies are renormalized by a factor of approximately 0.6 from those calculated by LSWT, almost irrespective of the momentum transfer. This inadequacy of LSWT, which is attributed to quantum fluctuations, provides evidence of negative quantum renormalization in the spin-$\frac{1}{2}$ kagome lattice antiferromagnet.


title: Mean Field Linear-Quadratic-Gaussian (LQG) Games of Forward-Backward Stochastic Differential Equations

abstract: This paper studies a new class of dynamic optimization problems of large-population (LP) system which consists of a large number of negligible and coupled agents. The most significant feature in our setup is the dynamics of individual agents follow the forward-backward stochastic differential equations (FBSDEs) in which the forward and backward states are coupled at the terminal time. This current paper is hence different to most existing large-population literature where the individual states are typically modeled by the SDEs including the forward state only. The associated mean-field linear-quadratic-Gaussian (LQG) game, in its forward-backward sense, is also formulated to seek the decentralized strategies. Unlike the forward case, the consistency conditions of our forward-backward mean-field games involve six Riccati and force rate equations. Moreover, their initial and terminal conditions are mixed thus some special decoupling technique is applied here. We also verify the $\epsilon$-Nash equilibrium property of the derived decentralized strategies. To this end, some estimates to backward stochastic system are employed. In addition, due to the adaptiveness requirement to forward-backward system, our arguments here are not parallel to those in its forward case.


title: The first, second and fourth Painlev\'{e} equations on weighted projective spaces

abstract: The first, second and fourth Painlev\'{e} equations are studied by means of dynamical systems theory and three dimensional weighted projective spaces $\C P^3(p,q,r,s)$ with suitable weights $(p,q,r,s)$ determined by the Newton diagrams of the equations or the versal deformations of vector fields. Singular normal forms of the equations, a simple proof of the Painlev\'{e} property and symplectic atlases of the spaces of initial conditions are given with the aid of the orbifold structure of $\C P^3(p,q,r,s)$. In particular, for the first Painlev\'{e} equation, a well known Painlev\'{e}'s transformation is geometrically derived, which proves to be the Darboux coordinates of a certain algebraic surface with a holomorphic symplectic form. The affine Weyl group, Dynkin diagram and the Boutroux coordinates are also studied from a view point of the weighted projective space.


title: The ground state and the character of the interaction between a colloidal particles in a liquid crystals

abstract: In this article is proposed the general approach to determination the character of the interaction between colloidal particles in a different liquid crystals. The main idea of this approach are in the presentation of the colloidal particle as sours of the possible deformation of the ground state of the director field. The ground state of liquid crystal imposes restrictions on the possible deformations and as result determine the character of the interaction between the colloidal particles. Based on this approach the Coulomb like interaction between dipole particles in a cholesteric liquid crystal and changing of the character of the interaction in a smactic liquid crystal has been predicted.


title: B-Machine Polarimeter: A Telescope to Measure the Polarization of the Cosmic Microwave Background

abstract: The B-Machine Telescope is the culmination of several years of development, construction, characterization and observation. The telescope is a departure from standard polarization chopping of correlation receivers to a half wave plate technique. Typical polarimeters use a correlation receiver to chop the polarization signal to overcome the $1/f$ noise inherent in HEMT amplifiers. B-Machine uses a room temperature half wave plate technology to chop between polarization states and measure the polarization signature of the CMB. The telescope has a demodulated $1/f$ knee of 5 mHz and an average sensitivity of 1.6 $\mathrm{mK}\sqrt{\mathrm{s}}$. This document examines the construction, characterization, observation of astronomical sources, and data set analysis of B-Machine. Preliminary power spectra and sky maps with large sky coverage for the first year data set are included.


title: Spectral Formulation of the Elastodynamic Boundary Integral Equations for Bi-material Interfaces

abstract: A spectral formulation of the plane-strain boundary integral equations for an interface between dissimilar elastic solids is presented. The boundary integral equations can be written in two equivalent forms: (a) The tractions can be written as a space-time convolution of the displacement continuities at the interface (Budiansky and Rice, 1979) (b) The displacement discontinuities can be written as a space-time convolution of the tractions at the interface (Kostrov, 1966). Prior work on spectral formulation of the boundary integral equations has adopted the former as the starting point. The present work has for its basis the latter form based on a space-time convolution of the tractions. Tractions and displacement components are given a spectral representation in the spatial coordinate along the interface. The radiation damping term is then explicitly extracted to avoid singularities in the convolution kernels. With the spectral forms introduced, the space-time convolutions reduce to convolutions in time for each Fourier mode. Due to continuity of tractions at the interface, this leads to a simpler formulation and form of the convolution kernels in comparison to the formulation involving convolutions over the slip and opening history at the bi-material interface. The convolution kernels are validated by studying some model problems to which analytical solutions are known. When coupled with a cohesive law or a friction law at the interface, the formulation proposed here is of wide applicability for studying spontaneous rupture propagation.


title: Spin wave nonreciprocity for logic device applications

abstract: The utilization of spin waves as eigenmodes of the magnetization dynamics for information processing and communication has been widely explored recently due to its high operational speed with low power consumption and possible applications for quantum computations. Previous proposals of spin wave Mach-Zehnder devices were based on the spin wave phase, a delicate entity which can be easily disrupted. Here, we propose a complete logic system based on the spin wave amplitude utilizing the nonreciprocal spin wave behavior excited by microstrip antennas. The experimental data reveal that the nonreciprocity of magnetostatic surface spin wave can be tuned by the bias magnetic field. Furthermore, engineering of the device structure could result in a high nonreciprocity factor for spin wave logic applications.


title: The complex singularity of a Stokes wave

abstract: Two-dimensional potential flow of the ideal incompressible fluid with free surface and infinite depth can be described by a conformal map of the fluid domain into the complex lower half-plane. Stokes wave is the fully nonlinear gravity wave propagating with the constant velocity. The increase of the scaled wave height $H/\lambda$ from the linear limit $H/\lambda=0$ to the critical value $H_{max}/\lambda$ marks the transition from the limit of almost linear wave to a strongly nonlinear limiting Stokes wave. Here $H$ is the wave height and $\lambda$ is the wavelength. We simulated fully nonlinear Euler equations, reformulated in terms of conformal variables, to find Stokes waves for different wave heights. Analyzing spectra of these solutions we found in conformal variables, at each Stokes wave height, the distance $v_c$ from the lowest singularity in the upper half-plane to the real line which corresponds to the fluid free surface. We also identified that this singularity is the square-root branch point. The limiting Stokes wave emerges as the singularity reaches the fluid surface. From the analysis of data for $v_c\to 0$ we suggest a new power law scaling $v_c\propto (H_{max}-H)^{3/2}$ as well as new estimate $H_{max}/\lambda \simeq 0.1410633$.

categories: physics.flu-dyn nlin.PS physics.ao-ph physics.comp-ph


title: Twist Accumulation and Topology Structure of a Solar Magnetic Flux Rope

abstract: To study the build up of a magnetic flux rope before a major flare and coronal mass ejection (CME), we compute the magnetic helicity injection, twist accumulation, and the topology structure of the three dimensional magnetic field, which is derived by the nonlinear force-free field model. The Extreme-ultraviolet Imaging Telescope on board the Solar and Heliospheric Observatory observed a series of confined flares without any CME before a major flare with a CME at 23:02 UT on 2005 January 15 in active region NOAA 10720. We derive the vector velocity at eight time points from 18:27 UT to 22:20 UT with the differential affine velocity estimator for vector magnetic fields, which were observed by the Digital Vector Magnetograph at Big Bear Solar Observatory. The injected magnetic helicity is computed with the vector magnetic and velocity fields. The helicity injection rate was (-16.47 \pm 3.52) \times 10^{40} Mx^2/hr. We find that only about 1.8% of the injected magnetic helicity became finally the internal helicity of the magnetic flux rope, whose twist increasing rate was -0.18 \pm 0.08 Turns/hr. The quasi-separatrix layers (QSLs) of the three dimensional magnetic field are computed by evaluating the squashing degree, Q. We find that the flux rope was wrapped by QSLs with large Q values, where the magnetic reconnection induced by the continuously injected magnetic helicity further produced the confined flares. We suggest that the flux rope was built up and heated by the magnetic reconnection in the QSLs.


title: A Parallel Simulated Annealing Approach for the Mirrored Traveling Tournament Problem

abstract: The Traveling Tournament Problem (TTP) is a benchmark problem in sports scheduling and has been extensively studied in recent years. The Mirrored Traveling Tournament Problem (mTTP) is variation of the TTP that represents certain types of sports scheduling problems where the main objective is to minimize the total distance traveled by all the participating teams. In this paper we test a parallel simulated annealing approach for solving the mTTP using OpenMP on shared memory systems and we found that this approach is superior especially with respect to the number of solution instances that are probed per second. We also see that there is significant speed up of 1.5x - 2.2x in terms of number of solutions explored per unit time.


title: Verifiable Control System Development for Gas Turbine Engines

abstract: A control software verification framework for gas turbine engines is developed. A stability proof is presented for gain scheduled closed-loop engine system based on global linearization and linear matrix inequality (LMI) techniques. Using convex optimization tools, a single quadratic Lyapunov function is computed for multiple linearizations near equilibrium points of the closed-loop system. With the computed stability matrices, ellipsoid invariant sets are constructed, which are used efficiently for DGEN turbofan engine control code stability analysis. Then a verifiable linear gain scheduled controller for DGEN engine is developed based on formal methods, and tested on the engine virtual test bench. Simulation results show that the developed verifiable gain scheduled controller is capable of regulating the engine in a stable fashion with proper tracking performance.


title: Thermodynamic and Logical Reversibilities Revisited

abstract: We review and investigate the general theory of thermodynamics of computation, and derive the fundamental inequalities that set the lower bounds of the work requirement and the heat emission during a computation. These inequalities constitute the generalized Landauer principle, where the information contents are involved in the second law of thermodynamics. We discuss in detail the relationship between the thermodynamic and logical reversibilities; the former is related to the entropy production in the total system including a heat bath, while the latter is related to the entropy change only in the logical states of the memory. In particular, we clarify that any logically irreversible computation can be performed in a thermodynamically reversible manner in the quasi-static limit, which does not contradict the conventional Landauer principle. Our arguments would serve as the theoretical foundation of thermodynamics of computation in terms of modern statistical physics.


title: Demand Side Management in Smart Grids using a Repeated Game Framework

abstract: Demand side management (DSM) is a key solution for reducing the peak-time power consumption in smart grids. To provide incentives for consumers to shift their consumption to off-peak times, the utility company charges consumers differential pricing for using power at different times of the day. Consumers take into account these differential prices when deciding when and how much power to consume daily. Importantly, while consumers enjoy lower billing costs when shifting their power usage to off-peak times, they also incur discomfort costs due to the altering of their power consumption patterns. Existing works propose stationary strategies for the myopic consumers to minimize their short-term billing and discomfort costs. In contrast, we model the interaction emerging among self-interested, foresighted consumers as a repeated energy scheduling game and prove that the stationary strategies are suboptimal in terms of long-term total billing and discomfort costs. Subsequently, we propose a novel framework for determining optimal nonstationary DSM strategies, in which consumers can choose different daily power consumption patterns depending on their preferences, routines, and needs. As a direct consequence of the nonstationary DSM policy, different subsets of consumers are allowed to use power in peak times at a low price. The subset of consumers that are selected daily to have their joint discomfort and billing costs minimized is determined based on the consumers' power consumption preferences as well as on the past history of which consumers have shifted their usage previously. Importantly, we show that the proposed strategies are incentive-compatible. Simulations confirm that, given the same peak-to-average ratio, the proposed strategy can reduce the total cost (billing and discomfort costs) by up to 50% compared to existing DSM strategies.


title: D$^3$PO - Denoising, Deconvolving, and Decomposing Photon Observations

abstract: The analysis of astronomical images is a non-trivial task. The D3PO algorithm addresses the inference problem of denoising, deconvolving, and decomposing photon observations. Its primary goal is the simultaneous but individual reconstruction of the diffuse and point-like photon flux given a single photon count image, where the fluxes are superimposed. In order to discriminate between these morphologically different signal components, a probabilistic algorithm is derived in the language of information field theory based on a hierarchical Bayesian parameter model. The signal inference exploits prior information on the spatial correlation structure of the diffuse component and the brightness distribution of the spatially uncorrelated point-like sources. A maximum a posteriori solution and a solution minimizing the Gibbs free energy of the inference problem using variational Bayesian methods are discussed. Since the derivation of the solution is not dependent on the underlying position space, the implementation of the D3PO algorithm uses the NIFTY package to ensure applicability to various spatial grids and at any resolution. The fidelity of the algorithm is validated by the analysis of simulated data, including a realistic high energy photon count image showing a 32 x 32 arcmin^2 observation with a spatial resolution of 0.1 arcmin. In all tests the D3PO algorithm successfully denoised, deconvolved, and decomposed the data into a diffuse and a point-like signal estimate for the respective photon flux components.

categories: astro-ph.IM cs.IT math.IT physics.data-an stat.CO


title: Configurable unitary transformations and linear logic gates using quantum memories

abstract: We show that a set of optical memories can act as a configurable linear optical network operating on frequency-multiplexed optical states. Our protocol is applicable to any quantum memories that employ off-resonant Raman transitions to store optical information in atomic spins. In addition to the configurability, the protocol also offers favourable scaling with an increasing number of modes where N memories can be configured to implement an arbitrary N-mode unitary operations during storage and readout. We demonstrate the versatility of this protocol by showing an example where cascaded memories are used to implement a conditional CZ gate.


title: Discrepancy estimates for variance bounding Markov chain quasi-Monte Carlo

abstract: Markov chain Monte Carlo (MCMC) simulations are modeled as driven by true random numbers. We consider variance bounding Markov chains driven by a deterministic sequence of numbers. The star-discrepancy provides a measure of efficiency of such Markov chain quasi-Monte Carlo methods. We define a pull-back discrepancy of the driver sequence and state a close relation to the star-discrepancy of the Markov chain-quasi Monte Carlo samples. We prove that there exists a deterministic driver sequence such that the discrepancies decrease almost with the Monte Carlo rate $n^{1/2}$. As for MCMC simulations, a burn-in period can also be taken into account for Markov chain quasi-Monte Carlo to reduce the influence of the initial state. In particular, our discrepancy bound leads to an estimate of the error for the computation of expectations. To illustrate our theory we provide an example for the Metropolis algorithm based on a ball walk. Furthermore, under additional assumptions we prove the existence of a driver sequence such that the discrepancy of the corresponding deterministic Markov chain sample decreases with order $n^{-1+\delta}$ for every $\delta>0$.


title: On cubic birational maps of $\mathbb{P}^3_\mathbb{C}$

abstract: We study the birational maps of $\mathbb{P}^3_\mathbb{C}$. More precisely we describe the irreducible components of the set of birational maps of bidegree $(3,3)$ (resp. $(3,4)$, resp. $(3,5)$).


title: Evolutionary genetic optimization of the injector beam dynamics for the ERL test facility at IHEP

abstract: The energy recovery linac test facility (ERL-TF), a compact ERL-FEL (free electron laser) two-purpose machine, was proposed at the Institute of High Energy Physics, Beijing. As one important component of the ERL-TF, the photo-injector started with a photocathode direct-current gun was designed and preliminarily optimized. In this paper an evolutionary genetic method, non-dominated sorting genetic algorithm II, is applied to optimize the injector beam dynamics, especially in the high-charge operation mode. Study shows that using an incident laser with rms transverse size of 1~1.2 mm, the normalized emittance of the electron beam can be kept below 1 mm.mrad at the end of the injector. This work, together with the previous optimization for the low-charge operation mode by using the iterative scan method, provides guidance and confidence for future constructing and commissioning of the ERL-TF injector.


title: Semiparametric testing of statistical functionals revisited

abstract: Along the lines of Janssen's and Pfanzagl's work the testing theory for statistical functionals is further developed for non-parametric one-sample problems. Efficient tests for the one-sided and two-sided problems are derived for nonparametric statistical functionals. The asymptotic power function is calculated under implicit alternatives and hypotheses, which are given by the functional itself, for the one- sided and two-sided cases. Under mild regularity assumptions it is shown that these tests are asymptotic most powerful. The combination of the modern theory of Le Cam and approximation in limit experiments provide a deep insight into the upper bounds for asymptotic power functions of tests for the one-sided and two-sided problems. As example tests based on the von Mises functional are treated in nonparametric context.


title: Optimality of Thompson Sampling for Gaussian Bandits Depends on Priors

abstract: In stochastic bandit problems, a Bayesian policy called Thompson sampling (TS) has recently attracted much attention for its excellent empirical performance. However, the theoretical analysis of this policy is difficult and its asymptotic optimality is only proved for one-parameter models. In this paper we discuss the optimality of TS for the model of normal distributions with unknown means and variances as one of the most fundamental example of multiparameter models. First we prove that the expected regret of TS with the uniform prior achieves the theoretical bound, which is the first result to show that the asymptotic bound is achievable for the normal distribution model. Next we prove that TS with Jeffreys prior and reference prior cannot achieve the theoretical bound. Therefore the choice of priors is important for TS and non-informative priors are sometimes risky in cases of multiparameter models.