title: A Semidefinite Programming Based Search Strategy for Feature Selection with Mutual Information Measure

abstract: Feature subset selection, as a special case of the general subset selection problem, has been the topic of a considerable number of studies due to the growing importance of data-mining applications. In the feature subset selection problem there are two main issues that need to be addressed: (i) Finding an appropriate measure function than can be fairly fast and robustly computed for high-dimensional data. (ii) A search strategy to optimize the measure over the subset space in a reasonable amount of time. In this article mutual information between features and class labels is considered to be the measure function. Two series expansions for mutual information are proposed, and it is shown that most heuristic criteria suggested in the literature are truncated approximations of these expansions. It is well-known that searching the whole subset space is an NP-hard problem. Here, instead of the conventional sequential search algorithms, we suggest a parallel search strategy based on semidefinite programming (SDP) that can search through the subset space in polynomial time. By exploiting the similarities between the proposed algorithm and an instance of the maximum-cut problem in graph theory, the approximation ratio of this algorithm is derived and is compared with the approximation ratio of the backward elimination method. The experiments show that it can be misleading to judge the quality of a measure solely based on the classification accuracy, without taking the effect of the non-optimum search strategy into account.


title: Newtonian View of General Relativistic Stars

abstract: Although general relativistic cosmological solutions, even in the presence of pressure, can be mimicked by using neo-Newtonian hydrodynamics, it is not clear whether there exists the same Newtonian correspondence for spherical static configurations. General relativity solutions for stars are known as the Tolman-Oppenheimer-Volkoff (TOV) equations. On the other hand, the Newtonian description does not take into account the total pressure effects and therefore can not be used in strong field regimes. We discuss how to incorporate pressure in the stellar equilibrium equations within the neo-Newtonian framework. We compare the Newtonian, neo-Newtonian and the full relativistic theory by solving the equilibrium equations for both three approaches and calculating the mass-radius diagrams for some simple neutron stars equation of state.


title: Performance of Stanford and Minipar Parser on Biomedical Texts

abstract: In this paper, the performance of two dependency parsers, namely Stanford and Minipar, on biomedical texts has been reported. The performance of te parsers to assignm dependencies between two biomedical concepts that are already proved to be connected is not satisfying. Both Stanford and Minipar, being statistical parsers, fail to assign dependency relation between two connected concepts if they are distant by at least one clause. Minipar's performance, in terms of precision, recall and the F-score of the attachment score (e.g., correctly identified head in a dependency), to parse biomedical text is also measured taking the Stanford's as a gold standard. The results suggest that Minipar is not suitable yet to parse biomedical texts. In addition, a qualitative investigation reveals that the difference between working principles of the parsers also play a vital role for Minipar's degraded performance.


title: Weather data analysis based on typical weather sequence analysis. Application: energy building simulation

abstract: In building studies dealing about energy efficiency and comfort, simulation software need relevant weather files with optimal time steps. Few tools generate extreme and mean values of simultaneous hourly data including correlation between the climatic parameters. This paper presents the C++ Runeole software based on typical weather sequences analysis. It runs an analysis process of a stochastic continuous multivariable phenomenon with frequencies properties applied to a climatic database. The database analysis associates basic statistics, PCA (Principal Component Analysis) and automatic classifications. Different ways of applying these methods will be presented. All the results are stored in the Runeole internal database that allows an easy selection of weather sequences. The extreme sequences are used for system and building sizing and the mean sequences are used for the determination of the annual cooling loads as proposed by Audrier-Cros (Audrier-Cros, 1984). This weather analysis was tested with the database of the French weather forecast utility Meteo France. Reunion Island experiences a lot of different micro-climates due to the high altitude (3069m), specific relief, and geographic situation (Tropic of Capricorn). Furthermore Reunion Island has the densest meteorological network in France and is an ideal place to validate the methodology with different climates. To test the efficiency of such analysis, simulations using the resulting weather sequences were carried out with the building simulation software CODYRUN. This analysis is the first step of a more global research concerning weather data generation. Future work will permit whole hourly typical meteorological year generation using neural networks.


title: The star formation history of galaxies in 3D: CALIFA perspective

abstract: We resolve spatially the star formation history of 300 nearby galaxies from the CALIFA integral field survey to investigate: a) the radial structure and gradients of the present stellar populations properties as a function of the Hubble type; and b) the role that plays the galaxy stellar mass and stellar mass surface density in governing the star formation history and metallicity enrichment of spheroids and the disks of galaxies. We apply the fossil record method based on spectral synthesis techniques to recover spatially and temporally resolved maps of stellar population properties of spheroids and spirals with galaxy mass from 10$^9$ to 7$\times$10$^{11}$ M$_{\odot}$. The individual radial profiles of the stellar mass surface density ($\mu_{*}$), stellar extinction (A$_{V}$), luminosity weighted ages ($\langle$ log age $\rangle_{L}$), and mass weighted metallicity ($\langle$ log Z/Z$_{\odot}$$\rangle_{M}$) are stacked in seven bins of galaxy morphology (E, S0, Sa, Sb, Sbc, Sc and Sd). All these properties show negative gradients as a sight of the inside-out growth of massive galaxies. However, the gradients depend on the Hubble type in different ways. For the same galaxy mass, E and S0 galaxies show the largest inner gradients in $\mu_{*}$; and Andromeda-like galaxies (Sb with log M$_{*}$(M$_{\odot}$) $\sim$ 11) show the largest inner age and metallicity gradients. In average, spiral galaxies have a stellar metallicity gradient $\sim$ -0.1 dex per half-light radius, in agreement with the value estimated for the ionized gas oxygen abundance gradient by CALIFA. A global (M$_{*}$-driven) and local ($\mu_{*}$- driven) stellar metallicity relation are derived. We find that in disks, the stellar mass surface density regulates the stellar metallicity; in spheroids, the galaxy stellar mass dominates the physics of star formation and chemical enrichment.


title: Modelling The Redshift-Space Three-Point Correlation Function in SDSS-III

abstract: We present the measurements of the redshift-space three-point correlation function (3PCF) for z~0.5 luminous red galaxies of the CMASS sample in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Data Release 11. The 3PCF measurements are interpreted within the halo occupation distribution framework using high-resolution N-body simulations, and the model successfully reproduces the 3PCF on scales larger than 1Mpc/h. As with the case for the redshift-space two-point correlation functions, we find that the redshift-space 3PCF measurements also favour the inclusion of galaxy velocity bias in the model. In particular, the central galaxy in a halo is on average in motion with respect to the core of the halo. We discuss the potential of the small-scale 3PCF to tighten the constraints on the relation between galaxies and dark matter haloes and on the phase-space distribution of galaxies.


title: A search for stellar tidal debris of defunct dwarf galaxies around globular clusters in the inner Galactic halo

abstract: In the hierarchical formation scenario in which the outer halo of the Milky Way is the result of the continuous accretion of low-mass galaxies, a fraction of the Galactic globular cluster system might have originated in and been accreted with already extinct dwarf galaxies. In this context, we expect that the remnants of these progenitor galaxies might be still populating the surroundings of those accreted globulars. In this work, we present wide-field photometry of a sample of 23 globular clusters in the Galactocentric distance range 10 < Rg < 40kpc, which we use to search for remnants of their hypothetical progenitor systems. Our deep photometry reveals the presence of underlying stellar populations along the line-of-sight of about half of the globulars included in our sample. Among the detections lying in the footprint of the Sagittarius tidal stream, which we identify via the comparison with its orbit derived from numerical simulations, only Whiting1 and NGC7492 seem to be inmersed in that remnant at a compatible heliocentric distance. We also confirm the existence of a subjacent Main-Sequence feature in the surroundings of NGC1851. A tentative detection of the vast Hercules-Aquila cloud is unveiled in the background of NGC7006.


title: The Unnatural Composite Higgs

abstract: Composite Higgs models can trivially satisfy precision-electroweak and flavour constraints by simply having a large spontaneous symmetry breaking scale, f > 10 TeV. This produces a 'split' spectrum, where the strong sector resonances have masses greater than 10 TeV and are separated from the pseudo Nambu-Goldstone bosons, which remain near the electroweak scale. Even though a tuning of order 10^{-4} is required to obtain the observed Higgs boson mass, the big hierarchy problem remains mostly solved. Intriguingly, models with a fully-composite right-handed top quark also exhibit improved gauge coupling unification. By restricting ourselves to models which preserve these features we find that the symmetry breaking scale cannot be arbitrarily raised, leading to an upper bound f < 100-1000 TeV. This implies that the resonances may be accessible at future colliders, or indirectly via rare-decay experiments. Dark matter is identified with a pseudo Nambu-Goldstone boson, and we show that the smallest coset space containing a stable, scalar singlet and an unbroken SU(5) symmetry is SU(7) / SU(6) x U(1). The colour-triplet pseudo Nambu-Goldstone boson also contained in this coset space is metastable due to a residual symmetry. It can decay via a displaced vertex when produced at colliders, leading to a distinctive signal of unnaturalness.


title: Soft Collinear Effective Theory for Heavy WIMP Annihilation

abstract: In a large class of models for Weakly Interacting Massive Particles (WIMPs), the WIMP mass $M$ lies far above the weak scale $m_W$. This work identifies universal Sudakov-type logarithms $\sim \alpha \log^2 (2\,M/m_W)$ that spoil the naive convergence of perturbation theory for annihilation processes. An effective field theory (EFT) framework is presented, allowing the systematic resummation of these logarithms. Another impact of the large separation of scales is that a long-distance wave-function distortion from electroweak boson exchange leads to observable modifications of the cross section. Careful accounting of momentum regions in the EFT allows the rigorously disentanglement of this so-called Sommerfeld enhancement from the short distance hard annihilation process. The WIMP is modeled as a heavy-particle field, while the light, energetic, final-state electroweak gauge bosons are treated as soft and collinear fields. Hard matching coefficients are computed at renormalization scale $\mu \sim 2\,M$, then evolved down to $\mu \sim m_W$, where electroweak symmetry breaking is incorporated and the matching onto the relevant quantum mechanical Hamiltonian is performed. The example of an $SU(2)_W$ triplet scalar dark matter candidate annihilating to line photons is used for concreteness, allowing the numerical exploration of the impact of next-to-leading order corrections and log resummation. For $M \simeq 3$ TeV, the resummed Sommerfeld enhanced cross section is reduced by a factor of $\sim 3$ with respect to the tree-level fixed order result.


title: Non-Thermal Electron Acceleration in Low Mach Number Collisionless Shocks. II. Firehose-Mediated Fermi Acceleration and its Dependence on Pre-Shock Conditions

abstract: Electron acceleration to non-thermal energies is known to occur in low Mach number (M<5) shocks in galaxy clusters and solar flares, but the electron acceleration mechanism remains poorly understood. Using two-dimensional (2D) particle-in-cell (PIC) plasma simulations, we showed in Paper I that electrons are efficiently accelerated in low Mach number (M=3) quasi-perpendicular shocks via a Fermi-like process. The electrons bounce between the upstream region and the shock front, with each reflection at the shock resulting in energy gain via shock drift acceleration. The upstream scattering is provided by oblique magnetic waves, that are self-generated by the electrons escaping ahead of the shock. In the present work, we employ additional 2D PIC simulations to address the nature of the upstream oblique waves. We find that the waves are generated by the shock-reflected electrons via the firehose instability, which is driven by an anisotropy in the electron velocity distribution. We systematically explore how the efficiency of wave generation and of electron acceleration depend on the magnetic field obliquity, the flow magnetization (or equivalently, the plasma beta), and the upstream electron temperature. We find that the mechanism works for shocks with high plasma beta (>20) at nearly all magnetic field obliquities, and for electron temperatures in the range relevant for galaxy clusters. Our findings offer a natural solution to the conflict between the bright radio synchrotron emission observed from the outskirts of galaxy clusters and the low electron acceleration efficiency usually expected in low Mach number shocks.


title: A resolved, au-scale gas disk around the B[e] star HD 50138

abstract: HD 50138 is a B[e] star surrounded by a large amount of circumstellar gas and dust. Its spectrum shows characteristics which may indicate either a pre- or a post-main-sequence system. Mapping the kinematics of the gas in the inner few au of the system contributes to a better understanding of its physical nature. We present the first high spatial and spectral resolution interferometric observations of the Br-gamma line of HD~50138, obtained with VLTI/AMBER. The line emission originates from a region more compact (up to 3 au) than the continuum-emitting region. Blue- and red-shifted emission originates from the two different hemispheres of an elongated structure perpendicular to the polarization angle. The velocity of the emitting medium decreases radially. An overall offset along the NW direction between the line- and continuum-emitting regions is observed. We compare the data with a geometric model of a thin Keplerian disk and a spherical halo on top of a Gaussian continuum. Most of the data are well reproduced by this model, except for the variability, the global offset and the visibility at the systemic velocity. The evolutionary state of the system is discussed; most diagnostics are ambiguous and may point either to a post-main-sequence or a pre-main-sequence nature.


title: GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods

abstract: We present two new Lagrangian methods for hydrodynamics, in a systematic comparison with moving-mesh, SPH, and stationary (non-moving) grid methods. The new methods are designed to simultaneously capture advantages of both smoothed-particle hydrodynamics (SPH) and grid-based/adaptive mesh refinement (AMR) schemes. They are based on a kernel discretization of the volume coupled to a high-order matrix gradient estimator and a Riemann solver acting over the volume 'overlap.' We implement and test a parallel, second-order version of the method with self-gravity & cosmological integration, in the code GIZMO: this maintains exact mass, energy and momentum conservation; exhibits superior angular momentum conservation compared to all other methods we study; does not require 'artificial diffusion' terms; and allows the fluid elements to move with the flow so resolution is automatically adaptive. We consider a large suite of test problems, and find that on all problems the new methods appear competitive with moving-mesh schemes, with some advantages (particularly in angular momentum conservation), at the cost of enhanced noise. The new methods have many advantages vs. SPH: proper convergence, good capturing of fluid-mixing instabilities, dramatically reduced 'particle noise' & numerical viscosity, more accurate sub-sonic flow evolution, & sharp shock-capturing. Advantages vs. non-moving meshes include: automatic adaptivity, dramatically reduced advection errors & numerical overmixing, velocity-independent errors, accurate coupling to gravity, good angular momentum conservation and elimination of 'grid alignment' effects. We can, for example, follow hundreds of orbits of gaseous disks, while AMR and SPH methods break down in a few orbits. However, fixed meshes minimize 'grid noise.' These differences are important for a range of astrophysical problems.

categories: astro-ph.CO astro-ph.GA astro-ph.IM physics.comp-ph physics.flu-dyn


title: Galaxy evolution across the optical emission-line diagnostic diagrams?

abstract: The discovery of the M-sigma relation, the local galaxy bimodality, and the link between black-hole and host-galaxy properties, have raised the question whether AGN play a role in galaxy evolution. Several theoretical models implement AGN feedback to explain the observed galaxy luminosity function, and possibly the color and morphological transformation of spiral galaxies into passive ellipticals. To understand the importance of AGN feedback, a study of the AGN populations in the radio-optical domain is crucial. A mass sequence linking star-forming galaxies and AGN has been already noted in previous works, and it is now investigated as possible evolutionary sequence. We observed a sample of 119 intermediate-redshift (0.04<z<0.4) SDSS-FIRST radio emitters with the Effelsberg 100-m telescope at 4.85 and 10.45 GHz and obtained spectral indices. We find indications of spectral index flattening in high-metallicity star-forming galaxies, composite galaxies, and Seyferts. This "flattening sequence" along the [NII]-based emission-line diagnostic diagram is consistent with the hardening of galaxy ionizing field, due to nuclear activity. After combining our data with FIRST measurements at 1.4 GHz, we find that the three-point radio spectra of Seyferts and LINERs show substantial differences, attributable to small radio core components and larger (arcsecond sized) jet/lobe components, respectively. A visual inspection of FIRST images seems to confirm this hypothesis. Galaxies along this sequence are hypothesized to be transitioning from the active star-forming galaxies (blue cloud) to the passive elliptical galaxies (red sequence). This supports the suggestion that AGN play a role in shutting down star-formation, and allow the transition from one galaxy class to the other.


title: Quantitative simulation of temperature dependent magnetization dynamics and equilibrium properties of elemental ferromagnets

abstract: Atomistic spin model simulations are immensely useful in determining temperature dependent magnetic prop- erties, but are known to give the incorrect dependence of the magnetization on temperature compared to exper- iment owing to their classical origin. We find a single parameter rescaling of thermal fluctuations which gives quantitative agreement of the temperature dependent magnetization between atomistic simulations and experi- ment for the elemental ferromagnets Ni, Fe, Co and Gd. Simulating the sub-picosecond magnetization dynam- ics of Ni under the action of a laser pulse we also find quantitative agreement with experiment in the ultrafast regime. This enables the quantitative determination of temperature dependent magnetic properties allowing for accurate simulations of magnetic materials at all temperatures.


title: The ancient heritage of water ice in the solar system

abstract: Identifying the source of Earth's water is central to understanding the origins of life-fostering environments and to assessing the prevalence of such environments in space. Water throughout the solar system exhibits deuterium-to-hydrogen enrichments, a fossil relic of low-temperature, ion-derived chemistry within either (i) the parent molecular cloud or (ii) the solar nebula protoplanetary disk. Utilizing a comprehensive treatment of disk ionization, we find that ion-driven deuterium pathways are inefficient, curtailing the disk's deuterated water formation and its viability as the sole source for the solar system's water. This finding implies that if the solar system's formation was typical, abundant interstellar ices are available to all nascent planetary systems.


title: Backreaction and Stochastic Effects in Single Field Inflation

abstract: The formalism of stochastic inflation is a powerful tool for analyzing the backreaction of cosmological perturbations, and making precise predictions for inflationary observables. We demonstrate this with the simple $m^2\phi^2$ model of inflation, wherein we obtain an effective field theory for IR modes of the inflaton, which remains coupled to UV modes through a classical noise. We compute slow-roll corrections to the evolution of UV modes (i.e. quantum fluctuations), and track this effect from the UV theory to the IR theory, where it manifests as a correction to the classical noise. We compute the stochastic correction to the spectral index of primordial perturbations, finding a small effect, and discuss models in which this effect can become large. We extend our analysis to tensor modes, and demonstrate that the stochastic approach allows us to recover the standard tensor tilt $n_T$, plus corrections.


title: Spectral synthesis of stellar populations in the 3D era: The CALIFA experience

abstract: Methods to recover the fossil record of galaxy evolution encoded in their optical spectra have been instrumental in processing the avalanche of data from mega-surveys along the last decade, effectively transforming observed spectra onto a long and rich list of physical properties: from stellar masses and mean ages to full star formation histories. This promoted progress in our understanding of galaxies as a whole. Yet, the lack of spatial resolution introduces undesirable aperture effects, and hampers advances on the internal physics of galaxies. This is now changing with 3D surveys. The mapping of stellar populations in data-cubes allows us to figure what comes from where, unscrambling information previously available only in integrated form. This contribution uses our starlight-based analysis of 300 CALIFA galaxies to illustrate the power of spectral synthesis applied to data-cubes. The selected results highlighted here include: (a) The evolution of the mass-metallicity and mass-density-metallicity relations, as traced by the mean stellar metallicity. (b) A comparison of star formation rates obtained from H{\alpha} to those derived from full spectral fits. (c) The relation between star formation rate and dust optical depth within galaxies, which turns out to mimic the Schmidt-Kennicutt law. (d) PCA tomography experiments.


title: Exotic Quantum Phase Transitions of $(2+1)d$ Dirac fermions

abstract: Using determinant quantum Monte Carlo (d-QMC) simulations, we demonstrate that an extended Hubbard model on a bilayer honeycomb lattice has two novel quantum phase transitions. The first is a quantum phase transition between the weakly interacting gapless Dirac fermion phase and a strongly interacting fully gapped and symmetric trivial phase, which cannot be described by the standard Gross-Neveu model. The second is a quantum critical point between a quantum spin Hall insulator with spin $S^z$ conservation and the previously mentioned strongly interacting fully gapped phase. At the latter quantum critical point the single particle excitations remain gapped, while spin and charge gap both close. We argue that the first quantum phase transition is related to the $\mathbb{Z}_{16}$ classification of the topological superconductor $^3\text{He-B}$ phase with interactions, while the second quantum phase transition is a topological phase transition described by a bosonic O(4) nonlinear sigma model field theory with a $\Theta$-term.


title: A Forecast for a South Heliopause Crossing by Voyager 2 in Late 2014 Using Intensity-time Features of Energetic Particles Observed by V1 and V2 in the North and South Heliosheaths

abstract: We have used corresponding marker features in the intensity vs. time profile of termination shock particles coming from within the heliosheath and also 5-12 MeV electrons which are coming from beyond the heliopause that are both observed at both V1 and V2 in the outer N and S heliosheaths, to estimate the thickness of the S heliosheath. These features in the N heliosheath have been observed on the CRS instrument on V1 in the N heliosheath beyond 110 AU or within 12 AU of the heliopause. Although the intensity-time features in the S heliosheath are quite different overall, certain features stand out at both V1 and V2. The correspondence of these features in the N and S heliosheath allows us to estimate that V2 will reach the S heliopause near the end of 2014.


title: Optimal high-level descriptions of dynamical systems

abstract: To analyze high-dimensional systems, many fields in science and engineering rely on high-level descriptions, sometimes called "macrostates," "coarse-grainings," or "effective theories". Examples of such descriptions include the thermodynamic properties of a large collection of point particles undergoing reversible dynamics, the variables in a macroeconomic model describing the individuals that participate in an economy, and the summary state of a cell composed of a large set of biochemical networks. Often these high-level descriptions are constructed without considering the ultimate reason for needing them in the first place. Here, we formalize and quantify one such purpose: the need to predict observables of interest concerning the high-dimensional system with as high accuracy as possible, while minimizing the computational cost of doing so. The resulting State Space Compression (SSC) framework provides a guide for how to solve for the {optimal} high-level description of a given dynamical system, rather than constructing it based on human intuition alone. In this preliminary report, we introduce SSC, and illustrate it with several information-theoretic quantifications of "accuracy", all with different implications for the optimal compression. We also discuss some other possible applications of SSC beyond the goal of accurate prediction. These include SSC as a measure of the complexity of a dynamical system, and as a way to quantify information flow between the scales of a system.

categories: cs.IT cond-mat.stat-mech cs.AI cs.CE math.IT q-bio.PE


title: Time-resolved infrared emission from radiation-driven central obscuring structures in Active Galactic Nuclei

abstract: The central engines of Seyfert galaxies are thought to be enshrouded by geometrically thick gas and dust structures. In this article, we derive observable properties for a self-consistent model of such toroidal gas and dust distributions, where the geometrical thickness is achieved and maintained with the help of X-ray heating and radiation pressure due to the central engine. Spectral energy distributions (SEDs) and images are obtained with the help of dust continuum radiative transfer calculations with RADMC-3D. For the first time, we are able to present time-resolved SEDs and images for a physical model of the central obscurer. Temporal changes are mostly visible at shorter wavelengths, close to the combined peak of the dust opacity as well as the central source spectrum and are caused by variations in the column densities of the generated outflow. Due to the three-component morphology of the hydrodynamical models -- a thin disc with high density filaments, a surrounding fluffy component (the obscurer) and a low density outflow along the rotation axis -- we find dramatic differences depending on wavelength: whereas the mid-infrared images are dominated by the elongated appearance of the outflow cone, the long wavelength emission is mainly given by the cold and dense disc component. Overall, we find good agreement with observed characteristics, especially for those models, which show clear outflow cones in combination with a geometrically thick distribution of gas and dust, as well as a geometrically thin, but high column density disc in the equatorial plane.


title: Black Hole Formation and Classicalization in Ultra-Planckian 2 -> N Scattering

abstract: We establish a connection between the ultra-Planckian scattering amplitudes in field and string theory and unitarization by black hole formation in these scattering processes. Using as a guideline an explicit microscopic theory in which the black hole represents a bound-state of many soft gravitons at the quantum critical point, we were able to identify and compute a set of perturbative amplitudes relevant for black hole formation. These are the tree-level N-graviton scattering S-matrix elements in a kinematical regime (called classicalization limit) where the two incoming ultra-Planckian gravitons produce a large number N of soft gravitons. We compute these amplitudes by using the Kawai-Lewellen-Tye relations, as well as scattering equations and string theory techniques. We discover that this limit reveals the key features of the microscopic corpuscular black hole N-portrait. In particular, the perturbative suppression factor of a N-graviton final state, derived from the amplitude, matches the non-perturbative black hole entropy when N reaches the quantum criticality value, whereas final states with different value of N are either suppressed or excluded by non-perturbative corpuscular physics. Thus we identify the microscopic reason behind the black hole dominance over other final states including non-black hole classical object. In the parameterization of the classicalization limit the scattering equations can be solved exactly allowing us to obtain closed expressions for the high-energy limit of the open and closed superstring tree-level scattering amplitudes for a generic number N of external legs. We demonstrate matching and complementarity between the string theory and field theory in different large-s and large-N regimes.


title: Supersymmetry and non-Abelian T-duality in type II supergravity

abstract: We study the effect of T-duality on supersymmetry in the context of type II supergravity. For both U(1) Abelian and SU(2) non-Abelian T-duality, we demonstrate that the supersymmetry variations after T-duality are related to the variations before T-duality through the Kosmann spinorial Lie derivative, which vanishes when the Killing spinors are independent of the T-duality directions. As a byproduct of our analysis, we present closed expressions for SU(2) T-duality in a class of spacetimes with diagonal Bianchi IX symmetry and comment on specific examples of T-dual geometries, including a novel AdS3 geometry with large N = (0,4) superconformal symmetry.


title: Dividing and weak quasidimensions in arbitary theories

abstract: We show that any countable model of a model complete theory has an elementary extension with a "pseudofinite-like" quasidimension that detects dividing.


title: LP-decodable multipermutation codes

abstract: In this paper, we introduce a new way of constructing and decoding multipermutation codes. Multipermutations are permutations of a multiset that may consist of duplicate entries. We first introduce a new class of matrices called multipermutation matrices. We characterize the convex hull of multipermutation matrices. Based on this characterization, we propose a new class of codes that we term LP-decodable multipermutation codes. Then, we derive two LP decoding algorithms. We first formulate an LP decoding problem for memoryless channels. We then derive an LP algorithm that minimizes the Chebyshev distance. Finally, we show a numerical example of our algorithm.


title: Generalized tight p-frames and spectral bounds for Laplace-like operators

abstract: We prove sharp upper bounds for sums of eigenvalues (and other spectral functionals) of Laplace-like operators, including bi-Laplacian and fractional Laplacian. We show that among linear images of a highly symmetric domain, our spectral functionals are maximal on the original domain. We exploit the symmetries of the domain, and the operator, avoiding necessity of finding good test functions for variational problems. This is especially important for fractional Laplacian, since exact solutions are not even known on intervals, making it hard to find good test functions. To achieve our goals we generalize tight $p$-fusion frames, to extract the best possible geometric results for domains with isometry groups admitting tight $p$-frames. Any such group generates a tight $p$-fusion frame via conjugation of arbitrary projection matrix. We show that generalized tight $p$-frames can also be obtained by conjugation of arbitrary rectangular matrix, with frame constant depending on the singular values of the matrix.


title: Revisiting Algebra and Complexity of Inference in Graphical Models

abstract: This paper studies the form and complexity of inference in graphical models using the abstraction offered by algebraic structures. In particular, we broadly formalize inference problems in graphical models by viewing them as a sequence of operations based on commutative semigroups. We then study the computational complexity of inference by organizing various problems into an "inference hierarchy". When the underlying structure of an inference problem is a commutative semiring -- i.e. a combination of two commutative semigroups with the distributive law -- a message passing procedure called belief propagation can leverage this distributive law to perform polynomial-time inference for certain problems. After establishing the NP-hardness of inference in any commutative semiring, we investigate the relation between algebraic properties in this setting and further show that polynomial-time inference using distributive law does not (trivially) extend to inference problems that are expressed using more than two commutative semigroups. We then extend the algebraic treatment of message passing procedures to survey propagation, providing a novel perspective using a combination of two commutative semirings. This formulation generalizes the application of survey propagation to new settings.


title: A Higher-order Framework for Decision Problems and Games

abstract: We introduce a new unified framework for modelling both decision problems and finite games based on quantifiers and selection functions. We show that the canonical utility maximisation is one special case of a quantifier and that our more abstract framework provides several additional degrees of freedom in modelling. In particular, incomplete preferences, non-maximising heuristics, and context-dependent motives can be taken into account when describing an agent's goal. We introduce a suitable generalisation of Nash equilibrium for games in terms of quantifiers and selection functions. Moreover, we introduce a refinement of Nash that captures context-dependency of goals. Modelling in our framework is compositional as the parts of the game are modular and can be easily exchanged. We provide an extended example where we illustrate concepts and highlight the benefits of our alternative modelling approach.


title: Electron Structure: Shape, Size and GPDs in QED

abstract: The shape of the electron is studied using lowest-order perturbation theory. Quantities used to probe the structure of the proton: form factors, generalized parton distributions, transverse densities, Wigner distributions and the angular momentum content are computed for the electron-photon component of the electron wave function. The influence of longitudinally polarized photons, demanded by the need for infrared regularization via a non-zero photon mass, is included. The appropriate value of the photon mass depends on experimental conditions, and consequently the size of the electron (as defined by the slope of its Dirac form factor) bound in a hydrogen atom is found to be about four times larger than when the electron is free. The shape of the electron, as determined from the transverse density and generalized parton distributions is shown to not be round, and the free electron is shown to be far less round than the bound electron. An electron distribution function (analogous to the quark distribution function) is defined, and that of the bound electron is shown to be suppressed compared to that of the free electron. If the relative transverse momentum of the virtual electron and photon is large compared with the electron mass, the virtual electron and photon each carry nearly the total angular momentum of the physical electron (1/2), with the orbital angular momentum being nearly (-1/2). Including the non-zero photon mass leads to the suppression of end-point contributions to form factors. Implications for proton structure and color transparency are discussed.


title: Staggered $\mathcal{PT}$-symmetric ladders with cubic nonlinearity

abstract: We introduce a ladder-shaped chain with each rung carrying a $\mathcal{PT}$ -symmetric gain-loss dipole. The polarity of the dipoles is staggered along the chain, meaning that a rung bearing gain-loss is followed by one bearing loss-gain. This renders the system $\mathcal{PT}$-symmetric in both horizontal and vertical directions. The system is governed by a pair of linearly coupled discrete nonlinear Schr\"{o}dinger (DNLS) equations with self-focusing or defocusing cubic onsite nonlinearity. Starting from the analytically tractable anti-continuum limit of uncoupled rungs and using the Newton's method for identifying solutions and parametric continuation in the inter-rung coupling for following the associated branches, we construct families of $\mathcal{PT}$-symmetric discrete solitons and identify their stability regions. Waveforms stemming from a single excited rung, as well as ones from multiple rungs are identified. Dynamics of unstable solitons is presented too.


title: Towards Quantum Simulating QCD

abstract: Quantum link models provide an alternative non-perturbative formulation of Abelian and non-Abelian lattice gauge theories. They are ideally suited for quantum simulation, for example, using ultracold atoms in an optical lattice. This holds the promise to address currently unsolvable problems, such as the real-time and high-density dynamics of strongly interacting matter, first in toy-model gauge theories, and ultimately in QCD.


title: CHARIS Science: Performance Simulations for the Subaru Telescope's Third-Generation of Exoplanet Imaging Instrumentation

abstract: We describe the expected scientific capabilities of CHARIS, a high-contrast integral-field spectrograph (IFS) currently under construction for the Subaru telescope. CHARIS is part of a new generation of instruments, enabled by extreme adaptive optics (AO) systems (including SCExAO at Subaru), that promise greatly improved contrasts at small angular separation thanks to their ability to use spectral information to distinguish planets from quasistatic speckles in the stellar point-spread function (PSF). CHARIS is similar in concept to GPI and SPHERE, on Gemini South and the Very Large Telescope, respectively, but will be unique in its ability to simultaneously cover the entire near-infrared $J$, $H$, and $K$ bands with a low-resolution mode. This extraordinarily broad wavelength coverage will enable spectral differential imaging down to angular separations of a few $\lambda/D$, corresponding to $\sim$$0.\!\!''1$. SCExAO will also offer contrast approaching $10^{-5}$ at similar separations, $\sim$$0.\!\!''1$--$0.\!\!''2$. The discovery yield of a CHARIS survey will depend on the exoplanet distribution function at around 10 AU. If the distribution of planets discovered by radial velocity surveys extends unchanged to $\sim$20 AU, observations of $\sim$200 mostly young, nearby stars targeted by existing high-contrast instruments might find $\sim$1--3 planets. Carefully optimizing the target sample could improve this yield by a factor of a few, while an upturn in frequency at a few AU could also increase the number of detections. CHARIS, with a higher spectral resolution mode of $R \sim 75$, will also be among the best instruments to characterize planets and brown dwarfs like HR 8799 cde and $\kappa$ And b.


title: A High-Resolution Study of Near-Infrared Diffuse Interstellar Bands

abstract: We present high-resolution echelle spectroscopic observations of the two near- infrared (NIR) Diffuse Interstellar Bands (DIBs) at 13175 A and 11797.5 A. The DIBs have been observed in a number of diffuse interstellar medium sightlines that exhibit a wide range of visual extinctions. Band profiles are similar to those seen in narrow DIBs, clearly asymmetric and can be closely fitted in most cases using two simple Gaussian components. Gaussian fits were generally found to be more successful than fits based on a multiple-cloud model using a template DIB profile. For a sample of 9 objects in which both bands are observed, the strength of both NIR DIBs generally increases with A(V), and we report a correlation between the two observed bands over a large A(V) range and widely-separated lines of sight. The strength of the two bands is also compared against those of two visual DIBs and the diffuse ISM aliphatic dust absorption feature at 3.4 microns previously detected in the same sightlines. We find that the NIR DIBs do not exhibit notable (anti)correlations with either. Implications of these observations on possible DIB carrier species are discussed.


title: Wavelength Spacing Tunable, Multiwavelength Q-switched Fiber Laser Mode-locked by Graphene Oxide

abstract: We demonstrate a wavelength spacing tunable, multiwavelength Q-switched mode-locked fiber laser (QML) based on a fiber taper deposited with graphene oxide. The operation of the laser can be understood in terms of the formation of bunches of QMLs which possess small temporal intervals, and multiwavelength spectra are generated due to the Fourier transformation. We find that the temporal spacing of the QMLs is highly sensitive to the pump power, and as a result, the wavelength spacing can be easily tuned by varying the pump power. Our experimental laser provides a wavelength spacing tuning range from ~0.001 nm to 0.145 nm with a pump power variation less than 10 mW. The laser could be developed into a low lost wavelength spacing tunable optical source for a wide range of applications, such as spectroscopy, microwave/terahertz signal generation, optical metrology, optical communications and sensing.


title: Modeling Charge-Sign Asymmetric Solvation Free Energies With Nonlinear Boundary Conditions

abstract: We show that charge-sign-dependent asymmetric hydration can be modeled accurately using linear Poisson theory but replacing the standard electric-displacement boundary condition with a simple nonlinear boundary condition. Using a single multiplicative scaling factor to determine atomic radii from molecular dynamics Lennard-Jones parameters, the new model accurately reproduces MD free-energy calculations of hydration asymmetries for (i) monatomic ions, (ii) titratable amino acids in both their protonated and unprotonated states, and (iii) the Mobley "bracelet" and "rod" test problems [J. Phys. Chem. B, v. 112:2408, 2008]. Remarkably, the model also justifies the use of linear response expressions for charging free energies. Our boundary-element method implementation demonstrates the ease with which other continuum-electrostatic solvers can be extended to include asymmetry.


title: Identifying the number of clusters in discrete mixture models

abstract: Research on cluster analysis for categorical data continues to develop, with new clustering algorithms being proposed. However, in this context, the determination of the number of clusters is rarely addressed. In this paper, we propose a new approach in which clustering of categorical data and the estimation of the number of clusters is carried out simultaneously. Assuming that the data originate from a finite mixture of multinomial distributions, we develop a method to select the number of mixture components based on a minimum message length (MML) criterion and implement a new expectation-maximization (EM) algorithm to estimate all the model parameters. The proposed EM-MML approach, rather than selecting one among a set of pre-estimated candidate models (which requires running EM several times), seamlessly integrates estimation and model selection in a single algorithm. The performance of the proposed approach is compared with other well-known criteria (such as the Bayesian information criterion-BIC), resorting to synthetic data and to two real applications from the European Social Survey. The EM-MML computation time is a clear advantage of the proposed method. Also, the real data solutions are much more parsimonious than the solutions provided by competing methods, which reduces the risk of model order overestimation and increases interpretability.


title: Hamiltonian commutators with large Hofer norm

abstract: We show that commutators of Hamiltonian diffeomorphisms may have arbitrarily large Hofer norm. The proposed technique is applicable to positive genus surfaces and their products. This gives partial answer to a question by McDuff and Polterovich.


title: Symmetry breaking for an elliptic equation involving the Fractional Laplacian

abstract: We study the symmetry breaking phenomenon for an elliptic equation involving the fractional Laplacian in a large ball. Our main tool is an extension of the Strauss radial lemma involving the fractional Laplacian, which might be of independent interest; and from which we derive compact embedding theorems for a Sobolev-type space of radial functions with power weights.


title: Homology of the braid group with coefficients in the reduced Burau representation

abstract: The reduced Burau representation $V_n$ of the braid group $B_n$ is obtained from the action of $B_n$ on the homology of an infinite cyclic cover of the $n$-punctured disc. In this note, we calculate $H_*(B_n;V_n)$ as a module over the Laurent polynomial ring $\mathbb{Q}[t,t^{-1}]$.


title: High-order boundary integral equation solution of high frequency wave scattering from obstacles in an unbounded linearly stratified medium

abstract: We apply boundary integral equations for the first time to the two-dimensional scattering of time-harmonic waves from a smooth obstacle embedded in a continuously-graded unbounded medium. In the case we solve the square of the wavenumber (refractive index) varies linearly in one coordinate, i.e. $(\Delta + E + x_2)u(x_1,x_2) = 0$ where $E$ is a constant; this models quantum particles of fixed energy in a uniform gravitational field, and has broader applications to stratified media in acoustics, optics and seismology. We evaluate the fundamental solution efficiently with exponential accuracy via numerical saddle-point integration, using the truncated trapezoid rule with typically 100 nodes, with an effort that is independent of the frequency parameter $E$. By combining with high-order Nystrom quadrature, we are able to solve the scattering from obstacles 50 wavelengths across to 11 digits of accuracy in under a minute on a desktop or laptop.


title: Eigenfunction Statistics for Anderson Model with H\"{o}lder continuous single site Potential

abstract: We consider random Schr\"{o}dinger operators on $\ell^2(\mathbb{Z}^d)$ when the distribution of single site potentials is $\alpha$-H\"{o}lder continuous ($0<\alpha\leq 1$). In localized regime we study the distribution of eigenfunctions simultaneously in space and energy. In a certain scaling limit we prove limits point are Poisson.


title: Net and Prune: A Linear Time Algorithm for Euclidean Distance Problems

abstract: We provide a general framework for getting expected linear time constant factor approximations (and in many cases FPTAS's) to several well known problems in Computational Geometry, such as $k$-center clustering and farthest nearest neighbor. The new approach is robust to variations in the input problem, and yet it is simple, elegant and practical. In particular, many of these well studied problems which fit easily into our framework, either previously had no linear time approximation algorithm, or required rather involved algorithms and analysis. A short list of the problems we consider include farthest nearest neighbor, $k$-center clustering, smallest disk enclosing $k$ points, $k$th largest distance, $k$th smallest $m$-nearest neighbor distance, $k$th heaviest edge in the MST and other spanning forest type problems, problems involving upward closed set systems, and more. Finally, we show how to extend our framework such that the linear running time bound holds with high probability.


title: Multi-Boson Correlation Sampling with Multi-mode Thermal Sources

abstract: The determination of the computational complexity of the boson sampling problem with single boson sources has opened a novel research direction in the quantum computing field. Some research effort has also been devoted towards the use of different input sources where the sampling from the output distribution of a linear interferometer is still a complex task not achievable by a classical device. In this letter, we discuss the use of multi-mode thermal sources, based on a generalization of the boson sampling problem which relies on time-correlation measurements. We derive, for any possible sample, two equivalent formulations for the respective multi-boson correlation functions in terms of matrix permanents. The properties of these permanents emerging from the physics of multi-boson interference raise essential questions in the field of complexity theory. Our results are a manifestation of the physics behind correlated measurements in either photonic or atomic interferometers with bosonic thermal sources. This fascinating physics paves the way towards important applications not only in quantum computing but also in quantum metrology and quantum communication, by taking advantage of different kinds of bosonic sources.


title: Anisotropy in Fracking: A Percolation Model for Observed Microseismicity

abstract: Hydraulic fracturing (fracking) using high pressures and a low viscosity fluid allow the extraction of large quantiles of oil and gas from very low permeability shale formations. The initial production of oil and gas at depth leads to high pressures and an extensive distribution of natural fractures which reduce the pressures. With time these fractures heal, sealing the remaining oil and gas in place. High volume fracking opens the healed fractures allowing the oil and gas to flow the horizontal productions wells. We model the injection process using invasion percolation. We utilize a 2D square lattice of bonds to model the sealed natural fractures. The bonds are assigned random strengths and the fluid, injected at a point, opens the weakest bond adjacent to the growing cluster of opened bonds. Our model exhibits burst dynamics in which the clusters extends rapidly into regions with weak bonds. We associate these bursts with the microseismic activity generated by fracking injections. A principal object of this paper is to study the role of anisotropic stress distributions. Bonds in the $y$-direction are assigned higher random strengths than bonds in the $x$-direction. We illustrate the spatial distribution of clusters and the spatial distribution of bursts (small earthquakes) for several degrees of anisotropy. The results are compared with observed distributions of microseismicity in a fracking injection. Both our bursts and the observed microseismicity satisfy Gutenberg-Richter frequency-size statistics.


title: Incidence of $q$-statistics in rank distributions

abstract: We show that size-rank distributions with power-law decay (often only over a limited extent) observed in a vast number of instances in a widespread family of systems obey Tsallis statistics. The theoretical framework for these distributions is analogous to that of a nonlinear iterated map near a tangent bifurcation for which the Lyapunov exponent is negligible or vanishes. The relevant statistical-mechanical expressions associated with these distributions are derived from a maximum entropy principle with the use of two different constraints, and the resulting duality of entropy indexes is seen to portray physically relevant information. While the value of the index $\alpha $ fixes the distribution's power-law exponent, that for the dual index $2-\alpha $ ensures the extensivity of the deformed entropy.


title: A planetary system and a highly eccentric brown dwarf around the giant stars HIP67851 and HIP97233

abstract: Context: So far more than 60 substellar companions have been discovered around giant stars. These systems present physical and orbital properties that contrast to those detected orbiting less evolved stars. Aims: We are conducting a radial velocity survey of 166 bright giant stars in the southern hemisphere. The main goals of our project are to detect and characterize planets in close-in orbits around giant stars in order to study the effects of the host star evolution on their orbital and physical properties. Methods: We have obtained precision radial velocities for the giant stars HIP67851 and HIP97233 that have revealed periodic signals, which are most likely induced by the presence of substellar companions. Results: We present the discovery of a planetary system and an eccentric brown dwarf orbiting the giant stars HIP67851 and HIP97233, respectively. The inner planet around HIP67851 has a period of 88.8 days, a projected mass of 1.4 Mjup and an eccentricity of 0.09. After Kepler 91b, HIP67851b is the closest-in known planet orbiting a giant star. Although the orbit of the outer object is not fully constrained, it is likely a super-Jupiter. The brown dwarf around HIP97233 has an orbital period of 1058.8 days, a minimum mass of 20.0 Mjup and an eccentricity of 0.61. This is the most eccentric known brown dwarf around a giant star.


title: How to repair tropicalizations of plane curves using modifications

abstract: Tropical geometry is sensitive to embeddings of algebraic varieties inside toric varieties. The purpose of this paper is to advertise tropical modifications as a tool to locally repair bad embeddings of plane curves, allowing the re-embedded tropical curve to better reflect the geometry of the input curve. Our approach is based on the close connection between analytic curves (in the sense of Berkovich) and tropical curves. We investigate the effect of these tropical modifications on the tropicalization map defined on the analytification of the given curve. Our study is motivated by the case of plane elliptic cubics, where good embeddings are characterized in terms of the j-invariant. Given a plane elliptic cubic whose tropicalization contains a cycle, we present an effective algorithm, based on non-Archimedean methods, to linearly re-embed the curve in dimension 4 so that its tropicalization reflects the j-invariant. We give an alternative elementary proof of this result by interpreting the initial terms of the A-discriminant of the defining equation as a local discriminant in the Newton subdivision.


title: Range-separated meta-GGA functional designed for noncovalent interactions

abstract: The accuracy of applying density functional theory to noncovalent interactions is hindered by errors arising from low-density regions of interaction-induced change in the density gradient, error compensation between correlation and exchange functionals, and dispersion double counting. A new exchange-correlation functional designed for noncovalent interactions is proposed to address these problems. The functional consists of the range-separated PBEsol exchange considered in two variants, pure and hybrid, and the semilocal correlation functional of Modrzejewski et al. [J. Chem. Phys. 137, 204121 (2012)] designed with the constraint satisfaction technique to smoothly connect with a dispersion term. Two variants of dispersion correction are appended to the correlation functional: the atom-atom pairwise additive DFT-D3 model and the density-dependent many-body dispersion with self-consistent screening (MBD-rsSCS). From these building blocks a set of four functionals is created to systematically examine the role of pure versus hybrid exchange and the underlying models for dispersion. The new functional is extensively tested on benchmark sets with diverse nature and size. The truly outstanding performance is demonstrated for water clusters of varying size, ionic hydrogen bonds, and thermochemistry of isodesmic n-alkane fragmentation reactions. The merits of each component of the new functional are discussed.


title: Approaching the Asymptotic Regime of Rapidly Rotating Convection: Boundary Layers vs Interior Dynamics

abstract: Rapidly rotating Rayleigh-B\'enard convection is studied by combining results from direct numerical simulations (DNS), laboratory experiments and asymptotic modeling. The asymptotic theory is shown to provide a good description of the bulk dynamics at low, but finite Rossby number. However, large deviations from the asymptotically predicted heat transfer scaling are found, with laboratory experiments and DNS consistently yielding much larger Nusselt numbers than expected. These deviations are traced down to dynamically active Ekman boundary layers, which are shown to play an integral part in controlling heat transfer even for Ekman numbers as small as $10^{-7}$. By adding an analytical parameterization of the Ekman transport to simulations using stress-free boundary conditions, we demonstrate that the heat transfer jumps from values broadly compatible with the asymptotic theory to states of strongly increased heat transfer, in good quantitative agreement with no-slip DNS and compatible with the experimental data. Finally, similarly to non-rotating convection, we find no single scaling behavior, but instead that multiple well-defined dynamical regimes exist in rapidly-rotating convection systems.


title: Throughput Analysis for Wireless Networks with Full-Duplex Radios

abstract: This paper investigates the throughput for wireless network with full-duplex radios using stochastic geometry. Full-duplex (FD) radios can exchange data simultaneously with each other. On the other hand, the downside of FD transmission is that it will inevitably cause extra interference to the network compared to half-duplex (HD) transmission. In this paper, we focus on a wireless network of nodes with both HD and FD capabilities and derive and optimize the throughput in such a network. Our analytical result shows that if the network is adapting an ALOHA protocol, the maximal throughput is always achieved by scheduling all concurrently transmitting nodes to work in FD mode instead of a mixed FD/HD mode or HD mode regardless of the network configurations. Moreover, the throughput gain of using FD transmission over HD transmission is analytically lower and upper bounded.


title: Wide Quantum Critical Region of Valence Fluctuations: Origin of Robust Quantum Criticality in Quasicrystal Yb15Al34Au51 under Pressure

abstract: The mechanism of the emergence of robust quantum criticality in the heavy-electron quasicrystal Yb15Al34Au51 is analyzed theoretically. By constructing a minimal model for the quasicrystal and its crystalline approximant, which contain concentric shell structures with Yb and Al-Au clusters, we show that a set of quantum critical points of the first-order valence transition of Yb appears as spots in the ground-state phase diagram. Their critical regions overlap each other, giving rise to a wide quantum critical region. This well explains the robust criticality observed in Yb15Al34Au51 under pressure, and predicts the emergence of the common criticality in the crystalline approximant under pressure. The wider critical region in the quasicrystal than that in the crystalline approximant in the T-P phase diagram and the field-induced valence-crossover "region" in the T-H phase diagram are predicted to appear.


title: D=4 supersymmetric Yang-Mills quantum mechanics with three colors

abstract: The spectrum of the $D=4$ supersymmetric Yang-Mills quantum mechanics with $SU(3)$ gauge group symmetry is computed in different channels with definite total angular momentum and the total number of fermions. In sectors with small number of fermions the spectrum is shown to be purely discrete. Continuous spectrum appears in channels with more fermions. Supersymmetry generators are used to identify supermultiplets and determine the level of restoration of supersymmetry for a finite cutoff.


title: Particle simulation of lower hybrid waves in tokamak plasmas

abstract: Global particle simulations of the lower hybrid waves have been carried out using fully kinetic ions and drift kinetic electrons with a realistic electron-to-ion mass ratio. The lower hybrid wave frequency, mode structure, and electron Landau damping from the electrostatic simulations agree very well with the analytic theory. Linear simulation of the propagation of a lower hybrid wave-packet in the toroidal geometry shows that the wave propagates faster in the high field side than the low field side, in agreement with a ray tracing calculation. Electromagnetic benchmarks of lower hybrid wave dispersion relation are also carried out. Electromagnetic mode conversion are observed in toroidal geometry, slow waves are launched at the plasma boundary and converts to fast waves at the mode conversion layer, which is consistent with linear theory.


title: Direct images of relative pluricanonical bundles

abstract: We discuss the local freeness and the numerical semipositivity of direct images of relative pluricanonical bundles for surjective morphisms between smooth projective varieties with connected fibers. We give a desirable semipositivity theorem under the assumption that the geometric generic fiber has a good minimal model.


title: QCD chiral symmetry restoration with a large number of quarks in a model with a confining propagator and dynamically massive gluons

abstract: Considering a QCD chiral symmetry breaking model where the gap equation contains an effective confining propagator and a dressed gluon propagator with a dynamically generated mass, we verify that the chiral symmetry is restored for a large number of quarks $n_{f}\approx 7-13$. We discuss the uncertainty in the results, that is related to the determination of the string tension ($K_{F}$), appearing in the confining propagator, and the effective gluon mass ($m_{g}$) at large $n_{f}$.


title: Quasi-exact-solvability of the $A_{2}/G_2$ Elliptic model: algebraic forms, $sl(3)/g^{(2)}$ hidden algebra, polynomial eigenfunctions

abstract: The potential of the $A_2$ quantum elliptic model (3-body Calogero-Moser elliptic model) is defined by the pairwise three-body interaction through Weierstrass $\wp$-function and has a single coupling constant. A change of variables has been found, which are $A_2$ elliptic invariants, such that the potential becomes a rational function, while the flat space metric as well as its associated vector are polynomials in two variables. It is shown that the model possesses the hidden $sl(3)$ algebra - the Hamiltonian is an element of the universal enveloping algebra $U_{sl(3)}$ for arbitrary coupling constant - thus, it is equivalent to $sl(3)$-quantum Euler-Arnold top. The integral, in a form of the third order differential operator with polynomial, is constructed explicitly, being also an element of $U_{sl(3)}$. It is shown that there exists a discrete sequence of the coupling constants for which a finite number of polynomial eigenfunctions, up to a (non-singular) gauge factor occur. The potential of the $G_2$ quantum elliptic model (3-body Wolfes elliptic model) is defined by the pairwise and three-body interactions through Weierstrass $\wp$-function and has two coupling constants. A change of variables has been found, which are $G_2$ elliptic invariants, such that the potential becomes a rational function, while the flat space metric as well as its associated vector are polynomials in two variables. It is shown the model possesses the hidden $g^{(2)}$ algebra. It is shown that there exists a discrete family of the coupling constants for which a finite number of polynomial eigenfunctions up to a (non-singular) gauge factor occur.


title: Quasinormal modes of charged fields around a Reissner-Nordstr\"om black hole

abstract: The quasinormal spectrum of a charged scalar field around a non-extremal Reissner-Nordstr\"om black hole, specially in the limit of large electromagnetic interaction, has been comprehensively studied only very recently. In this work, we extend the analysis to Dirac fields using the continued fraction method and compare the results with the scalar case. In particular, we study the behaviour of the fundamental quasinormal mode as a function of the black hole's charge and of the electromagnetic interaction parameter. We derive an analytical formula for the quasinormal frequencies in the limit of large electromagnetic interaction. As the extremal limit of black hole charge is approached, we show that, unlike the case of neutral fields, the imaginary part of the quasinormal frequencies approach zero for charged fields.


title: Model identification using the Efficient Determination Criterion

abstract: In the realm of the model selection context, Akaike's and Schwarz's information criteria, AIC and BIC, have been applied successfully for decades for model order identification. The Efficient Determination Criterion (EDC) is a generalization of these criteria, proposed originally to define a strongly consistent class of estimators for the dependency order of a multiple Markov chain. In this work, the EDC is generalized to partially nested models, which encompass many other order identification problems. Based on some assumptions, a class of strongly consistent estimators is established in this general environment. This framework is applied to BEKK multivariate GARCH models and, in particular, the strong consistency of the order estimator based on BIC is established for these models.


title: Filtering from Observations on Stiefel Manifolds

abstract: This paper considers the problem of optimal filtering for partially observed signals taking values on the rotation group. More precisely, one or more components are considered not to be available in the measurement of the attitude of a 3D rigid body. In such cases, the observed signal takes its values on a Stiefel manifold. It is demonstrated how to filter the observed signal through the anti-development built from observations. A particle filter implementation is proposed to perform the estimation of the signal partially observed and corrupted by noise. The sampling issue is also addressed and interpolation methods are introduced. Illustration of the proposed technique on synthetic data demonstrates the ability of the approach to estimate the angular velocity of a partially observed 3D system partially observed.


title: Ranking algorithms on directed configuration networks

abstract: This paper studies the distribution of a family of rankings, which includes Google's PageRank, on a directed configuration model. In particular, it is shown that the distribution of the rank of a randomly chosen node in the graph converges in distribution to a finite random variable $\mathcal{R}^*$ that can be written as a linear combination of i.i.d. copies of the endogenous solution to a stochastic fixed point equation of the form $$\mathcal{R} \stackrel{\mathcal{D}}{=} \sum_{i=1}^{\mathcal{N}} \mathcal{C}_i \mathcal{R}_i + \mathcal{Q},$$ where $(\mathcal{Q}, \mathcal{N}, \{ \mathcal{C}_i\})$ is a real-valued vector with $\mathcal{N} \in \{0,1,2,\dots\}$, $P(|\mathcal{Q}| > 0) > 0$, and the $\{\mathcal{R}_i\}$ are i.i.d. copies of $\mathcal{R}$, independent of $(\mathcal{Q}, \mathcal{N}, \{ \mathcal{C}_i\})$. Moreover, we provide precise asymptotics for the limit $\mathcal{R}^*$, which when the in-degree distribution in the directed configuration model has a power law imply a power law distribution for $\mathcal{R}^*$ with the same exponent.


title: Fifth Memoir or Letter on the Resistance of Air to the Movement of Pendulums

abstract: Translation of Denis Diderot's `Cinqui\`eme M\'emoire ou Lettre sur la R\'esitance de L'Air ou Mouvement des Pendules.'


title: The Theory of Witt Vectors

abstract: This is an introduction to the theory of Witt vectors. It includes a construction of the Witt rings, the Frobenius and Verschiebung endomorphisms, the canonical map from W to W^2 (its lambda-algebra structure), the relation to strict p-rings, and an account of the Artin-Hasse exponential.


title: The Mathematics and Physics of Diderot. I. On Pendulums and Air Resistance

abstract: In this article Denis Diderot's Fifth Memoir of 1748 on the problem of a pendulum damped by air resistance is discussed. Diderot wrote the Memoir in order to clarify an assumption Newton made without further justification in the first pages of the Principia in connection with an experiment to verify the Third Law of Motion using colliding pendulums. To explain the differences between experimental and theoretical values of momentum in the collision experiments he conducted Newton assumed that the bob was retarded by an air resistance $F_R$ proportional to the velocity $v$. By giving Newton's arguments a mathematical scaffolding and recasting his geometrical reasoning in the language of differential calculus, Diderot provides a step-by-step solution guide to the problem and proposes experiments to settle the question about the appropriate form of $F_R$, which for Diderot quadratic in $v$, that is $F_R \sim v^2$. The solution of Diderot is presented in full detail and his results are compared to those obtained from a Lindstedt-Poincare approximation for an oscillator with quadratic damping. It is shown that, up to a prefactor, both coincide. Some results that one can derive from his approach are presented and discussed for the first time. Experimental evidence to support Diderot's or Newton's claims is discussed together with the limitations of their solutions. Some misprints in the original memoir are pointed out.


title: Near invariance of the hypercube

abstract: We give an almost-complete description of orthogonal matrices $M$ of order $n$ that "rotate a non-negligible fraction of the Boolean hypercube $C_n=\{-1,1\}^n$ onto itself," in the sense that $$P_{x\in C_n}(Mx\in C_n) \ge n^{-C},\mbox{ for some positive constant } C,$$ where $x$ is sampled uniformly over $C_n$. In particular, we show that such matrices $M$ must be very close to products of permutation and reflection matrices. This result is a step toward characterizing those orthogonal and unitary matrices with large permanents, a question with applications to linear-optical quantum computing.


title: GRACES: Gemini remote access to CFHT ESPaDOnS Spectrograph through the longest astronomical fiber ever made (Experimental phase completed.)

abstract: The Gemini Remote Access to CFHT ESPaDONS Spectrograph has achieved first light of its experimental phase in May 2014. It successfully collected light from the Gemini North telescope and sent it through two 270 m optical fibers to the the ESPaDOnS spectrograph at CFHT to deliver high-resolution spectroscopy across the optical region. The fibers gave an average focal ratio degradation of 14% on sky, and a maximum transmittance of 85% at 800nm. GRACES achieved delivering spectra with a resolution power of R = 40,000 and R = 66,000 between 400 and 1,000 nm. It has a ~8% throughput and is sensitive to target fainter than 21st mag in 1 hour. The average acquisition time of a target is around 10 min. This project is a great example of a productive collaboration between two observatories on Maunakea that was successful due to the reciprocal involvement of the Gemini, CFHT, and NRC Herzberg teams, and all the staff involved closely or indirectly.


title: Scaling of distributions of sums of positions for chaotic dynamics at band-splitting points

abstract: The stationary distributions of sums of positions of trajectories generated by the logistic map have been found to follow a basic renormalization group (RG) structure: a nontrivial fixed-point multi-scale distribution at the period-doubling onset of chaos and a Gaussian trivial fixed-point distribution for all chaotic attractors. Here we describe in detail the crossover distributions that can be generated at chaotic band-splitting points that mediate between the aforementioned fixed-point distributions. Self affinity in the chaotic region imprints scaling features to the crossover distributions along the sequence of band splitting points. The trajectories that give rise to these distributions are governed first by the sequential formation of phase-space gaps when, initially uniformly-distributed, sets of trajectories evolve towards the chaotic band attractors. Subsequently, the summation of positions of trajectories already within the chaotic bands closes those gaps. The possible shapes of the resultant distributions depend crucially on the disposal of sets of early positions in the sums and the stoppage of the number of terms retained in them.


title: Two-stage Geometric Information Guided Image Reconstruction

abstract: In compressive sensing, it is challenging to reconstruct image of high quality from very few noisy linear projections. Existing methods mostly work well on piecewise constant images but not so well on piecewise smooth images such as natural images, medical images that contain a lot of details. We propose a two-stage method called GeoCS to recover images with rich geometric information from very limited amount of noisy measurements. The method adopts the shearlet transform that is mathematically proven to be optimal in sparsely representing images containing anisotropic features such as edges, corners, spikes etc. It also uses the weighted total variation (TV) sparsity with spatially variant weights to preserve sharp edges but to reduce the staircase effects of TV. Geometric information extracted from the results of stage I serves as an initial prior for stage II which alternates image reconstruction and geometric information update in a mutually beneficial way. GeoCS has been tested on incomplete spectral Fourier samples. It is applicable to other types of measurements as well. Experimental results on various complicated images show that GeoCS is efficient and generates high-quality images.


title: Anisotropic transverse magnetoresistivity in alpha-YbAlB4

abstract: We measured the transverse magnetoresistivity of the mixed valence compound $\alpha$-YbAlB$_4$. Two configurations were used where current was applied along [110] direction for both and magnetic field was applied along [-110] and $c$-axis. We found the transverse magnetoresistivity is highly anisotropic. In the weak field below 1 T, it is consistent with stronger $c$-$f$ hybridization in the $ab$ plane which was suggested from the previous zero field resistivity measurements. At the higher field above 3 T, we observed a negative transverse magnetoresistivity for the field applied along the $c$-axis. The temperature dependences of the resistivity measured at several different fields suggest the suppression of the heavy fermion behavior at the characteristic field of $\sim 5.5$ T.


title: An ultimate storage ring lattice with vertical emittance generated by damping wigglers

abstract: We discuss the approach of generating round beams for ultimate storage rings using vertical damping wigglers (with horizontal magnetic field). The vertical damping wigglers provide damping and excite vertical emittance. This eliminates the need to generate large linear coupling that is impractical with traditional off-axis injection. We use a PEP-X compatible lattice to demonstrate the approach. This lattice uses separate quadrupole and sextupole magnets with realistic gradient strengths. Intrabeam scattering effects are calculated. The horizontal and vertical emittances are 22.3 pm and 10.3 pm, respectively, for a 200 mA, 4.5 GeV beam, with a vertical damping wiggler of a total length of 90 meters, peak field of 1.5 T and wiggler period of 100 mm.


title: Non-divergent representation of non-Hermitian operator near the exceptional point with application to a quantum Lorentz gas

abstract: We propose a non-singular representation for a non-Hermitian operator even if the parameter space contains exceptional points (EPs), at which the operator cannot be diagonalized and the usual spectral representation ceases to exist. Our representation has a generalized Jordan block form and is written in terms of extended pseudo-eigenstates. Our method is free from a divergence in the spectral representation at EPs, at which multiple eigenvalues and eigenvectors coalesce and the eigenvectors cannot be normalized. Our representation improves the accuracy of numerical calculations of physical quantities near EPs. We also find that our method is applicable to various problems related to EPs in the parameter space of non-Hermitian operators. We demonstrate the usefulness of our representation by investigating Boltzmann's collision operator in a one-dimensional quantum Lorentz gas in the weak coupling approximation.


title: A Bayesian spatial temporal mixtures approach to kinetic parametric images in dynamic Positron Emission Tomography

abstract: We present a fully Bayesian statistical approach to the problem of compartmental modelling in the context of Positron Emission Tomography. We cluster homogeneous region of interest and perform kinetic parameter estimation simultaneously. A mixture modelling approach is adopted, incorporating both spatial and temporal information based on reconstructed dynamic PET image. Our modelling approach is flexible, and provides uncertainty estimates for the estimated kinetic parameters. Crucially, the proposed method allows us to determine the unknown number of clusters, which has a great impact on resulting estimated kinetic parameters. We demonstrate our method on simulated dynamic Myocardial PET data, and show that our method is superior to standard curve-fitting approach.


title: Point particles in 2+1 dimensions: general relativity and loop gravity descriptions

abstract: We develop a Hamiltonian description of point particles in (2+1)-dimensions using connection and frame-field variables for general relativity. The topology of each spatial hypersurface is that of a punctured two-sphere with particles residing at the punctures. We describe this topology with a CW complex (a collection of two-cells glued together along edges), and use this to fix a gauge and reduce the Hamiltonian. The equations of motion for the fields describe a dynamical triangulation where each vertex moves according to the equation of motion for a free relativistic particle. The evolution is continuous except for when triangles collapse (i.e. the edges become parallel) causing discrete, topological changes in the underlying CW complex. We then introduce the loop gravity phase space parameterized by holonomy-flux variables on a graph (a network of one-dimensional links). By embedding a graph within the CW complex, we find a description of this system in terms of loop variables. The resulting equations of motion describe the same dynamical triangulation as the connection and frame-field variables. In this framework, the collapse of a triangle causes a discrete change in the underlying graph, giving a concrete realization of the graph-changing moves that many expect to feature in full loop quantum gravity. The main result is a dynamical model of loop gravity which agrees with general relativity and is well-suited for quantization using existing methods.


title: Fiber-Optic Interferometry Using Narrowband Light Source and Electrical Spectrum Analyzer: Influence on Brillouin Measurement

abstract: We observe an interference pattern using a simple fiber-optic interferometer consisting of an electrical spectrum analyzer and a narrowband light source, which is commonly employed for observing the Brillouin gain spectrum. This interference pattern expands well beyond the frequency range corresponding to the Brillouin frequency shift in silica fibers (approximately 11 GHz at 1550 nm). Using both silica single-mode and polymer optical sensing fibers, we then experimentally prove that the distinctive noise in a self-heterodyne-based Brillouin measurement with an unoptimized polarization state originates from the interference between the reference light and the Fresnel-reflected light. This noise can be almost completely suppressed by employing a delay line that is longer than the coherence length of the light source and by artificially applying a high loss near the open end of the sensing fiber.


title: Icosahedra boron chain and sheets: new boron allotropic structures

abstract: The icosahedra boron chain and three icosahedra sheets (with {\alpha}, {\delta}4, and {\delta}6 symmetries), constructed by the icosahedra B12, have been obtained as new members of boron family using a highly efficient molecular dynamics scheme based on a transferable and reliable semi-empirical Hamiltonian. The icosahedral B12 in the icosahedra chain is slightly elongated along the china direction and directly bonded each other with the two-center covalent bonds. A deformation of the icosahedra B12 was also found in the two-dimensional icosahedra sheets. In addition to the three-center bonding nature inside the icosahedra B12, there are two types of directional inter-icosahedra bonds in the icosahedra sheet structures, one is the single strong covalent bond, and the other is a pair of the weak covalent {\delta} bonds. In contrast to the boron monolayer, there is no buckling found in these icosahedra sheets. The deformation of the icosahedra B12 and the special bonding nature in these new icosahedra structures induce the energy band gap of 0.74 eV in the icosahedra chain, 0.52 eV in the icosahedra {\delta}6 sheet, 0.39 eV in the icosahedra {\delta}4 sheet, and the gapless in the icosahedra {\alpha} sheet, respectively. The energy barrier per atom from the icosahedra {\delta}6 sheet to the icosahedra {\alpha} sheet is estimated to be 0.17 eV while it is estimated as 0.38 eV from the icosahedra {\delta}6 sheet to the icosahedra {\delta}4 sheet and 0.27 eV from the icosahedra {\alpha} sheet to the icosahedra {\delta}4 sheet, respectively. Such high energy barriers indicate that these icosahedra sheets are relatively stable.


title: Beyond Maximum Likelihood: from Theory to Practice

abstract: Maximum likelihood is the most widely used statistical estimation technique. Recent work by the authors introduced a general methodology for the construction of estimators for functionals in parametric models, and demonstrated improvements - both in theory and in practice - over the maximum likelihood estimator (MLE), particularly in high dimensional scenarios involving parameter dimension comparable to or larger than the number of samples. This approach to estimation, building on results from approximation theory, is shown to yield minimax rate-optimal estimators for a wide class of functionals, implementable with modest computational requirements. In a nutshell, a message of this recent work is that, for a wide class of functionals, the performance of these essentially optimal estimators with $n$ samples is comparable to that of the MLE with $n \ln n$ samples. In the present paper, we highlight the applicability of the aforementioned methodology to statistical problems beyond functional estimation, and show that it can yield substantial gains. For example, we demonstrate that for learning tree-structured graphical models, our approach achieves a significant reduction of the required data size compared with the classical Chow--Liu algorithm, which is an implementation of the MLE, to achieve the same accuracy. The key step in improving the Chow--Liu algorithm is to replace the empirical mutual information with the estimator for mutual information proposed by the authors. Further, applying the same replacement approach to classical Bayesian network classification, the resulting classifiers uniformly outperform the previous classifiers on 26 widely used datasets.


title: Orbital Angular Momentum and Spectral Flow in Two Dimensional Chiral Superfluids

abstract: We study the orbital angular momentum (OAM) $L_z$ in two dimensional chiral $(p_x+ip_y)^{\nu}$-wave superfluids (SF) of $N$ fermions on a disc at zero temperature, in terms of spectral asymmetry and spectral flow. It is shown that $L_z=\nu N/2$ for any integer $\nu$, in the BEC regime. In contrast, in the BCS limit, while the OAM is $L_z=N/2$ for the $p+ip$-wave SF, for chiral SF with $\nu\geq2$, the OAM is remarkably suppressed as $L_z=N\times O(\Delta_0/\varepsilon_F)\ll N$, where $\Delta_0$ is the gap amplitude and $\varepsilon_F$ is the Fermi energy. We demonstrate that the difference between the $p+ip$-wave SF and the other chiral SFs in the BCS regimes originates from the nature of edge modes and related depairing effects.


title: Projective symmetry of partons in Kitaev's honeycomb model

abstract: Low-energy states of quantum spin liquids are thought to involve partons living in a gauge-field background. We study the spectrum of Majorana fermions of Kitaev's honeycomb model on spherical clusters. The gauge field endows the partons with half-integer orbital angular momenta. As a consequence, the multiplicities reflect not the point-group symmetries of the cluster, but rather its projective symmetries, operations combining physical and gauge transformations. The projective symmetry group of the ground state is the double cover of the point group.


title: Autoencoder Trees

abstract: We discuss an autoencoder model in which the encoding and decoding functions are implemented by decision trees. We use the soft decision tree where internal nodes realize soft multivariate splits given by a gating function and the overall output is the average of all leaves weighted by the gating values on their path. The encoder tree takes the input and generates a lower dimensional representation in the leaves and the decoder tree takes this and reconstructs the original input. Exploiting the continuity of the trees, autoencoder trees are trained with stochastic gradient descent. On handwritten digit and news data, we see that the autoencoder trees yield good reconstruction error compared to traditional autoencoder perceptrons. We also see that the autoencoder tree captures hierarchical representations at different granularities of the data on its different levels and the leaves capture the localities in the input space.


title: The Early Universe with High-Scale Supersymmetry

abstract: A small tensor-to-scalar ratio $r$ may lead to distinctive phenomenology of high-scale supersymmetry. Assuming the same origin of SUSY breaking between the inflation and visible sector, we show model independent features. The simplest hybrid inflation, together with a new linear term for the inflaton field which is induced by large gravitino mass, is shown to be consistent with all experimental data for $r$ of order $10^{-5}$. For superpartner masses far above the weak scale we find that the reheating temperature after inflation might be beneath the value required by thermal leptogenesis if the inflaton decays to its products perturbatively, but above it if non-perturbatively instead. Remarkably, the gravitino overproduction can be evaded in such high-scale supersymmetry because of the kinematically blocking effect.


title: Revisiting random deposition with surface relaxation: approaches from growth rules to Edwards-Wilkinson equation

abstract: We present several approaches for deriving the coarse-grained continuous Langevin equation (or Edwards-Wilkinson equation) from a random deposition with surface relaxation (RDSR) model. First we introduce a novel procedure to divide the first transition moment into the three fundamental processes involved: deposition, diffusion and volume conservation. We show how the diffusion process is related to antisymmetric contribution and the volume conservation process is related to symmetric contribution, which renormalizes to zero in the coarse-grained limit. In another approach, we find the coefficients of the continuous Langevin equation, by regularizing the discrete Langevin equation. Finally, in a third approach, we derive these coefficients from the set of test functions supported by the stationary probability density function (SPDF) of the discrete model. The applicability of the used approaches to other discrete random deposition models with instantaneous relaxation to a neighboring site is discussed.


title: High-order Numerical Methods for Riesz Space Fractional Turbulent Diffusion Equation

abstract: Numerical methods for fractional calculus attract increasing interests due to its wide applications in various fields such as physics, mechanics, etc. In this paper, we focus on constructing high-order algorithms for Riesz derivatives, where the convergence orders cover from the second order to the sixth order. Then we apply the established schemes to the Riesz space fractional turbulent diffusion equation. Numerical experiments are displayed which support the theoretical analysis.


title: Error correcting codes and spatial coupling

abstract: These are notes from the lecture of R\"udiger Urbanke given at the autumn school "Statistical Physics, Optimization, Inference, and Message-Passing Algorithms", that took place in Les Houches, France from Monday September 30th, 2013, till Friday October 11th, 2013. The school was organized by Florent Krzakala from UPMC and ENS Paris, Federico Ricci-Tersenghi from La Sapienza Roma, Lenka Zdeborov\`a from CEA Saclay and CNRS, and Riccardo Zecchina from Politecnico Torino. The first three sections cover the basics of polar codes and low density parity check codes. In the last three sections, we see how the spatial coupling helps belief propagation decoding.

categories: cs.IT cond-mat.dis-nn cond-mat.stat-mech math.IT


title: Weierstrass points on the Drinfeld modular curve $X_0(\mathfrak{p})$

abstract: Consider the Drinfeld modular curve $X_0(\mathfrak{p})$ for $\mathfrak{p}$ a prime ideal of $\mathbb{F}_q[T]$. It was previously known that if $j$ is the $j$-invariant of a Weierstrass point of $X_0(\mathfrak{p})$, then the reduction of $j$ modulo $\mathfrak{p}$ is a supersingular $j$-invariant. In this paper we show the converse: Every supersingular $j$-invariant is the reduction modulo $\mathfrak{p}$ of the $j$-invariant of a Weierstrass point of $X_0(\mathfrak{p})$.


title: Trans-Planckian fluctuations and the stability of quantum mechanics

abstract: We present arguments suggesting that deviations from the Born probability rule could be generated for trans-Planckian field modes during inflation. Such deviations are theoretically possible in the de Broglie-Bohm pilot-wave formulation of quantum mechanics, according to which the Born rule describes a state of statistical equilibrium. We suggest that a stable equilibrium state can exist only in restricted conditions: on a classical background spacetime that is globally hyperbolic or in a mild quantum-gravity regime in which there is an effective Schr\"odinger equation with a well-defined time parameter. These arguments suggest that quantum equilibrium will be unstable at the Planck scale. We construct a model in which quantum nonequilibrium is generated by a time-dependent regulator for pilot-wave dynamics, where the regulator is introduced to eliminate phase singularities. Applying our model to trans-Planckian modes that exit the Planck radius, we calculate the corrected primordial power spectrum and show that it displays a power excess (above a critical wavenumber). We briefly consider how our proposals could be tested by measurements of the cosmic microwave background.


title: Intermittence and time fractional stochastic partial differential equations

abstract: We consider time fractional stochastic heat type equation $$\partial^\beta_tu(t,x)=-\nu(-\Delta)^{\alpha/2} u_t(x)+I^{1-\beta}_t[\sigma(u)\stackrel{\cdot}{W}(t,x)]$$ in $(d+1)$ dimensions, where $\nu>0$, $\beta\in (0,1)$, $\alpha\in (0,2]$, $d<\min\{2,\beta^{-1}\}\a$, $\partial^\beta_t$ is the Caputo fractional derivative, $-(-\Delta)^{\alpha/2} $ is the generator of an isotropic stable process, $\stackrel{\cdot}{W}(t,x)$ is space-time white noise, and $\sigma:\RR{R}\to\RR{R}$ is Lipschitz continuous. The time fractional stochastic heat type equations might be used to model phenomenon with random effects with thermal memory. We prove: (i) absolute moments of the solutions of this equation grows exponentially; and (ii) the distances to the origin of the farthest high peaks of those moments grow exactly linearly with time. These results extend the results of Foondun and Khoshnevisan \cite{foondun-khoshnevisan-09} %(Mohammud Foondun and Davar Khoshnevisan, Intermittence and nonlinear parabolic %stochastic partial differential equations, Electron. J. Probab. 14 (2009), no. 21, 548--568) and Conus and Khoshnevisan \cite{conus-khoshnevisan} % (On the existence and position of the farthest peaks of a family of stochastic %heat and wave equations, Probab. Theory Related Fields 152 (2012), no. 3-4, 681--701) on the parabolic stochastic heat equations.


title: Neutrino Oscillation Physics Potential of the T2K Experiment

abstract: The observation of the recent electron neutrino appearance in a muon neutrino beam and the high-precision measurement of the mixing angle $\theta_{13}$ have led to a re-evaluation of the physics potential of the T2K long-baseline neutrino oscillation experiment. Sensitivities are explored for CP violation in neutrinos, non-maximal $\sin^22\theta_{23}$, the octant of $\theta_{23}$, and the mass hierarchy, in addition to the measurements of $\delta_{CP}$, $\sin^2\theta_{23}$, and $\Delta m^2_{32}$, for various combinations of $\nu$-mode and \(\bar{\nu}\)-mode data-taking. With an exposure of $7.8\times10^{21}$~protons-on-target, T2K can achieve 1-$\sigma$ resolution of 0.050(0.054) on $\sin^2\theta_{23}$ and $0.040(0.045)\times10^{-3}~\rm{eV}^2$ on $\Delta m^2_{32}$ for 100\%(50\%) neutrino beam mode running assuming $\sin^2\theta_{23}=0.5$ and $\Delta m^2_{32} = 2.4\times10^{-3}$ eV$^2$. T2K will have sensitivity to the CP-violating phase $\delta_{\rm{CP}}$ at 90\% C.L. or better over a significant range. For example, if $\sin^22\theta_{23}$ is maximal (i.e $\theta_{23}$=$45^\circ$) the range is $-115^\circ<\delta_{\rm{CP}}<-60^\circ$ for normal hierarchy and $+50^\circ<\delta_{\rm{CP}}<+130^\circ$ for inverted hierarchy. When T2K data is combined with data from the NO$\nu$A experiment, the region of oscillation parameter space where there is sensitivity to observe a non-zero $\delta_{CP}$ is substantially increased compared to if each experiment is analyzed alone.


title: On a conjectural solution to open KdV and Virasoro

abstract: In this note, we present a recursive formula for the full partition function Z of descendent integrals over moduli spaces of open and closed Riemann surfaces, assuming the conjecture recently proposed by Pandharipande, Solomon and Tessler that Z satisfies the open KdV and Virasoro equations.


title: The Double Exponential Sinc Collocation Method for Singular Sturm-Liouville Problems

abstract: Sturm-Liouville problems are abundant in the numerical treatment of scientific and engineering problems. In the present contribution, we present an efficient and highly accurate method for computing eigenvalues of singular Sturm-Liouville boundary value problems. The proposed method uses the double exponential formula coupled with Sinc collocation method. This method produces a symmetric positive-definite generalized eigenvalue system and has exponential convergence rate. Numerical examples are presented and comparisons with single exponential Sinc collocation method clearly illustrate the advantage of using the double exponential formula.


title: The Expected Optimal Labeling Order Problem for Crowdsourced Joins and Entity Resolution

abstract: In the SIGMOD 2013 conference, we published a paper extending our earlier work on crowdsourced entity resolution to improve crowdsourced join processing by exploiting transitive relationships [Wang et al. 2013]. The VLDB 2014 conference has a paper that follows up on our previous work [Vesdapunt et al., 2014], which points out and corrects a mistake we made in our SIGMOD paper. Specifically, in Section 4.2 of our SIGMOD paper, we defined the "Expected Optimal Labeling Order" (EOLO) problem, and proposed an algorithm for solving it. We incorrectly claimed that our algorithm is optimal. In their paper, Vesdapunt et al. show that the problem is actually NP-Hard, and based on that observation, propose a new algorithm to solve it. In this note, we would like to put the Vesdapunt et al. results in context, something we believe that their paper does not adequately do.