title: Consistent picture of the octet-nodal gap and its evolution with doping in heavily overdoped Ba$_{1-x}$K$_{x}$Fe$_{2}$As$_{2}$

abstract: We investigate the pairing symmetry in heavily overdoped Ba$_{1-x}$K$_{x}$Fe$_{2}$As$_{2}$ based on the spin-fluctuation mechanism. The exotic octet nodes of the superconducting gap and the unusual evolution of the gap with doping observed by the recent experiments are well explained in a unified manner. We demonstrate that the scatterings of electrons on the Fermi patches is mainly responsible for the incommensurate spin fluctuations and consequently the Fermi-surface-dependent multi-gap structure, since the Fermi level is close to the flat band. In addition, we find that a $d$-wave pairing state will prevail over the s-wave pairing state around the Lifshitz transition point.


title: Improving Efficiency and Scalability of Formula-based Debugging

abstract: Formula-based debugging techniques are becoming increasingly popular, as they provide a principled way to identify potentially faulty statements together with information that can help fix such statements. Although effective, these approaches are computationally expensive, which limits their practical applicability. Moreover, they tend to focus on failing test cases alone, thus ignoring the wealth of information provided by passing tests. To mitigate these issues, we propose two techniques: on-demand formula computation (OFC) and clause weighting (CW). OFC improves the overall efficiency of formula-based debugging by exploring all and only the parts of a program that are relevant to a failure. CW improves the accuracy of formula-based debugging by leveraging statistical fault-localization information that accounts for passing tests. Our empirical results show that both techniques are effective and can improve the state of the art in formula-based debugging.


title: Properties of electrons scattered on a strong plane electromagnetic wave with a linear polarization: classical treatment

abstract: The relations among the components of the exit momenta of ultrarelativistic electrons scattered on a strong electromagnetic wave of a low (optical) frequency and linear polarization are established using the exact solutions to the equations of motion with radiation reaction included (the Landau-Lifshitz equation). It is found that the momentum components of the electrons traversed the electromagnetic wave depend weakly on the initial values of the momenta. These electrons are mostly scattered at the small angles to the direction of propagation of the electromagnetic wave. The maximum Lorentz factor of the electrons crossed the electromagnetic wave is proportional to the work done by the electromagnetic field and is independent of the initial momenta. The momentum component parallel to the electric field strength vector of the electromagnetic wave is determined only by the diameter of the laser beam measured in the units of the classical electron radius. As for the reflected electrons, they for the most part lose the energy, but remain relativistic. There is a reflection law for these electrons that relates the incident and the reflection angles and is independent of any parameters.


title: Compact spacelike surfaces whose mean curvature function satisfies a nonlinear inequality in a 3-dimensional Generalized Robertson-Walker spacetime

abstract: Spacelike surfaces in Generalized Robertson-Walker spacetimes whose mean curvature function satisfies a natural nonlinear inequality are analyzed. Several uniqueness and nonexistence results for such compact spacelike surfaces are proved. In the nonparametric case, new Calabi-Bernstein type problems are solved as a consequence.


title: Threats from space: 20 years of progress

abstract: It has been 20 years since planning began for the 1995 United Nations International Conference on Near-Earth Objects. The conference proceedings established the scientific basis for an international organizational framework to support research and collective actions to mitigate a potential near-Earth object (NEO) threat to the planet. Since that time, researchers have conducted telescope surveys that should, within the coming decade, answer many questions about the size, number, and Earth impact probability of these objects. Space explorations to asteroids and comets have been successfully carried out, including sample recovery. Laboratory experiments and computer simulations at Sandia National Laboratories have analyzed the effects of soft X-ray radiation on meteorites - which might help researchers develop a way to redirect an incoming asteroid by vaporizing a thin layer of its surface. An Action Team on NEOs, established in 2001 in response to recommendations of the Third United Nations Conference on the Exploration and Peaceful Uses of Outer Space, identified the primary components of NEO mitigation and emphasized the value of finding potentially hazardous NEOs as soon as possible. Recommendations from the Action Team are meant to ensure that all nations are aware of the NEO danger; and to coordinate mitigation activities among nations that could be affected by an impact, as well as those that might play an active role in any eventual deflection or disruption campaign.


title: Numerical Computations of Separability Probabilities

abstract: We compute the probability that a bipartite quantum state is separable by Monte Carlo sampling. This is carried out for rebits, qubits and quaterbits. We sampled $5\times 10^{11}$ points for each of these three cases. The results strongly support conjectures for certain rational values of these probabilities that have been found by other methods.


title: Risk-Limiting Dynamic Contracts for Direct Load Control

abstract: This paper proposes a novel continuous-time dynamic contract framework that has a risk-limiting capability. If a principal and an agent enter into such a contract, the principal can optimally manage its performance and risk with a guarantee that the agent's risk is less than or equal to a pre-specified level and that the agent's expected payoff is greater than or equal to another pre-specified threshold. We achieve such risk-management capabilities by formulating the contract design problem as mean-variance constrained risk-sensitive control. A dynamic programming-based method is developed to solve the problem. The key idea of our proposed solution method is to reformulate the inequality constraints on the mean and the variance of the agent's payoff as dynamical system constraints by introducing new state and control variables. The reformulations use the martingale representation theorem. The proposed contract method enables us to develop a new direct load control method that provides the load-serving entity with financial risk management solutions in real-time electricity markets. We also propose an approximate decomposition of the optimal contract design problem for multiple customers into multiple low-dimensional contract problems for one customer. This allows the direct load control program to work with a large number of customers without any scalability issues. Furthermore, the contract design procedure can be completely parallelized. The performance and usefulness of the proposed contract method and its application to direct load control are demonstrated using data on the electric energy consumption of customers in Austin, Texas as well as the Electricity Reliability Council of Texas' locational marginal price data.


title: Hypercontractivity for Stochastic Hamiltonian Systems

abstract: The hypercontractivity is proved for the Markov semigroup associated to a class of finite/infinite dimensional stochastic Hamiltonian systems. Consequently, the Markov semigroup is exponentially convergent to the invariant probability measure in entropy (thus, also in $L^2$), and is compact for large time. Since the log-Sobolev inequality is invalid for the associated Dirichlet form, we introduce a general result on the hypercontractivity using the Harnack inequality with power. The main results are illustrated by concrete examples.


title: Turbulent channel flow of dense suspensions of neutrally-buoyant spheres

abstract: Dense particle suspensions are widely encountered in many applications and in environmental flows. While many previous studies investigate their rheological properties in laminar flows, little is known on the behaviour of these suspensions in the turbulent/inertial regime. The present study aims to fill this gap by investigating the turbulent flow of a Newtonian fluid laden with solid neutrally-buoyant spheres at relatively high volume fractions in a plane channel. Direct Numerical Simulation are performed in the range of volume fractions Phi=0-0.2 with an Immersed Boundary Method used to account for the dispersed phase. The results show that the mean velocity profiles are significantly altered by the presence of a solid phase with a decrease of the von Karman constant in the log-law. The overall drag is found to increase with the volume fraction, more than one would expect just considering the increase of the system viscosity due to the presence of the particles. At the highest volume fraction here investigated, Phi=0.2, the velocity fluctuation intensities and the Reynolds shear stress are found to decrease. The analysis of the mean momentum balance shows that the particle-induced stresses govern the dynamics at high Phi and are the main responsible of the overall drag increase. In the dense limit, we therefore find a decrease of the turbulence activity and a growth of the particle induced stress, where the latter dominates for the Reynolds numbers considered here.


title: Dyadic shift randomization in classical discrepancy theory

abstract: Dyadic shifts of point distributions in the multi-dimensional unit cube are considered as a randomization. Explicit formulas for the discrepancies of such randomized distributions are given in the paper in terms of Rademacher functions. Relaying on the statistical independence of Rademacher functions, Khinchin's inequalities, and other related results, we obtain very sharp upper and lower bounds for the mean discrepancies.


title: Surface tension and the mechanics of liquid inclusions in compliant solids

abstract: Eshelby's theory of inclusions has wide-reaching implications across the mechanics of materials and structures including the theories of composites, fracture, and plasticity. However, it does not include the effects of surface stress, which has recently been shown to control many processes in soft materials such as gels, elastomers and biological tissue. To extend Eshelby's theory of inclusions to soft materials, we consider liquid inclusions within an isotropic, compressible, linear-elastic solid. We solve for the displacement and stress fields around individual stretched inclusions, accounting for the bulk elasticity of the solid and the surface tension (\textit{i.e.} isotropic strain-independent surface stress) of the solid-liquid interface. Surface tension significantly alters the inclusion's shape and stiffness as well as its near- and far-field stress fields. These phenomenon depend strongly on the ratio of inclusion radius, $R$, to an elastocapillary length, $L$. Surface tension is significant whenever inclusions are smaller than $100L$. While Eshelby theory predicts that liquid inclusions generically reduce the stiffness of an elastic solid, our results show that liquid inclusions can actually stiffen a solid when $R<3L/2$. Intriguingly, surface tension cloaks the far-field signature of liquid inclusions when $R=3L/2$. These results are have far-reaching applications from measuring local stresses in biological tissue, to determining the failure strength of soft composites.


title: A primal-simplex based Tardos' algorithm

abstract: In the mid-eighties Tardos proposed a strongly polynomial algorithm for solving linear programming problems for which the size of the coefficient matrix is polynomially bounded by the dimension. Combining Orlin's primal-based modification and Mizuno's use of the simplex method, we introduce a modification of Tardos' algorithm considering only the primal problem and using simplex method to solve the auxiliary problems. The proposed algorithm is strongly polynomial if the coefficient matrix is totally unimodular and the auxiliary problems are non-degenerate.


title: On quantum percolation in finite regular graphs

abstract: The aim of this paper is twofold. First, we study eigenvalues and eigenvectors of the adjacency matrix of a bond percolation graph when the base graph is finite and well approximated locally by an infinite regular graph. We relate quantitatively the empirical measure of the eigenvalues and the delocalization of the eigenvectors to the spectrum of the adjacency operator of the percolation on the infinite graph. Secondly, we prove that percolation on an infinite regular tree with degree at least $3$ preserves the existence of an absolutely continuous spectrum if the removal probability is small enough. These two results are notably relevant for bond percolation on a uniformly sampled regular graph or a Cayley graph with large girth.


title: Discrete constant mean curvature nets in space forms: Steiner's formula and Christoffel duality

abstract: We show that the discrete principal nets in quadrics of constant curvature that have constant mixed area mean curvature can be characterized by the existence of a K\"onigs dual in a concentric quadric.


title: A Martingale Approach to Detect Peak of News in Social Network

abstract: Nowadays, social medias such as Twitter, Memetracker and Blogs have become powerful tools to propagate information. They facilitate quick dissemination sequence of information such as news article, blog posts, user's interests and thoughts through large scale. Providing strong means to analyzing social networks structure and how information diffuse through them is essential. Many recent studies emphasize on modeling information diffusion and their patterns to gain some useful knowledge. In this paper, we propose a statistical approach to online detect peak points of news when spread over social networks, to the best of our knowledge has never investigated before. The proposed model use martingale approach to predict peak points when news reached the peak of its popularity. Experimental results on real datasets show good performance of our approach to online detect these peak points.


title: Complexity of a Problem Concerning Reset Words for Eulerian Binary Automata

abstract: A word is called a reset word for a deterministic finite automaton if it maps all the states of the automaton to a unique state. Deciding about the existence of a reset word of a given maximum length for a given automaton is known to be an NP-complete problem. We prove that it remains NP-complete even if restricted to Eulerian automata with binary alphabets, as it has been conjectured by Martyugin (2011).


title: Theory & observations of the PWN-SNR complex

abstract: In this work, we study theoretical and observational issues about pulsars (PSRs), pulsar wind nebulae (PWNe) and supernova remnants (SNRs). In particular, the spectral modeling of young PWNe and the X-ray analysis of SNRs with magnetars comparing their characteristics with those remnants surrounding canonical pulsars. The spectra of PWNe range from radio to $\gamma$-rays. They are the largest class of identified Galactic sources in $\gamma$-rays increasing the number from 1 to $\sim$30 during the last years. We have developed a detailed spectral code which reproduces the electromagnetic spectrum of PWNe in free expansion ($t_{age} \lesssim$10 kyr). We shed light and try to understand issues on time evolution of the spectra, the synchrotron self-Compton dominance in the Crab Nebula, the particle dominance in PWNe detected at TeV energies and how physical parameters constrain the detectability of PWNe at TeV. We make a systematic study of all Galactic, TeV-detected, young PWNe which allows to find correlations and trends between parameters. We also discuss about the spectrum of those PWNe not detected at TeV and if models with low magnetized nebulae can explain the lack of detection or, on the contrary, high-magnetization models are more favorable. Regarding the X-ray analysis of SNRs, we use X-ray spectroscopy in SNRs with magnetars to discuss about the formation mechanism of such extremely magnetized PSRs. The alpha-dynamo mechanism proposed in the 1990's produces an energy release that should have influence in the energy of the SN explosion. We extend the work done previously done by \cite{vink06} about the energetics of the SN explosion looking for this energy release and we look for the element ionization and the X-ray luminosity and we compare our results with other SNRs with an associated central source.


title: Three-level spin system under decoherence-minimizing driving fields: Application to nitrogen-vacancy spin dynamics

abstract: Within the framework of a general three-level problem, the dynamics of the nitrogen-vacancy (NV) spin is studied for the case of a special type of external driving consisting of a set of continuous fields with decreasing intensities. Such a set has been proposed for minimizing coherence losses. Each new driving field with smaller intensity is designed to protect against the fluctuations induced by the driving field at the preceding step with larger intensity. We show that indeed this particular type of external driving minimizes the loss of coherence, using purity and entropy as quantifiers for this purpose. As an illustration, we study the coherence loss of an NV spin due to a surrounding spin bath of $^{13}$C nuclei.


title: A new class of supermatrix algebras defined by transitive matrices

abstract: We give a natural definition for the transitivity of a matrix. Using an endomorphism d of a base ring R and a transitive nxn matrix over the center Z(R), we construct the subalgebra M_{n}(R,d,T) of the full nxn matrix algebra M_{n}(R) consisting of the so called nxn supermatrices. If the n-th power of d is the identity and T satisfies some extra conditions, then we exhibit an embedding of R into M_{n}(R,d,T). An other result is that M_{n}(R,d,T) is closed with respect to taking the (pre)adjoint. If R is Lie nilpotent and A is in M_{n}(R,d,T), then the use of the preadjoint and the corresponding determinants and characteristic polynomials yields a Cayley-Hamilton identity for A with right coefficients in the fixed ring Fix(d). The presence of a primitive n-th root of unity and the condition that the n-th power of d is the identity guarantee the right integrality of a Lie nilpotent R over Fix(d). We present essentially new supermatrix algebras over the Grassmann algebra.


title: A discussion on the most general torsion-gravity with electrodynamics for Dirac spinor matter fields

abstract: We consider the most general torsional completion of gravitation together with electrodynamics for the Dirac spinorial material fields, and we show that consistency arguments constrain torsion to be completely antisymmetric and the dynamics to be parity-invariant and described by actions that are either least-order derivative or renormalizable.


title: Computing the coefficients for the power series solution of the Lane-Emden equation with the Python library SymPy

abstract: It is shown how the Python library Sympy can be used to compute symbolically the coefficients of the power series solution of the Lane-Emden equation (LEE). Sympy is an open source Python library for symbolic mathematics. The power series solutions are compared to the numerically computed solutions using matplotlib. The results of a run time measurement of the implemented algorithm are discussed at the end.


title: Attractors of nonlinear Hamiltonian PDEs

abstract: This is a survey of results on long time behavior and attractors for nonlinear Hamiltonian partial differential equations, considering the global attraction to stationary states, stationary orbits, and solitons, the adiabatic effective dynamics of the solitons, and the asymptotic stability of the solitary manifolds. The corresponding numerical results and relations to quantum postulates are considered.


title: Medium-induced jet evolution: wave turbulence and energy loss

abstract: We study the gluon cascade generated via successive medium-induced branchings by an energetic parton propagating through a dense QCD medium. We focus on the high-energy regime where the energy $E$ of the leading particle is much larger than the characteristic medium scale $\omega_c=\hat q L^2/2$, with $\hat q$ the jet quenching parameter and $L$ the distance travelled through the medium. In this regime the leading particle loses only a small fraction $\sim\alpha_s(\omega_c/E)$ of its energy and can be treated as a steady source of radiation for gluons with energies $\omega\le\omega_c$. For this effective problem with a source, we obtain exact analytic solutions for the gluon spectrum and the energy flux. The solutions exhibit wave turbulence: the basic physical process is a continuing fragmentation which is `quasi-democratic' (i.e. quasi-local in energy) and which provides an energy transfer from the source to the medium at a rate (the energy flux $\mathcal{F}$) which is quasi-independent of $\omega$. The locality of the branching process implies a spectrum of the Kolmogorov-Obukhov type, i.e. a power-law spectrum which is a fixed point of the branching process and whose strength is proportional to the energy flux: $D(\omega)\sim\mathcal{F}/\sqrt\omega$ for $\omega\ll\omega_c$. Via this turbulent flow, the gluon cascade loses towards the medium an energy $\Delta E\sim\alpha_s^2\omega_c$, which is independent of the initial energy $E$ of the leading particle and of the details of the thermalization mechanism at the low-energy end of the cascade. This energy is carried away by very soft gluons, which propagate at very large angles with respect to the jet axis. Our predictions for the value of $\Delta E$ and for its angular distribution appear to agree quite well, qualitatively and even semi-quantitatively, with the phenomenology of di-jet asymmetry in nucleus-nucleus collisions at the LHC.


title: Broken axial symmetry as essential feature to predict radiative capture in heavy nuclei

abstract: Cross sections for neutron capture in the range of unresolved resonances are predicted for more than 140 spin-0 target nuclei with A > 50. Allowing the breaking of spherical and axial symmetry in nearly all these nuclei a combined parameterization for both, level density and photon strength is obtained which employs a surprisingly small number of parameters only. The strength functions used are based on a global fit to IVGDR shapes by the sum of three Lorentzians. They are based on theoretical predictions for the A-dependence of pole energies and spreading widths and add up to the TRK sum rule. For the small spins reached by capture resonance spacings are well described by a level density parameter close to the nuclear matter value; a significant collective enhancement is apparent due to the deviation from axial symmetry. Reliable predictions for compound nuclear reactions also outside the valley of stability (as important for nuclear astrophysics and for the transmutation of nuclear waste) are expected to result from the global parameterization presented.


title: From hollow gold nanoparticles to gold nanorings: Morphological tunability of the plasmonic response

abstract: The optical and morphological properties of hollow gold nanoparticles (HGNPs) can be finely modified by employing Poly-L-Lysine hydrobromide (PLL), an homo poly aminoacid of the L-lysine, used as reducer and stabilizer. We investigate locally the plasmonic response of these nanostructures by electron energy loss spectroscopy (EELS). The plasmon excitations are interpreted by discrete dipolar approximation (DDA) simulations. We demonstrate that this controlled top-down morphological modification a fine tuning of the optical response is possible. Unlike the traditional lithographic techniques, this has been achieved in a controlled manner using wet chemistry, enabling the potential use of these nanostructures for a broad range of plasmonic applications, including biomedicine, catalysis and quantum communications.


title: Statics and dynamics of selfish interactions in distributed service systems

abstract: We study a class of games which model the competition among agents to access some service provided by distributed service units and which exhibit congestion and frustration phenomena when service units have limited capacity. We propose a technique, based on the cavity method of statistical physics, to characterize the full spectrum of Nash equilibria of the game. The analysis reveals a large variety of equilibria, with very different statistical properties. Natural selfish dynamics, such as best-response, usually tend to large-utility equilibria, even though those of smaller utility are exponentially more numerous. Interestingly, the latter actually can be reached by selecting the initial conditions of the best-response dynamics close to the saturation limit of the service unit capacities. We also study a more realistic stochastic variant of the game by means of a simple and effective approximation of the average over the random parameters, showing that the properties of the average-case Nash equilibria are qualitatively similar to the deterministic ones.

categories: cs.GT cond-mat.dis-nn cs.MA physics.soc-ph


title: A novel nonlocal potential solver based on nonuniform FFT for efficient simulation of the Davey-Stewartson equations

abstract: We propose an efficient and accurate solver for the nonlocal potential in the Davey-Stewartson equation using nonuniform FFT (NUFFT). A discontinuity in the Fourier transform of the nonlocal potential causes accuracy locking if the potential is solved by standard FFT with periodic boundary conditions on a truncated domain. Using the fact that the discontinuity disappears in polar coordinates, we reformulate the potential integral and split it into high and low frequency parts. The high frequency part can be approximated by the standard FFT method, while the low frequency part is evaluated with a high order Gauss quadrature accelerated by nonuniform FFT. The NUFFT solver has $O(N\log N)$ complexity, where $N$ is the total number of discretization points, and achieves higher accuracy than standard FFT solver, which makes it a good alternative in simulation. Extensive numerical results show the efficiency and accuracy of the proposed method.


title: Actuator and sensor placement in linear advection PDE with building system applications

abstract: We study the problem of actuator and sensor placement in a linear advection partial differential equation (PDE). The problem is motivated by its application to actuator and sensor placement in building systems for the control and detection of a scalar quantity such as temperature and contaminants. We propose a gramian based approach to the problem of actuator and sensor placement. The special structure of the advection PDE is exploited to provide an explicit formula for the controllability and observability gramian in the form of a multiplication operator. The explicit formula for the gramian, as a function of actuator and sensor location, is used to provide test criteria for the suitability of a given sensor and actuator location. Furthermore, the solution obtained using gramian based criteria is interpreted in terms of the flow of the advective vector field. In particular, the almost everywhere stability property of the advective vector field is shown to play a crucial role in deciding the location of actuators and sensors. Simulation results are performed to support the main results of this paper.


title: Interacting particle systems at the edge of multilevel Dyson Brownian motions

abstract: We study the joint asymptotic behavior of spacings between particles at the edge of multilevel Dyson Brownian motions, when the number of levels tends to infinity. Despite the global interactions between particles in multilevel Dyson Brownian motions, we observe a decoupling phenomenon in the limit: the global interactions become negligible and only the local interactions remain. The resulting limiting objects are interacting particle systems which can be described as Brownian versions of certain totally asymmetric exclusion processes. This is the first appearance of a particle system with local interactions in the context of general $\beta$ random matrix models.


title: Consensus in multi-agent systems with non-periodic sampled-data exchange and uncertain network topology

abstract: In this paper consensus in second-order multi-agent systems with a non-periodic sampled-data exchange among agents is investigated. The sampling is random with bounded inter-sampling intervals. It is assumed that each agent has exact knowledge of its own state at any time instant. The considered local interaction rule is PD-type. Sufficient conditions for stability of the consensus protocol to a time-invariant value are derived based on LMIs. Such conditions only require the knowledge of the connectivity of the graph modeling the network topology. Numerical simulations are presented to corroborate the theoretical results.


title: Local strong maximal monotonicity and full stability for parametric variational systems

abstract: The paper introduces and characterizes new notions of Lipschitzian and H\"olderian full stability of solutions to general parametric variational systems described via partial subdifferential and normal cone mappings acting in Hilbert spaces. These notions, postulated certain quantitative properties of single-valued localizations of solution maps, are closely related to local strong maximal monotonicity of associated set-valued mappings. Based on advanced tools of variational analysis and generalized differentiation, we derive verifiable characterizations of the local strong maximal monotonicity and full stability notions under consideration via some positive-definiteness conditions involving second-order constructions of variational analysis. The general results obtained are specified for important classes of variational inequalities and variational conditions in both finite and infinite dimensions.


title: An ontology-based approach to the optimization of non-binary (2,v)-regular LDPC codes

abstract: A non-binary (2,v)-regular LDPC code is defined by a parity-check matrix with column weight 2 and row weight v. In this report, we give an ontology-based approach to the optimization for this class of codes. All possible inter-connected cycle patterns that lead to low symbol-weight codewords are identified to put together the ontology. The optimization goal is to improve the distance property of equivalent binary images. Using the proposed method, the estimation and optimization of bit-distance spectrum becomes easily handleable. Three codes in the CCSDS recommendation are analyzed and several codes with good minimum bit-distance are designed.


title: IST versus PDE, a comparative study

abstract: We survey and compare, mainly in the two-dimensional case, various results obtained by IST and PDE techniques for integrable equations. We also comment on what can be predicted from integrable equations on non integrable ones.


title: Multiwavelength campaign on Mrk 509 XIV. Chandra HETGS spectra

abstract: We present in this paper the results of a 270 ks Chandra HETGS observation in the context of a large multiwavelength campaign on the Seyfert galaxy Mrk 509. The HETGS spectrum allows us to study the high ionisation warm absorber and the Fe-K complex in Mrk 509. We search for variability in the spectral properties of the source with respect to previous observations in this campaign, as well as for evidence of ultra-fast outflow signatures. The Chandra HETGS X-ray spectrum of Mrk 509 was analysed using the SPEX fitting package. We confirm the basic structure of the warm absorber found in the 600 ks XMM-Newton RGS observation observed three years earlier, consisting of five distinct ionisation components in a multikinematic regime. We find little or no variability in the physical properties of the different warm absorber phases with respect to previous observations in this campaign, except for component D2 which has a higher column density at the expense of component C2 at the same outflow velocity (-240 km/s). Contrary to prior reports we find no -700 km/s outflow component. The O VIII absorption line profiles show an average covering factor of 0.81 +/- 0.08 for outflow velocities faster than -100 km/s, similar to those measured in the UV. This supports the idea of a patchy wind. The relative metal abundances in the outflow are close to proto-solar. The narrow component of the Fe Kalpha emission line shows no changes with respect to previous observations which confirms its origin in distant matter. The narrow line has a red wing that can be interpreted to be a weak relativistic emission line. We find no significant evidence of ultra-fast outflows in our new spectrum down to the sensitivity limit of our data.


title: Astrophysical Constraints on the Bumblebee Model

abstract: In this work the Bumblebee model for spontaneous Lorentz symmetry breaking is considered in the context of spherically symmetric astrophysical bodies. A discussion of the modified equations of motion is presented and constraints on the parameters of the model are perturbatively obtained.


title: Optimal investment with bounded above utilities in discrete time markets

abstract: We consider an arbitrage-free, discrete time and frictionless market. We prove that an investor maximising the expected utility of her terminal wealth can always find an optimal investment strategy provided that her dissatisfaction of infinite losses is infinite and her utility function is non-decreasing, continuous and bounded above. The same result is shown for cumulative prospect theory preferences, under additional assumptions.


title: Observed galaxy number counts on the lightcone up to second order: III. Magnification Bias

abstract: We study up to second order the galaxy number over-density that depends on magnification in redshift space on cosmological scales for a concordance model. The result contains all general relativistic effects up to second order which arise from observing on the past light cone, including all redshift and lensing distortions, contributions from velocities, Sachs-Wolfe, integrated SW and time-delay terms. We find several new terms and contributions that could be potentially important for an accurate calculation of the bias on estimates of non-Gaussianity and on precision parameter estimates.


title: Global branching laws by global Okounkov bodies

abstract: Let $G'$ be a complex semisimple group, and let $G \subseteq G'$ be a semisimple subgroup. We show that the branching cone of the pair $(G, G')$, which (asymptotically) parametrizes all pairs $(W, V)$ of irreducible finite-dimensional $G$-representations $W$ which occur as subrepresentations of a finite-dimensional irreducible $G'$-representation $V$, can be identified with the pseudo-effective cone, $\overline{\mbox{Eff}}(Y)$, of some GIT quotient $Y$ of the flag variety of the group $G \times G'$. Moreover, we prove that the quotient $Y$ is a Mori dream space. As a consequence, the global Okounkov body $\Delta(Y)$ of $Y$, with respect to some admissible flag of subvarieties of $Y$, is fibred over the branching cone of $(G, G')$, and the fibre $\Delta(Y)_{(W, V)}$ over a point $(W, V)$ carries information about (the asymptotics of) the multiplicity of $W$ in $V$. Using the global Okounkov body $\Delta(Y)$, we easily derive a multi-dimensional generalization of Okounkov's result about the log-concavity of asymptotic multiplicities.


title: Okounkov bodies for ample line bundles with applications to multiplicities for group representations

abstract: Let $\mathscr{L} \rightarrow X$ be an ample line bundle over a complex normal projective variety $X$. We construct a flag $X_0 \subseteq X_1 \subseteq \cdots \subseteq X_n=X$ of subvarieties for which the associated Okounkov body for $\mathscr{L}$ is a rational polytope. In the case when $X$ is a homogeneous surface, and the pseudoeffective cone of $X$ is rational polyhedral, we also show that the global Okounkov body is a rational polyhedral cone if the flag of subvarieties is suitably chosen. Finally, we provide an application to the asymptotic study of group representations.


title: Short-term Forecasting of Anomalous Load Using Rule-based Triple Seasonal Methods

abstract: Numerous methods have been proposed for forecasting load for normal days. Modeling of anomalous load, however, has often been ignored in the research literature. Occurring on special days, such as public holidays, anomalous load conditions pose considerable modeling challenges due to their infrequent occurrence and significant deviation from normal load. To overcome these limitations, we adopt a rule-based approach, which allows incorporation of prior expert knowledge of load profiles into the statistical model. We use triple seasonal Holt-Winters-Taylor (HWT) exponential smoothing, triple seasonal autoregressive moving average (ARMA), artificial neural networks (ANNs), and triple seasonal intraweek singular value decomposition (SVD) based exponential smoothing. These methods have been shown to be competitive for modeling load for normal days. The methodological contribution of this paper is to demonstrate how these methods can be adapted to model load for special days, when used in conjunction with a rule-based approach. The proposed rule-based method is able to model normal and anomalous load in a unified framework. Using nine years of half-hourly load for Great Britain, we evaluate point forecasts, for lead times from one half-hour up to a day ahead. A combination of two rule-based methods generated the most accurate forecasts.


title: Gabriel morphisms and the computability of Serre quotients with applications to coherent sheaves

abstract: The purpose of this paper is to develop an efficient computational model for Abelian categories of coherent sheaves over certain classes of varieties. These categories are naturally described as Serre quotient categories. Hence, our approach relies on describing general Serre quotient categories in a constructive way which leads to an efficient computer implementation.


title: Rationality of cycles over function fields of $SL_1(A)$-torsors

abstract: In this note we prove a result comparing rationality of algebraic cycles over the function field of a $SL_1(A)$-torsor for a central simple algebra $A$ and over the base field.


title: Algorithms and Polynomiography for Solving Quaternion Quadratic Equations

abstract: Solving a quadratic equation $P(x)=ax^2+bx+c=0$ with real coefficients is known to middle school students. Solving the equation over the quaternions is not straightforward. Huang and So \cite{Huang} give a complete set of formulas, breaking it into several cases depending on the coefficients. From a result of the second author in \cite{kalQ}, zeros of $P(x)$ can be expressed in terms of the zeros of a real quartic equation. This drastically simplifies solving a quadratic equation. Here we also consider solving $P(x)=0$ iteratively via Newton and Halley methods developed in \cite{kalQ}. We prove a property of the Jacobian of Newton and Halley methods and describe several 2D polynomiography based on these methods. The images not only encode the outcome of the iterative process, but by measuring the time taken to render them we find the relative speed of convergence for the methods.


title: Role of the host genetic variability in the influenza-A virus susceptibility; a review

abstract: The aftermath of influenza infection is determined by a complex set of host-pathogen interactions, where genomic variability on both viral and host sides influences the final outcome. Although there exists large body of literature describing influenza virus variability, only a very small fraction covers the issue of host variance. The goal of this review is to explore the variability of host genes responsible for host-pathogen interactions, paying particular attention to genes responsible for the presence of sialylated glycans in the host endothelial membrane, mucus, genes used by viral immune escape mechanisms, and genes particularly expressed after vaccination, since they are more likely to have a direct influence on the infection outcome.


title: Atomic electronic states: the L-S and j-j coupling schemes and their correlation

abstract: In the first part of this paper, we review the assumption of the L-S coupling scheme, with which we derive the electronic states arising from a given atomic configuration. Then, with the aid of the spectral data of Group 15 elements, it becomes clear that the assumption of the L-S coupling scheme is no longer valid as we go farther and farther down the Periodic Table. In the second part, we introduce the j-j coupling scheme, which is seldom covered in standard inorganic chemistry texts, and contrast the assumptions of the two schemes. Next, we use two worked examples to demonstrate the derivation of electronic states with the j-j coupling scheme. Finally, the correlation between the states derived by L-S and j-j schemes is pictorially shown. It is believed a student, by also studying j-j coupling schemes (by no means a difficult task) along with the L-S scheme, will gain a better understanding of the concept of atomic electronic states.


title: ESR studies on the spin-liquid candidate $\kappa$-(BEDT-TTF)$_2$Cu$_2$(CN)$_{3}$: anomalous response below $T=8$ K

abstract: The organic conductor $\kappa$-(BEDT-TTF)$_2$Cu$_2$(CN)$_{3}$ seems to form a quantum spin liquid, although at low temperatures unusual properties are seen in the charge, spin and lattice degrees of freedom. Here we report results of X-band ESR studies of $\kappa$-(BEDT-TTF)$_2$Cu$_2$(CN)$_{3}$ single crystals as a function of temperature and angle. We find indications of two anisotropic relaxation mechanisms at low temperatures and compare them to the spin-liquid behavior observed in other strongly correlated systems. In addition, we can recognize charge inhomogeneities in the copper ions of the anion layer. This disorder might be linked to the dielectric response measured in this compound.


title: Double roots of random Littlewood polynomials

abstract: We consider random polynomials whose coefficients are independent and uniform on {-1,1}. We prove that the probability that such a polynomial of degree n has a double root is o(n^{-2}) when n+1 is not divisible by 4 and asymptotic to $\frac{8\sqrt{3}}{\pi n^2}$ otherwise. This result is a corollary of a more general theorem that we prove concerning random polynomials with independent, identically distributed coefficients having a distribution which is supported on { -1, 0, 1} and whose largest atom is strictly less than 1/\sqrt{3}. In this general case, we prove that the probability of having a double root equals the probability that either -1, 0 or 1 are double roots up to an o(n^{-2}) factor and we find the asymptotics of the latter probability.


title: Gradient-based quantitative image reconstruction in ultrasound-modulated optical tomography: first harmonic measurement type in a linearised diffusion formulation

abstract: Ultrasound-modulated optical tomography is an emerging biomedical imaging modality which uses the spatially localised acoustically-driven modulation of coherent light as a probe of the structure and optical properties of biological tissues. In this work we begin by providing an overview of forward modelling methods, before deriving a linearised diffusion-style model which calculates the first-harmonic modulated flux measured on the boundary of a given domain. We derive and examine the correlation measurement density functions of the model which describe the sensitivity of the modality to perturbations in the optical parameters of interest. Finally, we employ said functions in the development of an adjoint-assisted gradient based image reconstruction method, which ameliorates the computational burden and memory requirements of a traditional Newton-based optimisation approach. We validate our work by performing reconstructions of optical absorption and scattering in two- and three-dimensions using simulated measurements with 1% proportional Gaussian noise, and demonstrate the successful recovery of the parameters to within +/-5% of their true values when the resolution of the ultrasound raster probing the domain is sufficient to delineate perturbing inclusions.


title: Interface fluctuations for deposition on enlarging flat substrates

abstract: We investigate solid-on-solid models that belong to the Kardar-Parisi-Zhang (KPZ) universality class on substrates that expand laterally at a constant rate by duplication of columns. Despite the null global curvature, we show that all investigated models have asymptotic height distributions and spatial covariances in agreement with those expected for the KPZ subclass for curved surfaces. In $1+1$ dimensions, the height distribution and covariance are given by the GUE Tracy-Widom distribution and the Airy$_2$ process, instead of the GOE and Airy$_1$ foreseen for flat interfaces. These results imply that, when the KPZ class splits into the curved and flat subclasses, as conventionally considered, the expanding substrate may play a role equivalent to, or perhaps more important than the global curvature. Moreover, the translational invariance of the interfaces evolving on growing domains allowed us to accurately determine, in $2+1$ dimensions, the analogue of the GUE Tracy-Widom distribution for height distribution and that of the Airy$_2$ process for spatial covariance. Temporal covariance is also calculated and shown to be universal in each dimension and in each of the two subclasses. A logarithmic correction associated to the duplication of column is observed and theoretically elucidated. Finally, crossover between regimes with fixed-size and enlarging substrates is also investigated.


title: A Simplified Hierarchical Dynamic Quantum Secret Sharing Protocol with Added Features

abstract: Generalizing the notion of dynamic quantum secret sharing (DQSS), a simplified protocol for hierarchical dynamic quantum secret sharing (HDQSS) is proposed and it is shown that the protocol can be implemented using any existing protocol of quantum key distribution, quantum key agreement or secure direct quantum communication. The security of this proposed protocol against eavesdropping and collusion attacks is discussed with specific attention towards the issues related to the composability of the subprotocols that constitute the proposed protocol. The security and qubit efficiency of the proposed protocol is also compared with that of other existing protocols of DQSS. Further, it is shown that it is possible to design a semi-quantum protocol of HDQSS and in principle, the protocols of HDQSS can be implemented using any quantum state. It is also noted that the completely orthogonal-state-based realization of HDQSS protocol is possible and that HDQSS can be experimentally realized using a large number of alternative approaches.


title: The Maximal Matching Energy of Tricyclic Graphs

abstract: Gutman and Wagner proposed the concept of the matching energy (ME) and pointed out that the chemical applications of ME go back to the 1970s. Let $G$ be a simple graph of order $n$ and $\mu_1,\mu_2,\ldots,\mu_n$ be the roots of its matching polynomial. The matching energy of $G$ is defined to be the sum of the absolute values of $\mu_{i}\ (i=1,2,\ldots,n)$. Gutman and Cvetkoi\'c determined the tricyclic graphs on $n$ vertices with maximal number of matchings by a computer search for small values of $n$ and by an induction argument for the rest. Based on this result, in this paper, we characterize the graphs with the maximal value of matching energy among all tricyclic graphs, and completely determine the tricyclic graphs with the maximal matching energy. We prove our result by using Coulson-type integral formula of matching energy, which is similar as the method to comparing the energies of two quasi-order incomparable graphs.


title: Matching Energy of Unicyclic and Bicyclic Graphs with a Given Diameter

abstract: Gutman and Wagner proposed the concept of matching energy (ME) and pointed out that the chemical applications of ME go back to the 1970s. Let $G$ be a simple graph of order $n$ and $\mu_1,\mu_2,\ldots,\mu_n$ be the roots of its matching polynomial. The matching energy of $G$ is defined to be the sum of the absolute values of $\mu_{i}\ (i=1,2,\ldots,n)$. In this paper, we characterize the graphs with minimal matching energy among all unicyclic and bicyclic graphs with a given diameter $d$.


title: Strain shielding from mechanically-activated covalent bond formation in nanoindentation

abstract: Mechanical failure of an ideal crystal is dictated either by an elastic instability or a soft-mode instability. We show that the ideal strength measurement of graphene based on nano-indentation experiments \cite{lee2008measurement, lee2013high}, however, indicates an anomaly: the inferred strain beneath the diamond indenter at the failure load is anomalously large compared to the fracture strain predicted by soft-mode analysis or acoustic analysis. Here we present a systematic investigation - based on multi-scale modeling combining the results of continuum, atomistic, and quantum calculations; and analysis of experiments - that identifies the operative mechanism responsible for the anomalous difference between the fracture strains. We suggest that a strain-shielding effect due to mechanically-activated covalent bond formation at graphene-indenter interface is responsible for this anomaly. Using Finite Elements Analysis (FEA) and MD simulations of the nanoindentation experiments, we explicitly show that the bonded interaction at the graphene-indenter interface substantially disperses (shields) the strain beneath the indenter, preventing the intensification of strain. Our calculations indicate that the extent of strain shielding depends upon the hydrogen coverage at the indenter surface; and at optimal hydrogen coverage, the strain-shielding effect can delay the onset of fracture to the experimentally-observed indentation load and depth.


title: The Buchberger resolution

abstract: We define the Buchberger resolution, which is a graded free resolution of a monomial ideal in a polynomial ring. Its construction uses a generalization of the Buchberger graph and encodes much of the combinatorics of the Buchberger algorithm. The Buchberger resolution is a cellular resolution that coincides with the Scarf resolution for generic monomial ideals, which is the case when it is minimal. The simplicial complex underlying the Buchberger resolution is of interest for its own sake and its combinatorics is not fully understood. We close with a conjecture on the clique complex of the Buchberger graph.


title: Recommendation Subgraphs for Web Discovery

abstract: Recommendations are central to the utility of many websites including YouTube, Quora as well as popular e-commerce stores. Such sites typically contain a set of recommendations on every product page that enables visitors to easily navigate the website. Choosing an appropriate set of recommendations at each page is one of the key features of backend engines that have been deployed at several e-commerce sites. Specifically at BloomReach, an engine consisting of several independent components analyzes and optimizes its clients' websites. This paper focuses on the structure optimizer component which improves the website navigation experience that enables the discovery of novel content. We begin by formalizing the concept of recommendations used for discovery. We formulate this as a natural graph optimization problem which in its simplest case, reduces to a bipartite matching problem. In practice, solving these matching problems requires superlinear time and is not scalable. Also, implementing simple algorithms is critical in practice because they are significantly easier to maintain in production. This motivated us to analyze three methods for solving the problem in increasing order of sophistication: a sampling algorithm, a greedy algorithm and a more involved partitioning based algorithm. We first theoretically analyze the performance of these three methods on random graph models characterizing when each method will yield a solution of sufficient quality and the parameter ranges when more sophistication is needed. We complement this by providing an empirical analysis of these algorithms on simulated and real-world production data. Our results confirm that it is not always necessary to implement complicated algorithms in the real-world and that very good practical results can be obtained by using heuristics that are backed by the confidence of concrete theoretical guarantees.


title: Isolation of Flow and Nonflow Correlations by Two- and Four-Particle Cumulant Measurements of Azimuthal Harmonics in $\sqrt{s_{_{\rm NN}}} =$ 200 GeV Au+Au Collisions

abstract: A data-driven method was applied to measurements of Au+Au collisions at $\sqrt{s_{_{\rm NN}}} =$ 200 GeV made with the STAR detector at RHIC to isolate pseudorapidity distance $\Delta\eta$-dependent and $\Delta\eta$-independent correlations by using two- and four-particle azimuthal cumulant measurements. We identified a component of the correlation that is $\Delta\eta$-independent, which is likely dominated by anisotropic flow and flow fluctuations. It was also found to be independent of $\eta$ within the measured range of pseudorapidity $|\eta|<1$. The relative flow fluctuation was found to be $34\% \pm 2\% (stat.) \pm 3\% (sys.)$ for particles of transverse momentum $p_{T}$ less than $2$ GeV/$c$. The $\Delta\eta$-dependent part may be attributed to nonflow correlations, and is found to be $5\% \pm 2\% (sys.)$ relative to the flow of the measured second harmonic cumulant at $|\Delta\eta| > 0.7$.


title: Equivariant Chern classes in Hopf cyclic cohomology

abstract: We present a geometric approach, in the spirit of the Chern-Weil theory, for constructing cocycles representing the classes of the Hopf cyclic cohomology of the Hopf algebra H(n) relative to GL(n, R). This provides an explicit description of the universal Hopf cyclic Chern classes, which complements our earlier geometric realization of the Hopf cyclic characteristic classes of foliations.


title: Global Convergence of Online Limited Memory BFGS

abstract: Global convergence of an online (stochastic) limited memory version of the Broyden-Fletcher- Goldfarb-Shanno (BFGS) quasi-Newton method for solving optimization problems with stochastic objectives that arise in large scale machine learning is established. Lower and upper bounds on the Hessian eigenvalues of the sample functions are shown to suffice to guarantee that the curvature approximation matrices have bounded determinants and traces, which, in turn, permits establishing convergence to optimal arguments with probability 1. Numerical experiments on support vector machines with synthetic data showcase reductions in convergence time relative to stochastic gradient descent algorithms as well as reductions in storage and computation relative to other online quasi-Newton methods. Experimental evaluation on a search engine advertising problem corroborates that these advantages also manifest in practical applications.


title: Correlation between amygdala BOLD activity and frontal EEG asymmetry during real-time fMRI neurofeedback training in patients with depression

abstract: Real-time fMRI neurofeedback (rtfMRI-nf) is an emerging approach for studies and novel treatments of major depressive disorder (MDD). EEG performed simultaneously with an rtfMRI-nf procedure allows an independent evaluation of rtfMRI-nf brain modulation effects. Frontal EEG asymmetry in the alpha band is a widely used measure of emotion and motivation that shows profound changes in depression. However, it has never been directly related to simultaneously acquired fMRI data. We report the first study investigating electrophysiological correlates of the rtfMRI-nf procedure, by combining rtfMRI-nf with simultaneous and passive EEG recordings. In this pilot study, MDD patients in the experimental group (n=13) learned to upregulate BOLD activity of the left amygdala using an rtfMRI-nf during a happy emotion induction task. MDD patients in the control group (n=11) were provided with a sham rtfMRI-nf. Correlations between frontal EEG asymmetry in the upper alpha band and BOLD activity across the brain were examined. Average individual changes in frontal EEG asymmetry during the rtfMRI-nf task for the experimental group showed a significant positive correlation with the MDD patients' depression severity ratings, consistent with an inverse correlation between the depression severity and frontal EEG asymmetry at rest. Temporal correlations between frontal EEG asymmetry and BOLD activity were significantly enhanced, during the rtfMRI-nf task, for the amygdala and many regions associated with emotion regulation. Our findings demonstrate an important link between amygdala BOLD activity and frontal EEG asymmetry. Our EEG asymmetry results suggest that the rtfMRI-nf training targeting the amygdala is beneficial to MDD patients, and that alpha-asymmetry EEG-nf would be compatible with the amygdala rtfMRI-nf. Combination of the two could enhance emotion regulation training and benefit MDD patients.


title: Untersuchungen zur Implementierung von Bildverarbeitungsalgorithmen mittels pulsgekoppelter neuronaler Netze

abstract: This thesis deals with the study of image processing algorithms which can be implemented by pulse-coupled neural nets. The inspiration for this choice is taken from biological image processing, which achieves with little computational effort in highly parallel processes image analysis tasks such as object recognition, image segmentation, velocity and distance estimation, etc. Conventional, serially implemented algorithms either cannot realize those tasks at all or will expend significantly more effort. Because the first stages of the visual system comprise a sensor interface, they are comparatively accessible with respect to defining their transfer or processing function. Some of those processing functions or principles are to be used in hardware implementations, with the focus on duplicating especially the highly parallel processing.


title: Quantum Belief Propagation Algorithm versus Suzuki-Trotter approach in the one-dimensional Heisenberg chains

abstract: Quantum systems are the future candidates for computers and information processing devices. Information about quantum states and processes may be incomplete and scattered in these systems. We use a quantum version of Belief Propagation(BP) Algorithm to integrate the distributed information. In this algorithm the distributed information, which is in the form of density matrix, can be approximated to local structures. The validity of this algorithm is measured in comparison with Suzuki-Trotter(ST) method, using simulated information. ST in 3-body Heisenberg example gives a more accurate answer, however Quantum Belief Propagation (QBP) runs faster based on complexity. In order to develop it in the future, we should be looking for ways to increase the accuracy of QBP.


title: Magnetochiral symmetry breaking in a M\"obius ring

abstract: We show that the interaction of the magnetic subsystem of a curved magnet with the magnet curvature results in coupling of a topologically nontrivial magnetization pattern and topology of the object. The mechanism of this coupling is explored and illustrated by an example of ferromagnetic M\"obius ring, where a topologically induced domain wall appears as a ground state in case of strong easy-normal anisotropy. For the M\"obius geometry the curvilinear form of the exchange interaction produces an additional effective Dzyaloshinskii-like term which leads to the coupling of the magnetochirality of the domain wall and chirality of the M\"obius ring. Two types of domain walls are found, transversal and longitudinal, which are oriented across and along the M\"obius ring, respectively. In both cases the effect of magnetochirality symmetry breaking is established. The dependence of the ground state of the M\"obius ring on its geometrical parameters and on the value of the easy-normal anisotropy is explored numerically.


title: Depth image hand tracking from an overhead perspective using partially labeled, unbalanced data: Development and real-world testing

abstract: We present the development and evaluation of a hand tracking algorithm based on single depth images captured from an overhead perspective for use in the COACH prompting system. We train a random decision forest body part classifier using approximately 5,000 manually labeled, unbalanced, partially labeled training images. The classifier represents a random subset of pixels in each depth image with a learned probability density function across all trained body parts. A local mode-find approach is used to search for clusters present in the underlying feature space sampled by the classified pixels. In each frame, body part positions are chosen as the mode with the highest confidence. User hand positions are translated into hand washing task actions based on proximity to environmental objects. We validate the performance of the classifier and task action proposals on a large set of approximately 24,000 manually labeled images.


title: Likelihood-based tree reconstruction on a concatenation of alignments can be positively misleading

abstract: The reconstruction of a species tree from genomic data faces a double hurdle. First, the (gene) tree describing the evolution of each gene may differ from the species tree, for instance, due to incomplete lineage sorting. Second, the aligned genetic sequences at the leaves of each gene tree provide merely an imperfect estimate of the topology of the gene tree. In this note, we demonstrate formally that a basic statistical problem arises if one tries to avoid accounting for these two processes and analyses the genetic data directly via a concatenation approach. More precisely, we show that, under the multi-species coalescent with a standard site substitution model, maximum likelihood estimation on sequence data that has been concatenated across genes and performed under the incorrect assumption that all sites have evolved independently and identically on a fixed tree is a statistically inconsistent estimator of the species tree. Our results provide a formal justification of simulation results described of Kubatko and Degnan (2007) and others, and complements recent theoretical results by DeGorgio and Degnan (2010) and Chifman and Kubtako (2014).


title: Gaudin subalgebras and wonderful models

abstract: Gaudin hamiltonians form families of r-dimensional abelian Lie subalgebras of the holonomy Lie algebra of the arrangement of reflection hyperplanes of a Coxeter group of rank r. We consider the set of principal Gaudin subalgebras, which is the closure in the appropriate Grassmannian of the set of spans of Gaudin hamiltonians. We show that principal Gaudin subalgebras form a smooth projective variety isomorphic to the De Concini-Procesi compactification of the projectivized complement of the arrangement of reflection hyperplanes.


title: Itinerant effects and enhanced magnetic interactions in Bi-based multilayer cuprates

abstract: The cuprate high temperature superconductors exhibit a pronounced trend in which the superconducting transition temperature, $T_{\rm c}$, increases with the number of CuO$_2$ planes, $n$, in the crystal structure. We compare the magnetic excitation spectrum of Bi$_{2+x}$Sr$_{2-x}$CuO$_{6+\delta}$ (Bi-2201) and Bi$_2$Sr$_2$Ca$_2$Cu$_3$O$_{10 + \delta}$ (Bi-2223), with $n=1$ and $n=3$ respectively, using Cu $L_3$-edge resonant inelastic x-ray scattering (RIXS). Near the anti-nodal zone boundary we find the paramagnon energy in Bi-2223 is substantially higher than that in Bi-2201, indicating that multilayer cuprates host stronger effective magnetic exchange interactions, providing a possible explanation for the $T_{\rm c}$ vs.\ $n$ scaling. In contrast, the nodal direction exhibits very strongly damped, almost non-dispersive excitations. We argue that this implies that the magnetism in the doped cuprates is partially itinerant in nature.


title: Another characterization of tilted algebras

abstract: We give a new characterization of tilted algebras by the existence of certain special subquivers in their Auslander-Reiten quiver. This result includes the existent characterizations of this kind and yields a way to obtain more tilted quotient algebras from a given algebra.


title: Zermelo Navigation in the Quantum Brachistochrone

abstract: We analyse the optimal times for implementing unitary quantum gates in a constrained finite dimensional controlled quantum system. The family of constraints studied is that the permitted set of (time dependent) Hamiltonians is the unit ball of a norm induced by an inner product on su(n). We also consider a generalisation of this to arbitrary norms. We construct a Randers metric, by applying a theorem of Shen on Zermelo navigation, the geodesics of which are the time optimal trajectories compatible with the prescribed constraint. We determine all geodesics and the corresponding time optimal Hamiltonian for a specific constraint on the control i.e. k (Tr(Hc(t)^2) = 1 for any given value of k > 0. Some of the results of Carlini et. al. are re-derived using alternative methods. A first order system of differential equations for the optimal Hamiltonian is obtained and shown to be of the form of the Euler Poincare equations. We illustrate that this method can form a methodology for determining which physical substrates are effective at supporting the implementation of fast quantum computation.


title: A One-Line Proof of the Fundamental Theorem of Algebra with Newton's Method as a Consequence

abstract: Many proofs of the fundamental theorem of algebra rely on the fact that the minimum of the modulus of a complex polynomial over the complex plane is attained at some complex number. The proof then follows by arguing the minimum value is zero. This can be done by proving that at any complex number that is not a zero of the polynomial we can exhibit a direction of descent for the modulus. In this note we present a very short and simple proof of the existence of such descent direction. In particular, our descent direction gives rise to Newton's method for solving a polynomial equation via modulus minimization and also makes the iterates definable at any critical point.


title: Erd\H{o}s-Ko-Rado for Perfect Matchings

abstract: A perfect matching of a complete graph $K_{2n}$ is a 1-regular subgraph that contains all the vertices. Two perfect matchings intersect if they share an edge. It is known that if $\mathcal{F}$ is family of intersecting perfect matchings of $K_{2n}$, then $|\mathcal{F}| \leq (2(n-1) - 1)!!$ and if equality holds, then $\mathcal{F} = \mathcal{F}_{ij}$ where $ \mathcal{F}_{ij}$ is the family of all perfect matchings of $K_{2n}$ that contain some fixed edge $ij$. We give a short algebraic proof of this result, resolving a question of Godsil and Meagher. Along the way, we show that if a family $\mathcal{F}$ is non-Hamiltonian, that is, $m \cup m' \not \cong C_{2n}$ for any $m,m' \in \mathcal{F}$, then $|\mathcal{F}| \leq (2(n-1) - 1)!!$ and this bound is met with equality if and only if $\mathcal{F} = \mathcal{F}_{ij}$. Our results make ample use of a somewhat understudied symmetric commutative association scheme arising from the Gelfand pair $(S_{2n},S_2 \wr S_n)$. We give an exposition of a few new interesting objects that live in this scheme as they pertain to our results.


title: Comment on chiral symmetry restoration at finite density in large-$N_c$ QCD

abstract: In the article "On chiral symmetry restoration at finite density in large-$N_c$ QCD" by Adhikari, Cohen, Ayyagari and Strother [Phys. Rev. C 83, 065201 (2011)] the description of dense nuclear matter by means of Skyrmions in hyperspherical unit cells is severely criticized. We point out that this criticism is based on invalid assumptions and therefore unwarranted.


title: Exactness of Second Order Ordinary Differential Equations and Integrating Factors

abstract: A new problem is studied, the concept of exactness of a second order nonlinear ordinary differential equations is established. A method is constructed to reduce this class into a first order equations. If the second order equation is not exact we introduce, under certain conditions, an integrating factor that transform it to an exact one.


title: Nonequilibrium stationary states of 3D self-gravitating systems

abstract: Three dimensional self-gravitating systems do not evolve to thermodynamic equilibrium, but become trapped in nonequilibrium quasistationary states. In this Letter we present a theory which allows us to a priori predict the particle distribution in a final quasistationary state to which a self-gravitating system will evolve from an initial condition which is isotropic in particle velocities and satisfies a virial constraint 2K=-U, where K is the total kinetic energy and U is the potential energy of the system.


title: Quantum Key Distribution without sending a Quantum Signal

abstract: Quantum Key Distribution is a quantum communication technique in which random numbers are encoded on quantum systems, usually photons, and sent from one party, Alice, to another, Bob. Using the data sent via the quantum signals, supplemented by classical communication, it is possible for Alice and Bob to share an unconditionally secure secret key. This is not possible if only classical signals are sent. Whilst this last statement is a long standing result from quantum information theory it turns out only to be true in a non-relativistic setting. If relativistic quantum field theory is considered we show it is possible to distribute an unconditionally secure secret key without sending a quantum signal, instead harnessing the intrinsic entanglement between different regions of space time. The protocol is practical in free space given horizon technology and might be testable in principle in the near term using microwave technology.


title: Surface salinity fields in the Arctic Ocean and statistical approaches to predicting anomalies and patterns

abstract: Significant salinity anomalies have been observed in the Arctic Ocean surface layer during the last decade. Using gridded data of winter salinity in the upper 50 m layer of the Arctic Ocean for the period 1950-1993 and 2007-2012, we investigated the inter-annual variability of the salinity fields, attempted to identify patterns and anomalies, and developed a statistical model for the prediction of surface layer salinity. The statistical model is based on linear regression equations linking the principal components with environmental factors, such as atmospheric circulation, river runoff, ice processes, and water exchange with neighboring oceans. Using this model, we obtained prognostic fields of the surface layer salinity for the winter period 2013-2014. The prognostic fields demonstrated the same tendencies of surface layer freshening that were observed previously. A phase portrait analysis involving the first two principal components exhibits a dramatic shift in behavior of the 2007-2012 data in comparison to earlier observations.


title: Focal points and sup-norms of eigenfunctions on analytic Riemannian manifolds II: the two-dimensional case

abstract: In the recent work arXiv:1311.3999, the authors proved that real analytic manifolds $(M, g)$ with maximal eigenfunction growth must have a self-focal point p whose first return map has an invariant L1 measure on $S^*_p M$. In this addendum we add a purely dynamical argument on circle maps to improve the conclusion to: all geodesics from p are smoothly closed.


title: Rethinking the Contention Resolution Mechanism in WiMAX Networks using Lattice Correlators for Improved Smart Grid Communication Performance

abstract: By using experimental measurements of smart grid applications compliant with IEC 61850 in trace-driven WiMAX simulations, we show that the WiMAX MAC protocol efficiency decreases as a function of the number of stations. To avoid this shortcoming, we propose and analyze a novel WiMAX MAC protocol for smart grid applications, which uses lattice correlators to improve the throughput-delay performance significantly. For the considered configurations, the obtained maximum throughput of the proposed MAC protocol outperforms the current WiMAX MAC protocol by up to 41%.


title: Convex Clustering: An Attractive Alternative to Hierarchical Clustering

abstract: The primary goal in cluster analysis is to discover natural groupings of objects. The field of cluster analysis is crowded with diverse methods that make special assumptions about data and address different scientific aims. Despite its shortcomings in accuracy, hierarchical clustering is the dominant clustering method in bioinformatics. Biologists find the trees constructed by hierarchical clustering visually appealing and in tune with their evolutionary perspective. Hierarchical clustering operates on multiple scales simultaneously. This is essential, for instance, in transcriptome data where one may be interested in making qualitative inferences about how lower-order relationships like gene modules lead to higher-order relationships like pathways or biological processes. The recently developed method of convex clustering preserves the visual appeal of hierarchical clustering while ameliorating its propensity to make false inferences in the presence of outliers and noise. The current paper exploits the proximal distance principle to construct a novel algorithm for solving the convex clustering problem. The solution paths generated by convex clustering reveal relationships between clusters that are hidden by static methods such as k-means clustering. Our convex clustering software separates parameters, accommodates missing data, and supports prior information on relationships. The software is implemented on ATI and nVidia graphics processing units (GPUs) for maximal speed. Several biological examples illustrate the strengths of convex clustering and the ability of the proximal distance algorithm to handle high-dimensional problems.


title: Crossover between Fermi liquid and non-Fermi liquid behavior in the non-centrosymmetric compound Yb$_2$Ni$_{12}$P$_7$

abstract: A crossover from a non-Fermi liquid to a Fermi liquid phase in Yb$_2$Ni$_{12}$P$_7$ is observed by analyzing electrical resistivity $\rho(T)$, magnetic susceptibility $\chi(T)$, specific heat $C(T)$, and thermoelectric power $S(T)$ measurements. The electronic contribution to specific heat, $C_{e}(T)$, behaves as $C_{e}(T)/T \sim -\ln(T)$ for 5 K $< T <$ 15 K, which is consistent with non-Fermi liquid behavior. Below $T \sim$ 4 K, the upturn in $C_{e}(T)/T$ begins to saturate, suggesting that the system crosses over into a Fermi-liquid ground state. This is consistent with robust $\rho(T) - \rho_0 = AT^2$ behavior below $T \sim$ 4 K, with the power-law exponent becoming sub-quadratic for $T >$ 4 K. A crossover between Fermi-liquid and non-Fermi liquid behavior suggests that Yb$_2$Ni$_{12}$P$_7$ is in close proximity to a quantum critical point, in agreement with results from recent measurements of this compound under applied pressure.


title: On the density of certain languages with $p^2$ letters

abstract: The sequence $(x_n)_{n\in\mathbb N} = (2,5,15,51,187,\dots)$ given by the rule $x_n=(2^n+1)(2^{n-1}+1)/3$ appears in several seemingly unrelated areas of mathematics. For example, $x_n$ is the density of a language of words of length $n$ with four different letters. It is also the cardinality of the quotient of $(\mathbb Z_2\times \mathbb Z_2)^n$ under the left action of the special linear group $\mathrm{SL}(2,\mathbb Z)$. In this paper we show how these two interpretations of $x_n$ are related to each other. More generally, for prime numbers $p$ we show a correspondence between a quotient of $(\mathbb Z_p\times\mathbb Z_p)^n$ and a language with $p^2$ letters and words of length $n$.


title: Quasi-Symmetries of Determinantal Point Processes

abstract: The main result of this paper is that determinantal point processes on the real line corresponding to projection operators with integrable kernels are quasi-invariant, in the continuous case, under the group of diffeomorphisms with compact support (Theorem 1.4); in the discrete case, under the group of all finite permutations of the phase space (Theorem 1.6). The Radon-Nikodym derivative is computed explicitly and is given by a regularized multiplicative functional. Theorem 1.4 applies, in particular, to the sine-process, as well as to determinantal point processes with the Bessel and the Airy kernels; Theorem 1.6 to the discrete sine-process and the Gamma kernel process. The paper answers a question of Grigori Olshanski.


title: Hyperbolic three manifolds of bounded volume and trace field degree II

abstract: In this paper, we prove the Bounded Height Conjecture which the author formulated in [2]. As a corollary, it follows that there are only a finite number of hyperbolic three manifolds of bounded volume and trace field degree.


title: Uniform disconnectedness and Quasi-Assouad Dimension

abstract: The uniform disconnectedness is an important invariant property under bi-Lipschitz mapping, and the Assouad dimension $\dim _{A}X<1$ implies the uniform disconnectedness of $X$. According to quasi-Lipschitz mapping, we introduce the quasi-Assouad dimension $\dim _{qA}$ such that $\dim _{qA}X<1$ implies its quasi uniform disconnectedness. We obtain $\overline{\dim } _{B}X\leq \dim _{qA}X\leq \dim _{A}X$ and compute the quasi-Assouad dimension of Moran set.


title: Untersuchungen zur Modellierung und Schaltungsrealisierung von synaptischer Plastizitaet

abstract: This manuscript deals with the analysis and VLSI implementation of adaptive information processing derived from biological measurements. Specifically, models for short term plasticity, long term plasticity and metaplasticity are derived from biological measurements and implemented in CMOS circuits.


title: The Mabuchi Geometry of Finite Energy Classes

abstract: We introduce different Finsler metrics on the space of smooth K\"ahler potentials that will induce a natural geometry on various finite energy classes $\mathcal E_{\tilde \chi}(X,\omega)$. Motivated by questions raised by R. Berman, V. Guedj and Y. Rubinstein, we characterize the underlying topology of these spaces in terms of convergence in energy and give applications of our results to existence of K\"ahler-Einstein metrics on Fano manifolds.


title: An NLP Assistant for Clide

abstract: This report describes an NLP assistant for the collaborative development environment Clide, that supports the development of NLP applications by providing easy access to some common NLP data structures. The assistant visualizes text fragments and their dependencies by displaying the semantic graph of a sentence, the coreference chain of a paragraph and mined triples that are extracted from a paragraph's semantic graphs and linked using its coreference chain. Using this information and a logic programming library, we create an NLP database which is used by a series of queries to mine the triples. The algorithm is tested by translating a natural language text describing a graph to an actual graph that is shown as an annotation in the text editor.


title: Bilateral zeta functions associated with the multiple sine functions

abstract: We introduce two types bilateral zeta functions, which are related to the primitive and normalized multiple sine functions respectively. Further, we establish their main properties, that is, Fourier expansions, analytic continuations, differential and difference equations, special values. By applying these results, we obtain not only some generalization of the primitive and normalized multiple sine functions but also simple construction of the multiple sine function theory.


title: Quantum Monte Carlo study of the $S_4$ symmetric microscopic model for iron-based superconductors

abstract: The $S_4$ symmetric microscopic model with two iso-spin components has been studied via constrained-path quantum Monte Carlo simulation. Our results demonstrate a stable $(\pi,0)$ or $(0,\pi)$ magnetic order which is significantly enhanced on increasing both the Coulomb repulsion $U$ and Hund's coupling strength $J$. Also, our simulation indicates that the magnetic order tends to be in an orthomagnetic one, in which the nearest-neighbour magnetic moment are orthogonal to each other, rather than in a collinear antiferromagnetic state. Interestingly, when the system is doped away from half filling, the magnetic order is obviously elevated in the low doping density, and then significantly suppressed when more electrons are introduced. Meanwhile, we find that an $A_{1g}$ $s_{\pm}$-wave pairing dominates all the singlet nearest-neighbour pairings, and is significantly enhanced via electron doping.


title: Electrically Tunable Optical Absorption in a Graphene-based Salisbury Screen

abstract: We demonstrate a graphene-based Salisbury screen consisting of a single layer of graphene placed in close proximity to a gold back reflector. The light absorption in the screen can be actively tuned by electrically gating the carrier density in the graphene layer with an ionic liquid/gel. The screen was designed to achieve maximum absorption at a target wavelength of 3.2 micrometer by using a 600 nm-thick, non-absorbing silica spacer layer. Spectroscopic reflectance measurements were performed in-situ as a function of gate bias. The changes in the reflectance spectra were analyzed using a Fresnel based transfer matrix model in which graphene was treated as an infinitesimally thin sheet with conductivity given by the Kubo formula. Temporal coupled mode theory was employed to analyze and intuitively understand the observed absorption changes in the Salisbury screen. We achieved ~ 6 % change in the optical absorption of graphene by tuning the applied gate bias from 0.8 V to 2.6 V, where 0.8 V corresponds to graphene's charge neutrality point.


title: Davenport constant for semigroups II

abstract: Let $\mathcal{S}$ be a finite commutative semigroup. The Davenport constant of $\mathcal{S}$, denoted ${\rm D}(\mathcal{S})$, is defined to be the least positive integer $\ell$ such that every sequence $T$ of elements in $\mathcal{S}$ of length at least $\ell$ contains a proper subsequence $T'$ ($T'\neq T$) with the sum of all terms from $T'$ equaling the sum of all terms from $T$. Let $q>2$ be a prime power, and let $\F_q[x]$ be the ring of polynomials over the finite field $\F_q$. Let $R$ be a quotient ring of $\F_q[x]$ with $0\neq R\neq \F_q[x]$. We prove that $${\rm D}(\mathcal{S}_R)={\rm D}(U(\mathcal{S}_R)),$$ where $\mathcal{S}_R$ denotes the multiplicative semigroup of the ring $R$, and $U(\mathcal{S}_R)$ denotes the group of units in $\mathcal{S}_R$.


title: Optimal prediction for positive self-similar Markov processes

abstract: This paper addresses the question of predicting when a positive self-similar Markov process X attains its pathwise global supremum or infimum before hitting zero for the first time (if it does at all). This problem has been studied in Glover et al. (2013) under the assumption that X is a positive transient diffusion. We extend their result to the class of positive self-similar Markov processes by establishing a link to Baurdoux and van Schaik (2013), where the same question is studied for a Levy process drifting to minus infinity. The connection to Baurdoux and van Schaik (2013) relies on the so-called Lamperti transformation which links the class of positive self-similar Markov processes with that of Levy processes. Our approach will reveal that the results in Glover et al. (2013) for Bessel processes can also be seen as a consequence of self-similarity.


title: Conjectured bounds for the sum of squares of positive eigenvalues of a graph

abstract: A well known upper bound for the spectral radius of a graph, due to Hong, is that $\mu_1^2 \le 2m - n + 1$. It is conjectured that for connected graphs $n - 1 \le s^+ \le 2m - n + 1$, where $s^+$ denotes the sum of the squares of the positive eigenvalues. The conjecture is proved for various classes of graphs, including bipartite, regular, complete $q$-partite, hyper-energetic, and barbell graphs. Various searches have found no counter-examples. The paper concludes with a brief discussion of the apparent difficulties of proving the conjecture in general.


title: Multiscale statistical testing for connectome-wide association studies in fMRI

abstract: Alterations in brain connectivity have been associated with a variety of clinical disorders using functional magnetic resonance imaging (fMRI). We investigated empirically how the number of brain parcels (or scale) impacted the results of a mass univariate general linear model (GLM) on connectomes. The brain parcels used as nodes in the connectome analysis were functionnally defined by a group cluster analysis. We first validated that a classic Benjamini-Hochberg procedure with parametric GLM tests did control appropriately the false-discovery rate (FDR) at a given scale. We then observed on realistic simulations that there was no substantial inflation of the FDR across scales, as long as the FDR was controlled independently within each scale, and the presence of true associations could be established using an omnibus permutation test combining all scales. Second, we observed both on simulations and on three real resting-state fMRI datasets (schizophrenia, congenital blindness, motor practice) that the rate of discovery varied markedly as a function of scales, and was relatively higher for low scales, below 25. Despite the differences in discovery rate, the statistical maps derived at different scales were generally very consistent in the three real datasets. Some seeds still showed effects better observed around 50, illustrating the potential benefits of multiscale analysis. On real data, the statistical maps agreed well with the existing literature. Overall, our results support that the multiscale GLM connectome analysis with FDR is statistically valid and can capture biologically meaningful effects in a variety of experimental conditions.


title: History-free Collision Response for Deformable Surfaces

abstract: Continuous collision detection (CCD) and response methods are widely adopted in dynamics simulation of deformable models. They are history-based, as their success is strictly based on an assumption of a collision-free state at the start of each time interval. On the other hand, in many applications surfaces have normals defined to designate their orientation (i.e. front- and back-face), yet CCD methods are totally blind to such orientation identification (thus are orientation-free). We notice that if such information is utilized, many penetrations can be untangled. In this paper we present a history-free method for separation of two penetrating meshes, where at least one of them has clarified surface orientation. This method first computes all edge-face (E-F) intersections with discrete collision detection (DCD), and then builds a number of penetration stencils. On response, the stencil vertices are relocated into a penetration-free state, via a global displacement minimizer. Our method is very effective for handling penetration between two meshes, being it an initial configuration or in the middle of physics simulation. The major limitation is that it is not applicable to self-collision within one mesh at the time being.


title: Kinetic mixing and symmetry breaking dependent interactions of the dark photon

abstract: We examine spontaneous symmetry breaking of a renormalisable U(1) x U(1) gauge theory coupled to fermions when kinetic mixing is present. We do not assume that the kinetic mixing parameter is small. A rotation plus scaling is used to remove the mixing and put the gauge kinetic terms in the canonical form. Fermion currents are also rotated in a non-orthogonal way by this basis transformation. Through suitable redefinitions the interaction is cast into a diagonal form. This framework, where mixing is absent, is used for subsequent analysis. The symmetry breaking determines the fermionic current which couples to the massless gauge boson. The strength of this coupling as well as the couplings of the massive gauge boson are extracted. This formulation is used to consider a gauged model for dark matter by identifying the massless gauge boson with the photon and the massive state to its dark counterpart. Matching the coupling of the residual symmetry with that of the photon sets a lower bound on the kinetic mixing parameter. We present analytical formulae of the couplings of the dark photon in this model and indicate some physics consequences.


title: Generalized Dirichlet to Neumann Maps for Linear Dispersive Equations on Half-Line

abstract: A large class of initial-boundary value problems of linear evolution partial differential equations formulated on the half-line is analyzed via the unified transform method. In particular, explicit formulae are presented for the generalized Dirichlet to Neumann maps. Namely, the determination of the unknown boundary values in terms of an essential set of given data.


title: Theory of Winds from Hot, Luminous Massive Stars

abstract: The high luminosities of massive stars drive strong stellar winds, through line scattering of the star's continuum radiation. This paper reviews the dynamics of such line driving, building first upon the standard CAK model for steady winds, and deriving the associated analytic scalings for the mass loss rate and wind velocity law. It next summarizes the origin and nature of the strong Line Deshadowing Instability (LDI) intrinsic to such line-driving, including also the role of a diffuse-line-drag effect that stabilizes the wind base, and then describes how both instability and drag are incorporated in the Smooth Source Function (SSF) method for time-dependent simulations of the nonlinear evolution of the resulting wind structure. The review concludes with a discussion of the effect of the resulting extensive structure in temperature, density and velocity for interpreting observational diagnostics. In addition to the usual clumping effect on density-squared diagnostics, the spatial porosity of optically thick clumps can reduce single-density continuum absorption, and a kind of velocity porosity, or vorocity, can reduce the absorption strength of spectral lines. An overall goal is to illuminate the rich physics of radiative driving and the challenges that lie ahead in developing dynamical models for the often complex structure and variability of hot-star winds.


title: Stochastic fields with paths in arbitrary rearrangement invariant spaces

abstract: We obtain sufficient conditions for belonging of almost all paths of a random process to some fixed rearrangement invariant (r.i.) Banach functional space, and to satisfying the Central Limit Theorem (CLT) in this space. We describe also some possible applications.


title: Simulations of Kinetic Electrostatic Electron Nonlinear (KEEN) Waves with Variable Velocity Resolution Grids and High-Order Time-Splitting

abstract: KEEN waves are nonlinear, non-stationary, self-organized asymptotic states in Vlasov plasmas outside the scope or purview of linear theory constructs such as electron plasma waves or ion acoustic waves. Nonlinear stationary mode theories such as those leading to BGK modes also do not apply. The range in velocity that is strongly perturbed by KEEN waves depends on the amplitude and duration of the ponderomotive force used to drive them. Smaller amplitude drives create highly localized structures attempting to coalesce into KEEN waves. These cases have much more chaotic and intricate time histories than strongly driven ones. The narrow range in which one must maintain adequate velocity resolution in the weakly driven cases challenges xed grid numerical schemes. What is missing there is the capability of resolving locally in velocity while maintaining a coarse grid outside the highly perturbed region of phase space. We here report on a new Semi-Lagrangian Vlasov-Poisson solver based on conservative non-uniform cubic splines in velocity that tackles this problem head on. An additional feature of our approach is the use of a new high-order time-splitting scheme which allows much longer simulations per computational e ort. This is needed for low amplitude runs which take a long time to set up KEEN waves, if they are able to do so at all. The new code's performance is compared to uniform grid simulations and the advantages quanti ed. The birth pains associated with KEEN waves which are weakly driven is captured in these simulations. These techniques allow the e cient simulation of KEEN waves in multiple dimensions which will be tackled next as well as generalizations to Vlasov-Maxwell codes which are essential to understanding the impact of KEEN waves in practice.


title: On the multiplier rules

abstract: We establish new results of first-order necessary conditions of optimality for finite-dimensional problems with inequality constraints and for problems with equality and inequality constraints, in the form of John's theorem and in the form of Karush-Kuhn-Tucker's theorem. In comparison with existing results we weaken assumptions of continuity and of differentiability.


title: Introducing Molly: Distributed Memory Parallelization with LLVM

abstract: Programming for distributed memory machines has always been a tedious task, but necessary because compilers have not been sufficiently able to optimize for such machines themselves. Molly is an extension to the LLVM compiler toolchain that is able to distribute and reorganize workload and data if the program is organized in statically determined loop control-flows. These are represented as polyhedral integer-point sets that allow program transformations applied on them. Memory distribution and layout can be declared by the programmer as needed and the necessary asynchronous MPI communication is generated automatically. The primary motivation is to run Lattice QCD simulations on IBM Blue Gene/Q supercomputers, but since the implementation is not yet completed, this paper shows the capabilities on Conway's Game of Life.


title: Perfrewrite -- Program Complexity Analysis via Source Code Instrumentation

abstract: Most program profiling methods output the execution time of one specific program execution, but not its computational complexity class in terms of the big-O notation. Perfrewrite is a tool based on LLVM's Clang compiler to rewrite a program such that it tracks semantic information while the program executes and uses it to guess memory usage, communication and computational complexity. While source code instrumentation is a standard technique for profiling, using it for deriving formulas is an uncommon approach.


title: On the asymptotics of random forests

abstract: The last decade has witnessed a growing interest in random forest models which are recognized to exhibit good practical performance, especially in high-dimensional settings. On the theoretical side, however, their predictive power remains largely unexplained, thereby creating a gap between theory and practice. The aim of this paper is twofold. Firstly, we provide theoretical guarantees to link finite forests used in practice (with a finite number M of trees) to their asymptotic counterparts. Using empirical process theory, we prove a uniform central limit theorem for a large class of random forest estimates, which holds in particular for Breiman's original forests. Secondly, we show that infinite forest consistency implies finite forest consistency and thus, we state the consistency of several infinite forests. In particular, we prove that q quantile forests---close in spirit to Breiman's forests but easier to study---are able to combine inconsistent trees to obtain a final consistent prediction, thus highlighting the benefits of random forests compared to single trees.


title: Probing the Accretion Scheme of the Dipping X-ray Binary 4U 1915-05 with Suzaku

abstract: The dipping low-mass X-ray binary 4U 1915-05 was observed by Suzaku on 2007 November 8 for a net exposure of 39 ksec. It was detected by the XIS with a 0.8-10 keV signal rate of 9.84+\-0.01 cts/s per camera, and HXD-PIN with a 12-45 keV signal rate of 0.29+/-0.01 cts/s. After removing the periodic dips and an X-ray burst, the 0.8 - 45 keV continuum was successfully described by an optically thick disk emission with an inner-disk temperature ~ 0.7 keV and a neutron-star blackbody emission with a temperature ~ 1.3 keV, on condition that the blackbody component, or possibly the disk emission too, is significantly Comptonized. This successful modeling is consistent with 4U 1915-05 being in a high-soft state in this observation, and implies that its broadband spectrum can be interpreted in the same scheme as for many non-dipping Low-mass X-ray binaries in the soft state. Its bolometric luminosity (~ 0.02 times the Eddington limit) is relatively low for the soft state, but within a tolerance, if considering the distance and inclination uncertainties. As a high-inclination binary, this source exhibited stronger Comptonization effect, with a larger Comptonizing y-parameter, compared to low and medium inclination binaries. This suggests that the Comptonizing coronae of these objects in the soft state is in an oblate (rather than spherical) shape, extending along the accretion disk plane, because the y-parameter would not depend on the inclination if the corona were spherical.