title: Global Spatio-temporal Patterns of Influenza in the Post-pandemic Era

abstract: We study the global spatio-temporal patterns of influenza dynamics. This is achieved by analysing and modelling weekly laboratory confirmed cases of influenza A and B from 138 countries between January 2006 and May 2014. The data were obtained from FluNet, the surveillance network compiled by the the World Health Organization. We report a pattern of {\it skip-and-resurgence} behavior between the years 2011 and 2013 for influenza H1N1/09, the strain responsible for the 2009 pandemic, in Europe and Eastern Asia. In particular, the expected H1N1/09 epidemic outbreak in 2011 failed to occur (or"skipped") in many countries across the globe, although an outbreak occurred in the following year. We also report a pattern of {\it well-synchronized} 2010 winter wave of H1N1/09 in the Northern Hemisphere countries, and a pattern of replacement of strain H1N1/77 by H1N1/09 between the 2009 and 2012 influenza seasons. Using both a statistical and a mechanistic mathematical model, and through fitting the data of 108 countries (108 countries in a statistical model and 10 large populations with a mechanistic model), we discuss the mechanisms that are likely to generate these events taking into account the role of multi-strain dynamics. A basic understanding of these patterns has important public health implications and scientific significance.


title: Supersymmetric field equations from momentum space

abstract: It is known that every irreducible unitary representation of positive energy of the Poincar\'e group can be realized as a subspace of tensor fields on Minkowski spacetime subjected to suitable partial differential equations. We first describe geometrically the general mechanism that produces, via Fourier transform, the invariant differential operators corresponding to those representations. Then, using a super-version of the Fourier transform, we show explicitly how a massive irreducible unitary representation of the super Poincar\'e group in dimension $(4|4)$ can be realized as a linear sub-supermanifold of suitably constrained superfunctions. In this way, we obtain supersymmetric equations in terms of ordinary (non-Grassmannian) fields. Finally, using the functor of points, we show how our equations can be related in a natural way to the Wess-Zumino equations for massive chiral superfields.


title: Nonzero Classical Discord

abstract: Quantum discord is the quantitative difference between two alternative expressions for bipartite mutual information, given respectively in terms of two distinct definitions for the conditional entropy. By constructing a stochastic model of shared states, classical discord can be similarly defined, quantifying the presence of some stochasticity in the measurement process. Therefore, discord can generally be understood as a quantification of the system's state disturbance due to local measurements, be it quantum or classical. We establish an operational meaning of classical discord in the context of state merging with noisy measurement and thereby show the quantum-classical separation in terms of a negative conditional entropy.


title: Traceable precision pA direct current measurements with the ULCA

abstract: A standard method for picoammeter calibrations is the capacitor charging technique, which allows generating traceable currents in the sub-nA range. However, its accuracy is limited by the ac-dc differences of the capacitances involved. The Ultrastable Low-noise Current Amplifier (ULCA) is a novel high-precision amperemeter for direct current measurements in the pA range, developed at PTB. Its amplifier stages, based on resistor networks and op-amps, can be calibrated traceably with a cryogenic current comparator (CCC) system. We compare the results from both independent calibration routes for two different ULCA prototypes. We find agreement between both methods at an uncertainty level below 10 microA/A, limited by the uncertainty of the currents generated with the capacitor charging method. The investigations confirm the superior performance of the new ULCA picoammeter.


title: Inconsistent treatment estimates from mis-specified logistic regression analyses of randomized trials

abstract: When the difference between treatments in a clinical trial is estimated by a difference in means, then it is well known that randomization ensures unbiassed estimation, even if no account is taken of important baseline covariates. However, when the treatment effect is assessed by other summaries, e.g. by an odds ratio if the outcome is binary, then bias can arise if some covariates are omitted, regardless of the use of randomization for treatment allocation or the size of the trial. We present accurate closed-form approximations for this asymptotic bias when important Normally distributed covariates are omitted from a logistic regression. We compare this approximation with ones in the literature and derive more convenient forms for some of these existing results. The expressions give insight into the form of the bias, which simulations show is usable for distributions other than the Normal. The key result applies even when there are additional binary covariates in the model.


title: Locally conformal symplectic nilmanifolds with no locally conformal K\"ahler metrics

abstract: We obtain an example of a compact locally conformal symplectic nilmanifold which admits no locally conformal K\"ahler metrics. This gives a new positive answer to a question raised by L. Ornea and M. Verbitsky.


title: A variational problem for curves on Riemann-Finsler surfaces

abstract: We study the variational problem for $N$-parallel curves on a Finslerian surface by means of Exterior Differential Systems using Griffiths' method. We obtain the conditions when these curves are extremals of a length functional and write the explicit form of Euler-Lagrange equations for this type of variational problem.


title: Charmed Pseudoscalar and Vector Mesons: a Comprehensive QCD Sum-Rule View of Their Decay Constants

abstract: In spite of undeniable similarities of the applied techniques, somewhat different challenges are encountered when extracting, from QCD sum rules derived from two-point correlators of appropriate interpolating heavy-light quark currents, the decay constants of charmed mesons of pseudoscalar nature, on the one hand, or of vector nature, on the other hand. This observation justifies a rather careful reassessment of the corresponding results.


title: Prime Coset Sum: A Systematic Method for Designing Multi-D Wavelet Filter Banks with Fast Algorithms

abstract: As constructing multi-D wavelets remains a challenging problem, we propose a new method called prime coset sum to construct multi-D wavelets. Our method provides a systematic way to construct multi-D non-separable wavelet filter banks from two 1-D lowpass filters, with one of whom being interpolatory. Our method has many important features including the following: 1) it works for any spatial dimension, and any prime scalar dilation, 2) the vanishing moments of the multi-D wavelet filter banks are guaranteed by certain properties of the initial 1-D lowpass filters, and furthermore, 3) the resulting multi-D wavelet filter banks are associated with fast algorithms that are faster than the existing fast tensor product algorithms.


title: Raking the Cocktail Party

abstract: We present the concept of an acoustic rake receiver---a microphone beamformer that uses echoes to improve the noise and interference suppression. The rake idea is well-known in wireless communications; it involves constructively combining different multipath components that arrive at the receiver antennas. Unlike spread-spectrum signals used in wireless communications, speech signals are not orthogonal to their shifts. Therefore, we focus on the spatial structure, rather than temporal. Instead of explicitly estimating the channel, we create correspondences between early echoes in time and image sources in space. These multiple sources of the desired and the interfering signal offer additional spatial diversity that we can exploit in the beamformer design. We present several "intuitive" and optimal formulations of acoustic rake receivers, and show theoretically and numerically that the rake formulation of the maximum signal-to-interference-and-noise beamformer offers significant performance boosts in terms of noise and interference suppression. Beyond signal-to-noise ratio, we observe gains in terms of the \emph{perceptual evaluation of speech quality} (PESQ) metric for the speech quality. We accompany the paper by the complete simulation and processing chain written in Python. The code and the sound samples are available online at \url{http://lcav.github.io/AcousticRakeReceiver/}.


title: A Factor-Adjusted Multiple Testing Procedure with Application to Mutual Fund Selection

abstract: In this article, we propose a factor-adjusted multiple testing (FAT) procedure based on factor-adjusted p-values in a linear factor model involving some observable and unobservable factors, for the purpose of selecting skilled funds in empirical finance. The factor-adjusted p-values were obtained after extracting the latent common factors by the principal component method. Under some mild conditions, the false discovery proportion can be consistently estimated even if the idiosyncratic errors are allowed to be weakly correlated across units. Furthermore, by appropriately setting a sequence of threshold values approaching zero, the proposed FAT procedure enjoys model selection consistency. Extensive simulation studies and a real data analysis for selecting skilled funds in the U.S. financial market are presented to illustrate the practical utility of the proposed method. Supplementary materials for this article are available online.


title: A DEIM Induced CUR Factorization

abstract: We derive a CUR matrix factorization based on the Discrete Empirical Interpolation Method (DEIM). For a given matrix $A$, such a factorization provides a low rank approximate decomposition of the form $A \approx C U R$, where $C$ and $R$ are subsets of the columns and rows of $A$, and $U$ is constructed to make $CUR$ a good approximation. Given a low-rank singular value decomposition $A \approx V S W^T$, the DEIM procedure uses $V$ and $W$ to select the columns and rows of $A$ that form $C$ and $R$. Through an error analysis applicable to a general class of CUR factorizations, we show that the accuracy tracks the optimal approximation error within a factor that depends on the conditioning of submatrices of $V$ and $W$. For large-scale problems, $V$ and $W$ can be approximated using an incremental QR algorithm that makes one pass through $A$. Numerical examples illustrate the favorable performance of the DEIM-CUR method, compared to CUR approximations based on leverage scores.


title: On the instability problem of a 3-D transonic oblique shock wave

abstract: In this paper, we are concerned with the instability problem of a 3-D transonic oblique shock wave for the steady supersonic flow past an infinitely long sharp wedge. The flow is assumed to be isentropic and irrotational. It was indicated in pages 317 of [9] that if a steady supersonic flow comes from minus infinity and hits a sharp symmetric wedge, then it follows from the Rankine-Hugoniot conditions and the physical entropy condition that there possibly appear a weak shock or a strong shock attached at the edge of the sharp wedge, which corresponds to a supersonic shock or a transonic shock, respectively. The question arises which of the two actually occurs. It has frequently been stated that the strong one is unstable and that, therefore, only the weak one could occur. However, a convincing proof of this instability has apparently never been given. The aim of this paper is to understand such a longstanding open question. We will show that the attached 3-D transonic oblique shock problem is overdetermined, which implies that the 3-D transonic shock is unstable in general.


title: Hardy inequalities for p-Laplacians with Robin boundary conditions

abstract: In this paper we study the best constant in a Hardy inequality for the p-Laplace operator on convex domains with Robin boundary conditions. We show, in particular, that the best constant equals $((p-1)/p)^p$ whenever Dirichlet boundary conditions are imposed on a subset of the boundary of non-zero measure. We also discuss some generalizations to non-convex domains.


title: A mathematical model for measurements in Quantum Mechanics

abstract: Let $V=\mathbb{C}^N$, and $H$ (an observable) a Hermitian linear operator on $V$. Let $v_1,..., v_n$ be an orthonormal basis for $V$. Let $\mathcal{M}$ be a measurement apparatus prepared to measure a state of an observed system and collapses the state to one of the $v_j$'s. Here we propose a simple model which explains the Born rule and is compatible with entanglement.


title: Continuous and Discontinuous Galerkin Time Stepping Methods for Nonlinear Initial Value Problems with Application to Finite Time Blow-Up

abstract: We consider continuous and discontinuous Galerkin time stepping methods of arbitrary order as applied to nonlinear initial value problems in real Hilbert spaces. Our only assumption is that the nonlinearities are continuous; in particular, we include the case of unbounded nonlinear operators. Specifically, we develop new techniques to prove general Peano-type existence results for discrete solutions. In particular, our results show that the existence of solutions is independent of the local approximation order, and only requires the local time steps to be sufficiently small (independent of the polynomial degree). The uniqueness of (local) solutions is addressed as well. In addition, our theory is applied to finite time blow-up problems with nonlinearities of algebraic growth. For such problems we develop a time step selection algorithm for the purpose of numerically computing the blow-up time, and provide a convergence result.


title: Spreading of intolerance under economic stress: results from a model with reputation

abstract: When a population is engaged in successive prisoner's dilemmas, indirect reciprocity through reputation fosters cooperation through the emergence of moral and action rules. A simplified model has recently been proposed where individuals choose between helping or not others, and are judged good or bad for it by the rest of the population. The reputation so acquired will condition future actions. In this model, eight strategies (referred to as 'leading eight') enforce a high level of cooperation, generate high payoffs and are therefore resistant to invasions by other strategies. Here we show that, by assigning each individual one out of two labels that peers can distinguish (e.g., political ideas, religion, skin colour...) and allowing moral and action rules to depend on the label, intolerant behaviours can emerge within minorities under sufficient economic stress. We analyse the sets of conditions where this can happen and also discuss the circumstances under which tolerance can be restored. Our results agree with empirical observations that correlate intolerance and economic stress, and predict a correlation between the degree of tolerance of a population and its composition and ethical stance.


title: Unidirectional magnetoelectric-field multiresonant tunneling

abstract: Unidirectional multi-resonant tunneling of the magneto-electric (ME)-field excitations through a sub wavelength (regarding the scales of regular electromagnetic radiation) vacuum or isotropic-dielectric regions has been observed in two-port microwave structures having a quasi-2D ferrite disk with magnetic-dipolar-mode (MDM) oscillations. The excitations manifest themselves as the Fano-resonance peaks in the scattering-matrix parameters at the stationary states of the MDM spectrum. The ME near-field excitations are quasi magnetostatic fields with non-zero helicity parameter. Topological-phase properties of ME fields are determined by edge chiral currents of MDM oscillations. We show that while for a given direction of a bias magnetic field (in other words, for a given direction of time), the ME-field excitations are considered as forward tunneling processes, in the opposite direction of a bias magnetic field (the opposite direction of time), there are backward tunneling processes. Unidirectional ME-field resonant tunneling is observed due to distinguishable topology of the forward and backward ME-field excitations. We establish a close connection between the Fano-resonance unidirectional tunneling and topology of ME fields in different microwave structures.


title: Water Dynamics at Rough Interfaces

abstract: We use molecular dynamics computer simulations and nuclear magnetic resonance experiments to investigate the dynamics of water at interfaces of molecular roughness and low mobility. We find that, when approaching such interfaces, the structural relaxation of water, i.e., the $\alpha$ process, slows down even when specific attractive interactions are absent. This prominent effect is accompanied by a smooth transition from Vogel to Arrhenius temperature dependence and by a growing importance of jump events. Consistently, at protein surfaces, deviations from Arrhenius behavior are weak when free water does not exist. Furthermore, in nanoporous silica, a dynamic crossover of liquid water occurs when a fraction of solid water forms near 225 K and, hence, the liquid dynamics changes from bulk-like to interface-dominated. At sufficiently low temperatures, water exhibits a quasi-universal $\beta$ process, which is characterized by an activation energy of $E_a\!=\!0.5$ eV and involves anisotropic reorientation about large angles. As a consequence of its large amplitude, the faster $\beta$ process destroys essentially all orientational correlation, rendering observation of a possible slower $\alpha$ process difficult in standard experiments. Nevertheless, we find indications for the existence of structural relaxation down to a glass transition of interfacial water near 185 K. Hydrated proteins show a highly restricted backbone motion with an amplitude, which decreases upon cooling and vanishes at comparable temperatures, providing evidence for a high relevance of water rearrangements in the hydration shell for secondary protein relaxations.


title: Process-Oriented Parallel Programming with an Application to Data-Intensive Computing

abstract: We introduce process-oriented programming as a natural extension of object-oriented programming for parallel computing. It is based on the observation that every class of an object-oriented language can be instantiated as a process, accessible via a remote pointer. The introduction of process pointers requires no syntax extension, identifies processes with programming objects, and enables processes to exchange information simply by executing remote methods. Process-oriented programming is a high-level language alternative to multithreading, MPI and many other languages, environments and tools currently used for parallel computations. It implements natural object-based parallelism using only minimal syntax extension of existing languages, such as C++ and Python, and has therefore the potential to lead to widespread adoption of parallel programming. We implemented a prototype system for running processes using C++ with MPI and used it to compute a large three-dimensional Fourier transform on a computer cluster built of commodity hardware components. Three-dimensional Fourier transform is a prototype of a data-intensive application with a complex data-access pattern. The process-oriented code is only a few hundred lines long, and attains very high data throughput by achieving massive parallelism and maximizing hardware utilization.


title: Hypothesis Testing For Network Data in Functional Neuroimaging

abstract: In recent years, it has become common practice in neuroscience to use networks to summarize relational information in a set of measurements, typically assumed to be reflective of either functional or structural relationships between regions of interest in the brain. One of the most basic tasks of interest in the analysis of such data is the testing of hypotheses, in answer to questions such as "Is there a difference between the networks of these two groups of subjects?" In the classical setting, where the unit of interest is a scalar or a vector, such questions are answered through the use of familiar two-sample testing strategies. Networks, however, are not Euclidean objects, and hence classical methods do not directly apply. We address this challenge by drawing on concepts and techniques from geometry, and high-dimensional statistical inference. Our work is based on a precise geometric characterization of the space of graph Laplacian matrices and a nonparametric notion of averaging due to Fr\'echet. We motivate and illustrate our resulting methodologies for testing in the context of networks derived from functional neuroimaging data on human subjects from the 1000 Functional Connectomes Project. In particular, we show that this global test is more statistical powerful, than a mass-univariate approach. In addition, we have also provided a method for visualizing the individual contribution of each edge to the overall test statistic.


title: Odd q-State Clock Spin-Glass Models in Three Dimensions, Asymmetric Phase Diagrams, and Multiple Algebraically Ordered Phases

abstract: Distinctive orderings and phase diagram structures are found, from renormalization-group theory, for odd q-state clock spin-glass models in d=3 dimensions. These models exhibit asymmetric phase diagrams, as is also the case for quantum Heisenberg spin-glass models. No finite-temperature spin-glass phase occurs. For all odd $q\geqslant 5$, algebraically ordered antiferromagnetic phases occur. One such phase is dominant and occurs for all $q\geqslant 5$. Other such phases occupy small low-temperature portions of the phase diagrams and occur for $5 \leqslant q \leqslant 15$. All algebraically ordered phases have the same structure, determined by an attractive finite-temperature sink fixed point where a dominant and a subdominant pair states have the only non-zero Boltzmann weights. The phase transition critical exponents quickly saturate to the high q value.


title: Quality of Experience (QoE) beyond Quality of Service (QoS) as its baseline: QoE at the Interface of Experience Domains

abstract: In this work, a new approach to the definition of the quality of experience is presented. By considering the quality of service as a baseline, that portion of the QoE that can be inferred from the QoS is excluded, and then the rest of the QoE is approached with the notion of QoE at a Boundary (QoEaaB). With the QoEaaB as the core of the proposed approach, various potential boundaries, and their associated unseen opportunities to improve the QoE are discussed. In particular, property, contract, SLA, and content are explored in terms of their boundaries and also their associated QoEaaB. With an interest in online video delivery, management of resource sharing and isolation associated with multi-tenant operations is considered. It is concluded that the proposed QoEaaB can bring a new perspective in QoE modeling and assessment toward a more enriched approach to improving the experience based on innovation and deep connectivity among actors.


title: Arbitrage-free prediction of the implied volatility smile

abstract: This paper gives an arbitrage-free prediction for future prices of an arbitrary co-terminal set of options with a given maturity, based on the observed time series of these option prices. The statistical analysis of such a multi-dimensional time series of option prices corresponding to $n$ strikes (with $n$ large, e.g. $n\geq 40$) and the same maturity, is a difficult task due to the fact that option prices at any moment in time satisfy non-linear and non-explicit no-arbitrage restrictions. Hence any $n$-dimensional time series model also has to satisfy these implicit restrictions at each time step, a condition that is impossible to meet since the model innovations can take arbitrary values. We solve this problem for any $n\in\NN$ in the context of Foreign Exchange (FX) by first encoding the option prices at each time step in terms of the parameters of the corresponding risk-neutral measure and then performing the time series analysis in the parameter space. The option price predictions are obtained from the predicted risk-neutral measure by effectively integrating it against the corresponding option payoffs. The non-linear transformation between option prices and the risk-neutral parameters applied here is \textit{not} arbitrary: it is the standard mapping used by market makers in the FX option markets (the SABR parameterisation) and is given explicitly in closed form. Our method is not restricted to the FX asset class nor does it depend on the type of parameterisation used. Statistical analysis of FX market data illustrates that our arbitrage-free predictions outperform the naive random walk forecasts, suggesting a potential for building management strategies for portfolios of derivative products, akin to the ones widely used in the underlying equity and futures markets.


title: Route to chaos in optomechanics

abstract: We establish the emergence of chaotic motion in optomechanical systems. Chaos appears at negative detuning for experimentally accessible values of the pump power and other system parameters. We describe the sequence of period doubling bifurcations that leads to chaos, and state the experimentally observable signatures in the optical spectrum. In addition to the semi-classical dynamics we analyze the possibility of chaotic motion in the quantum regime. We find that quantum mechanics protects the optomechanical system against irregular dynamics, such that simple periodic orbits reappear and replace the classically chaotic motion. In this way observation of the dynamical signatures makes it possible to pin down the crossover from quantum to classical mechanics.


title: Event-by-event mean $p_{\rm T}$ fluctuations in pp and Pb-Pb collisions at the LHC

abstract: Event-by-event fluctuations of the mean transverse momentum of charged particles produced in pp collisions at $\sqrt{s}$ = 0.9, 2.76 and 7 TeV, and Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV are studied as a function of the charged-particle multiplicity using the ALICE detector at the LHC. Dynamical fluctuations indicative of correlated particle emission are observed in all systems. The results in pp collisions show little dependence on collision energy. The Monte Carlo event generators PYTHIA and PHOJET are in qualitative agreement with the data. Peripheral Pb-Pb data exhibit a similar multiplicity dependence as that observed in pp. In central Pb-Pb, the results deviate from this trend, featuring a significant reduction of the fluctuation strength. The results in Pb--Pb are in qualitative agreement with previous measurements in Au-Au at lower collision energies and with expectations from models that incorporate collective phenomena.


title: Modeling the dynamics of single-bubble sonoluminescence

abstract: Sonoluminescence (SL) is the phenomenon in which acoustic energy is (partially) transformed into light. It may occur by means of many or just one bubble of gas inside a liquid medium, giving rise to the terms multi-bubble- and single-bubble sonoluminescence (MBSL and SBSL). In the last years some models have been proposed to explain this phenomenon, but there is still no complete theory for the light emission mechanism (especially in the case of SBSL). In this work, we will not address this more complicated particular issue, but only present a simple model describing the dynamical behaviour of the sonoluminescent bubble, in the SBSL case. Using simple numerical techniques within the software Matlab, we discuss solutions considering various possibilities for some of the parameters involved: liquid compressibility, superficial tension, viscosity, and type of gas. The model may be used as an introductory study of sonoluminescence in physics courses at undergraduate or graduate levels, as well as a quite clarifying example of a physical system exhibiting large nonlinearity.


title: Measurement of the production cross-section of $\psi(2S)\to J/\psi(\to\mu^+\mu^-)\pi^+\pi^-$ in $pp$ collisions at $\sqrt{s}=7$ TeV at ATLAS

abstract: The prompt and non-prompt production cross-sections for $\psi(2S)$ mesons are measured using 2.1 fb$^{-1}$ of $pp$ collision data at a centre-of-mass energy of 7 TeV recorded by the ATLAS experiment at the LHC. The measurement exploits the $\psi(2S)\to J/\psi(\to\mu^+\mu^-)\pi^+\pi^-$ decay mode, and probes $\psi(2S)$ mesons with transverse momenta in the range $10\leq p_T<100$ GeV and rapidity $|y|<2.0$. The results are compared to other measurements of $\psi(2S)$ production at the LHC and to various theoretical models for prompt and non-prompt quarkonium production.


title: Resource Letter TTSM-1: Teaching thermodynamics and statistical mechanics in introductory physics, chemistry, and biology

abstract: This Resource Letter draws on discipline-based education research from physics, chemistry, and biology to collect literature on the teaching of thermodynamics and statistical mechanics in the three disciplines. While the overlap among the disciplinary literatures is limited at present, we hope this Resource Letter will spark more interdisciplinary interaction.


title: Spin-dependent recombination at arsenic donors in ion-implanted silicon

abstract: Spin-dependent transport processes in thin near-surface doping regions created by low energy ion implantation of arsenic in silicon are detected by two methods, spin-dependent recombination (SDR) using microwave photoconductivity and electrically detected magnetic resonance (EDMR) monitoring the DC current through the sample. The high sensitivity of these techniques allows the observation of the magnetic resonance in particular of As in weak magnetic fields and at low resonance frequencies (40-1200 MHz), where high-field-forbidden transitions between the magnetic substates can be observed due to the mixing of electron and nuclear spin states. Several implantation-induced defects are present in the samples studied and act as spin readout partner. We explicitly demonstrate this by electrically detected electron double resonance experiments and identify a pair recombination of close pairs formed by As donors and oxygen-vacancy centers in an excited triplet state (SL1) as the dominant spin-dependent process in As-implanted Czochralski-grown Si.


title: Macroscopic quantum effects for classical light

abstract: Macroscopic quantum optical effects (Schrodinger cat states, squeezing, collapse and revival) for light beams propagating in an inhomogeneous linear medium are demonstrated using exact analytical solutions of wave equation. It is shown that the coherent superposition of macroscopically distinguishable states is generated via mode interference from an initial off-axis single wave-packet. Squeezed cat states with a fidelity > 99% arise periodically and disappear rapidly within limited intervals of a propagation distance. Collapse and revival of wave packets at long-term non-paraxial evolution due to mode interference is demonstrated. Oscillations of the beam trajectory with extremely small amplitude of the order of 10-19 m which is typical of the estimated displacement caused by cosmic gravitational waves in gravity-wave detectors occur.


title: Multichannel Compressive Sensing MRI Using Noiselet Encoding

abstract: The incoherence between measurement and sparsifying transform matrices and the restricted isometry property (RIP) of measurement matrix are two of the key factors in determining the performance of compressive sensing (CS). In CS-MRI, the randomly under-sampled Fourier matrix is used as the measurement matrix and the wavelet transform is usually used as sparsifying transform matrix. However, the incoherence between the randomly under-sampled Fourier matrix and the wavelet matrix is not optimal, which can deteriorate the performance of CS-MRI. Using the mathematical result that noiselets are maximally incoherent with wavelets, this paper introduces the noiselet unitary bases as the measurement matrix to improve the incoherence and RIP in CS-MRI, and presents a method to design the pulse sequence for the noiselet encoding. This novel encoding scheme is combined with the multichannel compressive sensing (MCS) framework to take the advantage of multichannel data acquisition used in MRI scanners. An empirical RIP analysis is presented to compare the multichannel noiselet and multichannel Fourier measurement matrices in MCS. Simulations are presented in the MCS framework to compare the performance of noiselet encoding reconstructions and Fourier encoding reconstructions at different acceleration factors. The comparisons indicate that multichannel noiselet measurement matrix has better RIP than that of its Fourier counterpart, and that noiselet encoded MCS-MRI outperforms Fourier encoded MCS-MRI in preserving image resolution and can achieve higher acceleration factors. To demonstrate the feasibility of the proposed noiselet encoding scheme, two pulse sequences with tailored spatially selective RF excitation pulses was designed and implemented on a 3T scanner to acquire the data in the noiselet domain from a phantom and a human brain.


title: Tractable Model for Rate in Self-Backhauled Millimeter Wave Cellular Networks

abstract: Millimeter wave (mmW) cellular systems will require high gain directional antennas and dense base station (BS) deployments to overcome high near field path loss and poor diffraction. As a desirable side effect, high gain antennas provide interference isolation, providing an opportunity to incorporate self-backhauling--BSs backhauling among themselves in a mesh architecture without significant loss in throughput--to enable the requisite large BS densities. The use of directional antennas and resource sharing between access and backhaul links leads to coverage and rate trends that differ significantly from conventional microwave ($\mu$W) cellular systems. In this paper, we propose a general and tractable mmW cellular model capturing these key trends and characterize the associated rate distribution. The developed model and analysis is validated using actual building locations from dense urban settings and empirically-derived path loss models. The analysis shows that in sharp contrast to the interference limited nature of $\mu$W cellular networks, the spectral efficiency of mmW networks (besides total rate) also increases with BS density particularly at the cell edge. Increasing the system bandwidth, although boosting median and peak rates, does not significantly influence the cell edge rate. With self-backhauling, different combinations of the wired backhaul fraction (i.e. the faction of BSs with a wired connection) and BS density are shown to guarantee the same median rate (QoS).


title: Supernatural numbers and a new topology on the arithmetic site

abstract: In arXiv:1405.4527 Connes and Consani introduced and studied the arithmetic site and showed that the isomorphism classes of points are in canonical bijection with the finite adele classes $\mathbb{Q}^*_+ \backslash \mathbb{A}^f_{\mathbb{Q}} / \widehat{\mathbb{Z}}^*$. The induced topology of $\mathbb{A}^f_{\mathbb{Q}}$ on this set is trivial, whence this space is usually studied via noncommutative geometry. However, we can define another topology on this set of points, which shares several properties one might expect of the mythical object $\overline{\mathbf{Spec}(\mathbb{Z})}/\mathbb{F}_1$: it is compact, has an uncountable basis of opens, each non-empty open being dense, and it satisfies the $T_1$ separation property for incomparable points.


title: An eigenvalue inequality for Schr\"odinger operators with $\delta$ and $\delta'$-interactions supported on hypersurfaces

abstract: We consider self-adjoint Schr\"odinger operators in $L^2 (\mathbb{R}^d)$ with a $\delta$-interaction of strength $\alpha$ and a $\delta'$-interaction of strength $\beta$, respectively, supported on a hypersurface, where $\alpha$ and $\beta^{-1}$ are bounded, real-valued functions. It is known that the inequality $0 < \beta \leq 4/\alpha$ implies inequality of the eigenvalues of these two operators below the bottoms of the essential spectra. We show that this eigenvalue inequality is strict whenever $\beta < 4 / \alpha$ on a nonempty, open subset of the hypersurface. Moreover, we point out special geometries of the interaction support, such as broken lines or infinite cones, for which strict inequality of the eigenvalues even holds in the borderline case $\beta = 4 / \alpha$.


title: Electron Spin Dephasing and Optical Pumping of Nuclear Spins in GaN

abstract: We have measured the donor-bound electron spin dynamics in cubic GaN by time-resolved Kerr rotation experiments. The ensemble electron spin dephasing time in this quantum dot like system characterized by a Bohr radius of 2.5 nm is of the order of 1.5 ns as a result of the interaction with the fluctuating nuclear spins. It increases drastically when an external magnetic field as small as 10 mT is applied. We extract a dispersion of the nuclear hyperfine field {\delta}Bn $\sim$ 4 mT, in agreement with calculations. We also demonstrate for the first time in GaN based systems the optical pumping of nuclear spin yielding the build-up of a significant nuclear polarization.


title: Canonical decomposition of linear differential operators with selected differential Galois groups

abstract: We revisit an order-six linear differential operator having a solution which is a diagonal of a rational function of three variables. Its exterior square has a rational solution, indicating that it has a selected differential Galois group, and is actually homomorphic to its adjoint. We obtain the two corresponding intertwiners giving this homomorphism to the adjoint. We show that these intertwiners are also homomorphic to their adjoint and have a simple decomposition, already underlined in a previous paper, in terms of order-two self-adjoint operators. From these results, we deduce a new form of decomposition of operators for this selected order-six linear differential operator in terms of three order-two self-adjoint operators. We then generalize the previous decomposition to decompositions in terms of an arbitrary number of self-adjoint operators of the same parity order. This yields an infinite family of linear differential operators homomorphic to their adjoint, and, thus, with a selected differential Galois group. We show that the equivalence of such operators is compatible with these canonical decompositions. The rational solutions of the symmetric, or exterior, squares of these selected operators are, noticeably, seen to depend only on the rightmost self-adjoint operator in the decomposition. These results, and tools, are applied on operators of large orders. For instance, it is seen that a large set of (quite massive) operators, associated with reflexive 4-polytopes defining Calabi-Yau 3-folds, obtained recently by P. Lairez, correspond to a particular form of the decomposition detailed in this paper.


title: Cyclic homogeneous Riemannian manifolds

abstract: In spin geometry, traceless cyclic homogeneous Riemannian manifolds equipped with a homogeneous spin structure can be viewed as the simplest manifolds after Riemannian symmetric spin spaces. In this paper, we give some characterizations and properties of cyclic and traceless cyclic homogeneous Riemannian manifolds and we obtain the classification of simply-connected cyclic homogeneous Riemannian manifolds of dimension less than or equal to four. We also present a wide list of examples of non-compact irreducible Riemannian $3$-symmetric spaces admitting cyclic metrics and give the expression of these metrics.


title: Spin pumping efficiency in room-temperature CdSe nanocrystal quantum dots

abstract: To understand and optimize optical spin initialization in room temperature CdSe nanocrystal quantum dots (NCQDs) we studied the dependence of the time-resolved Faraday rotation signal on pump energy $E_p$ in a series of NCQD samples with different sizes. In larger NCQDs, we observe two peaks in the spin signal vs. $E_p$, whereas in smaller NQCDs, only a single peak is observed before the signal falls to a low, broad plateau at higher energies. We calculate the spin-dependent oscillator strengths of optical transitions using a simple effective mass model to understand these results. The observed $E_p$ dependence of the spin pumping efficiency (SPE) arises from the competition between the heavy hole (hh), light hole (lh) and split-off (so) band contributions to transitions to the conduction band. The two latter contributions lead to an electron spin polarization in the opposite direction from the former. At lower $E_p$ the transitions are dominated by the hh band, giving rise to the low energy peaks. At higher $E_p$, the increasing contributions from the lh and so bands lead to a reduction in SPE. The different number of peaks in larger and smaller NCQDs is attributed to size-dependence of the ordering of the valence band states.


title: Countable families of solutions of a limit stationary semilinear fourth-order Cahn--Hilliard equation I. Mountain pass and Lusternik--Schnirel'man patterns in R^N

abstract: Solutions of the stationary semilinear Cahn--Hilliard equation -\Delta^2 u - u -\Delta(|u|^{p-1}u)=0 in R^N, with p>1, which are exponentially decaying at infinity, are studied. Using the Mounting Pass Lemma allows us the determination of two different solutions. On the other hand, the application of Lusternik--Schnirel'man (L--S) Category Theory shows the existence of, at least, a countable family of solutions. However, through numerical methods it is shown that the whole set of solutions, even in 1D, is much wider. This suggests that, actually, there exists, at least, a countable set of countable families of solutions, in which only the first one can be obtained by the L--S min-max approach.


title: Spectromicroscopy of electronic phase separation in K$_x$Fe$_{2-y}$Se$_2$ superconductor

abstract: Structural phase separation in A$_x$Fe$_{2-y}$Se$_2$ system has been studied by different experimental techniques, however, it should be important to know how the electronic uniformity is influenced, on which length scale the electronic phases coexist, and what is their spatial distribution. Here, we have used novel scanning photoelectron microscopy (SPEM) to study the electronic phase separation in K$_x$Fe$_{2-y}$Se$_2$, providing a direct measurement of the topological spatial distribution of the different electronic phases. The SPEM results reveal a peculiar interconnected conducting filamentary phase that is embedded in the insulating texture. The filamentary structure with a particular topological geometry could be important for the high T$_c$ superconductivity in the presence of a phase with a large magnetic moment in A$_x$Fe$_{2-y}$Se$_2$ materials.


title: Analysis of orbit accumulation points and the Greene-Krantz conjecture

abstract: In $\mathbb{C}^2$, we classify the domains for which $\rm Aut(\Omega)$ is noncompact and describe these domains by their defining functions. This note is based on the technique of the scaling method introduced by Frankel \cite{Fr86} and Kim \cite{Ki90}. One feature of this article is that we are able to analyze the defining functions of infinite type boundary. As a corollary, we also prove a result that under some conditions, $\rm Aut(\Omega)$ contains $\mathbb{R}$, which is an extension of \cite{Fr86}.


title: Reading the Source Code of Social Ties

abstract: Though online social network research has exploded during the past years, not much thought has been given to the exploration of the nature of social links. Online interactions have been interpreted as indicative of one social process or another (e.g., status exchange or trust), often with little systematic justification regarding the relation between observed data and theoretical concept. Our research aims to breach this gap in computational social science by proposing an unsupervised, parameter-free method to discover, with high accuracy, the fundamental domains of interaction occurring in social networks. By applying this method on two online datasets different by scope and type of interaction (aNobii and Flickr) we observe the spontaneous emergence of three domains of interaction representing the exchange of status, knowledge and social support. By finding significant relations between the domains of interaction and classic social network analysis issues (e.g., tie strength, dyadic interaction over time) we show how the network of interactions induced by the extracted domains can be used as a starting point for more nuanced analysis of online social data that may one day incorporate the normative grammar of social interaction. Our methods finds applications in online social media services ranging from recommendation to visual link summarization.


title: Homogeneous and inhomogeneous sources of optical transition broadening in room temperature CdSe/ZnS nanocrystal quantum dots

abstract: We perform photoluminescence excitation measurements on individual CdSe/ZnS nanocrystal quantum dots (NCQDs) at room temperature to study optical transition energies and broadening. The observed features in the spectra are identified and compared to calculated transition energies using an effective mass model. The observed broadening is attributed to phonon broadening, spectral diffusion and size and shape inhomogeneity. The former two contribute the broadening transitions in individual QDs while the latter contributes to the QD-to-QD variation. We find that phonon broadening is often not the dominant contribution to transition line widths, even at room temperature, and that broadening does not necessarily increase with transition energy. This may be explained by differing magnitude of spectral diffusion for different quantum-confined states.


title: Comment on 'Counterintuitive Dispersion Violating Kramers-Kronig Relations in Gain Slabs'

abstract: In this Comment, we refute conclusions made in Phys. Rev. Lett. 112, 233601 (2014) by L.-G. Wang, L. Wang, M. Al-Amri, S.-Y. Zhu, and M. S. Zubairy. These conclusions stem from the use of the linear theory, which is not applicable to systems having poles of the reflection and transmission coefficients in the upper half of the complex frequency plane.


title: S-wave ground state charmonium decays of $B_c$ mesons in the perturbative QCD approach

abstract: We make a systematic investigation on the two-body nonleptonic decays $B_c\rightarrow J/\Psi(\eta_c),M$ by employing the perturbative QCD approach based on $k_T$ factorization, where $M$ is a light pseudoscalar or vector or a heavy charmed meson. We predict the branching ratios and direct CP asymmetries of these $B_c$ decays and also the transverse polarization fractions of $B_c\rightarrow J/\Psi V,J/\psi D^{*}_{(s)}$ decays. It is found that these decays have a large branching ratios of the order of $10^{-4} - 10^{-2}$ and could be measured by the future LHC-b experiment. Our predictions for the ratios of branching fractions $\frac{\mathcal {BR}(B_c^+\rightarrow J/\Psi D_s^+)} {\mathcal {BR}(B_c^+\rightarrow J/\Psi \pi^+)}$ ,$\frac{\mathcal {BR}(B_c^+\rightarrow J/\Psi D_s^{*+})} {\mathcal {BR}(B_c^+\rightarrow J/\Psi D_s^+)}$ and $\frac{\mathcal {BR}(B_c^+\rightarrow J/\Psi K^+)} {\mathcal {BR}(B_c^+\rightarrow J/\Psi \pi^+)}$ are in good agreement with the data. A large transverse polarization fraction which can reach $48\%$ is predicted in $B_c^+\rightarrow J/\Psi D_s^{*+}$ decay, which is consistent with the data. We find a possible direct CP violation in $B_c\rightarrow J/\psi D^{*}$ decays, which are helpful to test the CP violating effects in $B_c$ decays.


title: Equivalence between formulations in Cosmological Perturbation Theory: The primordial magnetic fields as an example

abstract: Nowdays, Cosmological Perturbation Theory is a standard and useful tool in theoretical cosmology. In this work, we compare the 1+3 covariant formalism in perturbation theory (Ellis et al.) to the gauge invariant approach (Bruni et al.), and we show the equivalence of these formalisms to fix the choice of the perturbed variables (gauge choice) in magnetogenesis. We analyze the evolution of primordial magnetic fields through perturbation theory and we discuss the similarities and differences between these two approaches. We get the Maxwell's equations and show a cosmic dynamo like equation written in Poisson gauge, computing the evolution of primordial magnetic fields. Finally, prospects around these formalisms in the study of magnetogenesis are discussed.


title: Inequalities involving the generating function for the number of partitions into odd parts

abstract: Fibonacci numbers can be expressed in terms of multinomial coefficients as sums over integer partitions into odd parts. We use this fact to introduce a family of double inequalities involving the generating function for the number of partitions into odd parts and the generating function for the number of odd divisors.


title: Privacy-Preserving Filtering for Event Streams

abstract: Many large-scale information systems such as intelligent transportation systems, smart grids or smart buildings collect data about the activities of their users to optimize their operations. To encourage participation and adoption of these systems, it is becoming increasingly important that the design process take privacy issues into consideration. In a typical scenario, signals originate from many sensors capturing events involving the users, and several statistics of interest need to be continuously published in real-time. This paper considers the problem of providing differential privacy guarantees for such multi-input multi-output systems processing event streams. We show how to construct and optimize various extensions of the zero-forcing equalization mechanism, which we previously proposed for single-input single-output systems. Some of these extensions can take a model of the input signals into account. We illustrate our privacy-preserving filter design methodology through the problem of privately monitoring and forecasting occupancy in a building equipped with multiple motion detection sensors.


title: The Cauchy problem for a tenth-order thin film equation II. Oscillatory source-type and fundamental similarity solutions

abstract: Fundamental global similarity solutions of the standard form u_\g(x,t)=t^{-\a_\g} f_\g(y), with the rescaled variable y= x/{t^{\b_\g}}, \b_\g= \frac {1-n \a_\g}{10}, where \a_\g>0 are real nonlinear eigenvalues (\g is a multiindex in R^N) of the tenth-order thin film equation (TFE-10) u_{t} = \nabla \cdot(|u|^{n} \n \D^4 u) in R^N \times R_+, n>0, are studied. The present paper continues the study began by the authors in the previous paper P. Alvarez-Caudevilla, J.D.Evans, and V.A. Galaktionov, The Cauchy problem for a tenth-order thin film equation I. Bifurcation of self-similar oscillatory fundamental solutions, Mediterranean Journal of Mathematics, No. 4, Vol. 10 (2013), 1759-1790. Thus, the following questions are also under scrutiny: (I) Further study of the limit n \to 0, where the behaviour of finite interfaces and solutions as y \to infinity are described. In particular, for N=1, the interfaces are shown to diverge as follows: |x_0(t)| \sim 10 \left( \frac{1}{n}\sec\left( \frac{4\pi}{9} \right) \right)^{\frac 9{10}} t^{\frac 1{10}} \to \infty as n \to 0^+. (II) For a fixed n \in (0, \frac 98), oscillatory structures of solutions near interfaces. (III) Again, for a fixed n \in (0, \frac 98), global structures of some nonlinear eigenfunctions \{f_\g\}_{|\g| \ge 0} by a combination of numerical and analytical methods.


title: Global behavior of N competing species with strong diffusion: diffusion leads to exclusion

abstract: It is known that the competitive exclusion principle holds for a large kind of models involving several species competing for a single resource in an homogeneous environment. Various works indicate that the coexistence is possible in an heterogeneous environment. We propose a spatially heterogeneous system modeling the competition of several species for a single resource. If spatial movements are fast enough, we show that our system can be well approximated by a spatially homogeneous system, called aggregated model, which can be explicitly computed. Moreover, we show that if the competitive exclusion principle holds for the aggregated model, it holds for the spatially heterogeneous model too.


title: Variational approach for a class of cooperative systems

abstract: The aim of this work is to ascertain the characterization of the existence of coexistence states for a class of cooperative systems supported by the study of an associated non--local equation through classical variational methods. Thanks to those results we are able to obtain the blow--up behaviour of the solutions in the whole domain for certain values of the main continuation parameter.


title: The Cauchy problem for tenth-order thin film equation I. Bifurcation of oscillatory fundamental solutions

abstract: Fundamental global similarity solutions of the tenth-order thin film equation u_{t} = \nabla \cdot(|u|^{n} \n \D^4 u) in R^N \times R_+, where n>0 are studied. The main approach consists in passing to the limit n \to 0^+ by using Hermitian non-self-adjoint spectral theory corresponding to the rescaled linear poly-harmonic equation u_t= \D^5 u in R^N \times \re_+.


title: Second order QCD corrections to gluonic jet production at hadron colliders

abstract: We report on the calculation of the next-to-next-to-leading order (NNLO) QCD corrections to the production of two gluonic jets at hadron colliders. In previous work, we discussed gluonic dijet production in the gluon-gluon channel. Here, for the first time, we update our numerical results to include the leading colour contribution to the production of two gluonic jets via quark-antiquark scattering.


title: Some Liouville theorems for the fractional Laplacian

abstract: In this paper, we prove the following result. Let $\alpha$ be any real number between $0$ and $2$. Assume that $u$ is a solution of $$ \left\{\begin{array}{ll} (-\Delta)^{\alpha/2} u(x) = 0 , \;\; x \in \mathbb{R}^n ,\\ \displaystyle\underset{|x| \to \infty}{\underline{\lim}} \frac{u(x)}{|x|^{\gamma}} \geq 0 , \end{array} \right. $$ for some $0 \leq \gamma \leq 1$ and $\gamma < \alpha$. Then $u$ must be constant throughout $\mathbb{R}^n$. This is a Liouville Theorem for $\alpha$-harmonic functions under a much weaker condition. For this theorem we have two different proofs by using two different methods: One is a direct approach using potential theory. The other is by Fourier analysis as a corollary of the fact that the only $\alpha$-harmonic functions are affine.


title: Efficient boundary integral solution for acoustic wave scattering by irregular surfaces

abstract: The left-right operator splitting method is studied for the efficient calculation of acoustic fields scattered by arbitrary rough surfaces. Here the governing boundary integral is written as a sum of left- and right-going components, and the solution expressed as an iterative series, expanding about the predominant direction of propagation. Calculation of each term is computationally inexpensive both in time and memory, and the field is often accurately captured using one or two terms. The convergence and accuracy are examined by comparison with exact solution for smaller problems, and a series of much larger problems are tackled. The method is also immediately applicable to other scatterers such as waveguides, of which examples are given.


title: Washing out of the 0-pi transition in Josephson junctions

abstract: We consider a Josephson junction formed by a quantum dot connected to two bulk superconductors in presence of Coulomb interaction and coupling to both an electromagnetic environment and a finite density of electronic quasi-particles. In the limit of large superconducting gap we obtain a Born-Markov description of the system dynamics. We calculate the current-phase relation and we find that the experimentally unavoidable presence of quasi-particles can dramatically modify the 0-pi standard transition picture. We show that photon-assisted quasi-particles absorption allows the dynamic switching from the 0- to the pi-state and vice-versa, washing out the 0-pi transition predicted by purely thermodynamic arguments.


title: A hybrid variational principle for the Keller-Segel system in $\mathbb R^2$

abstract: We construct weak global in time solutions to the classical Keller-Segel system cell movement by chemotaxis in two dimensions when the total mass is below the well-known critical value. Our construction takes advantage of the fact that the Keller-Segel system can be realized as a gradient flow in a suitable functional product space. This allows us to employ a hybrid variational principle which is a generalisation of the minimising implicit scheme for Wasserstein distances introduced by Jordan, Kinderlehrer and Otto (1998).


title: Uniform Hausdorff measure of the level sets of the Brownian tree

abstract: Let $(\mathcal{T},d)$ be the random real tree with root $\rho$ coded by a Brownian excursion. So $(\mathcal{T},d)$ is (up to normalisation) Aldous CRT \cite{AldousI} (see Le Gall \cite{LG91}). The $a$-level set of $\mathcal{T}$ is the set $\mathcal{T}(a)$ of all points in $\mathcal{T}$ that are at distance $a$ from the root. We know from Duquesne and Le Gall \cite{DuLG06} that for any fixed $a\in (0, \infty)$, the measure $\ell^a$ that is induced on $\mathcal{T}(a)$ by the local time at $a$ of the Brownian excursion, is equal, up to a multiplicative constant, to the Hausdorff measure in $\mathcal{T}$ with gauge function $g(r)= r \log\log1/r$, restricted to $\mathcal{T}(a)$. As suggested by a result due to Perkins \cite{Per88,Per89} for super-Brownian motion, we prove in this paper a more precise statement that holds almost surely uniformly in $a$, and we specify the multiplicative constant. Namely, we prove that almost surely for any $a\in (0, \infty)$, $\ell^a(\cdot) = \frac{1}{2} \mathscr{H}_g (\, \cdot \, \cap \mathcal{T}(a))$, where $\mathscr{H}_g$ stands for the $g$-Hausdorff measure.


title: A. Stern's analysis of the nodal sets of some families of spherical harmonics revisited

abstract: In this paper, we revisit the analyses of Antonie Stern (1925) and Hans Lewy (1977) devoted to the construction of spherical harmonics with two or three nodal domains. Our method yields sharp quantitative results and a better understanding of the occurrence of bifurcations in the families of nodal sets.This paper is a natural continuation of our critical reading of A. Stern's results for Dirichlet eigenfunctions in the square, see arXiv:14026054.


title: On the consistency of Sobol indices with respect to stochastic ordering of model parameters

abstract: In the past decade, Sobol's variance decomposition have been used as a tool - among others - in risk management. We show some links between global sensitivity analysis and stochastic ordering theories. This gives an argument in favor of using Sobol's indices in uncertainty quantification, as one indicator among others.


title: Inverse problem on a tree-shaped network

abstract: In this article, we prove a uniqueness result for a coefficient inverse problems regarding a wave, a heat or a Schr\"odinger equation set on a tree-shaped network, as well as the corresponding stability result of the inverse problem for the wave equation. The objective is the determination of the potential on each edge of the network from the additional measurement of the solution at all but one external end-points. Our idea for proving the uniqueness is to use a traditional approach in coefficient inverse problem by Carleman estimate. Afterwards, using an observability estimate on the whole network, we apply a compactness-uniqueness argument and prove the stability for the wave inverse problem.


title: A New Effective Asymptotic Formula for the Stieltjes Constants

abstract: We derive a new integral formula for the Stieltjes constants. The new formula permits easy computations as well as an exact approximate asymptotic formula. Both the sign oscillations and the leading order of growth are provided. The formula can also be easily extended to generalized Euler constants.


title: On Form Factors in N=4 SYM Theory and Polytopes

abstract: In this paper we discuss different recursion relations (BCFW and all-line shift) for the form factors of the operators from the $\mathcal{N}=4$ SYM stress-tensor current supermultiplet $T^{AB}$ in momentum twistor space. We show that cancelations of spurious poles and equivalence between different types of recursion relations can be naturally understood using geometrical interpretation of the form factors as special limit of the volumes of polytopes in $\mathbb{C}\mathbb{P}^4$ in close analogy with the amplitude case. We also show how different relations for the IR pole coefficients can be easily derived using momentum twistor representation. This opens an intriguing question - which of powerful on-shell methods and ideas can survive off-shell ?


title: Optical Properties of (162173) 1999 JU3: In Preparation for the JAXA Hayabusa 2 Sample Return Mission

abstract: We investigated the magnitude-phase relation of (162173) 1999 JU3, a target asteroid for the JAXA Hayabusa 2 sample return mission. We initially employed the international Astronomical Union's H-G formalism but found that it fits less well using a single set of parameters. To improve the inadequate fit, we employed two photometric functions, the Shevchenko and Hapke functions. With the Shevchenko function, we found that the magnitude-phase relation exhibits linear behavior in a wide phase angle range (alpha = 5-75 deg) and shows weak nonlinear opposition brightening at alpha< 5 deg, providing a more reliable absolute magnitude of Hv = 19.25 +- 0.03. The phase slope (0.039 +- 0.001 mag/deg) and opposition effect amplitude (parameterized by the ratio of intensity at alpha=0.3 deg to that at alpha=5 deg, I(0.3)/I(5)=1.31+-0.05) are consistent with those of typical C-type asteroids. We also attempted to determine the parameters for the Hapke model, which are applicable for constructing the surface reflectance map with the Hayabusa 2 onboard cameras. Although we could not constrain the full set of Hapke parameters, we obtained possible values, w=0.041, g=-0.38, B0=1.43, and h=0.050, assuming a surface roughness parameter theta=20 deg. By combining our photometric study with a thermal model of the asteroid (Mueller et al. in preparation), we obtained a geometric albedo of pv = 0.047 +- 0.003, phase integral q = 0.32 +- 0.03, and Bond albedo AB = 0.014 +- 0.002, which are commensurate with the values for common C-type asteroids.


title: A WISE view of novae. I. The data

abstract: We present the results of data-mining the Wide-field Infrared Survey Explorer (WISE) archive for data on classical and recurrent novae. We find that the detections are consistent with dust emission, line emission, emission by a stellar photosphere, or a combination of these. Of the 36 novae detected in one or more WISE bands, 16 are detected in all four; thirty-one known novae are not detected by WISE. We also searched for WISE data on post-WISE novae, to gain information about nova progenitors. In this first paper we consider only the WISE data. In future papers we will provide a more detailed modelling of the WISE data, and discuss WISE data on post-WISE novae - including their variability - and will complement the WISE data with data from other IR surveys.


title: Arithmetical rank of strings and cycles

abstract: Let $R$ be a polynomial ring over a field $K$. To a given squarefree monomial ideal $I \subset R$, one can associate a hypergraph $H(I)$. In this article, we prove that the arithmetical rank of $I$ is equal to the projective dimension of $R/I$ when $H(I)$ is a string or a cycle hypergraph.


title: Secrecy Capacity Region of Some Classes of Wiretap Broadcast Channels

abstract: This work investigates the secrecy capacity of the Wiretap Broadcast Channel (WBC) with an external eavesdropper where a source wishes to communicate two private messages over a Broadcast Channel (BC) while keeping them secret from the eavesdropper. We derive a non-trivial outer bound on the secrecy capacity region of this channel which, in absence of security constraints, reduces to the best known outer bound to the capacity of the standard BC. An inner bound is also derived which follows the behavior of both the best known inner bound for the BC and the Wiretap Channel. These bounds are shown to be tight for the deterministic BC with a general eavesdropper, the semi-deterministic BC with a more-noisy eavesdropper and the Wiretap BC where users exhibit a less-noisiness order between them. Finally, by rewriting our outer bound to encompass the characteristics of parallel channels, we also derive the secrecy capacity region of the product of two inversely less-noisy BCs with a more-noisy eavesdropper. We illustrate our results by studying the impact of security constraints on the capacity of the WBC with binary erasure (BEC) and binary symmetric (BSC) components.


title: Disentangling phase transitions and critical points in the proton-neutron interacting boson model by catastrophe theory

abstract: We introduce the basic concepts of catastrophe theory needed to derive analytically the phase diagram of the proton-neutron interacting boson model (IBM-2). Previous studies [1,2,3] were based on numerical solutions. We here explain the whole IBM-2 phase diagram including the precise order of the phase transitions in terms of the cusp catastrophe.


title: A Novel Hybrid Crossover based Artificial Bee Colony Algorithm for Optimization Problem

abstract: Artificial bee colony (ABC) algorithm has proved its importance in solving a number of problems including engineering optimization problems. ABC algorithm is one of the most popular and youngest member of the family of population based nature inspired meta-heuristic swarm intelligence method. ABC has been proved its superiority over some other Nature Inspired Algorithms (NIA) when applied for both benchmark functions and real world problems. The performance of search process of ABC depends on a random value which tries to balance exploration and exploitation phase. In order to increase the performance it is required to balance the exploration of search space and exploitation of optimal solution of the ABC. This paper outlines a new hybrid of ABC algorithm with Genetic Algorithm. The proposed method integrates crossover operation from Genetic Algorithm (GA) with original ABC algorithm. The proposed method is named as Crossover based ABC (CbABC). The CbABC strengthens the exploitation phase of ABC as crossover enhances exploration of search space. The CbABC tested over four standard benchmark functions and a popular continuous optimization problem.


title: Long fully commutative elements in affine Coxeter groups

abstract: An element of a Coxeter group $W$ is called fully commutative if any two of its reduced decompositions can be related by a series of transpositions of adjacent commuting generators. In the preprint "Fully commutative elements in finite and affine Coxeter groups" (arXiv:1402.2166), R. Biagioli and the authors proved among other things that, for each irreducible affine Coxeter group, the sequence counting fully commutative elements with respect to length is ultimately periodic. In the present work, we study this sequence in its periodic part for each of these groups, and in particular we determine the minimal period. We also observe that in type $A$ affine we get an instance of the cyclic sieving phenomenon.


title: Non-geometric orbifolds and wrapping rules

abstract: We show that the number of half-supersymmetric p-branes in the Type II theories compactified on orbifolds is determined by the wrapping rules recently introduced, provided that one accounts correctly for both geometric and non-geometric T-dual configurations. Starting from the Type II theories compactified on K3, we analyze their toroidal dimensional reductions, showing how the resulting half-supersymmetric p-branes satisfy the wrapping rules only by taking into account all the possible higher-dimensional origins. We then consider Type II theories compactified on the orbifold T^6/(Z_2 \times Z_2 ), whose massless four-dimensional theory is an N=2 supergravity. Again, the wrapping rules are obeyed only if one includes the complete orbit of the T-duality group, namely either Type IIA or Type IIB theories compactified on either the geometric or the non-geometric T-dual orbifold. Finally, we comment on the interpretation of our results in the framework of the duality between the Heterotic string compactified on K3 \times T^2 and the Type II string compactified on a Calabi-Yau threefold.


title: Detection efficiency loophole in Pusey-Barrett-Rudolph theorem

abstract: Detection efficiency loophole poses a significant problem for experimental tests of Bell inequalities. Recently discovered Pusey-Barrett-Rudolph (PBR) theorem suffers from the same vulnerability. In this paper we calculate the critical detection efficiency, below which the PBR argument for the ontic nature of quantum state is inconclusive. This is done for the maximally $\psi$-epistemic models. We use two different definitions of this property. The optimal number of parties, for which the critical detection efficiency is the lowest is given. We also approach the problem from the opposite direction. We provide a function which enables us to specify which epistemic models are ruled out by the results of an experiment with a given detection efficiency.


title: A special point problem of Andr\'e-Pink-Zannier in the universal family of abelian varieties

abstract: The Andr\'e-Pink-Zannier conjecture concerns the intersection of subvarieties and the generalized Hecke orbit of a given point in mixed Shimura varieties. It is part of the Zilber-Pink conjecture. In this paper we focus on the universal family of principally polarized abelian varieties. We explain the moduli interpretation of the conjecture in this case and prove several different cases for this conjecture: its overlap with the Andr\'e-Oort conjecture; when the subvariety is contained in an abelian scheme over a curve and the point is a torsion point on its fiber; when the subvariety is a curve.


title: Improved Reduced Models for Single-Pass and Reflective Semiconductor Optical Amplifiers

abstract: We present highly accurate and easy to implement, improved lumped semiconductor optical amplifier (SOA) models for both single-pass and reflective semiconductor optical amplifiers (RSOA). The key feature of the model is the inclusion of the internal losses and we show that a few subdivisions are required to achieve an accuracy of 0.12 dB. For the case of RSOAs, we generalize a recently published model to account for the internal losses that are vital to replicate observed RSOA behavior. The results of the improved reduced RSOA model show large overlap when compared to a full bidirectional travelling wave model over a 40 dB dynamic range of input powers and a 20 dB dynamic range of reflectivity values. The models would be useful for the rapid system simulation of signals in communication systems, i.e. passive optical networks that employ RSOAs, signal processing using SOAs and for implementing digital back propagation to undo amplifier induced signal distortions.


title: Non-local formulation of ghost-free bigravity theory

abstract: We study the ghost-free bimetric theory of Hassan and Rosen, with parameters $\beta_i$ such that a flat Minkowski solution exists for both metrics. We show that, expanding around this solution and eliminating one of the two metrics with its own equation of motion, the remaining metric is governed by the Einstein-Hilbert action plus a non-local term proportional to $W_{\mu\nu\rho\sigma} (\Box-m^2)^{-1}W^{\mu\nu\rho\sigma}$, where $W_{\mu\nu\rho\sigma}$ is the Weyl tensor. The result is valid to quadratic order in the metric perturbation and to all orders in the derivative expansion. This example shows, in a simple setting, how such non-local extensions of GR can emerge from an underlying consistent theory, at the purely classical level.


title: Wind and Wave Extremes over the World Oceans from Very Large Ensembles

abstract: Global return values of marine wind speed and significant wave height are estimated from very large aggregates of archived ensemble forecasts at +240-h lead time. Long lead time ensures that the forecasts represent independent draws from the model climate. Compared with ERA-Interim, a reanalysis, the ensemble yields higher return estimates for both wind speed and significant wave height. Confidence intervals are much tighter due to the large size of the dataset. The period (9 yrs) is short enough to be considered stationary even with climate change. Furthermore, the ensemble is large enough for non-parametric 100-yr return estimates to be made from order statistics. These direct return estimates compare well with extreme value estimates outside areas with tropical cyclones. Like any method employing modeled fields, it is sensitive to tail biases in the numerical model, but we find that the biases are moderate outside areas with tropical cyclones.


title: Electronic Instabilities of the AA-Honeycomb Bilayer

abstract: We use a functional renormalization group approach to study the instabilities due to electron-electron interactions in a bilayer honeycomb lattice model with AA stacking, as it might be relevant for layered graphene with this structure. Starting with a tight- binding description for the four $\pi$-bands, we integrate out the modes of the dispersion by successively lowering an infrared cutoff and determine the leading tendencies in the effective interactions. The antiferromagnetic spin-density wave is an expected instability for dominant local repulsion among the electrons, but for nonlocal interaction terms also other instabilities occur. We discuss the phase diagrams depending on the model parameters. We compare our results to single-layer graphene and the more common AB-stacked bilayer, both qualitatively and quantitatively.


title: A New Determination of the Spectra and Luminosity Function of Gamma-Ray Millisecond Pulsars

abstract: In this article, we revisit the gamma-ray emission observed from millisecond pulsars and globular clusters. Based on 5.6 years of data from the Fermi Gamma-Ray Space Telescope, we report gamma-ray spectra for 61 millisecond pulsars, finding most to be well fit by a power-law with an exponential cutoff, producing to a spectral peak near ~1-2 GeV (in $E^2 dN/dE$ units). Additionally, while most globular clusters exhibit a similar spectral shape, we identify a few with significantly softer spectra. We also determine the gamma-ray luminosity function of millisecond pulsars using the population found in the nearby field of the Milky Way, and within the globular cluster 47 Tucanae. We find that the gamma-ray emission observed from globular clusters is dominated by a relatively small number of bright millisecond pulsars, and that low-luminosity pulsars account for only a small fraction of the total flux. Our results also suggest that the gamma-ray emission from millisecond pulsars is more isotropic and less strongly beamed than the emission at X-ray wavelengths. Furthermore, the observed distribution of apparent gamma-ray efficiencies provides support for the slot gap or the outer gap models over those in which the gamma-ray emission originates from regions close to the neutron star's magnetic poles (polar cap models).


title: D-modules on rigid analytic spaces

abstract: We give an overview of the theory of $\wideparen{\mathcal{D}}$-modules on rigid analytic spaces and its applications to admissible locally analytic representations of $p$-adic Lie groups.


title: Potential Theoretic Hyperbolicity and $\mathbf{L^2}$ extension. Part I: Stein manifolds

abstract: We establish a new generalization of an $L^2$ extension theorem of Ohsawa-Takegoshi type. The improvement in the theorem is that it allows the usual curvature assumptions to be significantly weakened in certain favorable settings. The favorable settings come about, for example, when the underlying structure (e.g., the underlying manifold, or the holomorphic Hermitian line bundle whose sections are being extended) has certain potential theoretic positivity. The simplest case occurs when the underlying manifold supports a function with self-bounded gradient in the sense of McNeal, but there are other cases. We demonstrate the improvement over the usual Ohsawa-Takegoshi-type extension theorems through a number of examples.


title: Variation in Microbial Growth under Hypergravity

abstract: We report bacterial growth under hypergravitational stress. Cultures of E. coli and B. subtilis were subjected to the gravitational stress (38g) and their growth curves were measured using UV-VIS spectrophotometer. Experiments were also carried out to investigate nutrient consumption under hypergravitational conditions. Our results show considerable difference between samples subjected to hypergravity and normal conditions. This study has importance to understand bacterial response to external stress factors like gravity and changes in bacterial system in order to adapt with stress conditions for its survival.


title: Weihrauch degrees of finding equilibria in sequential games

abstract: We consider the degrees of non-computability (Weihrauch degrees) of finding winning strategies (or more generally, Nash equilibria) in infinite sequential games with certain winning sets (or more generally, outcome sets). In particular, we show that as the complexity of the winning sets increases in the difference hierarchy, the complexity of constructing winning strategies increases in the effective Borel hierarchy.


title: Study of nonclassicality in tripartite correlations beyond standard Bell nonlocality paradigm

abstract: In the nonsignaling framework, nonclassicality in correlation arising from two spatially separated input-output devices gets manifested, solely, through its \emph{nonlocal} behavior. Study of correlations based on this said feature is commonly known as local-nonlocal paradigm. While in two-party scenario correlations can be of only two types either local or nonlocal, the situation gets more involved for multi-party scenario, \emph{e.g.,} for tripartite scenario, correlations can be of three types: fully local, two-way local, and genuinely nonlocal. Nonsignaling correlations having quantum realization are termed physical. Fully local and certain quantum realizable two-way local tripartite correlations always have a quantum realization with tripartite biseparable states if there is no restriction on the local Hilbert-space dimensions. In this work, we study the quantum simulation of fully local and two-way local tripartite correlations with restricted local Hilbert-space dimensions, in particular we consider $\mathbb{C}^2\otimes\mathbb{C}^2\otimes\mathbb{C}^2$ quantum systems. Interestingly, in this restricted simulation scenario we find that simulation of certain fully local and two-way local correlations necessarily requires \emph{genuine quantumness} in the three qubit states. This, going beyond the standard nonlocality paradigm, captures a new notion of genuine nonclassicality even in the fully local and two-way local correlations. To explore this newly introduced notion of genuine nonclassicality, we propose two quantities of interest, called \emph{Svetlichny strength} and \emph{Mermin strength}, and extensively study their properties.


title: Quantum Correlations in Cavity QED Networks

abstract: We present a Review of the dynamical features such as generation, propagation, distribution, sudden transition and freezing of the various quantum correlation measures, as Concurrence, Entanglement of Formation, Quantum Discord, as well their geometrical measure counterparts within the models of Cavity Quantum Electrodynamics Networks. Dissipation and thermal effects are discussed both in the generation of quantum correlations as well as their effect on the sudden changes and freezing of the classical and quantum correlations in a cavity quantum electrodynamical network. For certain initial conditions, double transitions in the Bures geometrical discord are found. One of these transitions tends to disappear at a critical temperature


title: Learning about Learning: Human Brain Sub-Network Biomarkers in fMRI Data

abstract: It has become increasingly popular to study the brain as a network due to the realization that functionality cannot be explained exclusively by independent activation of specialized regions. Instead, across a large spectrum of behaviors, function arises due to the dynamic interactions between brain regions. The existing literature on functional brain networks focuses mainly on a battery of network properties characterizing the "resting state" using for example the modularity, clustering, or path length among regions. In contrast, we seek to uncover subgraphs of functional connectivity that predict or drive individual differences in sensorimotor learning across subjects. We employ a principled approach for the discovery of significant subgraphs of functional connectivity, induced by brain activity (measured via fMRI imaging) while subjects perform a motor learning task. Our aim is to uncover patterns of functional connectivity that discriminate between high and low rates of learning among subjects. The discovery of such significant discriminative subgraphs promises a better data-driven understanding of the dynamic brain processes associated with brain plasticity.


title: Autonomous models on a Cayley tree

abstract: The most general single species autonomous reaction-diffusion model on a Cayley tree with nearest-neighbor interactions is introduced. The stationary solutions of such models, as well as their dynamics, are discussed. To study dynamics of the system, directionally-symmetric Green function for evolution equation of average number density is obtained. In some limiting cases the Green function is studied. Some examples are worked out in more detail.


title: Generalization of the Fedorova-Schmidt method for determining particle size distributions

abstract: One reports the integral transform that determines the particle size distribution of a given sample from the small-angle scattering intensity under the assumption that the particle correlation function is a polynomial of degree M. The Fedorova-Schmidt solution [J. Appl. Cryst. 11, 405, (1978)] corresponds to the case M = 3. The procedure for obtaining a polynomial approximation to a particle correlation function is discussed and applied to the cases of polidisperse particles of tetrahedral or octahedral or cubical shape.


title: Deterministic Versus Randomized Kaczmarz Iterative Projection

abstract: Kaczmarz's alternating projection method has been widely used for solving a consistent (mostly over-determined) linear system of equations Ax=b. Because of its simple iterative nature with light computation, this method was successfully applied in computerized tomography. Since tomography generates a matrix A with highly coherent rows, randomized Kaczmarz algorithm is expected to provide faster convergence as it picks a row for each iteration at random, based on a certain probability distribution. It was recently shown that picking a row at random, proportional with its norm, makes the iteration converge exponentially in expectation with a decay constant that depends on the scaled condition number of A and not the number of equations. Since Kaczmarz's method is a subspace projection method, the convergence rate for simple Kaczmarz algorithm was developed in terms of subspace angles. This paper provides analyses of simple and randomized Kaczmarz algorithms and explain the link between them. It also propose new versions of randomization that may speed up convergence.


title: Fermi-Level Pinning, Charge Transfer, and Relaxation of Spin-Momentum Locking at Metal Contacts to Topological Insulators

abstract: Topological insulators are of interest for many applications in electronics and optoelectronics, but harnessing their unique properties requires detailed understanding and control of charge injection at electrical contacts. Here we present large-scale ab initio calculations of the electronic properties of Au, Ni, Pt, Pd, and graphene contacts to Bi2Se3. We show that regardless of the metal, the Fermi level is located in the conduction band, leading to n-type Ohmic contact to the first quintuplet. Furthermore, we find strong charge transfer and band-bending in the first few quintuplets, with no Schottky barrier for charge injection even when the topoplogical insulator is undoped. Our calculations indicate that Au and graphene leave the spin-momentum locking mostly unaltered, but on the other hand, Ni, Pd, and Pt strongly hybridize with Bi2Se3 and relax spin-momentum locking. Our results indicate that judicious choice of the contact metal is essential to reveal the unique surface features of topological insulators.


title: Dissipation in Josephson tunneling junctions at low temperatures

abstract: It is important to know the decoherence mechanism of a qubit based on Josephson junctions. At low temperatures, as quasiparticle concentration becomes exponentially small, one needs to consider energy transfer from tunneling electrons to other degrees of freedom to find dissipation in Josephson junctions and decoherence in qubits. Here we discuss the energy transfer to two-level systems, i.e. the transitions between two different configurations of ions inside insulating layer separated by a potential barrier. We derive a general equation of motion for the phase difference between two superconducting electrodes and we find a retarded dissipation term due to electromagnetic mechanism and also contribution due to electron tunneling mechanism. Using the equation of motion we calculate the decay of Rabi oscillations and frequency shift in qubits due to the presence of the two-level systems. In the long time limit our results coincide with those obtained by Martinis et al. [J. M. Martinis et al.. Phys. Rev. Lett. 95, 210503 (2005)] within the Fermi's Golden rule approach up to a numerical factor.


title: Impact of quantized vibrations on the efficiency of interfacial charge separation in photovoltaic devices

abstract: We demonstrate that charge separation at donor-acceptor interfaces is a complex process that is controlled by the combined action of Coulomb binding for electron-hole pairs and partial relaxation due to quantized phonons. A joint electron-vibration quantum dynamical study reveals that high energy vibrations sensitively tune the charge transfer probability as a function of time and injection energy, due to polaron formation. These results have bearings for the optimization of energy transfer both in organic and quantum dot photovoltaics, as well as in biological light harvesting complexes.

categories: cond-mat.mtrl-sci cond-mat.mes-hall physics.chem-ph quant-ph


title: Causality Constraints on Corrections to the Graviton Three-Point Coupling

abstract: We consider higher derivative corrections to the graviton three-point coupling within a weakly coupled theory of gravity. Lorentz invariance allows further structures beyond the one present in the Einstein theory. We argue that these are constrained by causality. We devise a thought experiment involving a high energy scattering process which leads to causality violation if the graviton three-point vertex contains the additional structures. This violation cannot be fixed by adding conventional particles with spins $J \leq 2$. But, it can be fixed by adding an infinite tower of extra massive particles with higher spins, $J > 2$. In AdS theories this implies a constraint on the conformal anomaly coefficients $\left|{a - c \over c} \right| \lesssim {1 \over \Delta_{gap}^2}$ in terms of $\Delta_{gap}$, the dimension of the lightest single particle operator with spin $J > 2$. For inflation, or de Sitter-like solutions, it indicates the existence of massive higher spin particles if the gravity wave non-gaussianity deviates significantly from the one computed in the Einstein theory.


title: Fractional Gaussian fields: a survey

abstract: We discuss a family of random fields indexed by a parameter $s\in \mathbb{R}$ which we call the fractional Gaussian fields, given by \[ \mathrm{FGF}_s(\mathbb{R}^d)=(-\Delta)^{-s/2} W, \] where $W$ is a white noise on $\mathbb{R}^d$ and $(-\Delta)^{-s/2}$ is the fractional Laplacian. These fields can also be parameterized by their Hurst parameter $H = s-d/2$. In one dimension, examples of $\mathrm{FGF}_s$ processes include Brownian motion ($s = 1$) and fractional Brownian motion ($1/2 < s < 3/2$). Examples in arbitrary dimension include white noise ($s = 0$), the Gaussian free field ($s = 1$), the bi-Laplacian Gaussian field ($s = 2$), the log-correlated Gaussian field ($s = d/2$), L\'evy's Brownian motion ($s = d/2 + 1/2$), and multidimensional fractional Brownian motion ($d/2 < s < d/2 + 1$). These fields have applications to statistical physics, early-universe cosmology, finance, quantum field theory, image processing, and other disciplines. We present an overview of fractional Gaussian fields including covariance formulas, Gibbs properties, spherical coordinate decompositions, restrictions to linear subspaces, local set theorems, and other basic results. We also define a discrete fractional Gaussian field and explain how the $\mathrm{FGF}_s$ with $s \in (0,1)$ can be understood as a long range Gaussian free field in which the potential theory of Brownian motion is replaced by that of an isotropic $2s$-stable L\'evy process.


title: Scalable Kernel Methods via Doubly Stochastic Gradients

abstract: The general perception is that kernel methods are not scalable, and neural nets are the methods of choice for nonlinear learning problems. Or have we simply not tried hard enough for kernel methods? Here we propose an approach that scales up kernel methods using a novel concept called "doubly stochastic functional gradients". Our approach relies on the fact that many kernel methods can be expressed as convex optimization problems, and we solve the problems by making two unbiased stochastic approximations to the functional gradient, one using random training points and another using random functions associated with the kernel, and then descending using this noisy functional gradient. We show that a function produced by this procedure after $t$ iterations converges to the optimal function in the reproducing kernel Hilbert space in rate $O(1/t)$, and achieves a generalization performance of $O(1/\sqrt{t})$. This doubly stochasticity also allows us to avoid keeping the support vectors and to implement the algorithm in a small memory footprint, which is linear in number of iterations and independent of data dimension. Our approach can readily scale kernel methods up to the regimes which are dominated by neural nets. We show that our method can achieve competitive performance to neural nets in datasets such as 8 million handwritten digits from MNIST, 2.3 million energy materials from MolecularSpace, and 1 million photos from ImageNet.


title: A multiblob approach to colloidal hydrodynamics with inherent lubrication

abstract: This work presents an intermediate resolution model of the hydrodynamics of colloidal particles based on a mixed Eulerian-Lagrangian formulation. The particle is constructed with a small set of overlapping Peskin's Immersed Boundary kernels (blobs) which are held together by springs to build up a particle impenetrable core. Here, we used 12 blobs placed in the vertexes of an icosahedron with an extra one in its center. Although the particle surface is not explicitly resolved, we show that the short-distance hydrodynamic responses (flow profiles, translational and rotational mobilities, lubrication, etc) agree with spherical colloids and provide consistent effective radii. A remarkable property of the present multiblob model is that it naturally presents a "divergent" lubrication force at finite inter-particle distance. This permits to resolve the large viscosity increase at dense colloidal volume fractions. The intermediate resolution model is able to recover highly non-trivial (many-body) hydrodynamics using small particles whose radii are similar to the grid size $h$ (in the range $[1.6-3.2]\,h$). Considering that the cost of the embedding fluid phase scales like the cube of the particle radius, this result brings about a significant computational speed-up. Our code Fluam works in Graphics Processor Units (GPU's) and uses Fast Fourier Transform for the Poisson solver, which further improves its efficiency.


title: The Inner Life of the Kondo Ground State: An Answer to Kenneth Wilson's Question

abstract: The Kondo ground state has been investigated by numerical and exact methods, but the physics behind these results remains veiled. Nobel prize winner Wilson, who engineered the break through in his numerical renormalization group theory, commented in his review article "the author has no simple explanation ...for the crossover from weak to strong coupling". In this article a graphical interpretation is given for the extraordinary properties of the Kondo ground state. At the crossover all electron states in the low energy range of k_{B}T_{K} are synchronized. An internal orthogonality catastrophe is averted.


title: Predictive support recovery with TV-Elastic Net penalty and logistic regression: an application to structural MRI

abstract: The use of machine-learning in neuroimaging offers new perspectives in early diagnosis and prognosis of brain diseases. Although such multivariate methods can capture complex relationships in the data, traditional approaches provide irregular (l2 penalty) or scattered (l1 penalty) predictive pattern with a very limited relevance. A penalty like Total Variation (TV) that exploits the natural 3D structure of the images can increase the spatial coherence of the weight map. However, TV penalization leads to non-smooth optimization problems that are hard to minimize. We propose an optimization framework that minimizes any combination of l1, l2, and TV penalties while preserving the exact l1 penalty. This algorithm uses Nesterov's smoothing technique to approximate the TV penalty with a smooth function such that the loss and the penalties are minimized with an exact accelerated proximal gradient algorithm. We propose an original continuation algorithm that uses successively smaller values of the smoothing parameter to reach a prescribed precision while achieving the best possible convergence rate. This algorithm can be used with other losses or penalties. The algorithm is applied on a classification problem on the ADNI dataset. We observe that the TV penalty does not necessarily improve the prediction but provides a major breakthrough in terms of support recovery of the predictive brain regions.


title: Formation of granular structures in trapped Bose-Einstein condensates under oscillatory excitations

abstract: We present experimental observations and numerical simulations of nonequilibrium spatial structures in a trapped Bose-Einstein condensate subject to oscillatory perturbations. In experiment, first, there appear collective excitations, followed by quantum vortices. Increasing the amount of the injected energy leads to the formation of vortex tangles representing quantum turbulence. We study what happens after the regime of quantum turbulence, with increasing further the amount of injected energy. In such a strongly nonequilibrium Bose-condensed system of trapped atoms, vortices become destroyed and there develops a new kind of spatial structure exhibiting essentially heterogeneous spatial density. The structure reminds fog consisting of high-density droplets, or grains, surrounded by the regions of low density. The grains are randomly distributed in space, where they move. They live sufficiently long time to be treated as a type of metastable objects. Such structures have been observed in nonequilibrium trapped Bose gases of $^{87}$Rb, subject to the action of alternating fields. Here we present experimental results and support them by numerical simulations. The granular, or fog structure is essentially different from the state of wave turbulence that develops after increasing further the amount of injected energy.


title: Weak Galerkin method for the coupled Darcy-Stokes flow

abstract: A family of weak Galerkin finite element discretization is developed for solving the coupled Darcy-Stokes equation. The equation in consideration admits the Beaver-Joseph-Saffman condition on the interface. By using the weak Galerkin approach, in the discrete space we are able to impose the normal continuity of velocity explicitly. Or in other words, strong coupling is achieved in the discrete space. Different choices of weak Galerkin finite element spaces are discussed, and error estimates are given.


title: Log-correlated Gaussian fields: an overview

abstract: We survey the properties of the log-correlated Gaussian field (LGF), which is a centered Gaussian random distribution (generalized function) $h$ on $\mathbb R^d$, defined up to a global additive constant. Its law is determined by the covariance formula $$\mathrm{Cov}\bigl[ (h, \phi_1), (h, \phi_2) \bigr] = \int_{\mathbb R^d \times \mathbb R^d} -\log|y-z| \phi_1(y) \phi_2(z)dydz$$ which holds for mean-zero test functions $\phi_1, \phi_2$. The LGF belongs to the larger family of fractional Gaussian fields obtained by applying fractional powers of the Laplacian to a white noise $W$ on $\mathbb R^d$. It takes the form $h = (-\Delta)^{-d/4} W$. By comparison, the Gaussian free field (GFF) takes the form $(-\Delta)^{-1/2} W$ in any dimension. The LGFs with $d \in \{2,1\}$ coincide with the 2D GFF and its restriction to a line. These objects arise in the study of conformal field theory and SLE, random surfaces, random matrices, Liouville quantum gravity, and (when $d=1$) finance. Higher dimensional LGFs appear in models of turbulence and early-universe cosmology. LGFs are closely related to cascade models and Gaussian branching random walks. We review LGF approximation schemes, restriction properties, Markov properties, conformal symmetries, and multiplicative chaos applications.