title: Protoplanetary disk lifetimes vs stellar mass and possible implications for giant planet populations

abstract: We study the dependence of protoplanetary disk evolution on stellar mass using a large sample of young stellar objects in nearby young star-forming regions. We update the protoplanetary disk fractions presented in our recent work (paper I of this series) derived for 22 nearby (< 500 pc) associations between 1 and 100 Myr. We use a subsample of 1 428 spectroscopically confirmed members to study the impact of stellar mass on protoplanetary disk evolution. We divide this sample into two stellar mass bins (2 M$_{\odot}$ boundary) and two age bins (3 Myr boundary), and use infrared excesses over the photospheric emission to classify objects in three groups: protoplanetary disks, evolved disks, and diskless. The homogeneous analysis and bias corrections allow for a statistically significant inter-comparison of the obtained results. We find robust statistical evidence of disk evolution dependence with stellar mass. Our results, combined with previous studies on disk evolution, confirm that protoplanetary disks evolve faster and/or earlier around high-mass (> 2 M$_{\odot}$) stars. We also find a roughly constant level of evolved disks throughout the whole age and stellar mass spectra. We conclude that protoplanetary disk evolution depends on stellar mass. Such a dependence could have important implications for gas giant planet formation and migration, and could contribute to explaining the apparent paucity of hot Jupiters around high-mass stars.


title: Hunting for Supermassive Black Holes in Nearby Galaxies with the Hobby-Eberly Telescope

abstract: We have conducted an optical long-slit spectroscopic survey of 1022 galaxies using the 10m Hobby-Eberly Telescope (HET) at McDonald Observatory. The main goal of the HET Massive Galaxy Survey (HETMGS) is to find nearby galaxies that are suitable for black hole mass measurements. In order to measure accurately the black hole mass, one should kinematically resolve the region where the black hole dominates the gravitational potential. For most galaxies, this region is much less than an arcsecond. Thus, black hole masses are best measured in nearby galaxies with telescopes that obtain high-spatial resolution. The HETMGS focuses on those galaxies predicted to have the largest sphere-of-influence, based on published stellar velocity dispersions or the galaxy fundamental plane. To ensure coverage over galaxy types, the survey targets those galaxies across a face-on projection of the fundamental plane. We present the sample selection and resulting data products from the long-slit observations, including central stellar kinematics and emission line ratios. The full dataset, including spectra and resolved kinematics, is available online. Additionally, we show that the current crop of black hole masses are highly biased towards dense galaxies and that especially large disks and low dispersion galaxies are under-represented. This survey provides the necessary groundwork for future systematic black hole mass measurement campaigns.


title: Made-to-Measure models of the Galactic Box/Peanut bulge: stellar and total mass in the bulge region

abstract: We construct dynamical models of the Milky Way's Box/Peanut (B/P) bulge, using the recently measured 3D density of Red Clump Giants (RCGs) as well as kinematic data from the BRAVA survey. We match these data using the NMAGIC Made-to-Measure method, starting with N-body models for barred discs in different dark matter haloes. We determine the total mass in the bulge volume of the RCGs measurement (+-2.2 x +- 1.4 x +- 1.2 kpc) with unprecedented accuracy and robustness to be 1.84 +- 0.07 x10^10 Msun. The stellar mass in this volume varies between 1.25-1.6 x10^10 Msun, depending on the amount of dark matter in the bulge. We evaluate the mass-to-light and mass-to-clump ratios in the bulge and compare them to theoretical predictions from population synthesis models. We find a mass-to-light ratio in the K-band in the range 0.8-1.1. The models are consistent with a Kroupa or Chabrier IMF, but a Salpeter IMF is ruled out for stellar ages of 10 Gyr. To match predictions from the Zoccali IMF derived from the bulge stellar luminosity function requires about 40% or 0.7 x10^10 Msun dark matter in the bulge region. The BRAVA data together with the RCGs 3D density imply a low pattern speed for the Galactic B/P bulge of 25-30 km.s-1.kpc-1. This would place the Galaxy among the slow rotators (R >= 1.5). Finally, we show that the Milky Way's B/P bulge has an off-centred X structure, and that the stellar mass involved in the peanut shape accounts for at least 20% of the stellar mass of the bulge, significantly larger than previously thought.


title: A New Channel for the Formation of Kinematically Decoupled Cores in Early-type galaxies

abstract: We present the formation of a Kinematically Decoupled Core (KDC) in an elliptical galaxy, resulting from a major merger simulation of two disk galaxies. We show that although the two progenitor galaxies are initially following a prograde orbit, strong reactive forces during the merger can cause a short-lived change of their orbital spin; the two progenitors follow a retrograde orbit right before their final coalescence. This results in a central kinematic decoupling and the formation of a large-scale (~2 kpc radius) counter-rotating core (CRC) at the center of the final elliptical-like merger remnant (M*=1.3x10^11 Msun), while its outer parts keep the rotation direction of the initial orbital spin. The stellar velocity dispersion distribution of the merger remnant galaxy exhibits two symmetrical off-centered peaks, comparable to the observed "2-sigma galaxies". The KDC/CRC consists mainly of old, pre-merger population stars (older than 5 Gyr), remaining prominent in the center of the galaxy for more than 2 Gyr after the coalescence of its progenitors. Its properties are consistent with KDCs observed in massive elliptical galaxies. This new channel for the formation of KDCs from prograde mergers is in addition to previously known formation scenarios from retrograde mergers and can help towards explaining the substantial fraction of KDCs observed in early-type galaxies.


title: Optimal analysis of the CMB trispectrum

abstract: We develop a general framework for data analysis and phenomenology of the CMB four-point function or trispectrum. To lowest order in the derivative expansion, the inflationary action admits three quartic operators consistent with symmetry: $\dot\sigma^4$, $\dot\sigma^2 (\partial\sigma^2)$, and $(\partial\sigma)^4$. In single field inflation, only the first of these operators can be the leading non-Gaussian signal. A Fisher matrix analysis shows that there is one near-degeneracy among the three CMB trispectra, so we parameterize the trispectrum with two coefficients $g_{NL}^{\dot\sigma^4}$ and $g_{NL}^{(\partial\sigma)^4}$, in addition to the coefficient $g_{NL}^{\rm loc}$ of $\zeta^3$-type local non-Gaussianity. This three-parameter space is analogous to the parameter space $(f_{NL}^{\rm loc}, f_{NL}^{\rm equil}, f_{NL}^{\rm orth})$ commonly used to parameterize the CMB three-point function. We next turn to data analysis and show how to represent these trispectra in a factorizable form which leads to computationally fast operations such as evaluating a CMB estimator or simulating a non-Gaussian CMB. We discuss practical issues in CMB analysis pipelines, and perform an optimal analysis of WMAP data. Our minimum-variance estimates are $g_{NL}^{\rm loc} = (-3.80 \pm 2.19) \times 10^5$, $g_{NL}^{\dot\sigma^4} = (-3.20 \pm 3.09) \times 10^6$, and $g_{NL}^{(\partial\sigma)^4} = (-10.8 \pm 6.33) \times 10^5$ after correcting for the effects of CMB lensing. No evidence of a nonzero inflationary four-point function is seen.


title: Adiabatic hydrodynamics: The eightfold way to dissipation

abstract: We provide a complete solution to hydrodynamic transport at all orders in the gradient expansion compatible with the second law constraint. The key new ingredient we introduce is the notion of adiabaticity, which allows us to take hydrodynamics off-shell. Adiabatic fluids are such that off-shell dynamics of the fluid compensates for entropy production. The space of adiabatic fluids is quite rich, and admits a decomposition into seven distinct classes. Together with the dissipative class this establishes the eightfold way of hydrodynamic transport. Furthermore, recent results guarantee that dissipative terms beyond leading order in the gradient expansion are agnostic of the second law. While this completes a transport taxonomy, we go on to argue for a new symmetry principle, an Abelian gauge invariance that guarantees adiabaticity in hydrodynamics. We suggest that this symmetry is the macroscopic manifestation of the microscopic KMS invariance. We demonstrate its utility by explicitly constructing effective actions for adiabatic transport. The theory of adiabatic fluids, we speculate, provides a useful starting point for a new framework to describe non-equilibrium dynamics, wherein dissipative effects arise by Higgsing the Abelian symmetry.

categories: hep-th cond-mat.stat-mech gr-qc physics.flu-dyn


title: The Metagalactic Ionizing Background: A Crisis in UV Photon Production or Incorrect Galaxy Escape Fractions?

abstract: Recent suggestions of a "photon underproduction crisis" (Kollmeier \etal\ 2014) have generated concern over the intensity and spectrum of ionizing photons in the metagalactic ultraviolet background (UVB). The balance of hydrogen photoionization and recombination determines the opacity of the low-redshift intergalactic medium (IGM). We calibrate the hydrogen photoionization rate ($\Gamma_{\rm H}$) by comparing {\it Hubble Space Telescope} spectroscopic surveys of the low-redshift column density distribution of \HI\ absorbers and the observed ($z < 0.4$) mean \Lya\ flux decrement, $D_A = (0.014)(1+z)^{2.2}$, to new cosmological simulations. The distribution, $f(N_{\rm HI}, z) \equiv d^2 {\cal N} / d(\log N_{\rm HI}) dz$, is consistent with an increased UVB that includes contributions from both quasars and galaxies. Our recommended fit, $\Gamma_{\rm H}(z) = (4.6 \times 10^{-14}$ s$^{-1})(1+z)^{4.4}$ for $0 < z < 0.47$, corresponds to unidirectional LyC photon flux $\Phi_0 \approx 5700$~cm$^{-2}$~s$^{-1}$ at $z = 0$. This flux agrees with observed IGM metal ionization ratios (\CIII/\CIV\ and \SiIII/\SiIV) and suggests a 25-30\% contribution of \Lya\ absorbers to the cosmic baryon inventory. The primary uncertainties in the low-redshift UVB are the contribution from massive stars in galaxies and the LyC escape fraction ($f_{\rm esc}$), a highly directional quantity that is difficult to constrain statistically. We suggest that both quasars and low-mass starburst galaxies are important contributors to the ionizing UVB at $z < 2$. Their additional ionizing flux would resolve any crisis in photon underproduction.


title: ANIR : Atacama Near-Infrared Camera for the 1.0-m miniTAO Telescope

abstract: We have developed a near-infrared camera called ANIR (Atacama Near-InfraRed camera) for the University of Tokyo Atacama Observatory 1.0m telescope (miniTAO) installed at the summit of Cerro Chajnantor (5640 m above sea level) in northern Chile. The camera provides a field of view of 5'.1 $\times$ 5'.1 with a spatial resolution of 0".298 /pixel in the wavelength range of 0.95 to 2.4 $\mu$m. Taking advantage of the dry site, the camera is capable of hydrogen Paschen-$\alpha$ (Pa$\alpha$, $\lambda=$1.8751 $\mu$m in air) narrow-band imaging observations, at which wavelength ground-based observations have been quite difficult due to deep atmospheric absorption mainly from water vapor. We have been successfully obtaining Pa$\alpha$ images of Galactic objects and nearby galaxies since the first-light observation in 2009 with ANIR. The throughputs at the narrow-band filters ($N1875$, $N191$) including the atmospheric absorption show larger dispersion (~10%) than those at broad-band filters (a few %), indicating that they are affected by temporal fluctuations in Precipitable Water Vapor (PWV) above the site. We evaluate the PWV content via the atmospheric transmittance at the narrow-band filters, and derive that the median and the dispersion of the distribution of the PWV are 0.40+/-0.30 mm for $N1875$ and 0.37+/-0.21 mm for $N191$, which are remarkably smaller (49+/-38% for $N1875$ and 59+/-26% for $N191$) than radiometry measurements at the base of Cerro Chajnantor (5100 m alt.). The decrease in PWV can be explained by the altitude of the site when we assume that the vertical distribution of the water vapor is approximated at an exponential profile with scale heights within 0.3-1.9 km (previously observed values at night). We thus conclude that miniTAO/ANIR at the summit of Cerro Chajnantor indeed provides us an excellent capability for a "ground-based" Pa$\alpha$ observation.


title: Reminiscences of my work with Richard Lewis Arnowitt

abstract: This article contains reminiscences of the collaborative work that Richard Arnowitt and I did together which stretched over many years and encompasses several areas of particle theory. The article is an extended version of my talk at the Memorial Symposium in honor of Richard Arnowitt at Texas A\&M, College Station, Texas, September 19-20, 2014.


title: ALMA constraints on the faint millimetre source number counts and their contribution to the cosmic infrared background

abstract: We have analysed 18 ALMA continuum maps in Bands 6 and 7, with rms down to 7.8$\mu$Jy, to derive differential number counts down to 60$\mu$Jy and 100$\mu$Jy at $\lambda=$1.3 mm and $\lambda=$1.1 mm, respectively. The area covered by the combined fields is $\rm 9.5\times10^{-4}deg^2$ at 1.1mm and $\rm 6.6\times10^{-4}deg^{2}$ at 1.3mm. We improved the source extraction method by requiring that the dimension of the detected sources be consistent with the beam size. This method enabled us to remove spurious detections that have plagued the purity of the catalogues in previous studies. We detected 50 faint sources with S/N$>$3.5 down to 60$\mu$Jy, hence improving the statistics by a factor of four relative to previous studies. The inferred differential number counts are $\rm dN/d(Log_{10}S)=1\times10^5~deg^2$ at a 1.1 mm flux $S_{\lambda = 1.1~mm} = 130~\mu$Jy, and $\rm dN/d(Log_{10}S)=1.1\times10^5~deg^2$ at a 1.3 mm flux $\rm S_{\lambda = 1.3~mm} = 60~\mu$Jy. At the faintest flux limits, i.e. 30$\mu$Jy and 40$\mu$Jy, we obtain upper limits on the differential number counts of $\rm dN/d(Log_{10}S) < 7\times10^5~deg^2$ and $\rm dN/d(Log_{10}S)<3\times10^5~deg^2$, respectively. Our results provide a new lower limit to CIB intensity of 17.2${\rm Jy\ deg^{-2}}$ at 1.1mm and of 12.9${\rm Jy\ deg^{-2}}$ at 1.3mm. Moreover, the flattening of the integrated number counts at faint fluxes strongly suggests that we are probably close to the CIB intensity. Our data imply that galaxies with SFR$<40~M_{\odot}/yr$ certainly contribute less than 50% to the CIB while more than 50% of the CIB must be produced by galaxies with $\rm SFR>40~M_{\odot}/yr$. The differential number counts are in nice agreement with recent semi-analytical models of galaxy formation even as low as our faint fluxes. Consequently, this supports the galaxy evolutionary scenarios and assumptions made in these models.


title: A discretized Chern-Simons gauge theory on arbitrary graphs

abstract: In this paper, we show how to discretize the abelian Chern-Simons gauge theory on generic planar lattices/graphs (with or without translational symmetries) embedded in arbitrary 2D closed orientable manifolds. We find that, as long as a one-to-one correspondence between vertices and faces can be defined on the graph such that each face is paired up with a neighboring vertex (and vice versa), a discretized Chern-Simons theory can be constructed consistently. We further verify that all the essential properties of the Chern-Simons gauge theory are preserved in the discretized setup. In addition, we find that the existence of such a one-to-one correspondence is not only a sufficient condition for discretizing a Chern-Simons gauge theory but, for the discretized theory to be nonsingular and to preserve some key properties of the topological field theory, this correspondence is also a necessary one. A specific example will then be provided, in which we discretize the Chern-Simons gauge theory on a tetrahedron.


title: A Chandra X-ray study of the mixed-morphology supernova remnant 3C400.2

abstract: We present an analysis of archival Chandra observations of the mixed-morphology remnant 3C400.2. We analysed spectra of different parts of the remnant to observe if the plasma properties provide hints on the origin of the mixed-morphology class. These remnants often show overionization, which is a sign of rapid cooling of the thermal plasma, and super-solar abundances of elements which is a sign of ejecta emission. Our analysis shows that the thermal emission of 3C400.2 can be well explained by a two component non-equilibrium ionization model, of which one component is underionized, has a high temperature ($kT \approx 3.9$ keV) and super-solar abundances, while the other component has a much lower temperature ($kT \approx 0.14$ keV), solar abundances and shows signs of overionization. The temperature structure, abundance values and density contrast between the different model components suggest that the hot component comes from ejecta plasma, while the cooler component has an interstellar matter origin. This seems to be the first instance of an overionized plasma found in the outer regions of a supernova remnant, whereas the ejecta component of the inner region is underionized. In addition, the non-ionization equilibrium plasma component associated with the ejecta is confined to the central, brighter parts of the remnant, whereas the cooler component is present mostly in the outer regions. Therefore our data can most naturally be explained by an evolutionary scenario in which the outer parts of the remnant are cooling rapidly due to having swept up high density ISM, while the inner parts are very hot and cooling inefficiently due to low density of the plasma. This is also known as the relic X-ray scenario.


title: BICEP2 / Keck Array V: Measurements of B-mode Polarization at Degree Angular Scales and 150 GHz by the Keck Array

abstract: The Keck Array is a system of cosmic microwave background (CMB) polarimeters, each similar to the BICEP2 experiment. In this paper we report results from the 2012 and 2013 observing seasons, during which the Keck Array consisted of five receivers all operating in the same (150 GHz) frequency band and observing field as BICEP2. We again find an excess of B-mode power over the lensed-$\Lambda$CDM expectation of $> 5 \sigma$ in the range $30 < \ell < 150$ and confirm that this is not due to systematics using jackknife tests and simulations based on detailed calibration measurements. In map difference and spectral difference tests these new data are shown to be consistent with BICEP2. Finally, we combine the maps from the two experiments to produce final Q and U maps which have a depth of 57 nK deg (3.4 $\mu$K arcmin) over an effective area of 400 deg$^2$ for an equivalent survey weight of 250,000 $\mu$K$^{-2}$. The final BB band powers have noise uncertainty a factor of 2.3 times better than the previous results, and a significance of detection of excess power of $> 6\sigma$.


title: Gamma-ray emission from young supernova remnants: hadronic or leptonic?

abstract: The debate on the nature of the gamma-ray emission from young supernova remnants is still open. Ascribing such emission to hadronic rather than leptonic processes would provide an evidence for the acceleration of protons and nuclei, and this fact would fit with the very popular (but not proven) paradigm that supernova remnants are the sources of Galactic cosmic rays. Here, we discuss this issue with a particular focus on the best studied gamma-ray-bright supernova remnant: RX~J1713.7-3946.


title: Evidence for jet launching close to the black hole in GRB 101219B - a Fermi GRB dominated by thermal emission

abstract: We present observations by the Fermi Gamma-Ray Space Telescope Gamma-Ray Burst Monitor (GBM) of the nearby (z=0.55) GRB 101219B. This burst is a long GRB, with an associated supernova and with a blackbody component detected in the early afterglow observed by the Swift X-ray Telescope (XRT). Here we show that the prompt gamma-ray emission has a blackbody spectrum, making this the second such burst observed by Fermi GBM. The properties of the blackbody, together with the redshift and our estimate of the radiative efficiency, makes it possible to calculate the absolute values of the properties of the outflow. We obtain an initial Lorentz factor Gamma=138\pm 8, a photospheric radius r_phot=4.4\pm 1.9 \times 10^{11} cm and a launch radius r_0=2.7\pm 1.6 \times 10^{7} cm. The latter value is close to the black hole and suggests that the jet has a relatively unobstructed path through the star. There is no smooth connection between the blackbody components seen by GBM and XRT, ruling out the scenario that the late emission is due to high-latitude effects. In the interpretation that the XRT blackbody is prompt emission due to late central engine activity, the jet either has to be very wide or have a clumpy structure where the emission originates from a small patch. Other explanations for this component, such as emission from a cocoon surrounding the jet, are also possible.


title: Constraining the Progenitor Companion of the Nearby Type Ia SN 2011fe with a Nebular Spectrum at +981 Days

abstract: We present an optical nebular spectrum of the nearby Type Ia supernova 2011fe, obtained 981 days after explosion. SN 2011fe exhibits little evolution since the +593 day optical spectrum, but there are several curious aspects in this new extremely late-time regime. We suggest that the persistence of the $\sim5800$~\AA\ feature is due to Na I D, and that a new emission feature at $\sim7300$~\AA\ may be [Ca II]. Also, we discuss whether the new emission feature at $\sim6400$~\AA\ might be [Fe I] or the high-velocity hydrogen predicted by Mazzali et al. The nebular feature at 5200~\AA\ exhibits a linear velocity evolution of $\sim350$ $\rm km\ s^{-1}$ per 100 days from at least +220 to +980 days, but the line's shape also changes in this time, suggesting that line blending contributes to the evolution. At $\sim 1000$ days after explosion, flux from the SN has declined to a point where contribution from a luminous secondary could be detected. In this work we make the first observational tests for a post-impact remnant star and constrain its temperature and luminosity to $T \gtrsim 10^4$ $\rm K$ and $L \lesssim 10^4$ $\rm L_{\odot}$. Additionally, we do not see any evidence for narrow H$\alpha$ emission in our spectrum. We conclude that observations continue to strongly exclude many single-degenerate scenarios for SN 2011fe.


title: Minimax Robust Hypothesis Testing

abstract: The minimax robust hypothesis testing problem for the case where the nominal probability distributions are subject to both modeling errors and outliers is studied in twofold. First, a robust hypothesis testing scheme based on a relative entropy distance is designed. This approach provides robustness with respect to modeling errors and is a generalization of a previous work proposed by Levy. Then, it is shown that this scheme can be combined with Huber's robust test through a composite uncertainty class, for which the existence of a saddle value condition is also proven. The composite version of the robust hypothesis testing scheme as well as the individual robust tests are extended to fixed sample size and sequential probability ratio tests. The composite model is shown to extend to robust estimation problems as well. Simulation results are provided to validate the proposed assertions.


title: Sign Structure of Susceptibilities of Conserved Charges in the (2+1) Polyakov Quark Meson Model

abstract: The sign structure of correlations of conserved charges are investigated in a QCD like model: the (2+1) flavor Polyakov Quark Meson model. We compute all susceptibilities of the conserved charges on the $(\mu_{B}-T)$ plane up to fourth order and a few at higher order as well. By varying the mass of the sigma meson, we are able to study and compare scenarios with as well as without a critical point. In the hadron-quark transition regime we identify certain correlations that turn negative unlike expectation from ideal hadron resonance gas calculations. The striking feature being that these remain negative deep into the hadronic side and thus could be measured in experiments. Measurement of such quantities in heavy ion collision experiments can elucidate the location of the QCD transition curve and possibly the critical point.


title: Notophs of N = 8 supergravity

abstract: We study the tensor gauge fields ("notophs") of ungauged N=8 D=4 supergravity in superspace. These are described by 2-form potentials in the adjoint representation of G=E(7(+7)). The consistency of the natural candidates for the superspace constraints for their field strengths fixes the form of the generalized Bianchi identities and also requires the potentials with indices of G/H=E(7(+7))/SU(8) generators to be dual to the scalars of the N=8,D=4 supergravity multiplet. In contrast, the field strengths of the 2-form potentials corresponding to the SU(8) generators are dual to fermionic bilinears, so that these potentials are auxiliary rather than physical fields. Their presence, however, is essential to formulate a tensor hierarchy of N=8,D=4 supergravity consistent with its U-duality group E(7(+7)).


title: Rigid inner forms vs isocrystals

abstract: We compare two statements of the refined local Langlands correspondence for connected reductive groups defined over a p-adic field -- one involving Kottwitz's set B(G) of isocrystals with additional structure, and one involving the cohomology set H^1(u -> W,Z -> G) introduced in arXiv:1304.3292. We show that if either statement is valid for all connected reductive groups, then so is the other. We also discuss how the second statement depends on the choice of element of H^1(u -> W,Z -> G).


title: Efficient Electro-Optical Modulation Based on Indium Tin Oxide

abstract: We experimentally demonstrate several electro-optical modulators based on transparent conducting oxides. Our previous work demonstrated the modulator structure on glass substrate with broadband bias polarity-dependent modulation. Further exploration shows similar modulation effect of the modulator on quartz and silicon substrate.


title: Learning the Matching Function

abstract: The matching function for the problem of stereo reconstruction or optical flow has been traditionally designed as a function of the distance between the features describing matched pixels. This approach works under assumption, that the appearance of pixels in two stereo cameras or in two consecutive video frames does not change dramatically. However, this might not be the case, if we try to match pixels over a large interval of time. In this paper we propose a method, which learns the matching function, that automatically finds the space of allowed changes in visual appearance, such as due to the motion blur, chromatic distortions, different colour calibration or seasonal changes. Furthermore, it automatically learns the importance of matching scores of contextual features at different relative locations and scales. Proposed classifier gives reliable estimations of pixel disparities already without any form of regularization. We evaluated our method on two standard problems - stereo matching on KITTI outdoor dataset, optical flow on Sintel data set, and on newly introduced TimeLapse change detection dataset. Our algorithm obtained very promising results comparable to the state-of-the-art.


title: Derivations and Observations of Prominence Bulk Motions and Mass

abstract: In this chapter we review observations and techniques for measuring both bulk flows in prominences and prominence mass. Measuring these quantities is essential to development and testing of models discussed throughout this book. Prominence flows are complex and various, ranging from the relatively linear flows along prominence spines to the complex, turbulent patterns exhibited by hedgerow prominences. Techniques for measuring flows include time slice and optical flow techniques used for motions in the plane of the sky and the use of spectral line profiles to determine Doppler velocities along the line of sight. Prominence mass measurement is chiefly done via continuum absorption measurements, but mass has also been estimated using cloud modeling and white light measurements.


title: Quantifying three dimensional reconnection in fragmented current layers

abstract: There is growing evidence that when magnetic reconnection occurs in high Lundquist number plasmas such as in the Solar Corona or the Earth's Magnetosphere it does so within a fragmented, rather than a smooth current layer. Within the extent of these fragmented current regions the associated magnetic flux transfer and energy release occurs simultaneously in many different places. This investigation focusses on how best to quantify the rate at which reconnection occurs in such layers. An analytical theory is developed which describes the manner in which new connections form within fragmented current layers in the absence of magnetic nulls. It is shown that the collective rate at which new connections form can be characterized by two measures; a total rate which measures the true rate at which new connections are formed and a net rate which measures the net change of connection associated with the largest value of the integral of $E_{\|}$ through all of the non-ideal regions. Two simple analytical models are presented which demonstrate how each should be applied and what they quantify.


title: Graviton Amplitudes from Collinear Limits of Gauge Amplitudes

abstract: We express all tree-level graviton amplitudes in Einstein's gravity as the collinear limits of a linear combination of pure Yang-Mills amplitudes in which each graviton is represented by two gauge bosons, each of them carrying exactly one half of graviton's momentum and helicity.


title: Alignment based Network Coding for Two-Unicast-Z Networks

abstract: In this paper, we study the wireline two-unicast-Z communication network over directed acyclic graphs. The two-unicast-Z network is a two-unicast network where the destination intending to decode the second message has apriori side information of the first message. We make three contributions in this paper: 1. We describe a new linear network coding algorithm for two-unicast-Z networks over directed acyclic graphs. Our approach includes the idea of interference alignment as one of its key ingredients. For graphs of a bounded degree, our algorithm has linear complexity in terms of the number of vertices, and polynomial complexity in terms of the number of edges. 2. We prove that our algorithm achieves the rate-pair (1, 1) whenever it is feasible in the network. Our proof serves as an alternative, albeit restricted to two-unicast-Z networks over directed acyclic graphs, to an earlier result of Wang et al. which studied necessary and sufficient conditions for feasibility of the rate pair (1, 1) in two-unicast networks. 3. We provide a new proof of the classical max-flow min-cut theorem for directed acyclic graphs.


title: The Structure of Dark Molecular Gas in the Galaxy - I: A Pilot Survey for 18-cm OH Emission Towards $l \approx 105^{\deg}, b \approx +1^{\deg}$

abstract: We report the first results from a survey for 1665, 1667, and 1720 MHz OH emission over a small region of the Outer Galaxy centered at $l \approx 105.0\deg , b \approx +1.0\deg$ . This sparse, high-sensitivity survey ($\Delta Ta \approx \Delta Tmb \approx 3.0 - 3.5$ mK rms in 0.55 km/s channels), was carried out as a pilot project with the Green Bank Telescope (GBT, FWHM $\approx 7.6'$) on a 3 X 9 grid at $0.5\deg$ spacing. The pointings chosen correspond with those of the existing $^{12}$CO(1-0) CfA survey of the Galaxy (FWHM $\approx 8.4'$). With 2-hr integrations, 1667 MHz OH emission was detected with the GBT at $\gtrsim 21$ of the 27 survey positions ($\geq 78\%$ ), confirming the ubiquity of molecular gas in the ISM as traced by this spectral line. With few exceptions, the main OH lines at 1665 and 1667 MHz appear in the ratio of 5:9 characteristic of LTE at our sensitivity levels. No OH absorption features are recorded in the area of the present survey, in agreement with the low levels of continuum background emission in this direction. At each pointing the OH emission appears in several components extending over a range of radial velocity and coinciding with well-known features of Galactic structure such as the Local Arm and the Perseus Arm. In contrast, little CO emission is seen in the survey area; less than half of the $\gtrsim 50$ identified OH components show detectable CO at the CfA sensitivity levels, and these are generally faint. There are no CO profiles without OH emission. With few exceptions, peaks in the OH profiles coincide with peaks in the GBT HI spectra (obtained concurrently, FWHM $8.9'$), although the converse is not true. We conclude that main-line OH emission is a promising tracer for the "dark molecular gas" in the Galaxy discovered earlier in Far-IR and gamma-ray emission. Further work is needed to establish the quantitative details of this connection.


title: Throwing Out the Baby with the Bathwater: The Undesirable Effects of National Research Assessment Exercises on Research

abstract: The evaluation of the quality of research at a national level has become increasingly common. The UK has been at the forefront of this trend having undertaken many assessments since 1986, the latest being the Research Excellence Framework in 2014. The argument of this paper is that, whatever the intended results in terms of evaluating and improving research, there have been many, presumably unintended, results that are highly undesirable for research and the university community more generally. We situate our analysis using Bourdieu's theory of cultural reproduction and then focus on the peculiarities of the 2008 RAE and the 2014 REF the rules of which allowed for, and indeed encouraged, significant game-playing on the part of striving universities. We conclude with practical recommendations to maintain the general intention of research assessment without the undesirable side-effects.


title: A new $F(R)$-gravity model

abstract: We propose a new model of modified $F(R)$ gravity theory with the function $F(R) = (1/\beta) \arcsin(\beta R)$. Constant curvature solutions corresponding to the flat and de Sitter spacetime are obtained. The Jordan and Einstein frames are considered; the potential and the mass of the scalar degree of freedom are found. We show that the flat spacetime is stable and the de Sitter spacetime is unstable. The slow-roll parameters $\epsilon$, $\eta$, and the $e$-fold number of the model are evaluated in the Einstein frame. The index of the scalar spectrum power-law $n_s$ and the tensor-to-scalar ratio $r$ are calculated. Critical points of autonomous equations for the de Sitter phase and the matter dominated epoch are found and studied. We obtain the approximate solution of equations of motion which is the deviation from the de Sitter phase in the Jordan frame. It is demonstrated that the model passes the matter stability test.


title: Machian space quanta

abstract: A new model for space and matter is obtained by joining every pair of point charges in the observable universe by an ethereal string. Positive gravitational potential energy in each string gives an attractive gravitational force due to the action of an energy conservation constraint. Newton's laws of motion are derived and inertia is explained in accordance with Mach's principle. The Machian string model gives a surprisingly simple way to understand the expansion history of the Universe. The decelerating expansion in the radiation era and the matter era is explained without using General Relativity and the transition from deceleration to acceleration is explained without the need to introduce a separate 'dark energy' component. The interaction between Machian strings gives a physical model for modified Newtonian dynamics (MOND) and is therefore an alternative to 'dark matter'.


title: Dynamics of Spatial Heterogeneity in Landfill - A Stochastic Analysis

abstract: A landfill represents a complex and dynamically evolving structure that can be stochastically perturbed by exogenous factors. Both thermodynamic (equilibrium) and time varying (non-steady state) properties of a landfill are affected by spatially heterogenous and nonlinear subprocesses that combine with constraining initial and boundary conditions arising from the associated surroundings. While multiple approaches have been made to model landfill statistics by incorporating spatially dependent parameters on the one hand (data based approach) and continuum dynamical mass-balance equations on the other (equation based modelling), practically no attempt has been made to amalgamate these two approaches while also incorporating inherent stochastically induced fluctuations affecting the process overall. In this article, we will implement a minimalist scheme of modelling the time evolution of a realistic three dimensional landfill through a reaction-diffusion based approach, focusing on the coupled interactions of four key variables - solid mass density, hydrolysed mass density, acetogenic mass density and methanogenic mass density. In a marked departure from previous predictions, our results indicate that close to the linearly stable limit, the large time steady state properties, arising out of a series of complex coupled interactions between the stochastically driven variables, are scarcely affected by the biochemical growth-decay statistics. Our results clearly show that an equilibrium landfill structure relates to plant production times of approximately 20-30 years instead of the previous (incorrect) deterministic model predictions of 50 years and above.


title: Candidate Quantum Spin Liquid in the Ce\textsuperscript{3+} Pyrochlore Stannate Ce$_2$Sn$_2$O$_7$

abstract: We report the low temperature magnetic properties of Ce$_2$Sn$_2$O$_7$, a rare-earth pyrochlore. Our susceptibility and magnetization measurements show that due to the thermal isolation of a Kramers doublet ground state, Ce$_2$Sn$_2$O$_7$ has Ising-like magnetic moments of $\sim1.18$ $\mu_\mathrm{B}$. The magnetic moments are confined to the local trigonal axes, as in a spin ice, but the exchange interactions are antiferromagnetic. Below 1 K the system enters a regime with antiferromagnetic correlations. In contrast to predictions for classical $\langle 111 \rangle$-Ising spins on the pyrochlore lattice, there is no sign of long-range ordering down to 0.02 K. Our results suggest that Ce$_2$Sn$_2$O$_7$ features an antiferromagnetic liquid ground state with strong quantum fluctuations.


title: Stability of Observations of Partial Differential Equations under Uncertain Perturbations

abstract: We analyse stability of observability estimates for solutions to wave and Scr\" odinger equations subjected to additive perturbations. The paper generalises the recent averaged observability/control result by allowing for systems consisting of operators of different types. The method also applies to the simultaneous observability problem by which one tries to estimate the energy of each component of a system under consideration. The analysis relies on microlocal defect tools; in particular on standard H-measures, when the main dynamic of the system is governed by the wave operator, while parabolic H-measures are explored in the case of the Schr\" odinger operator.


title: Analytical description of the structure of chaos

abstract: We consider analytical formulae that describe the chaotic regions around the main periodic orbit $(x=y=0)$ of the H\'{e}non map. Following our previous paper (Efthymiopoulos, Contopoulos, Katsanikas $2014$) we introduce new variables $(\xi, \eta)$ in which the product $\xi\eta=c$ (constant) gives hyperbolic invariant curves. These hyperbolae are mapped by a canonical transformation $\Phi$ to the plane $(x,y)$, giving "Moser invariant curves". We find that the series $\Phi$ are convergent up to a maximum value of $c=c_{max}$. We give estimates of the errors due to the finite truncation of the series and discuss how these errors affect the applicability of analytical computations. For values of the basic parameter $\kappa$ of the H\'{e}non map smaller than a critical value, there is an island of stability, around a stable periodic orbit $S$, containing KAM invariant curves. The Moser curves for $c \leq 0.32$ are completely outside the last KAM curve around $S$, the curves with $0.32<c<0.41$ intersect the last KAM curve and the curves with $0.41\leq c< c_{max} \simeq 0.49$ are completely inside the last KAM curve. All orbits in the chaotic region around the periodic orbit $(x=y=0)$, although they seem random, belong to Moser invariant curves, which, therefore define a "structure of chaos". Orbits starting close and outside the last KAM curve remain close to it for a stickiness time that is estimated analytically using the series $\Phi$. We finally calculate the periodic orbits that accumulate close to the homoclinic points, i.e. the points of intersection of the asymptotic curves from $x=y=0$, exploiting a method based on the self-intersections of the invariant Moser curves. We find that all the computed periodic orbits are generated from the stable orbit $S$ for smaller values of the H\'{e}non parameter $\kappa$, i.e. they are all regular periodic orbits.


title: Minimax estimation of linear and quadratic functionals on sparsity classes

abstract: For the Gaussian sequence model, we obtain non-asymptotic minimax rates of estimation of the linear, quadratic and the L2-norm functionals on classes of sparse vectors and construct optimal estimators that attain these rates. The main object of interest is the class s-sparse vectors for which we also provide completely adaptive estimators (independent of s and of the noise variance) having only logarithmically slower rates than the minimax ones. Furthermore, we obtain the minimax rates on the Lq-balls where 0 < q < 2. This analysis shows that there are, in general, three zones in the rates of convergence that we call the sparse zone, the dense zone and the degenerate zone, while a fourth zone appears for estimation of the quadratic functional. We show that, as opposed to estimation of the vector, the correct logarithmic terms in the optimal rates for the sparse zone scale as log(d/s^2) and not as log(d/s). For the sparse class, the rates of estimation of the linear functional and of the L2-norm have a simple elbow at s = sqrt(d) (boundary between the sparse and the dense zones) and exhibit similar performances, whereas the estimation of the quadratic functional reveals more complex effects and is not possible only on the basis of sparsity described by the sparsity condition on the vector. Finally, we apply our results on estimation of the L2-norm to the problem of testing against sparse alternatives. In particular, we obtain a non-asymptotic analog of the Ingster-Donoho-Jin theory revealing some effects that were not captured by the previous asymptotic analysis.


title: Negative probability

abstract: This article was written for the Logic in Computer Science column in the February 2015 issue of the Bulletin of the European Association for Theoretical Computer Science. The intended audience is general computer science audience. The uncertainty principle asserts a limit to the precision with which position x and momentum p of a particle can be known simultaneously. You may know the probability distributions of x and p individually but the joint distribution makes no physical sense. Yet Wigner exhibited such a joint distribution f(x,p). There was, however, a little trouble with it: some of its values were negative. Nevertheless Wigner's discovery attracted attention and found applications. There are other joint distribution, all with negative values, which produce the correct marginal distributions of x and p. But only Wigner's distribution produces the correct marginal distributions for all linear combinations of position and momentum. We offer a simple proof of the uniqueness and discuss related issues.


title: Galilean invariance at quantum Hall edge

abstract: We construct the theory of a chiral Luttinger liquid that lives on the boundary of a Galilean invariant quantum Hall fluid. In contrast to previous studies, Galilean invariance of the total (bulk plus edge) theory is guaranteed. We consider electromagnetic response at the edge and calculate momentum- and frequency-dependent electric conductivity and argue that its experimental measurement can provide a new means to determine the shift and bulk Hall viscosity.


title: Energy gap for Yang-Mills connections, II: Arbitrary closed Riemannian manifolds

abstract: In this sequel to [arXiv:1412.4114], we prove an $L^{d/2}$ energy gap result for Yang-Mills connections on principal $G$-bundles, $P$, over arbitrary, closed, Riemannian, smooth manifolds of dimension $d\geq 2$. We apply our version of the Lojasiewicz-Simon gradient inequality [arXiv:1409.1525, arXiv:1510.03815] to remove a positivity constraint on a combination of the Ricci and Riemannian curvatures in a previous $L^{d/2}$-energy gap result due to Gerhardt (2010) and a previous $L^\infty$-energy gap result due to Bourguignon, Lawson, and Simons (1981, 1979), as well as an $L^2$-energy gap result due to Nakajima (1987) for a Yang-Mills connection over the sphere, $S^d$, but with an arbitrary Riemannian metric. The main correction in this version involves replacement of the role of Corollary 4.3 due to Uhlenbeck (1985) and Theorem 5.1 due to the author in the published version of this article at http://dx.doi.org/10.1016/j.aim.2017.03.023 by that of Theorems 1 and 9 due to the author in arXiv:1906.03954.


title: On quantum computation, anyons, and categories

abstract: We explain the use of category theory in describing certain sorts of anyons. Yoneda's lemma leads to a simplification of that description. For the particular case of Fibonacci anyons, we also exhibit some calculations that seem to be known to the experts but not explicit in the literature.


title: Operator Positivity and Analytic Models of Commuting Tuples of Operators

abstract: We study analytic models of operators of class $C_{\cdot 0}$ with natural positivity assumptions. In particular, we prove that for an $m$-hypercontraction $T \in C_{\cdot 0}$ on a Hilbert space $\mathcal{H}$, there exists a Hilbert space $\mathcal{E}$ and a partially isometric multiplier $\theta \in \mathcal{M}(H^2(\mathcal{E}), A^2_m(\mathcal{H}))$ such that \[\mathcal{H} \cong \mathcal{Q}_{\theta} = A^2_m(\mathcal{H}) \ominus \theta H^2(\mathcal{E}), \quad \quad \mbox{and} \quad \quad T \cong P_{\mathcal{Q}_{\theta}} M_z|_{\mathcal{Q}_{\theta}},\]where $A^2_m$ is the weighted Bergman space and $H^2$ is the Hardy space over the unit disc $\mathbb{D}$. We then proceed to study and develop analytic models for doubly commuting $n$-tuples of operators and investigate their applications to joint shift co-invariant subspaces of reproducing kernel Hilbert spaces over polydisc. In particular, we completely analyze doubly commuting quotient modules of a large class of reproducing kernel Hilbert modules, in the sense of Arazy and Englis, over the unit polydisc $\mathbb{D}^n$.


title: The Radical Pair Mechanism and the Avian Chemical Compass: Quantum Coherence and Entanglement

abstract: We review the spin radical pair mechanism which is a promising explanation of avian navigation. This mechanism is based on the dependence of product yields on (1) the hyperfine interaction involving electron spins and neighboring nuclear spins and (2) the intensity and orientation of the geomagnetic field. One surprising result is that even at ambient conditions quantum entanglement of electron spins can play an important role in avian magnetoreception. This review describes the general scheme of chemical reactions involving radical pairs generated from singlet and triplet precursors; the spin dynamics of the radical pairs; and the magnetic field dependence of product yields caused by the radical pair mechanism. The main part of the review includes a description of the chemical compass in birds. We review: the general properties of the avian compass; the basic scheme of the radical pair mechanism; the reaction kinetics in cryptochrome; quantum coherence and entanglement in the avian compass; and the effects of noise. We believe that the "quantum avian compass" can play an important role in avian navigation and can also provide the foundation for a new generation of sensitive and selective magnetic-sensing nano-devices.


title: Conceptual design of first toroidal electron cyclotron resonance ion source and modeling of ion extraction from it

abstract: Electron Cyclotron Resonance Ion Sources (ECRIS) progressed to higher and higher ion currents and charge states by adopting stronger magnetic fields (beneficial for confinement) and proportionally higher ECR frequencies. Further improvements would require the attainment of "triple products" of density, temperature and confinement time comparable with major fusion experiments. For this, we propose a new, toroidal rather than linear, ECRIS geometry, which would at the same time improve confinement and make better use of the magnetic field. Ion extraction is more complicated than from a linear device, but feasible, as our modelling suggests: single-particle tracings showed successful extraction by at least two techniques, making use respectively of a magnetic extractor and of ExB drifts. Additional techniques are briefly discussed.


title: Mesoscale simulations of shockwave energy dissipation via chemical reactions

abstract: We use a particle-based mesoscale model that incorporates chemical reactions at a coarse-grained level to study the response of materials that undergo volume-reducing chemical reactions under shockwave-loading conditions. We find that such chemical reactions can attenuate the shockwave and characterize how the parameters of the chemical model affect this behavior. The simulations show that the magnitude of the volume collapse and velocity at which the chemistry propagates are critical to weaken the shock, whereas the energetics in the reactions play only a minor role. Shock loading results in transient states where the material is away from local equilibrium and, interestingly, chemical reactions can nucleate under such non-equilibrium states. Thus, the timescales for equilibration between the various degrees of freedom in the material affect the shock-induced chemistry and its ability to attenuate the propagating shock.


title: An equilibrium model for spot and forward prices of commodities

abstract: We consider a market model that consists of financial investors and producers of a commodity. Producers optionally store some production for future sale and go short on forward contracts to hedge the uncertainty of the future commodity price. Financial investors take positions in these contracts in order to diversify their portfolios. The spot and forward equilibrium commodity prices are endogenously derived as the outcome of the interaction between producers and investors. Assuming that both are utility maximizers, we first prove the existence of an equilibrium in an abstract setting. Then, in a framework where the consumers' demand and the exogenously priced financial market are correlated, we provide semi-explicit expressions for the equilibrium prices and analyze their dependence on the model parameters. The model can explain why increased investors' participation in forward commodity markets and higher correlation between the commodity and the stock market could result in higher spot prices and lower forward premia.


title: Heavy hadrons on $N_f=2$ and $2+1$ improved clover-Wilson lattices

abstract: We present the masses of singly ($B$, $B_s$, $\Lambda_b$, $\Sigma_b$, etc.), doubly ($B_c$, $\eta_b$, $\Upsilon$, $\Xi_{bc}$, $\Xi_{bb}$, etc.), and triply ($\Omega_{bcc}$, $\Omega_{bbc}$, $\Omega_{bbb}$, etc.) heavy hadrons arising from (QCDSF-UKQCD) lattices with improved clover-Wilson light quarks. For the bottom quark, we use an $O(a,v^4)$-improved version of lattice NRQCD. Part of the bottomonia spectrum is used to provide an alternative scale and to determine the physical quark mass and radiative corrections used in the heavy-quark action. Results for spin splittings, opposite parities, and, in some cases, excited states are presented. Higher lying states and baryons with two light quarks appear to be especially affected by the relatively small volumes of this (initially) initial study. This and other systematics are briefly discussed.


title: The XMM-Newton survey in the H-ATLAS field

abstract: Wide area X-ray and far infrared surveys are a fundamental tool to investigate the link between AGN growth and star formation, especially in the low-redshift universe (z<1). The Herschel Terahertz Large Area survey (H-ATLAS) has covered 550 deg^2 in five far-infrared and sub-mm bands, 16 deg^2 of which have been presented in the Science Demonstration Phase (SDP) catalogue. Here we introduce the XMM-Newton observations in H-ATLAS SDP area, covering 7.1 deg^2 with flux limits of 2e-15, 6e-15 and 9e-15 erg/s/cm^2 in the 0.5--2, 0.5--8 and 2--8 keV bands, respectively. We present the source detection and the catalogue, which includes 1700, 1582 and 814 sources detected by Emldetect in the 0.5--8, 0.5--2 and 2--8 keV bands, respectively; the number of unique sources is 1816. We extract spectra and derive fluxes from power-law fits for 398 sources with more than 40 counts in the 0.5--8 keV band. We compare the best-fit fluxes with the catalogue ones, obtained by assuming a common photon index of Gamma=1.7; we find no bulk difference between the fluxes, and a moderate dispersion of s=0.33 dex. Using wherever possible the fluxes from the spectral fits, we derive the 2--10 keV LogN-LogS, which is consistent with a Euclidean distribution. Finally, we release computer code for the tools developed for this project.


title: The Crab Pulsar at Centimeter Wavelengths: I. Ensemble Characteristics

abstract: We have observed the pulsar in the Crab Nebula at high radio frequencies and high time resolution. We present continuously sampled data at 640-ns time resolution, and individual bright pulses recorded at down to 0.25-ns time resolution. Combining our new data with previous data from our group and from the literature shows the dramatic changes in the pulsar's radio emission between low and high radio frequencies. Below about 5 GHz the mean profile is dominated by the bright Main Pulse and Low-Frequency Interpulse. Everything changes, however, above about 5 GHz; the Main Pulse disappears, the mean profile of the Crab pulsar is dominated by the High-Frequency Interpulse (which is quite different from its low-frequency counterpart) and the two High-Frequency Components. We present detailed observational characteristics of these different components which future models of the pulsar's magnetosphere must explain.


title: Geometry of the intersection ring and vanishing relations in the cohomology of the moduli space of parabolic bundles on a curve

abstract: We study the ring generated by the Chern classes of tautological line bundles on the moduli space of parabolic bundles of arbitrary rank on a Riemann surface. We show the Poincar\'e duals to these Chern classes have simple geometric representatives. We use this construction to show that the ring generated by these Chern classes vanishes below the dimension of the moduli space, in analogy with the Newstead-Ramanan conjecture for stable bundles.


title: Competition and Coalition Formation of Renewable Power Producers

abstract: We investigate group formations and strategic behaviors of renewable power producers in electricity markets. These producers currently bid into the day-ahead market in a conservative fashion because of the real-time risk associated with not meeting their bid amount. It has been suggested in the literature that producers would bid less conservatively if they can form large groups to take advantages of spatial diversity to reduce the uncertainty in their aggregate output. We show that large groups of renewable producers would act strategically to lower the aggregate output because of market power. To maximize renewable power production, we characterize the trade-off between market power and generation uncertainty as a function of the size of the groups. We show there is a sweet spot in the sense that there exists groups that are large enough to achieve the uncertainty reduction of the grand coalition, but are small enough such that they have no significant market power.We consider both independent and correlated forecast errors under a fixed real-time penalty. We also consider a real-time market where both selling and buying of energy are allowed. We validate our claims using PJM and NREL data.


title: Quasi-Centralized Limit Order Books

abstract: A quasi-centralized limit order book (QCLOB) is a limit order book (LOB) in which financial institutions can only access the trading opportunities offered by counterparties with whom they possess sufficient bilateral credit. We perform an empirical analysis of a recent, high-quality data set from a large electronic trading platform that utilizes QCLOBs to facilitate trade. We find many significant differences between our results and those widely reported for other LOBs. We also uncover a remarkable empirical universality: although the distributions describing order flow and market state vary considerably across days, a simple, linear rescaling causes them to collapse onto a single curve. Motivated by this finding, we propose a semi-parametric model of order flow and market state in a QCLOB on a single trading day. Our model provides similar performance to that of parametric curve-fitting techniques, while being simpler to compute and faster to implement.

categories: q-fin.TR nlin.AO physics.data-an physics.soc-ph stat.AP


title: Quantum enhanced estimation of optical detector efficiencies

abstract: Quantum mechanics establishes the ultimate limit to the scaling of the precision on any parameter, by iden- tifying optimal probe states and measurements. While this paradigm is, at least in principle, adequate for the metrology of quantum channels involving the estimation of phase and loss parameters, we show that estimat- ing the loss parameters associated with a quantum channel and a realistic quantum detector are fundamentally different. While Fock states are provably optimal for the former, we identify a crossover in the nature of the optimal probe state for estimating detector imperfections as a function of the loss parameter. We provide explicit results for on-off and homodyne detectors, the most widely used detectors in quantum photonics technologies.


title: Stern-Brocot Trees from Weighted Mediants

abstract: In this paper we discuss a natural generalization of the Stern Brocot tree which comes from the introduction of weighted mediants. We focus our attention on the case $k = 3$, in which $(2a + c)/(2b + d)$ and $(a + 2c)/(b + 2d)$ are the two mediants inserted between $a/b$ and $c/d$. Our main result is a determination of which rational numbers between the starting terms appear in the tree. We extend this result to arbitrary reduction schemes as well.


title: Critical String from Non-Abelian Vortex in Four Dimensions

abstract: In a class of non-Abelian solitonic vortex strings supported in certain N=2 super-Yang-Mills theories we search for the vortex which can behave as a critical fundamental string. We use the Polchinski-Strominger criterion of the ultraviolet completeness. We identify an appropriate four-dimensional bulk theory: it has the U(2) gauge group, the Fayet-Iliopoulos term and four flavor hypermultiplets. It supports semilocal vortices with the world-sheet theory for orientational (size) moduli described by the weighted CP(2,2) model. The latter is superconformal. Its target space is six-dimensional. The overall Virasoro central charge is critical. We show that the world-sheet theory on the vortex supported in this bulk model is the bona fide critical string.


title: Average Transmission Decrease in One-Dimensional Photonic Structures by Widening the Random Layer Thickness Distribution

abstract: In this work we have studied the optical properties of disordered photonic structures, in which we have controlled the distribution of the random layer thickness. Such structures are characterized by an alternation of high and low refractive index layers, but the layer thicknesses follow the aforementioned distributions. We have used two types of distribution: a distribution in which each thickness has the same probability to occur and one in which the thickness follows a Gaussian function. We have simulated the average transmission all over the spectrum for photonic structure characterized by a different width of the distribution. We have found that the choice of the distribution of the layer thickness is a control of the average transmission of a random photonic structure.


title: Disentangling defects and sound modes in disordered solids

abstract: We develop a new method to isolate localized defects from extended vibrational modes in disordered solids. This method augments particle interactions with an artificial potential that acts as a high-pass filter: it preserves small-scale structures while pushing extended vibrational modes to higher frequencies. The low-frequency modes that remain are "bare" defects; they are exponentially localized without the quadrupolar tails associated with elastic interactions. We demonstrate that these localized excitations are excellent predictors of plastic rearrangements in the solid. We characterize several of the properties of these defects that appear in mesoscopic theory of plasticity, including their distribution of energy barriers, number density, and size, which is a first step in testing and revising continuum models for plasticity in disordered solids.


title: Renormalization of spectra by phase competition in the half-filled Hubbard-Holstein model

abstract: We present electron and phonon spectral functions calculated from determinant quantum Monte Carlo simulations of the half-filled two-dimensional Hubbard-Holstein model on a square lattice. By tuning the relative electron-electron ($e$-$e$) and electron-phonon ($e$-$ph$) interaction strengths, we show the electron spectral function evolving between antiferromagnetic insulating, metallic, and charge density wave insulating phases. The phonon spectra concurrently gain a strong momentum dependence and soften in energy upon approaching the charge density wave phase. In particular, we study how the $e$-$e$ and $e$-$ph$ interactions renormalize the spectra, and analyze how the interplay of these interactions influence the spectral renormalizations. We find that the presence of both interactions suppresses the amount of renormalization at low energy, thus allowing the emergence of a metallic phase. These findings demonstrate the importance of considering the influence of multiple interactions in spectroscopically determining any one interaction strength in strongly correlated materials.


title: Global infinite energy solutions for the 2D gravity water waves system

abstract: We prove global existence and modified scattering property for the solutions of the $2D$ gravity water waves system in the infinite depth setting for a class of initial data, which is only required to be small above the level $\dot{H}^{1/5}\times \dot{H}^{1/5+1/2}$. No assumption is assumed below this level, therefore, it allows to have infinite energy. As a direct consequence, the momentum condition assumed on the physical velocity in all previous small energy results by Ionescu-Pusateri, Alazard-Delort and Ifrim-Tataru is removed.


title: Charged black holes in expanding Einstein-de Sitter universes

abstract: Inspired in a previous work by McClure and Dyer (Classical Quantum Gravity 23, 1971 (2006)), we analyze some solutions of the Einstein-Maxwell equations which were originally written to describe charged black holes in cosmological backgrounds. A detailed analysis of the electromagnetic sources for a sufficiently general metric is performed, and then we focus on deriving the electromagnetic four-current as well as the conserved electric charge of each metric. The charged McVittie solution is revisited and a brief study of its causal structure is performed, showing that it may represent a charged black hole in an expanding universe, with the black hole horizon being formed at infinite late times. Charged versions of solutions originally put forward by Vaidya (Vd) and Sultana and Dyer (SD) are also analyzed. It is shown that the charged Sultana-Dyer metric requires a global electric current, besides a central (pointlike) electric charge. With the aim of comparing to the charged McVittie metric, new charged solutions of Vd and SD type are considered. In these cases, the original mass and charge parameters are replaced by particular functions of the cosmological time. In the new generalized charged Vaidya metric the black hole horizon never forms, whereas in the new generalized Sultana-Dyer case both the Cauchy and the black hole horizons develop at infinite late times. A charged version of the Thakurta metric is also studied here. It is also a new solution. As in the charged Sultana-Dyer case, the natural source of the electromagnetic field is a central electric charge with an additional global electric current. The global structure is briefly studied and it is verified that the corresponding spacetime may represent a charged black hole in a cosmological background. All the solutions present initial singularities as found in the McVittie metric.


title: Finite cyclicity of some graphics through a nilpotent point of saddle type inside quadratic systems

abstract: In this paper we show the finite cyclicity of the two graphics $(I_{12}^1)$ and $(I_{13}^1)$ through a triple nilpotent point of saddle type inside quadratic vector fields. These results contribute to the program launched in 1994 by Dumortier, Roussarie and Rousseau (DRR program) to show the existence of a uniform upper bound for the number of limit cycles for planar quadratic vector fields.


title: Coupling Human Mobility and Social Ties

abstract: Studies using massive, passively data collected from communication technologies have revealed many ubiquitous aspects of social networks, helping us understand and model social media, information diffusion, and organizational dynamics. More recently, these data have come tagged with geographic information, enabling studies of human mobility patterns and the science of cities. We combine these two pursuits and uncover reproducible mobility patterns amongst social contacts. First, we introduce measures of mobility similarity and predictability and measure them for populations of users in three large urban areas. We find individuals' visitations patterns are far more similar to and predictable by social contacts than strangers and that these measures are positively correlated with tie strength. Unsupervised clustering of hourly variations in mobility similarity identifies three categories of social ties and suggests geography is an important feature to contextualize social relationships. We find that the composition of a user's ego network in terms of the type of contacts they keep is correlated with mobility behavior. Finally, we extend a popular mobility model to include movement choices based on social contacts and compare it's ability to reproduce empirical measurements with two additional models of mobility.


title: Massless Neutrino Oscillations via Quantum Tunneling

abstract: In the present theory, the masses of neutrinos are required in order to realize oscillations. By analogy with the oscillation of quantum two-state system, we assume that neutrino oscillations may be regarded as quantum tunneling process. The difference of the quantum numbers between two particles may build a barrier. Thus neutrinos with vanishing mass will oscillate. The hypothesis can also be applied to the lepton flavor violation and quark mixing as well as the oscillation of neutral kaons. We point out that quark mixing is the consequence of quark oscillations via barrier tunneling in hadrons as bound states, while neutrinos oscillate as free particles. We expect that the recent experimental results on anomalous magnetic moment of the muon and the breaking of lepton universality in beauty-quark decays will be explained with flavor-changing via barrier tunneling instead of introducing new interactions and particles.


title: Stable cheapest nonconforming finite elements for the Stokes equations

abstract: We introduce two pairs of stable cheapest nonconforming finite element space pairs to approximate the Stokes equations. One pair has each component of its velocity field to be approximated by the $P_1$ nonconforming quadrilateral element while the pressure field is approximated by the piecewise constant function with globally two-dimensional subspaces removed: one removed space is due to the integral mean--zero property and the other space consists of global checker--board patterns. The other pair consists of the velocity space as the $P_1$ nonconforming quadrilateral element enriched by a globally one--dimensional macro bubble function space based on $DSSY$ (Douglas-Santos-Sheen-Ye) nonconforming finite element space; the pressure field is approximated by the piecewise constant function with mean--zero space eliminated. We show that two element pairs satisfy the discrete inf-sup condition uniformly. And we investigate the relationship between them. Several numerical examples are shown to confirm the efficiency and reliability of the proposed methods.


title: Deformation Classification of Typical Configurations of 7 Points in the Real Projective Plane

abstract: A configuration of 7 points in RP2 is called typical if it has no collinear triples and no coconic sextuples of points. We show that there exist 14 deformation classes of such configurations. This yields classification of real Aronhold sets.


title: Lusin-type theorems for Cheeger derivatives on metric measure spaces

abstract: A theorem of Lusin states that every Borel function on $R$ is equal almost everywhere to the derivative of a continuous function. This result was later generalized to $R^n$ in works of Alberti and Moonens-Pfeffer. In this note, we prove direct analogs of these results on a large class of metric measure spaces, those with doubling measures and Poincar\'e inequalities, which admit a form of differentiation by a famous theorem of Cheeger.


title: Assessing the Fault Proneness Degree (DFP) by Estimating the Impact of Change Request Artifacts Correlation

abstract: Exploring the impact of change requests applied to a software maintenance project helps to assess the fault-proneness of the change request to be handled further, which is perhaps a bug fix or even a new feature demand. In practice, the major development community stores change requests and related data using bug tracking systems such as Bugzilla. These data, together with the data stored in a versioning system, such as Concurrent Versioning Systems, are a valuable source of information to create descriptions and also can perform useful analyzes. In our earlier work, we proposed a novel statistical bipartite weighted graph-based approach to assessing the degree of fault-proneness of the change request and Change Request artifacts. With the motivation gained from this model, here we propose a novel strategy that estimates the degree of fault-proneness of a change request by assessing the impact of a change request artifact towards fault-proneness that considers the correlation between change requests artifact as another factor, which is in addition to our earlier strategy. The proposed model can be titled as Assessing the Fault Proneness Degree of Change Request Artifacts by estimating the impact of change requests correlation (DFP-CRC). As stated in our earlier model, the method DFP-CRC also makes use of information retrieval methods to identify the change request artifacts of the devised change request. And further evaluates the degree of fault-proneness of the Change Requests by estimating the correlation between change requests. The proposed method is evaluated by applying on concurrent versioning and Change request logs of the production level maintenance project.


title: Lebesgue spaces norm estimates for fractional integrals and derivatives

abstract: We study the problem estimation of classical Lebesgue-Riesz and Grand Lebesgue Norm for the fractional integrals and derivatives for the functions from the classical Lebesgue-Riesz spaces as well as from the modified Besov's spaces.


title: Gap Eigenvalues and Asymptotic Dynamics of Geometric Wave Equations on Hyperbolic Space

abstract: In this paper we study $k$-equivariant wave maps from the hyperbolic plane into the $2$-sphere as well as the energy critical equivariant $SU(2)$ Yang-Mills problem on $4$-dimensional hyperbolic space. The latter problem bears many similarities to a $2$-equivariant wave map into a surface of revolution. As in the case of $1$-equivariant wave maps considered in~\cite{LOS1}, both problems admit a family of stationary solutions indexed by a parameter that determines how far the image of the map wraps around the target manifold. Here we show that if the image of a stationary solution is contained in a geodesically convex subset of the target, then it is asymptotically stable in the energy space. However, for a stationary solution that covers a large enough portion of the target, we prove that the Schr\"odinger operator obtained by linearizing about such a harmonic map admits a simple positive eigenvalue in the spectral gap. As there is no a priori nonlinear obstruction to asymptotic stability, this gives evidence for the existence of metastable states (i.e., solutions with anomalously slow decay rates) in these simple geometric models.


title: Large-angle quasi-self-collimation effect in a rod-type photonic crystal

abstract: A rod-type photonic crystal (PC) with a rectangular lattice shows a large-angle quasi-self-collimation (quasi-SC) effect by changing the symmetry of its rectangular lattice to straighten one of the isofrequency contours. To investigate the straightness of the isofrequency contour as well as the quasi-SC effect, we propose a straightness factor L based on the method of least squares. With L smaller than L0 (L0 = 0.01 is the critical value), the isofrequency contour is sufficiently straight to induce quasi-SC effect with the beam quasi-collimating in the structure. Furthermore, the efficiency of light coupling to the quasi-SC PC is studied, and can be greatly improved by applying a carefully designed antireflection layer.


title: On random subgraphs of Kneser and Schrijver graphs

abstract: A Kneser graph $KG_{n,k}$ is a graph whose vertices are in one-to-one correspondence with $k$-element subsets of $[n],$ with two vertices connected if and only if the corresponding sets do not intersect. A famous result due to Lov\'asz states that the chromatic number of a Kneser graph $KG_{n,k}$ is equal to $n-2k+2$. In this paper we study the chromatic number of a random subgraph of a Kneser graph $KG_{n,k}$ as $n$ grows. A random subgraph $KG_{n,k}(p)$ is obtained by including each edge of $KG_{n,k}$ with probability $p$. For a wide range of parameters $k = k(n), p = p(n)$ we show that $\chi(KG_{n,k}(p))$ is very close to $\chi(KG_{n,k}),$ a.a.s. differing by at most 4 in many cases. Moreover, we obtain the same bounds on the chromatic numbers for the so-called Schrijver graphs, which are known to be vertex-critical induced subgraphs of Kneser graphs.


title: Ab initio effective interactions for sd-shell valence nucleons

abstract: We perform \textit{ab initio} no-core shell-model calculations for $A=18$ and $19$ nuclei in a $4\hbar\Omega$, or $N_{\rm max}=4$, model space by using the effective JISP16 and chiral N3LO nucleon-nucleon potentials and transform the many-body effective Hamiltonians into the $0\hbar\Omega$ model space to construct the $A$-body effective Hamiltonians in the $sd$-shell. We separate the $A$-body effective Hamiltonians with $A=18$ and $A=19$ into inert core, one-, and two-body components. Then, we use these core, one-, and two-body components to perform standard shell-model calculations for the $A=18$ and $A=19$ systems with valence nucleons restricted to the $sd$ shell. Finally, we compare the standard shell-model results in the $0\hbar\Omega$ model space with the exact no-core shell model results in the $4\hbar\Omega$ model space for the $A=18$ and $A=19$ systems and find good agreement.


title: Phonon and Raman scattering of two-dimensional transition metal dichalcogenides from monolayer, multilayer to bulk material

abstract: Two-dimensional (2D) transition metal dichalcogenide (TMD) nanosheets exhibit remarkable electronic and optical properties. The 2D features, sizable bandgaps, and recent advances in the synthesis, characterization, and device fabrication of the representative MoS$_2$, WS$_2$, WSe$_2$, and MoSe$_2$ TMDs make TMDs very attractive in nanoelectronics and optoelectronics. Similar to graphite and graphene, the atoms within each layer in 2D TMDs are joined together by covalent bonds, while van der Waals interactions keep the layers together. This makes the physical and chemical properties of 2D TMDs layer dependent. In this review, we discuss the basic lattice vibrations of monolayer, multilayer, and bulk TMDs, including high-frequency optical phonons, interlayer shear and layer breathing phonons, the Raman selection rule, layer-number evolution of phonons, multiple phonon replica, and phonons at the edge of the Brillouin zone. The extensive capabilities of Raman spectroscopy in investigating the properties of TMDs are discussed, such as interlayer coupling, spin--orbit splitting, and external perturbations. The interlayer vibrational modes are used in rapid and substrate-free characterization of the layer number of multilayer TMDs and in probing interface coupling in TMD heterostructures. The success of Raman spectroscopy in investigating TMD nanosheets paves the way for experiments on other 2D crystals and related van der Waals heterostructures.


title: Hybrid Orthogonal Projection and Estimation (HOPE): A New Framework to Probe and Learn Neural Networks

abstract: In this paper, we propose a novel model for high-dimensional data, called the Hybrid Orthogonal Projection and Estimation (HOPE) model, which combines a linear orthogonal projection and a finite mixture model under a unified generative modeling framework. The HOPE model itself can be learned unsupervised from unlabelled data based on the maximum likelihood estimation as well as discriminatively from labelled data. More interestingly, we have shown the proposed HOPE models are closely related to neural networks (NNs) in a sense that each hidden layer can be reformulated as a HOPE model. As a result, the HOPE framework can be used as a novel tool to probe why and how NNs work, more importantly, to learn NNs in either supervised or unsupervised ways. In this work, we have investigated the HOPE framework to learn NNs for several standard tasks, including image recognition on MNIST and speech recognition on TIMIT. Experimental results have shown that the HOPE framework yields significant performance gains over the current state-of-the-art methods in various types of NN learning problems, including unsupervised feature learning, supervised or semi-supervised learning.


title: Study of in-medium mass modification at J-PARC

abstract: Study of in-medium mass modification has attracted interest in terms of the restoration of the spontaneously broken chiral symmetry, which is responsible for the generation of hadron mass. Many experiments were performed to measure in-medium property of hadrons but there is no consensus yet. J-PARC E16 has been proposed to study in-medium property of vector mesons via dilepton decay channel. The status of spectrometer R \& D is explained. Other related experiments planned at J-PARC are also introduced.


title: Nonlinear X-ray Compton Scattering

abstract: X-ray scattering is a weak linear probe of matter. It is primarily sensitive to the position of electrons and their momentum distribution. Elastic X-ray scattering forms the basis of atomic structural determination while inelastic Compton scattering is often used as a spectroscopic probe of both single-particle excitations and collective modes. X-ray free-electron lasers (XFELs) are unique tools for studying matter on its natural time and length scales due to their bright and coherent ultrashort pulses. However, in the focus of an XFEL the assumption of a weak linear probe breaks down, and nonlinear light-matter interactions can become ubiquitous. The field can be sufficiently high that even non-resonant multiphoton interactions at hard X-rays wavelengths become relevant. Here we report the observation of one of the most fundamental nonlinear X-ray-matter interactions, the simultaneous Compton scattering of two identical photons producing a single photon at nearly twice the photon energy. We measure scattered photons with an energy near 18 keV generated from solid beryllium irradiated by 8.8-9.75 keV XFEL pulses. The intensity in the X-ray focus reaches up to 4x20 W/cm2, which corresponds to a peak electric field two orders of magnitude higher than the atomic unit of field-strength and within four orders of magnitude of the quantum electrodynamic critical field. The observed signal scales quadratically in intensity and is emitted into a non-dipolar pattern, consistent with the simultaneous two-photon scattering from free electrons. However, the energy of the generated photons shows an anomalously large redshift only present at high intensities. This indicates that the instantaneous high-intensity scattering effectively interacts with a different electron momentum distribution than linear Compton scattering, with implications for the study of atomic-scale structure and dynamics of matter


title: Recovery of Piecewise Smooth Images from Few Fourier Samples

abstract: We introduce a Prony-like method to recover a continuous domain 2-D piecewise smooth image from few of its Fourier samples. Assuming the discontinuity set of the image is localized to the zero level-set of a trigonometric polynomial, we show the Fourier transform coefficients of partial derivatives of the signal satisfy an annihilation relation. We present necessary and sufficient conditions for unique recovery of piecewise constant images using the above annihilation relation. We pose the recovery of the Fourier coefficients of the signal from the measurements as a convex matrix completion algorithm, which relies on the lifting of the Fourier data to a structured low-rank matrix; this approach jointly estimates the signal and the annihilating filter. Finally, we demonstrate our algorithm on the recovery of MRI phantoms from few low-resolution Fourier samples.


title: Towards age/rotation/magnetic activity relation with seismology

abstract: The knowledge of stellar ages directly impacts the characterization of a planetary system as it puts strong constraints on the moment when the system was born. Unfortunately, the determination of precise stellar ages is a very difficult task. Different methods can be used to do so (based on isochrones or chemical element abundances) but they usually provide large uncertainties. During its evolution a star goes through processes leading to loss of angular momentum but also changes in its magnetic activity. Building rotation, magnetic, age relations would be an asset to infer stellar ages model independently. Several attempts to build empirical relations between rotation and age (namely gyrochronology) were made with a focus on cluster stars where the age determination is easier and for young stars on the main sequence. For field stars, we can now take advantage of high-precision photometric observations where we can perform asteroseismic analyses to improve the accuracy of stellar ages. Furthermore, the variability in the light curves allow us to put strong constraints on the stellar rotation and magnetic activity. By combining these precise measurements, we are on the way of understanding and improving relations between magnetic activity, rotation, and age, in particular at different stages of stellar evolution. I will review the status on gyrochronology relationships based on observations of young cluster stars. Then I will focus on solar-like stars and describe the inferences on stellar ages, rotation, and magnetism that can be provided by high-quality photometric observations such as the ones of the Kepler mission, in particular through asteroseismic analyses.


title: Searching for quantum optimal controls under severe constraints

abstract: The success of quantum optimal control for both experimental and theoretical objectives is connected to the topology of the corresponding control landscapes, which are free from local traps if three conditions are met: (1) the quantum system is controllable, (2) the Jacobian of the map from the control field to the evolution operator is of full rank, and (3) there are no constraints on the control field. This paper investigates how the violation of assumption (3) affects gradient searches for globally optimal control fields. The satisfaction of assumptions (1) and (2) ensures that the control landscape lacks fundamental traps, but certain control constraints can still introduce artificial traps. Proper management of these constraints is an issue of great practical importance for numerical simulations as well as optimization in the laboratory. Using optimal control simulations, we show that constraints on quantities such as the number of control variables, the control duration, and the field strength are potentially severe enough to prevent successful optimization of the objective. For each such constraint, we show that exceeding quantifiable limits can prevent gradient searches from reaching a globally optimal solution. These results demonstrate that careful choice of relevant control parameters helps to eliminate artificial traps and facilitate successful optimization.


title: Vizing's Conjecture for Almost All Pairs of Graphs

abstract: For any graph $G=(V,E)$, a subset $S\subseteq V$ $dominates$ $G$ if all vertices are contained in the closed neighborhood of $S$, that is $N[S]=V$. The minimum cardinality over all such $S$ is called the domination number, written $\gamma(G)$. In 1963, V.G. Vizing conjectured that $\gamma(G \square H) \geq \gamma(G)\gamma(H)$ where $\square$ stands for the Cartesian product of graphs. In this note, we prove that if $\left|G\right|\geq \gamma(G)\gamma(H)$ and $\left|H\right|\geq \gamma(G)\gamma(H)$, then the conjecture holds. This result quickly implies Vizing's conjecture for almost all pairs of graphs $G,H$ with $\left|G\right|\geq \left|H\right|$, satisfying $\left|G\right|\leq q^{\frac{\left|H\right|}{\log_q\left|H\right|}}$ for $q=\frac{1}{1-p}$ and $p$ the edge probability of the Erd\H{o}s-R\'enyi random graph.


title: How to collect matches that will catch fire

abstract: How can we select a cohort of promising astrophysicists before they have made their discoveries? This is the fundamental challenge of academic planning. I argue that science can only blossom if young researchers are rewarded for acquired skills and growth rather than inherited academic ancestry.


title: Day's fixed point theorem, Group cohomology and Quasi-isometric rigidity

abstract: In this note we explain how Day's fixed point theorem can be used to conjugate certain groups of biLipschitz maps of a metric space into special subgroups like similarity groups. In particular, we use Day's theorem to establish Tukia-type theorems and to give new proofs of quasi-isometric rigidity results.


title: Spectral Asymptotics for $V$-variable Sierpinski Gaskets

abstract: The family of $V$-variable fractals provides a means of interpolating between two families of random fractals previously considered in the literature; scale irregular fractals ($V=1$) and random recursive fractals ($V=\infty$). We consider a class of $V$-variable affine nested fractals based on the Sierpinski gasket with a general class of measures. We calculate the spectral exponent for a general measure and find the spectral dimension for these fractals. We show that the spectral properties and on-diagonal heat kernel estimates for $V$-variable fractals are closer to those of scale irregular fractals, in that it is the fluctuations in scale that determine their behaviour but that there are also effects of the spatial variability.


title: Deep Boosting: Layered Feature Mining for General Image Classification

abstract: Constructing effective representations is a critical but challenging problem in multimedia understanding. The traditional handcraft features often rely on domain knowledge, limiting the performances of exiting methods. This paper discusses a novel computational architecture for general image feature mining, which assembles the primitive filters (i.e. Gabor wavelets) into compositional features in a layer-wise manner. In each layer, we produce a number of base classifiers (i.e. regression stumps) associated with the generated features, and discover informative compositions by using the boosting algorithm. The output compositional features of each layer are treated as the base components to build up the next layer. Our framework is able to generate expressive image representations while inducing very discriminate functions for image classification. The experiments are conducted on several public datasets, and we demonstrate superior performances over state-of-the-art approaches.


title: Computing Supply Function Equilibria via Spline Approximations

abstract: The supply function equilibrium (SFE) is a model for competition in markets where each firm offers a schedule of prices and quantities to face demand uncertainty, and has been successfully applied to wholesale electricity markets. However, characterizing the SFE is difficult, both analytically and numerically. In this paper, we first present a specialized algorithm for capacity constrained asymmetric duopoly markets with affine costs. We show that solving the first order conditions (a system of differential equations) using spline approximations is equivalent to solving a least squares problem, which makes the algorithm highly efficient. We also propose using splines as a way to improve a recently introduced general algorithm, so that the equilibrium can be found more easily and faster with less user intervention. We show asymptotic convergence of the approximations to the true equilibria for both algorithms, and illustrate their performance with numerical examples.


title: Exploiting the Preferred Domain of FDD Massive MIMO Systems with Uniform Planar Arrays

abstract: Massive multiple-input multiple-output (MIMO) systems hold the potential to be an enabling technology for 5G cellular. Uniform planar array (UPA) antenna structures are a focus of much commercial discussion because of their ability to enable a large number of antennas in a relatively small area. With UPA antenna structures, the base station can control the beam direction in both the horizontal and vertical domains simultaneously. However, channel conditions may dictate that one dimension requires higher channel state information (CSI) accuracy than the other. We propose the use of an additional one bit of feedback information sent from the user to the base station to indicate the preferred domain on top of the feedback overhead of CSI quantization in frequency division duplexing (FDD) massive MIMO systems. Combined with variable-rate CSI quantization schemes, the numerical studies show that the additional one bit of feedback can increase the quality of CSI significantly for UPA antenna structures.


title: IceCube potential for detecting Q-ball dark matter in gauge mediation

abstract: We study Q-ball dark matter in gauge-mediated supersymmetry breaking, and seek the possibility of detection in the IceCube experiment. We find that the Q balls would be the dark matter in the parameter region different from that for gravitino dark matter. In particular, the Q ball is a good dark matter candidate for low reheating temperature, which may be suitable for the Affleck-Dine baryogenesis and/or nonthermal leptogenesis. Dark matter Q balls are detectable by IceCube-like experiments in the future, which is the peculiar feature compared to the case of gravitino dark matter.


title: Optimal dynamic program for r-domination problems over tree decompositions

abstract: There has been recent progress in showing that the exponential dependence on treewidth in dynamic programming algorithms for solving NP-hard problems are optimal under the Strong Exponential Time Hypothesis (SETH). We extend this work to $r$-domination problems. In $r$-dominating set, one wished to find a minimum subset $S$ of vertices such that every vertex of $G$ is within $r$ hops of some vertex in $S$. In connected $r$-dominating set, one additionally requires that the set induces a connected subgraph of $G$. We give a $O((2r+1)^{\mathrm{tw}} n)$ time algorithm for $r$-dominating set and a $O((2r+2)^{\mathrm{tw}} n^{O(1)})$ time algorithm for connected $r$-dominating set in $n$-vertex graphs of treewidth $\mathrm{tw}$. We show that the running time dependence on $r$ and $\mathrm{tw}$ is the best possible under SETH. This adds to earlier observations that a "+1" in the denominator is required for connectivity constraints.


title: Beyond Pixels: A Comprehensive Survey from Bottom-up to Semantic Image Segmentation and Cosegmentation

abstract: Image segmentation refers to the process to divide an image into nonoverlapping meaningful regions according to human perception, which has become a classic topic since the early ages of computer vision. A lot of research has been conducted and has resulted in many applications. However, while many segmentation algorithms exist, yet there are only a few sparse and outdated summarizations available, an overview of the recent achievements and issues is lacking. We aim to provide a comprehensive review of the recent progress in this field. Covering 180 publications, we give an overview of broad areas of segmentation topics including not only the classic bottom-up approaches, but also the recent development in superpixel, interactive methods, object proposals, semantic image parsing and image cosegmentation. In addition, we also review the existing influential datasets and evaluation metrics. Finally, we suggest some design flavors and research directions for future research in image segmentation.


title: Product Reservoir Computing: Time-Series Computation with Multiplicative Neurons

abstract: Echo state networks (ESN), a type of reservoir computing (RC) architecture, are efficient and accurate artificial neural systems for time series processing and learning. An ESN consists of a core of recurrent neural networks, called a reservoir, with a small number of tunable parameters to generate a high-dimensional representation of an input, and a readout layer which is easily trained using regression to produce a desired output from the reservoir states. Certain computational tasks involve real-time calculation of high-order time correlations, which requires nonlinear transformation either in the reservoir or the readout layer. Traditional ESN employs a reservoir with sigmoid or tanh function neurons. In contrast, some types of biological neurons obey response curves that can be described as a product unit rather than a sum and threshold. Inspired by this class of neurons, we introduce a RC architecture with a reservoir of product nodes for time series computation. We find that the product RC shows many properties of standard ESN such as short-term memory and nonlinear capacity. On standard benchmarks for chaotic prediction tasks, the product RC maintains the performance of a standard nonlinear ESN while being more amenable to mathematical analysis. Our study provides evidence that such networks are powerful in highly nonlinear tasks owing to high-order statistics generated by the recurrent product node reservoir.


title: A Second Order Approximation for the Caputo Fractional Derivative

abstract: When $0<\alpha<1$, the approximation for the Caputo derivative $$y^{(\alpha)}(x) = \frac{1}{\Gamma(2-\alpha)h^\alpha}\sum_{k=0}^n \sigma_k^{(\alpha)} y(x-kh)+O\bigl(h^{2-\alpha}\bigr),$$ where $\sigma_0^{(\alpha)} = 1, \sigma_n^{(\alpha)} = (n-1)^{1-a}-n^{1-a}$ and $$\sigma_k^{(\alpha)} = (k-1)^{1-\alpha}-2k^{1-a}+(k+1)^{1-\alpha},\quad (k=1...,n-1),$$ has accuracy $O\bigl(h^{2-\alpha}\bigr)$. We use the expansion of $\sum_{k=0}^n k^\alpha$ to determine an approximation for the fractional integral of order $2-\alpha$ and the second order approximation for the Caputo derivative $$y^{(\alpha)}(x) = \frac{1}{\Gamma(2-\alpha)h^\alpha}\sum_{k=0}^n \delta_k^{(\alpha)} y(x-kh)+O\bigl(h^{2}\bigr),$$ where $\delta_k^{(\alpha)} = \sigma_k^{(\alpha)}$ for $2\leq k\leq n$, $$\delta_0^{(\alpha)} = \sigma_0^{(\alpha)}-\zeta(\alpha-1), \delta_1^{(\alpha)} = \sigma_1^{(\alpha)}+2\zeta(\alpha-1),\delta_2^{(\alpha)} = \sigma_2^{(\alpha)}-\zeta(\alpha-1),$$ and $\zeta(s)$ is the Riemann zeta function. The numerical solutions of the fractional relaxation and subdiffusion equations are computed.