title: Intuition, iteration, induction

abstract: In Mathematical Thought and Its Objects, Charles Parsons argues that our knowledge of the iterability of functions on the natural numbers and of the validity of complete induction is not intuitive knowledge; Brouwer disagrees on both counts. I will compare Parsons' argument with Brouwer's and defend the latter. I will not argue that Parsons is wrong once his own conception of intuition is granted, as I do not think that that is the case. But I will try to make two points: (1) Using elements from Husserl and from Brouwer, Brouwer's claims can be justified in more detail than he has done; (2) There are certain elements in Parsons' discussion that, when developed further, would lead to Brouwer's notion thus analysed, or at least something relevantly similar to it. (This version contains a postscript of May, 2015.)


title: Evolution of initially contracting Bianchi Class A models in the presence of an ultra-stiff anisotropic pressure fluid

abstract: We study the behaviour of Bianchi class A universes containing an ultra-stiff isotropic ghost field and a fluid with anisotropic pressures which is also ultra-stiff on the average. This allows us to investigate whether cyclic universe scenarios, like the ekpyrotic model, do indeed lead to isotropisation on approach to a singularity (or bounce) in the presence of dominant ultra-stiff pressure anisotropies. We specialise to consider the closed Bianchi type IX universe and show that when the anisotropic pressures are stiffer on average than any isotropic ultra-stiff fluid then, if they dominate on approach to the singularity, it will be anisotropic. We include an isotropic ultra-stiff ghost fluid with negative energy density in order to create a cosmological bounce at finite volume in the absence of the anisotropic fluid. When the dominant anisotropic fluid is present it leads to an anisotropic cosmological singularity rather than an isotropic bounce. The inclusion of anisotropic stresses generated by collisionless particles in an anisotropically expanding universe is therefore essential for a full analysis of the consequences of a cosmological bounce or singularity in cyclic universes.


title: Low-mass right-handed gauge bosons from minimal grand unified theories

abstract: Prediction of low-mass $W_R$ and $Z_R$ gauge bosons in popular grand unified theories has been the subject of considerable attention over the last three decades. In this work we show that when gravity induced corrections due to dim.5 operator are included the minimal symmetry breaking chain of $SO(10)$ and $E_6$ GUTs can yield $W_R^{\pm}$ and $Z_R$ bosons with masses in the range $(3-10)$ TeV which are accessible to experimental tests at the Large Hordan Collider. The RH neutrinos turn out to be heavy pseudo-Dirac fermions. The model can fit all fermion masses and manifest in rich structure of lepton flavor violation while proton life time is predicted to be much longer than the accessible limit of Super-Kamiokande or planned Hyper-Kamiokande collaborations.


title: Observation of mixed singlet-triplet Rb$_2$ Rydberg molecules

abstract: We present high-resolution spectroscopy of Rb$_\text{2}$ ultralong-range Rydberg molecules bound by mixed singlet-triplet electron-neutral atom scattering. The mixing of the scattering channels is a consequence of the hyperfine interaction in the ground-state atom, as predicted recently by Anderson et al. \cite{Anderson2014b}. Our experimental data enables the determination of the effective zero-energy singlet $s$-wave scattering length for Rb. We show that an external magnetic field can tune the contributions of the singlet and the triplet scattering channels and therefore the binding energies of the observed molecules. This mixing of molecular states via the magnetic field results in observed shifts of the molecular line which differ from the Zeeman shift of the asymptotic atomic states. Finally, we calculate molecular potentials using a full diagonalization approach including the $p$-wave contribution and all orders in the relative momentum $k$, and compare the obtained molecular binding energies to the experimental data.


title: Intensity-only optical compressive imaging using a multiply scattering material and a double phase retrieval approach

abstract: In this paper, the problem of compressive imaging is addressed using natural randomization by means of a multiply scattering medium. To utilize the medium in this way, its corresponding transmission matrix must be estimated. To calibrate the imager, we use a digital micromirror device (DMD) as a simple, cheap, and high-resolution binary intensity modulator. We propose a phase retrieval algorithm which is well adapted to intensity-only measurements on the camera, and to the input binary intensity patterns, both to estimate the complex transmission matrix as well as image reconstruction. We demonstrate promising experimental results for the proposed algorithm using the MNIST dataset of handwritten digits as example images.


title: Dynamical Domain Wall and Localization

abstract: Based on the previous works (arXiv:1202.5375 and 1402.1346), we investigate the localization of the fields on the dynamical domain wall, where the four dimensional FRW universe is realized on the domain wall in the five dimensional space-time. Especially we show that the chiral spinor can localize on the domain wall, which has not been succeeded in the past works as the seminal work in arXiv:0810.3746.


title: SWS promoting MHb-IPN-MRN circuit opposes the theta promoting circuit, active wake and REM sleep

abstract: Hippocampus was by multisynaptic axonal tracts linked to medial habenula in my DTI study. These tracts linked hippocampus to septum, and amygdala to BNST. Septal and BNST axons passed by anteromedial thalamic nucleus to MHb, which projected to pineal gland. Question is what is the MHb doing and why it receives information from hippocampus via septum? This study explores the MHb connectivity, predicts its functional role and how it is linked to memory. Combination of known findings about the septum and MHb connectivity and function led to this circuit-based idea that posterior septum activates MHb, MHb activates IPN,and IPN stimulates MRN and its serotonin release. Proposed idea is that this MHb-IPN-MRN circuit promotes slow wave sleep (SWS), high serotonin and low acetylcholine state. My prediction is that this SWS-promoting circuit reciprocally suppresses the theta oscillations promoting circuit, linked to high acetylcholine levels in brain, and formed by supramamillary area projections to the medial septum (MS) that induces theta rhythm in hippocampus and other theta-coupled regions. The MHb-IPN-MRN pathway likely inhibits, possibly reciprocally,some regions that have stimulating input to the theta-generating SUM and MS, such as the wake promoting nucleus incertus, posterior and lateral hypothalamus and LDT, but also the REM sleep inducing neurons in LDT and PnO. Thus, proposed SWS- promoting circuit attenuates the output of theta-promoting regions, both the active wake-on and REM-on regions. The theta rhythm in wake state is linked to recording and binding information with their spatio-temporal and relational context by hippocampus, while the SWS supports rest and replay of hippocampaly stored information and their cortical reactivations, e. g. in retrosplenial cortex linked to autobiographic memory or in prefrontal cortex that can combine information from any source.


title: Stationary shapes of deformable particles moving at low Reynolds numbers

abstract: Lecture Notes of the Summer School ``Microswimmers -- From Single Particle Motion to Collective Behaviour'', organised by the DFG Priority Programme SPP 1726 (Forschungszentrum J{\"{u}}lich, 2015).


title: Describing the response of saturated SiPMs

abstract: We have developed a function which describes SiPM response in both small signal and highly saturated regimes. The function includes the reactivation of SiPM pixels during a single input light pulse, and results in an approximately linear increase of SiPM response in the highly saturated regime, as observed in real SiPMs. This article shows that the function can accurately describe the measured response of real SiPM devices over a wide range of signal intensities.


title: Blocking estimators and inference under the Neyman-Rubin model

abstract: We derive the variances of estimators for sample average treatment effects under the Neyman-Rubin potential outcomes model for arbitrary blocking assignments and an arbitrary number of treatments.


title: Delocalization of two interacting particles in the two-dimensional Harper model

abstract: We study the problem of two interacting particles in a two-dimensional quasiperiodic potential of the Harper model. We consider an amplitude of the quasiperiodic potential such that in absence of interactions all eigenstates are exponentially localized while the two interacting particles are delocalized showing anomalous subdiffusive spreading over the lattice with the spreading exponent $b \approx 0.5$ instead of a usual diffusion with $b=1$. This spreading is stronger than in the case of a correlated disorder potential with a one particle localization length as for the quasiperiodic potential. At the same time we do not find signatures of ballistic FIKS pairs existing for two interacting particles in the one-dimensional Harper model.


title: Evidence of Excitonic Optical Tamm States using Molecular Materials

abstract: We report the first experimental observation of an Excitonic Optical Tamm State supported at the interface between a periodic multilayer dielectric structure and an organic dye-doped polymer layer. The existence of such states is enabled by the metal-like optical properties of the excitonic layer based on aggregated dye molecules. Experimentally determined dispersion curves, together with simulated data, including field profiles, allow us to identify the nature of these new modes. Our results demonstrate the potential of organic excitonic materials as a powerful means to control light at the nanoscale, offering the prospect of a new alternative type of nanophotonics based on molecular materials.

categories: physics.optics cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.soft


title: Fundamental Speed Limits to the Generation of Quantumness

abstract: Quantum physics dictates fundamental speed limits during time evolution. We present a quantum speed limit governing the generation of nonclassicality and the mutual incompatibility of two states connected by time evolution. This result is used to characterize the timescale required to generate a given amount of quantumness under an arbitrary physical process. The bound is shown to be tight under pure dephasing dynamics. More generally, our analysis reveals the dependence on the initial state and non-Markovian effects.


title: Drag reduction in numerical two-phase Taylor-Couette turbulence using an Euler-Lagrange approach

abstract: Two-phase turbulent Taylor-Couette (TC) flow is simulated using an Euler-Lagrange approach to study the effects of a secondary phase dispersed into a turbulent carrier phase (here bubbles dispersed into water). The dynamics of the carrier phase is computed using Direct Numerical Simulations (DNS) in an Eulerian framework, while the bubbles are tracked in a Lagrangian manner by modelling the effective drag, lift, added mass and buoyancy force acting on them. Two-way coupling is implemented between the dispersed phase and the carrier phase which allows for momentum exchange among both phases and to study the effect of the dispersed phase on the carrier phase dynamics. The radius ratio of the TC setup is fixed to $\eta=0.833$, and a maximum inner cylinder Reynolds number of $Re_i=8000$ is reached. We vary the Froude number ($Fr$), which is the ratio of the centripetal to the gravitational acceleration of the dispersed phase and study its effect on the net torque required to drive the TC system. For the two-phase TC system, we observe drag reduction, i.e., the torque required to drive the inner cylinder is less compared to that of the single phase system. The net drag reduction decreases with increasing Reynolds number $Re_i$, which is consistent with previous experimental findings (Murai et al. 2005, 2008). The drag reduction is strongly related to the Froude number: for fixed Reynolds number we observe higher drag reduction when $Fr < 1$ than for with$ Fr > 1$. This buoyancy effect is more prominent in low $Re_i$ systems and decreases with increasing Reynolds number $Re_i$. We trace the drag reduction back to the weakening of the angular momentum carrying Taylor rolls by the rising bubbles.


title: Stability of rotating magnetized jets in the solar atmosphere. I. Kelvin-Helmholtz instability

abstract: Observations show various jets in the solar atmosphere with significant rotational motions, which may undergo instabilities leading to heat ambient plasma. We study the Kelvin-Helmholtz (KH) instability of twisted and rotating jets caused by the velocity jumps near the jet surface. We derive a dispersion equation with appropriate boundary condition for total pressure (including centrifugal force of tube rotation), which governs the dynamics of incompressible jets. Then, we obtain analytical instability criteria of Kelvin-Helmholtz instability in various cases, which were verified by numerical solutions to the dispersion equation. We find that twisted and rotating jets are unstable to KH instability when the kinetic energy of rotation is more than the magnetic energy of the twist. Our analysis shows that the azimuthal magnetic field of 1-5 G can stabilize observed rotations in spicule/macrospicules and X-ray/EUV jets. On the other hand, non-twisted jets are always unstable to KH instability. In this case, the instability growth time is several seconds for spicule/macrospicules and few minutes (or less) for EUV/X-ray jets. We also find that standing kink and torsional Alfven waves are always unstable near the antinodes due to the jump of azimuthal velocity at the surface, while the propagating waves are generally stable. KH vortices may lead to enhanced turbulence development and heating of surrounding plasma, therefore rotating jets may provide energy for chromospheric and coronal heating.


title: On type IIn/Ia-CSM supernovae as exemplified by SN 2012ca

abstract: We present the complete set of ultra-violet, optical and near-infrared photometry and spectroscopy for SN 2012ca, covering the period from 6~days prior to maximum light, until 531 days after maximum. The spectroscopic time series for SN 2012ca is essentially unchanged over 1.5 years, and appear to be dominated at all epochs by signatures of interaction with a dense circumstellar medium rather than the underlying supernova (SN). At late phases, we see a near infrared excess in flux which is possibly associated with dust formation, although without any signs of accompanying line shifts. SN 2012ca is a member of the set of type of the ambiguous IIn/Ia-CSM SNe, the nature of which have been debated extensively in the literature. The two leading scenarios are either a type Ia SN exploding within a dense CSM from a non-degenerate, evolved companion, or a core-collapse SN from a massive star. While some members of the population have been unequivocally associated with type Ia SNe, in other cases the association is less certain. While it is possible that SN 2012ca does arise from a thermonuclear SN, this would require a relatively high (between 20 and 70 per cent) efficiency in converting kinetic energy to optical luminosity, and a massive (~2.3-2.6 Msun) circumstellar medium. On the basis of energetics, and the results of simple modelling, we suggest that SN 2012ca is more likely associated with a core-collapse SN. This would imply that the observed set of similar SNe to SN 2012ca is in fact originated by two populations, and while these are drawn from physically distinct channels, they can have observationally similar properties.


title: Spectral decomposition of the Mueller matrix: A geometrical model

abstract: An arbitrary Mueller matrix can be decomposed into a sum of up to four deterministic Mueller-Jones matrices, with strengths given by the eigenvalues of an associated Hermitian matrix. A geometrical representation of the eigenvalues in terms of the matrix invariants is presented.


title: Quasi periodic oscillations in active galactic nuclei

abstract: Quasi-periodic oscillations (QPOs) are coherent peaks of variability power observed in the X-ray power spectra (PSDs) of stellar mass X-ray binaries (XRBs). A scale invariance of the accretion process implies they should be present in the active galactic nuclei. The first robust detection was a $\sim 1$ hr periodicity in the Seyfert galaxy RE J1034+396 from a $\sim 90$ ks XMM-Newton observation, however, subsequent observations failed to detect the QPO in the $0.3-10.0$ keV band. In this talk we present the recent detection of the $\sim 1$ hr periodicity in the $1.0-4.0$ keV band of 4 further low-flux/spectrally-harder observations of RE J1034+396 (see Alston et al 2014). We also present recent work on the discovery of a QPO in the Seyfert galaxy, MS 22549-3712, which again is only detected in energy bands associated with the primary power-law continuum emission (Alston et al 2015). We conclude these features are most likely analogous to the high-frequency QPOs observed in XRBs. In both sources, we also see evidence for X-ray reverberation at the QPO frequency, where soft X-ray bands and Iron K$\alpha$ emission lag the primary X-ray continuum. These time delays may provide another diagnostic for understanding the underlying QPO mechanism observed in accreting black holes.


title: First Principles Investigation of Hydrogen Physical Adsorption on Graphynes' layers

abstract: Graphynes are 2D porous structures deriving from graphene featuring triangular and regularly distributed subnanometer pores, which may be exploited to host small gaseous species. First principles adsorption energies of molecular hydrogen (H2) on graphene, graphdiyne and graphtriyne molecular prototypes are obtained at the MP2C level of theory. First, a single layer is investigated and it is found that graphynes are more suited than graphene for H2 physical adsorption since they provide larger binding energies at equilibrium distances much closer to the 2D plane. In particular, for graphtriyne a flat minimum located right in the geometric center of the pore is identified. A novel graphite composed of graphtriyne stacked sheets is then proposed and an estimation of its 3D arrangement is obtained at the DFT level of theory. In contrast to pristine graphite this new carbon material allow both H2 intercalation and out-of-plane diffusion by exploiting the larger volume provided by its nanopores. Related H2 binding energies for intercalation and in-pore adsorption are around 0.1 eV and they could lead to high storage capacities. The proposed carbon-based layered material could represent a safer and potentially cheaper alternative for hydrogen on-board storage than conventional solutions based on cryogenic liquefaction and/or high compression.


title: RAID: A Relation-Augmented Image Descriptor

abstract: As humans, we regularly interpret images based on the relations between image regions. For example, a person riding object X, or a plank bridging two objects. Current methods provide limited support to search for images based on such relations. We present RAID, a relation-augmented image descriptor that supports queries based on inter-region relations. The key idea of our descriptor is to capture the spatial distribution of simple point-to-region relationships to describe more complex relationships between two image regions. We evaluate the proposed descriptor by querying into a large subset of the Microsoft COCO database and successfully extract nontrivial images demonstrating complex inter-region relations, which are easily missed or erroneously classified by existing methods.


title: A Piecewise Deterministic Markov Toy Model for Traffic/Maintenance and Associated Hamilton-Jacobi Integrodifferential Systems on Networks

abstract: We study optimal control problems in infinite horizon when the dynamics belong to a specific class of piecewise deterministic Markov processes constrained to star-shaped networks (inspired by traffic models). We adapt the results in [H. M. Soner. Optimal control with state-space constraint. II. SIAM J. Control Optim., 24(6):1110.1122, 1986] to prove the regularity of the value function and the dynamic programming principle. Extending the networks and Krylov's ''shaking the coefficients'' method, we prove that the value function can be seen as the solution to a linearized optimization problem set on a convenient set of probability measures. The approach relies entirely on viscosity arguments. As a by-product, the dual formulation guarantees that the value function is the pointwise supremum over regular subsolutions of the associated Hamilton-Jacobi integrodifferential system. This ensures that the value function satisfies Perron's preconization for the (unique) candidate to viscosity solution. Finally, we prove that the same kind of linearization can be obtained by combining linearization for classical (unconstrained) problems and cost penalization. The latter method works for very general near-viable systems (possibly without further controllability) and discontinuous costs.


title: Linking the Structural Properties of Galaxies and their Star Formation Histories with STAGES

abstract: We study the links between star formation history and structure for a large mass-selected galaxy sample at 0.05 < z_phot < 0.30. The galaxies inhabit a very broad range of environments, from cluster cores to the field. Using HST images, we quantify their structure following Hoyos et al. (2012), and divide them into disturbed and undisturbed. We also visually identify mergers. Additionally, we provide a quantitative measure of the degree of disturbance for each galaxy ("roughness"). The majority of elliptical and lenticular galaxies have relaxed structure, showing no signs of ongoing star formation. Structurally-disturbed galaxies, which tend to avoid the lowest-density regions, have higher star-formation activity and younger stellar populations than undisturbed systems. Cluster spirals with reduced/quenched star formation have somewhat less disturbed morphologies than spirals with "normal" star-formation activity, suggesting that these "passive" spirals have started their morphological transformation into S0s. Visually identified mergers and galaxies not identified as mergers but with similar roughness have similar specific star formation rates and stellar ages. The degree of enhanced star formation is thus linked to the degree of structural disturbance, regardless of whether it is caused by major mergers or not. This suggests that merging galaxies are not special in terms of their higher-than-normal star-formation activity. Any physical process that produces "roughness", or regions of enhanced luminosity density, will increase the star-formation activity in a galaxy with similar efficiency. An alternative explanation is that star formation episodes increase the galaxies' roughness similarly, regardless of whether they are merger-induced or not.


title: Centrality metrics and localization in core-periphery networks

abstract: Two concepts of centrality have been defined in complex networks. The first considers the centrality of a node and many different metrics for it has been defined (e.g. eigenvector centrality, PageRank, non-backtracking centrality, etc). The second is related to a large scale organization of the network, the core-periphery structure, composed by a dense core plus an outlying and loosely-connected periphery. In this paper we investigate the relation between these two concepts. We consider networks generated via the Stochastic Block Model, or its degree corrected version, with a strong core-periphery structure and we investigate the centrality properties of the core nodes and the ability of several centrality metrics to identify them. We find that the three measures with the best performance are marginals obtained with belief propagation, PageRank, and degree centrality, while non-backtracking and eigenvector centrality (or MINRES}, showed to be equivalent to the latter in the large network limit) perform worse in the investigated networks.


title: On the escape rate of unique beta-expansions

abstract: Let $1<\beta \leq 2$. It is well-known that the set of points in $% [0,1/(\beta -1)]$ having unique $\beta $-expansion, in other words, those points whose orbits under greedy $\beta $-transformation escape a hole depending on $\beta $, is of zero Lebesgue measure. The corresponding escape rate is investigated in this paper. A formula which links the Hausdorff dimension of univoque set and escape rate is established in this study. Then we also proved that such rate forms a devil's staircase function with respect to $\beta $.


title: Illustration of iterative linear solver behavior on simple 1D and 2D problems

abstract: In geometry processing, numerical optimization methods often involve solving sparse linear systems of equations. These linear systems have a structure that strongly resembles to adjacency graphs of the underlying mesh. We observe how classic linear solvers behave on this specific type of problems. For the sake of simplicity, we minimise either the squared gradient or the squared Laplacian, evaluated by finite differences on a regular 1D or 2D grid. We observed the evolution of the solution for both energies, in 1D and 2D, and with different solvers: Jacobi, Gauss-Seidel, SSOR (Symmetric successive over-relaxation) and CG (conjugate gradient [She94]). Plotting results at different iterations allows to have an intuition of the behavior of these classic solvers.


title: Amplitude equations for weakly nonlinear surface waves in variational problems

abstract: Among hyperbolic Initial Boundary Value Problems (IBVP), those coming from a variational principle 'generically' admit linear surface waves, as was shown by Serre [J. Funct. Anal. 2006]. At the weakly nonlinear level, the behavior of surface waves is expected to be governed by an amplitude equation that can be derived by means of a formal asymptotic expansion. Amplitude equations for weakly nonlinear surface waves were introduced by Lardner [Int. J. Engng Sci. 1983], Parker and co-workers [J. Elasticity 1985] in the framework of elasticity, and by Hunter [Contemp. Math. 1989] for abstract hyperbolic problems. They consist of nonlocal evolution equations involving a complicated, bilinear Fourier multiplier in the direction of propagation along the boundary. It was shown by the authors in an earlier work [Arch. Ration. Mech. Anal. 2012] that this multiplier, or kernel, inherits some algebraic properties from the original IBVP. These properties are crucial for the (local) well-posedness of the amplitude equation, as shown together with Tzvetkov [Adv. Math., 2011]. Properties of amplitude equations are revisited here in a somehow simpler way, for surface waves in a variational setting. Applications include various physical models, from elasticity of course to the director-field system for liquid crystals introduced by Saxton [Contemp. Math. 1989] and studied by Austria and Hunter [Commun. Inf. Syst. 2013]. Similar properties are eventually shown for the amplitude equation associated with surface waves at reversible phase boundaries in compressible fluids, thus completing a work initiated by Benzoni-Gavage and Rosini [Comput. Math. Appl. 2009].


title: Modelling the synchrotron light curves in recurrent novae V745 Scorpii and RS Ophiuchi

abstract: In this paper, we present the synchrotron light curve at 610 MHz from the recurrent nova V745 Sco following its outburst on 6 February 2014. The system has been detected and periodically monitored with the Giant Metrewave Radio Telescope (GMRT) since 9 February 2014 as part of the Galactic Nova with GMRT (GNovaG) project. The light curves are well fit by a model of synchrotron emitting region obscured by foreground thermal gas which eventually becomes optically thin to the low GMRT frequencies. We present the model fit to the 2014 data on V745 Sco and discuss it alongwith the model fit to the 1.4 GHz data of the recurrent nova RS Ophiuchi following its outburst in 1985.


title: Range and critical generations of a random walk on Galton-Watson trees

abstract: In this paper we consider a random walk in random environment on a tree and focus on the boundary case for the underlying branching potential. We study the range $R\_n$ of this walk up to time $n$ and obtain its correct asymptotic in probability which is of order $n/\log n$. Thisresult is a consequence of the asymptotical behavior of the number of visited sites at generations of order $(\log n)^2$,which turn out to be the most visited generations. Our proof which involves a quenched analysisgives a description of the typical environments responsible for the behavior of $R\_n$.


title: On The Evolution Of User Support Topics in Computational Science and Engineering Software

abstract: We investigate ten years of user support emails in the large-scale solver library PETSc in order to identify changes in user requests. For this purpose we assign each email thread to one or several categories describing the type of support request. We find that despite several changes in hardware architecture as well programming models, the relative share of emails for the individual categories does not show a notable change over time. This is particularly remarkable as the total communication volume has increased four-fold in the considered time frame, indicating a considerable growth of the user base. Our data also demonstrates that user support cannot be substituted with what is often referred to as 'better documentation' and that the involvement of core developers in user support is essential.


title: From a Kac-like particle system to the Landau equation for hard potentials and Maxwell molecules

abstract: We prove a quantitative result of convergence of a conservative stochastic particle system to the solution of the homogeneous Landau equation for hard potentials. There are two main difficulties: (i) the known stability results for this class of Landau equations concern regular solutions and seem difficult to extend to study the rate of convergence of some empirical measures; (ii) the conservativeness of the particle system is an obstacle for (approximate) independence. To overcome (i), we prove a new stability result for the Landau equation for hard potentials concerning very general measure solutions. Due to (ii), we have to couple, our particle system with some non independent nonlinear processes, of which the law solves, in some sense, the Landau equation. We then prove that these nonlinear processes are not so far from being independent. Using finally some ideas of Rousset [25], we show that in the case of Maxwell molecules, the convergence of the particle system is uniform in time.


title: The semi-classical limit of large fermionic systems

abstract: We study a system of $N$ fermions in the regime where the intensity of the interaction scales as $1/N$ and with an effective semi-classical parameter $\hbar=N^{-1/d}$ where $d$ is the space dimension. For a large class of interaction potentials and of external electromagnetic fields, we prove the convergence to the Thomas-Fermi minimizers in the limit $N\to\infty$. The limit is expressed using many-particle coherent states and Wigner functions. The method of proof is based on a fermionic de Finetti-Hewitt-Savage theorem in phase space and on a careful analysis of the possible lack of compactness at infinity.


title: An empirical study of Unfairness and Oscillation in ETSI DCC

abstract: -Performance of Vehicular Adhoc Networks (VANETs) in high node density situation has long been a major field of studies. Particular attention has been paid to the frequent exchange of Cooperative Awareness Messages (CAMs) on which many road safety applications rely. In the present paper, se focus on the European Telecommunications Standard Institute (ETSI) Decentralized Congestion Control (DCC) mechanism, particularly on the evaluation of its facility layers component when applied in the context of dense networks. For this purpose, a set of simulations has been conducted over several scenarios, considering rural highway and urban mobility in order to investigate unfairness and oscillation issues, and analyze the triggering factors. The experimental results show that the latest technical specification of the ETSI DCC presents a significant enhancement in terms of fairness. In contrast, the stability criterion leaves room for improvement as channel load measurement presents (i) considerable fluctuations when only the facility layer control is applied and (i.i) severe state oscillation when different DCC control methods are combined.


title: Elementary Superconductivity in Nonlinear Electrodynamics Coupled to Gravity

abstract: Source-free equations of nonlinear electrodynamics minimally coupled to gravity admit regular axially symmetric asymptotically Kerr-Newman solutions which describe charged rotating black holes and electromagnetic spinning solitons (lumps). Asymptotic analysis of solutions shows, for both black holes and solitons, the existence of de Sitter vacuum interior which has the properties of a perfect conductor and ideal diamagnetic and displays superconducting behaviour which can be responsible for practically unlimited life time of an object. Superconducting current flows on the equatorial ring replacing the Kerr ring singularity of the Kerr-Newman geometry. Interior de Sitter vacuum supplies the electron with the finite positive electromagnetic mass related the interior de Sitter vacuum of the electroweak scale and to breaking of space-time symmetry, which allows to explain the mass-square differences for neutrino and the appearance of the minimal length scale in the annihilation reaction $e^{+}e^{-}\rightarrow\gamma\gamma(\gamma)$.


title: Liaison Linkages

abstract: The complete classification of hexapods - also known as Stewart Gough platforms - of mobility one is still open. To tackle this problem, we can associate to each hexapod of mobility one an algebraic curve, called the configuration curve. In this paper we establish an upper bound for the degree of this curve, assuming the hexapod is general enough. Moreover, we provide a construction of hexapods with curves of maximal degree, which is based on liaison, a technique used in the theory of algebraic curves.


title: The age of Taurus - environmental effects on disc lifetimes

abstract: Using semi-empirical isochrones, we find the age of the Taurus star-forming region to be 3-4 Myr. Comparing the disc fraction in Taurus to young massive clusters suggests discs survive longer in this low density environment. We also present a method of photometrically de-reddening young stars using $iZJH$ data.


title: Ray mapping approach for the efficient design of continuous freeform surfaces

abstract: The efficient design of continuous freeform surfaces, which transform a given source into an arbitrary target intensity, remains a challenging problem. A popular approach are ray-mapping methods, where first a ray mapping between the intensities is calculated and in a subsequent step the surface is constructed. The challenging part hereby is the to find an integrable mapping ensuring a continuous surface. Based on the law of reflection/refraction and the well-known integrability condition, we derive a general condition for the surface and ray mapping for a collimated input beam. It is shown that in a small-angle approximation a proper mapping can be calculated via optimal mass transport. We show that the surface can be constructed by solving a linear advection equation with appropriate boundary conditions. The results imply that the optimal mass transport mapping is approximately integrable over a wide range of distances between the freeform and the target plane and offer an efficient way to construct the surface by solving standard integrals. The efficiency is demonstrated by applying it to two challenging design examples.


title: Bregman Iteration for Correspondence Problems: A Study of Optical Flow

abstract: Bregman iterations are known to yield excellent results for denoising, deblurring and compressed sensing tasks, but so far this technique has rarely been used for other image processing problems. In this paper we give a thorough description of the Bregman iteration, unifying thereby results of different authors within a common framework. Then we show how to adapt the split Bregman iteration, originally developed by Goldstein and Osher for image restoration purposes, to optical flow which is a fundamental correspondence problem in computer vision. We consider some classic and modern optical flow models and present detailed algorithms that exhibit the benefits of the Bregman iteration. By making use of the results of the Bregman framework, we address the issues of convergence and error estimation for the algorithms. Numerical examples complement the theoretical part.


title: Observation of top quark pairs produced in association with a vector boson in pp collisions at sqrt(s) = 8 TeV

abstract: Measurements of the cross sections for top quark pairs produced in association with a W or Z boson are presented, using 8 TeV pp collision data corresponding to an integrated luminosity of 19.5 inverse femtobarns, collected by the CMS experiment at the LHC. Final states are selected in which the associated W boson decays to a charged lepton and a neutrino or the Z boson decays to two charged leptons. Signal events are identified by matching reconstructed objects in the detector to specific final state particles from ttW or ttZ decays. The ttW cross section is measured to be 382 +117 -102 fb with a significance of 4.8 standard deviations from the background-only hypothesis. The ttZ cross section is measured to be 242 +65 -55 fb with a significance of 6.4 standard deviations from the background-only hypothesis. These measurements are used to set bounds on five anomalous dimension-six operators that would affect the ttW and ttZ cross sections.


title: A logarithmic Schr\"{o}dinger equation with asymptotic conditions on the potential

abstract: In this paper we consider a class of logarithmic Schr\"{o}dinger equations with a potential which may change sign. When the potential is coercive, we obtain infinitely many solutions by adapting some arguments of the Fountain theorem, and in the case of bounded potential we obtain a ground state solution, i.e. a nontrivial solution with least possible energy. The functional corresponding to the problem is the sum of a smooth and a convex lower semicontinuous term.


title: A New Method for Compensation and Rematch of Cavity Failure in C-ADS Linac

abstract: For proton linear accelerators used in applications such as accelerator-driven systems, due to the nature of the operation, it is essential for the beam failure rate to be several orders of magnitude lower than usual performance of similar accelerators. A fault-tolerant mechanism should be mandatorily imposed in order to maintain short recovery time, high uptime and extremely low frequency of beam loss. This paper proposes an innovative and challenging way for compensation and rematch of cavity failure using fast electronic devices and Field Programmable Gate Arrays (FPGAs) instead of embedded computers to complete the computation of beam dynamics. A method of building an equivalent model for the FPGA, with optimization using a genetic algorithm, is shown. Results based on the model and algorithm are compared with TRACEWIN simulation to show the precision and correctness of the mechanism.


title: A System for Precise End-to-End Delay Measurements in Video Communication

abstract: Low delay video transmission is becoming increasingly important. Delay critical, video enabled applications range from teleoperation scenarios such as controlling drones or telesurgery to autonomous control through computer vision algorithms applied on real-time video. To judge the quality of the video transmission in such a system, it is important to be able to precisely measure the end-to-end (E2E) delay of the transmitted video. We present a low-complexity system that automatically takes pairwise independent measurements of E2E delay. The precision can be far below the millisecond order, mainly limited by the sampling rate of the measurement system. In our implementation, we achieve a precision of 0.5 milliseconds with a sampling rate of 2kHz.


title: Feedback and Time are Essential for the Optimal Control of Computing Systems

abstract: The performance, reliability, cost, size and energy usage of computing systems can be improved by one or more orders of magnitude by the systematic use of modern control and optimization methods. Computing systems rely on the use of feedback algorithms to schedule tasks, data and resources, but the models that are used to design these algorithms are validated using open-loop metrics. By using closed-loop metrics instead, such as the gap metric developed in the control community, it should be possible to develop improved scheduling algorithms and computing systems that have not been over-engineered. Furthermore, scheduling problems are most naturally formulated as constraint satisfaction or mathematical optimization problems, but these are seldom implemented using state of the art numerical methods, nor do they explicitly take into account the fact that the scheduling problem itself takes time to solve. This paper makes the case that recent results in real-time model predictive control, where optimization problems are solved in order to control a process that evolves in time, are likely to form the basis of scheduling algorithms of the future. We therefore outline some of the research problems and opportunities that could arise by explicitly considering feedback and time when designing optimal scheduling algorithms for computing systems.


title: A study of a sample of high rotation measure AGNs through multifrequency single dish observations

abstract: We characterised and studied, in the radio band, a sample of candidates of high Rotation Measure (RM). These point-like objects show a strong depolarisation at 21cm. This feature suggests the presence of a very dense medium surrounding them in a combination of a strong magnetic field. This work aims at selecting and studying a sample of radio sources with high RM, thus to study their physical conditions and their status with respect to their surrounding medium. We want to understand if any connection is present between the AGN hosting galaxy medium with some evolutionary track and/or some restarting phase of the AGN itself. Multifrequency single-dish observations were performed with the 100-m Effelsberg telescope to define the initial sample, to characterise the SED of the final sample (30 targets) and to determine their RM in the 11 to 2 cm wavelength range. From the observations, the SED together with polarisation information, i.e. the fractional polarisation and the polarisation angle, have been determined. Three different object types were revealed from the SEDs analysis: Older, GPS-like and Mixed. For each of the targets, the rotation measure has been found and the depolarisation has been modelled. No significant correlation have been found between the depolarisation behaviours and the SEDs, while a correlation has been found between sources with mixed SED (with an old component at low frequency and compact components at high frequencies) and high values of the rotation measure (with values in the rest frame larger than 1000 rad/m^2). This work helps us to define and identify a sample of sources with high RM. From the analysis we can conclude that the sources showing a restarting phase at high frequency (with a Mixed SED), are characterised by a really dense and/or a magnetised medium that strongly rotates the polarisation angle at the different frequencies, leading to a high RM.


title: The dearth of nuclear star clusters in bright galaxies

abstract: We investigate the interaction of a massive globular cluster (GC) with a super massive black hole (SMBH), located at the centre of its host galaxy, by means of direct N-body simulations. The results show that tidal distortions induced by the stellar background and the SMBH act on a time shorter than that of dynamical friction decay for a $10^6$ M$_\odot$ GC whenever the SMBH mass exceeds $\sim 10^8$ M$_\odot$ . This implies an almost complete dissolution of the infalling GC before it reaches the inner region ($\lesssim 5$ pc) of the parent galaxy. The generalization of this result to a larger sample of infalling GCs shows that such destructive process may prevent the formation and growth of a bright galactic nucleus. Another interesting, serendipitous, result we obtained is that the close interaction between the SMBH and the GC produces a "wave" of stars that escape from the cluster and, in a fraction, even from the whole galaxy.


title: The Coulomb potential V(r)=1/r and other radial problems on the Bethe lattice

abstract: We study the problem of a particle hopping on the Bethe lattice in the presence of a Coulomb potential. We obtain an exact solution to the particle's Green's function along with the full energy spectrum. In addition, we present a mapping of a generalized radial potential problem defined on the Bethe lattice to an infinite number of one dimensional problems that are easily accessible numerically. The latter method is particularly useful when the problem admits no analytical solution.

categories: cond-mat.stat-mech cond-mat.dis-nn math-ph math.MP


title: Performance Peculiarities of Viterbi Decoder in Mathworks Simulink, GNU Radio and Other Systems with Likewise Implementation

abstract: The performance of convolutional codes decoding by the Viterbi algorithm should not depend on the particular distribution of zeros and ones in the input messages, as they are linear. However, it was identified that specific implementations of Add-Compare-Select unit for the Viterbi Algorithm demonstrate the decoding performance that depends on proportion of elements in the input message. It is conjectured that the modern commercial hard- and software defined communication equipment may also feature similar implementation and as such their decoding performance could also vary.


title: Deconstructing Bataknese Gorga Computationally

abstract: The carved and painted decorations in traditional Batak houses and buildings, gorga, are the source of their exoticism. There are no identical patterns of the ornaments within Batak houses and the drawings are closely related to the way ancient Batak capture the dynamicity of the growing 'tree of life', one of central things within their cosmology and mythology. The survey of ornaments of Batak houses and buildings in Northern Sumatera Indonesia has made us possible to observe the complex pattern. The fractal dimensions of the geometrical shapes in gorga are calculated and they are conjectured into 1.5-1.6, between the dimensional of a line and a plane. The way gorga is drawn is captured by using some modification to the turtle geometry of L-System model, a popular model to model the dynamics of growing plants. The result is a proposal to see Bataknese gorga as one of traditional heritage that may enrich the studies to the generative art.


title: On the algebraicity of some products of special values of Barnes' multiple gamma function

abstract: We consider partial zeta functions $\zeta(s,c)$ associated with ray classes $c$'s of a totally real field. Stark's conjecture implies that an appropriate product of $\exp(\zeta'(0,c))$'s is an algebraic number which is called a Stark unit. Shintani gave an explicit formula for $\exp(\zeta'(0,c))$ in terms of Barnes' multiple gamma function. Yoshida ``decomposed'' Shintani's formula: he defined the symbol $X(c,\iota)$ satisfying that $\exp(\zeta'(0,c))=\prod_{\iota} \exp(X(c,\iota))$ where $\iota$ runs over all real embeddings of $F$. Hence we can decompose a Stark unit into a product of $[F:\mathbb Q]$ terms. The main result is to show that $([F:\mathbb Q]-1)$ of them are algebraic numbers. We also study a relation between Yoshida's conjecture on CM-periods and Stark's conjecture.


title: Comprehensive Census and Complete Characterization of Nearby Debris Disk Stars

abstract: Debris disks are intimately linked to planetary system evolution since the rocky material surrounding the host stars is due to secondary generation from the collisions of planetesimals. With the conclusion and lack of future large scale infrared excess survey missions, it is time to summarize the history of using excess emission in the infrared as a tracer of debris. We have compiled a catalog of infrared excess stars from peer reviewed articles and performed an extensive search for new debris disks by cross correlating the Tycho 2 and AllWISE catalogs. This study will examine each debris disk stars parameters obtained through high resolution spectroscopy at various facilities. We will maintain a webpage devoted to these infrared excess sources and provide various resources related to our catalog creation, SED fitting, and data reduction.


title: Photocurrents in a Single InAs Nanowire/ Silicon Heterojunction

abstract: We investigate the optoelectronic properties of single indium arsenide nanowires, which are grown vertically on p-doped silicon substrates. We apply a scanning photocurrent microscopy to study the optoelectronic properties of the single heterojunctions. The measured photocurrent characteristics are consistent with an excess charge carrier transport through mid-gap trap states, which form at the Si/InAs heterojunctions. Namely, the trap states add an additional transport path across a heterojunction, and the charge of the defects changes the band bending at the junction. The bending gives rise to a photovoltaic effect at a small bias voltage. In addition, we observe a photoconductance effect within the InAs nanowires at large biases.


title: Minimal number of edges in hypergraph guaranteeing perfect fractional matching and MMS conjecture, complete version

abstract: In this paper we prove Ahlswede- Khachatrian conjecture~\cite{1} up to finite number of cases, which can be checked using modern computers. From this conjecture follows conjecture from~\cite{2} and Manickam-Mikl\'{o}s-Singhi conjecture. This paper is combined from papers~\cite{09},~\cite{08}.


title: Reduced Precision Checking to Detect Errors in Floating Point Arithmetic

abstract: In this paper, we use reduced precision checking (RPC) to detect errors in floating point arithmetic. Prior work explored RPC for addition and multiplication. In this work, we extend RPC to a complete floating point unit (FPU), including division and square root, and we present precise analyses of the errors undetectable with RPC that show bounds that are smaller than prior work. We implement RPC for a complete FPU in RTL and experimentally evaluate its error coverage and cost.


title: Universal aspects in the equation of state for Yang-Mills theories

abstract: We present high-precision lattice calculations of the thermodynamics of Yang-Mills theories with different gauge groups. In the confining phase, we show that the equation of state is described remarkably well by a gas of massive, non-interacting glueballs, provided that an effective bosonic closed-string model is used to derive an exponentially growing Hagedorn spectrum for the heavy states. In particular, this model describes very accurately the results for the SU(3) theory reported by Bors\'anyi et al. in JHEP 07 (2012) 056, as well as a novel set of lattice data for the SU(2) theory. In addition, we also also show that the equation of state in the deconfined phase exhibits a near perfect proportionality to the number of gluon degrees of freedom, including for the Yang-Mills theory based on the exceptional, center-less gauge group $G_2$.


title: Local spectral statistics of Gaussian matrices with correlated entries

abstract: We prove optimal local law, bulk universality and non-trivial decay for the off-diagonal elements of the resolvent for a class of translation invariant Gaussian random matrix ensembles with correlated entries.


title: Visual Tracking via Nonnegative Regularization Multiple Locality Coding

abstract: This paper presents a novel object tracking method based on approximated Locality-constrained Linear Coding (LLC). Rather than using a non-negativity constraint on encoding coefficients to guarantee these elements nonnegative, in this paper, the non-negativity constraint is substituted for a conventional $\ell_2$ norm regularization term in approximated LLC to obtain the similar nonnegative effect. And we provide a detailed and adequate explanation in theoretical analysis to clarify the rationality of this replacement. Instead of specifying fixed K nearest neighbors to construct the local dictionary, a series of different dictionaries with pre-defined numbers of nearest neighbors are selected. Weights of these various dictionaries are also learned from approximated LLC in the similar framework. In order to alleviate tracking drifts, we propose a simple and efficient occlusion detection method. The occlusion detection criterion mainly depends on whether negative templates are selected to represent the severe occluded target. Both qualitative and quantitative evaluations on several challenging sequences show that the proposed tracking algorithm achieves favorable performance compared with other state-of-the-art methods.


title: Operator-Theoretic Characterization of Eventually Monotone Systems

abstract: Monotone systems are dynamical systems whose solutions preserve a partial order in the initial condition for all positive times. It stands to reason that some systems may preserve a partial order only after some initial transient. These systems are usually called eventually monotone. While monotone systems have a characterization in terms of their vector fields (i.e. Kamke-Muller condition), eventually monotone systems have not been characterized in such an explicit manner. In order to provide a characterization, we drew inspiration from the results for linear systems, where eventually monotone (positive) systems are studied using the spectral properties of the system (i.e. Perron-Frobenius property). In the case of nonlinear systems, this spectral characterization is not straightforward, a fact that explains why the class of eventually monotone systems has received little attention to date. In this paper, we show that a spectral characterization of nonlinear eventually monotone systems can be obtained through the Koopman operator framework. We consider a number of biologically inspired examples to illustrate the potential applicability of eventual monotonicity.


title: Finite Dimensional Fokker-Planck Equations for Continuous Time Random Walks

abstract: Continuous Time Random Walk(CTRW) is a model where particle's jumps in space are coupled with waiting times before each jump. A Continuous Time Random Walk Limit(CTRWL) is obtained by a limit procedure on a CTRW and can be used to model anomalous diffusion. The distribution $p\left(dx,t\right)$ of a CTRWL $X_{t}$ satisfies a Fractional Fokker-Planck Equation(FFPE). Since CTRWLs are usually not Markovian, their one dimensional FFPE is not enough to completely define them. In this paper we find the FFPEs of the distribution of $X_{t}$ at multiple times , i.e. the distribution of the random vector $\left(X_{t_{1}},...,X_{t_{n}}\right)$ for $t_{1}<...<t_{n}$ for a large class of CTRWLs. This allows us to define CTRWLs by their finite dimensional FFPEs.


title: Stellar oscillations. II The non-adiabatic case

abstract: A leap forward has been performed due to the space-borne missions, MOST, CoRoT and Kepler. They provided a wealth of observational data, and more precisely oscillation spectra, which have been (and are still) exploited to infer the internal structure of stars. While an adiabatic approach is often sufficient to get information on the stellar equilibrium structures it is not sufficient to get a full understanding of the physics of the oscillation. Indeed, it does not permit one to answer some fundamental questions about the oscillations, such as: What are the physical mechanisms responsible for the pulsations inside stars? What determines the amplitudes? To what extent the adiabatic approximation is valid? All these questions can only be addressed by considering the energy exchanges between the oscillations and the surrounding medium. This lecture therefore aims at considering the energetical aspects of stellar pulsations with particular emphasis on the driving and damping mechanisms. To this end, the full non-adiabatic equations are introduced and thoroughly discussed. Two types of pulsation are distinguished, namely the self-excited oscillations that result from an instability and the solar-like oscillations that result from a balance between driving and damping by turbulent convection. For each type, the main physical principles are presented and illustrated using recent observations obtained with the ultra-high precision photometry space-borne missions (MOST, CoRoT and Kepler). Finally, we consider in detail the physics of scaling relations, which relates the seismic global indices with the global stellar parameters and gave birth to the development of statistical (or ensemble) asteroseismology. Indeed, several of these relations rely on the same cause: the physics of non-adiabatic oscillations.


title: Ages of records in random walks

abstract: We consider random walks with continuous and symmetric step distributions. We prove universal asymptotics for the average proportion of the age of the kth longest lasting record for k=1,2,... and for the probability that the record of the kth longest age is broken at step n. Furthermore, we show that the ranked sequence of proportions of ages converges to the Poisson-Dirichlet distribution.


title: Properties of Isostables and Basins of Attraction of Monotone Systems

abstract: In this paper, we investigate geometric properties of monotone systems by studying their isostables and basins of attraction. Isostables are boundaries of specific forward-invariant sets defined by the so-called Koopman operator, which provides a linear infinite-dimensional description of a nonlinear system. First, we study the spectral properties of the Koopman operator and the associated semigroup in the context of monotone systems. Our results generalize the celebrated Perron-Frobenius theorem to the nonlinear case and allow us to derive geometric properties of isostables and basins of attraction. Additionally, we show that under certain conditions we can characterize the bounds on the basins of attraction under parametric uncertainty in the vector field. We discuss computational approaches to estimate isostables and basins of attraction and illustrate the results on two and four state monotone systems.


title: Finite system scheme for infinite rate mutually catalytic branching with infinite branching rate

abstract: For many stochastic diffusion processes with mean field interaction, convergence of the rescaled total mass processes towards a diffusion process is known. Here we show convergence of the so-called finite system scheme for interacting jump-type processes known as mutually catalytic branching processes with infinite branching rate. Due to the lack of second moments the rescaling of time is different from the finite rate mutually catalytic case. The limit of rescaled total mass processes is identified as the finite rate mutually catalytic branching diffusion. The convergence of rescaled processes holds jointly with convergence of coordinate processes, whereas the latter converge at a different time scale.


title: Balancing the Communication Load of Asynchronously Parallelized Machine Learning Algorithms

abstract: Stochastic Gradient Descent (SGD) is the standard numerical method used to solve the core optimization problem for the vast majority of machine learning (ML) algorithms. In the context of large scale learning, as utilized by many Big Data applications, efficient parallelization of SGD is in the focus of active research. Recently, we were able to show that the asynchronous communication paradigm can be applied to achieve a fast and scalable parallelization of SGD. Asynchronous Stochastic Gradient Descent (ASGD) outperforms other, mostly MapReduce based, parallel algorithms solving large scale machine learning problems. In this paper, we investigate the impact of asynchronous communication frequency and message size on the performance of ASGD applied to large scale ML on HTC cluster and cloud environments. We introduce a novel algorithm for the automatic balancing of the asynchronous communication load, which allows to adapt ASGD to changing network bandwidths and latencies.


title: Clustered Low-Rank Tensor Format: Introduction and Application to Fast Construction of Hartree-Fock Exchange

abstract: Clustered Low Rank (CLR) framework for block-sparse and block-low-rank tensor representation and computation is described. The CLR framework depends on 2 parameters that control precision: one controlling the CLR block rank truncation and another that controls screening of small contributions in arithmetic operations on CLR tensors. As these parameters approach zero CLR representation and arithmetic become exact. There are no other ad-hoc heuristics, such as domains. Use of the CLR format for the order-2 and order-3 tensors that appear in the context of density fitting (DF) evaluation of the Hartree-Fock (exact) exchange significantly reduced the storage and computational complexities below their standard $\mathcal{O}(N^3)$ and $\mathcal{O}(N^4)$ figures. Even for relatively small systems and realistic basis sets CLR-based DF HF becomes more efficient than the standard DF approach, and significantly more efficient than the conventional non-DF HF, while negligibly affecting molecular energies and properties.


title: A Note on Jointly Modeling Edges and Node Attributes of a Network

abstract: We are interested in modeling networks in which the connectivity among the nodes and node attributes are random variables and interact with each other. We propose a probabilistic model that allows one to formulate jointly a probability distribution for these variables. This model can be described as a combination of a latent space model and a Gaussian graphical model: given the node variables, the edges will follow independent logistic distributions, with the node variables as covariates; given edges, the node variables will be distributed jointly as multivariate Gaussian, with their conditional covariance matrix depending on the graph induced by the edges. We will present some basic properties of this model, including a connection between this model and a dynamical network process involving both edges and node variables, the marginal distribution of the model for edges as a random graph model, its one-edge conditional distributions, the FKG inequality, and the existence of a limiting distribution for the edges in a infinite graph.


title: Design and Management of Vehicle Sharing Systems: A Survey of Algorithmic Approaches

abstract: Vehicle (bike or car) sharing represents an emerging transportation scheme which may comprise an important link in the green mobility chain of smart city environments. This chapter offers a comprehensive review of algorithmic approaches for the design and management of vehicle sharing systems. Our focus is on one-way vehicle sharing systems (wherein customers are allowed to pick-up a vehicle at any location and return it to any other station) which best suits typical urban journey requirements. Along this line, we present methods dealing with the so-called asymmetric demand-offer problem (i.e. the unbalanced offer and demand of vehicles) typically experienced in one-way sharing systems which severely affects their economic viability as it implies that considerable human (and financial) resources should be engaged in relocating vehicles to satisfy customer demand. The chapter covers all planning aspects that affect the effectiveness and viability of vehicle sharing systems: the actual system design (e.g. number and location of vehicle station facilities, vehicle fleet size, vehicles distribution among stations); customer incentivisation schemes to motivate customer-based distribution of bicycles/cars (such schemes offer meaningful incentives to users so as to leave their vehicle to a station different to that originally intended and satisfy future user demand); cost-effective solutions to schedule operator-based repositioning of bicycles/cars (by employees explicitly enrolled in vehicle relocation) based on the current and future (predicted) demand patterns (operator-based and customer-based relocation may be thought as complementary methods to achieve the intended distribution of vehicles among stations).


title: An MCMC study of general squark flavour mixing in the MSSM

abstract: We present an extensive study of non-minimally flavour violating (NMFV) terms in the Lagrangian of the Minimal Supersymmetric Standard Model (MSSM). We impose a variety of theoretical and experimental constraints and perform a detailed scan of the parameter space by means of a Markov Chain Monte-Carlo (MCMC) setup. This represents the first study of several non-zero flavour-violating elements within the MSSM. We present the results of the MCMC scan with a special focus on the flavour-violating parameters. Based on these results, we define benchmark scenarios for future studies of NMFV effects at the LHC.


title: Completeness of Sums of Subspace of Bounded Functions and Applications

abstract: We give a new proof of a characterization of the closeness of the range of a continuous linear operator and of the closeness of the sum of two closed vector subspaces of a Banach space. Then we state sufficient conditions for the closeness of the sum of two closed subspaces of the Banach space of bounded functions and apply this result on various pseudo almost periodic spaces and pseudo almost automorphic spaces.


title: eROSITA - Nearby Young Stars in X-rays

abstract: X-ray surveys are well suited to detect, identify and study young stars based on their high levels of magnetic activity and thus X-ray brightness. The eROSITA instrument onboard the Spectrum-Roentgen-Gamma (SRG) satellite will perform an X-ray all-sky survey that surpasses existing data by a sensitivity increase of more than an order of magnitude. The 4 yr survey is expected to detect more than half a million stars and stellar systems in X-rays.


title: Multiplons in the two-hole excitation spectra of the one-dimensional Hubbard model

abstract: Using the density-matrix renormalization group in combination with the Chebyshev polynomial expansion technique, we study the two-hole excitation spectrum of the one-dimensional Hubbard model in the entire filling range from the completely occupied band (n=2) down to half-filling (n=1). For strong interactions, the spectra reveal multiplon physics, i.e., relevant final states are characterized by two (doublon), three (triplon), four (quadruplon) and more holes, potentially forming stable compound objects or resonances with finite lifetime. These give rise to several satellites in the spectra with largely different spectral weights as well as to different continua. The complex multiplon phenomenology is analyzed by interpreting not only local and k-resolved two-hole spectra but also three- and four-hole spectra for the Hubbard model and by referring to effective low-energy models. In addition, a filter-operator technique is presented which allows to extract specific information on the final states at a given excitation energy. While multiplons composed of an odd number of holes do neither form stable compounds nor well-defined resonances unless a nearest-neighbor density interaction V is added to the Hamiltonian, the doublon and the quadruplon are well-defined resonances. The k-resolved four-hole spectrum at n=2 represents an interesting special case where a completely stable quadruplon turns into a resonance by merging with the doublon-doublon continuum at a critical wave vector. For all fillings with n > 1, the doublon lifetime is strongly k-dependent and is even infinite at the Brillouin zone edges as demonstrated by k-resolved two-hole spectra. This can be traced back to the "hidden" charge-SU(2) symmetry of the model which is explicitly broken off half-filling and gives rise to a massive collective excitation, even for arbitrary higher-dimensional but bipartite lattices.


title: On the convergence rate of grid search for polynomial optimization over the simplex

abstract: We consider the approximate minimization of a given polynomial on the standard simplex, obtained by taking the minimum value over all rational grid points with given denominator ${r} \in \mathbb{N}$. It was shown in [De Klerk, E., Laurent, M., Sun, Z.: An error analysis for polynomial optimization over the simplex based on the multivariate hypergeometric distribution. {\em SIAM J. Optim.} 25(3) 1498--1514 (2015)] that the relative accuracy of this approximation depends on $r$ as $O(1/r^2)$ if there exists a rational global minimizer. In this note we show that the rational minimizer condition is not necessary to obtain the $O(1/r^2)$ bound.


title: Cross-phase modulation of a probe stored in a waveguide for non-destructive detection of photonic qubits

abstract: Non-destructive detection of photonic qubits is an enabling technology for quantum information processing and quantum communication. For practical applications such as quantum repeaters and networks, it is desirable to implement such detection in a way that allows some form of multiplexing as well as easy integration with other components such as solid-state quantum memories. Here we propose an approach to non-destructive photonic qubit detection that promises to have all the mentioned features. Mediated by an impurity-doped crystal, a signal photon in an arbitrary time-bin qubit state modulates the phase of an intense probe pulse that is stored during the interaction. Using a thulium-doped waveguide in LiNbO$_3$, we perform a proof-of-principle experiment with macroscopic signal pulses, demonstrating the expected cross-phase modulation as well as the ability to preserve the coherence between temporal modes. Our findings open the path to a new key component of quantum photonics based on rare-earth-ion doped crystals.


title: Tight lower bounds on the number of bicliques in false-twin-free graphs

abstract: A \emph{biclique} is a maximal bipartite complete induced subgraph of $G$. Bicliques have been studied in the last years motivated by the large number of applications. In particular, enumeration of the maximal bicliques has been of interest in data analysis. Associated with this issue, bounds on the maximum number of bicliques were given. In this paper we study bounds on the minimun number of bicliques of a graph. Since adding false-twin vertices to $G$ does not change the number of bicliques, we restrict to false-twin-free graphs. We give a tight lower bound on the minimum number bicliques for a subclass of $\{C_4$,false-twin$\}$-free graphs and for the class of $\{K_3$,false-twin$\}$-free graphs. Finally we discuss the problem for general graphs.


title: Pentagonal Monolayer Crystals of Carbon, Boron Nitride, and Silver Azide

abstract: In this study we present a theoretical investigation of structural, electronic and mechanical properties of pentagonal monolayers of carbon (p-graphene), boron nitride (p-B$_{2}$N$_{4}$ and p-B$_{4}$N$_{2}$) and silver azide (p-AgN$_{3}$) by performing state-of-the-art first principles calculations. Our total energy calculations suggest feasible formation of monolayer crystal structures composed entirely of pentagons. In addition, electronic band dispersion calculations indicate that while p-graphene and p-AgN$_{3}$ are semiconductors with indirect bandgaps, p-BN structures display metallic behavior. We also investigate the mechanical properties (in-plane stiffness and the Poisson's ratio) of four different pentagonal structures under uniaxial strain. p-graphene is found to have the highest stiffness value and the corresponding Poisson's ratio is found to be negative. Similarly, p-B$_{2}$N$_{4}$ and p-B$_{4}$N$_{2}$ have negative Poisson's ratio values. On the other hand, the p-AgN$_{3}$ has a large and positive Poisson's ratio. In dynamical stability tests based on calculated phonon spectra of these pentagonal monolayers, we find that only p-graphene and p-B$_{2}$N$_{4}$ are stable, but p-AgN$_{3}$ and p-B$_{4}$N$_{2}$ are vulnerable against vibrational excitations.


title: On the number of unary-binary tree-like structures with restrictions on the unary height

abstract: We consider various classes of Motzkin trees as well as lambda-terms for which we derive asymptotic enumeration results. These classes are defined through various restrictions concerning the unary nodes or abstractions, respectively: We either bound their number or the allowed levels of nesting. The enumeration is done by means of a generating function approach and singularity analysis. The generating functions are composed of nested square roots and exhibit unexpected phenomena in some of the cases. Furthermore, we present some observations obtained from generating such terms randomly and explain why usually powerful tools for random generation, such as Boltzmann samplers, face serious difficulties in generating lambda-terms.


title: Note on the stability of viscous roll-waves

abstract: In this note, we announce a complete classification of stability of periodic roll-wave solutions of the viscous shallow-water equations, from their onset at Froude number $F\approx 2$ up to the infinite-Froude limit. For intermediate Froude numbers, we obtain numerically a particularly simple power-law relation between $F$ and the boundaries of the region of stable periods, that appears potentially useful in hydraulic engineering applications. In the asymptotic regime $F\to 2$ (onset), we provide an analytic expression of the stability boundaries whereas in the limit $F\to\infty$, we show that roll-waves are always unstable.


title: Gravity and nonabelian gauge fields in noncommutative space-time

abstract: Noncommutative geometric construction of gravity in the two sheeted spacetime can be viewed as a discretized version of a Kaluza-Klein theory. In this paper, we show that it is possible to incorporate the nonabelian gauge fields in the same framework. The generalized Hilbert-Einstein action is gauge invariant only in two cases. In the first case, the gauge group must be abelian on one sheet of spacetime and nonabelian on the other one. In the second case, the gauge group must be the same on two sheets of spacetime. Actually, the theories of electroweak and strong interactions are exactly these two cases.


title: Quasinormal mode theory and modelling of electron energy loss spectroscopy for plasmonic nanostructures

abstract: Understanding light-matter interactions using localized surface plasmons (LSPs) is of fundamental interest in classical and quantum plasmonics and has a wide range of applications. In order to understand the spatial properties of LSPs, electron energy loss spectroscopy (EELS) is a common and powerful method of spatially resolving the extreme localized fields that can be obtained with metal resonators. However, modelling EELS for general shaped resonators presents a major challenge in computational electrodynamics, requiring the full photon Green function as a function of two space points and frequency. Here we present an intuitive and computationally simple method for computing EELS maps of plasmonic resonators using a quasinormal mode (QNM) expansion technique. By separating the contribution of the QNM and the bulk material, we give closed-form analytical formulas for the plasmonic QNM contribution to the EELS maps. We exemplify our technique for a split ring resonator, a gold nanorod, and a nanorod dimer structure. The method is accurate, intuitive, and gives orders of magnitude improvements over direct dipole simulations that numerically solve the full 3D Maxwell equations. We also show how the same QNM Green function can be used to obtain the Purcell factor (and projected local density of optical states) from quantum dipole emitters or two level atoms, and we demonstrate how the spectral features differ in general to the EELS spectrum.


title: On the Online Frank-Wolfe Algorithms for Convex and Non-convex Optimizations

abstract: In this paper, the online variants of the classical Frank-Wolfe algorithm are considered. We consider minimizing the regret with a stochastic cost. The online algorithms only require simple iterative updates and a non-adaptive step size rule, in contrast to the hybrid schemes commonly considered in the literature. Several new results are derived for convex and non-convex losses. With a strongly convex stochastic cost and when the optimal solution lies in the interior of the constraint set or the constraint set is a polytope, the regret bound and anytime optimality are shown to be ${\cal O}( \log^3 T / T )$ and ${\cal O}( \log^2 T / T)$, respectively, where $T$ is the number of rounds played. These results are based on an improved analysis on the stochastic Frank-Wolfe algorithms. Moreover, the online algorithms are shown to converge even when the loss is non-convex, i.e., the algorithms find a stationary point to the time-varying/stochastic loss at a rate of ${\cal O}(\sqrt{1/T})$. Numerical experiments on realistic data sets are presented to support our theoretical claims.


title: Consistent Pricing of VIX and Equity Derivatives with the 4/2 Stochastic Volatility Plus Jumps Model

abstract: In this paper, we develop a 4/2 stochastic volatility plus jumps model, namely, a new stochastic volatility model including the Heston model and 3/2 model as special cases. Our model is highly tractable by applying the Lie symmetries theory for PDEs, which means that the pricing procedure can be performed efficiently. In fact, we obtain a closed-form solution for the joint Fourier-Laplace transform so that equity and realized-variance derivatives can be priced. We also employ our model to consistently price equity and VIX derivatives. In this process, the quasi-closed-form solutions for future and option prices are derived. Furthermore, through adopting data on daily VIX future and option prices, we investigate our model along with the Heston model and 3/2 model and compare their different performance in practice. Our result illustrates that the 4/2 model with an instantaneous volatility of the form $(a\sqrt{V_t}+b/\sqrt{V_t})$ for some constants $a, b$ presents considerable advantages in pricing VIX derivatives.


title: Lagrangian and Hamiltonian structures in an integrable hierarchy and space-time duality

abstract: We define and illustrate the novel notion of dual integrable hierarchies, on the example of the nonlinear Schr\"odinger (NLS) hierarchy. For each integrable nonlinear evolution equation (NLEE) in the hierarchy, dual integrable structures are characterized by the fact that the zero-curvature representation of the NLEE can be realized by two Hamiltonian formulations stemming from two distinct choices of the configuration space, yielding two inequivalent Poisson structures on the corresponding phase space and two distinct Hamiltonians. This is fundamentally different from the standard bi-Hamiltonian or generally multitime structure. The first formulation chooses purely space-dependent fields as configuration space; it yields the standard Poisson structure for NLS. The other one is new: it chooses purely time-dependent fields as configuration space and yields a different Poisson structure at each level of the hierarchy. The corresponding NLEE becomes a {\it space} evolution equation. We emphasize the role of the Lagrangian formulation as a unifying framework for deriving both Poisson structures, using ideas from covariant field theory. One of our main results is to show that the two matrices of the Lax pair satisfy the same form of ultralocal Poisson algebra (up to a sign) characterized by an $r$-matrix structure, whereas traditionally only one of them is involved in the classical $r$-matrix method. We construct explicit dual hierarchies of Hamiltonians, and Lax representations of the triggered dynamics, from the monodromy matrices of either Lax matrix. An appealing procedure to build a multi-dimensional lattice of Lax pair, through successive uses of the dual Poisson structures, is briefly introduced.


title: Size, vacancy and temperature effects on Young's modulus of silicene nanoribbons

abstract: We report results on the Young's modulus (YM) of defect-free and defective silicene nanoribbons (SNRs) as a function of length and temperature. In this study, we perform molecular dynamics simulations using the Environment-Dependent Interatomic Potential (EDIP) to describe the interaction of the Si atoms. We find that the Young's modulus of pristine and defective SNRs increases with the ribbon length in both chirality directions. It is shown that the Young's modulus of defective SNRs exhibit a complex dependence on the combinations of vacancies. With respect to temperature, we find that YM for SNRs with and without vacancy defects shows a nonlinear behavior and it could be tailoring for a given length and chirality.


title: Cross-Device Tracking: Matching Devices and Cookies

abstract: The number of computers, tablets and smartphones is increasing rapidly, which entails the ownership and use of multiple devices to perform online tasks. As people move across devices to complete these tasks, their identities becomes fragmented. Understanding the usage and transition between those devices is essential to develop efficient applications in a multi-device world. In this paper we present a solution to deal with the cross-device identification of users based on semi-supervised machine learning methods to identify which cookies belong to an individual using a device. The method proposed in this paper scored third in the ICDM 2015 Drawbridge Cross-Device Connections challenge proving its good performance.


title: Energy-Efficient Data Transmission with Non-FIFO Packets

abstract: This paper investigates the problem of energy-efficient packet transmission with arbitrary arrival instants and deadline constraints over a point-to-point Additive White Gaussian Noise (AWGN) channel. This is different from previous work where it is assumed that the packets follow a First-In-First-Out (FIFO) order in that the packets that arrive earlier will have a deadline that is also earlier. We first investigate the necessary and sufficient conditions of the optimal transmission scheduler. We then propose an algorithm which finds the transmission schedule of each packet in the order of the packets with the largest transmission rate to the packets with the smallest transmission rate. Finally, we show that our algorithm satisfies the sufficient conditions of the optimal transmission scheduler and thus, is optimal.


title: Powerful, Rotating Disk Winds from Stellar-mass Black Holes

abstract: We present an analysis of ionized X-ray disk winds observed in the Fe K band of four stellar-mass black holes observed with Chandra, including 4U 1630-47, GRO J1655-40, H 1743-322, and GRS 1915+105. High-resolution photoionization grids were generated in order to model the data. Third-order gratings spectra were used to resolve complex absorption profiles into atomic effects and multiple velocity components. The Fe XXV line is found to be shaped by contributions from the intercombination line (in absorption), and the Fe XXVI line is detected as a spin-orbit doublet. The data require 2-3 absorption zones, depending on the source. The fastest components have velocities approaching or exceeding 0.01c, increasing mass outflow rates and wind kinetic power by orders of magnitude over prior single-zone models. The first-order spectra require re-emission from the wind, broadened by a degree that is loosely consistent with Keplerian orbital velocities at the photoionization radius. This suggests that disk winds are rotating with the orbital velocity of the underlying disk, and provides a new means of estimating launching radii -- crucial to understanding wind driving mechanisms. Some aspects of the wind velocities and radii correspond well to the broad-line region (BLR) in active galactic nuclei, suggesting a physical connection. We discuss these results in terms of prevalent models for disk wind production and disk accretion itself, and implications for massive black holes in active galactic nuclei.


title: Analytic functions on tubes of non-Archimedean analytic spaces

abstract: Let $k$ be a discretely valued non-Archimedean field. We give an explicit description of analytic functions whose norm is bounded by a given real number $r$ on tubes of reduced $k$-analytic spaces associated to special formal schemes (those include $k$-affinoid spaces as well as open polydiscs). As an application we study the connectedness of these tubes. This generalizes (in the discretely valued case) a result of Siegfried Bosch. We use as a main tool a result of A.J. de Jong relating formal and analytic functions on special formal schemes and a generalization of de Jong's result which is proved in the joint appendix with Christian Kappen.


title: Evolvable Autonomic Management

abstract: Autonomic management is aimed at adapting to uncertainty. Hence, it is devised as m-connected k-dominating set problem, resembled by dominator and dominate, such that dominators are resilient up to m-1 uncertainty among them and dominate are resilient up to k-1 uncertainty on their way to dominators. Therefore, an evolutionary algorithm GENESIS is proposed, which resolves uncertainty by evolving population of solutions, while considering uncertain constraints as sub-problems, started by initial populations by a greedy algorithm AVIDO. Theoretical analysis first justifies original problem as NP-hard problem. Eventually, the absence of polynomial time approximation scheme necessitates justification of original problem as multiobjective optimization problem. Furthermore, approximation to Pareto front is verified to be decomposed into scalar optimization sub-problems, which lays out the theoretical foundation for decomposition based evolutionary solution. Finally, case-study, feasibility analysis and exemplary implication are presented for evolvable autonomic management in combined cancer treatment with in-vivo sensor networks.


title: Modeling Aftershocks as a Stretched Exponential Relaxation

abstract: The decay rate of aftershocks has been modeled as a power law since the pioneering work of Omori in the late nineteenth century. Considered the second most fundamental empirical law after the Gutenberg-Richter relationship, the power law paradigm has rarely been challenged by the seismological community. By taking a view of aftershock research not biased by prior conceptions of Omori power law decay and by applying statistical methods recommended in applied mathematics, I show that all aftershock sequences tested in three regional earthquake catalogs (Southern and Northern California, Taiwan) and with three declustering techniques (nearest-neighbor, second-order moment, window methods) follow a stretched exponential instead of a power law. These results infer that aftershocks are due to a simpler relaxation process than originally thought, in accordance with most other relaxation processes observed in Nature.


title: Searches for a heavy scalar boson H decaying to a pair of 125 GeV Higgs bosons hh or for a heavy pseudoscalar boson A decaying to Zh, in the final states with h to tau tau

abstract: A search for a heavy scalar boson H decaying into a pair of lighter standard-model-like 125 GeV Higgs bosons h and a search for a heavy pseudoscalar boson A decaying into a Z and an h boson are presented. The searches are performed on a data set corresponding to an integrated luminosity of 19.7 inverse femtobarns of pp collision data at a centre-of-mass energy of 8 TeV, collected by CMS in 2012. A final state consisting of two tau leptons and two b jets is used to search for the H to hh decay. A final state consisting of two tau leptons from the h boson decay, and two additional leptons from the Z boson decay, is used to search for the decay A to Zh. The results are interpreted in the context of the two-Higgs-doublet models. No excess is found above the standard model expectation and upper limits are set on the heavy boson production cross sections in the mass ranges 260 < m[H] < 350 GeV and 220 < m[A] < 350 GeV.


title: Resonant-state expansion of dispersive open optical systems

abstract: A resonant-state expansion (RSE) for open optical systems with a general frequency dispersion of the relative permittivity, described by a finite number of simple poles, is presented. As in the non-dispersive case, the RSE of dispersive systems converts Maxwell's wave equation into a linear matrix eigenvalue problem in the basis of unperturbed resonant states, in this way numerically exactly determining all relevant eigenmodes of the optical system. This dispersive RSE is verified by application to the analytically solvable system of a sphere in vacuum, with a dispersion of the dielectric constant described by the Drude and Drude-Lorentz models. We calculate the change of the optical modes when converting the sphere material from gold to non-dispersive silica and back to gold, and evaluate the accuracy using the exact solutions.


title: Robinson--Trautman solution with nonlinear electrodynamics

abstract: Explicit Robinson--Trautman solutions with electromagnetic field satisfying nonlinear field equations are derived and analyzed. The solutions are generated from the spherically symmetric ones. In all cases the electromagnetic field singularity is removed while the gravitational one persists. The models resolving curvature singularity were not possible to generalize to Robinson--Trautman geometry indicating that the removal of singularity in associated spherically symmetric case is just a consequence of high symmetry. We show that the solutions are generally of algebraic type II but reduce to type D in spherical symmetry. Asymptotically they tend to the spherically symmetric case as well.


title: Born-Infeld Gravity with a Unique Vacuum and a Massless Graviton

abstract: We construct an n-dimensional Born-Infeld type gravity theory that has the same properties as Einstein's gravity in terms of the vacuum and particle content: Namely, the theory has a unique viable vacuum (maximally symmetric solution) and a single massless unitary spin-2 graviton about this vacuum. The BI gravity, in some sense, is the most natural, minimal generalization of Einstein's gravity with a better UV behavior, and hence, is a potentially viable proposal for low energy quantum gravity. The Gauss-Bonnet combination plays a non-trivial role in the construction of the theory. As an extreme example, we consider the infinite dimensional limit where an interesting exponential gravity arises.


title: Surface induced magnetization reversal of MnP nanoclusters embedded in GaP

abstract: We investigate the quasi-static magnetic behavior of ensembles of non-interacting ferromagnetic nanoparticles consisting of MnP nanoclusters embedded in GaP(001) epilayers grown at 600, 650 and 700{\deg}C. We use a phenomenological model, in which surface effects are included, to reproduce the experimental hysteresis curves measured as a function of temperature (120-260 K) and direction of the applied field. The slope of the hysteresis curve during magnetization reversal is determined by the MnP nanoclusters size distribution, which is a function of the growth temperature. Our results show that the coercive field is very sensitive to the strength of the surface anisotropy, which reduces the energy barrier between the two states of opposite magnetization. Notably, this reduction in the energy barrier increases by a factor of 3 as the sample temperature is lowered from 260 to 120 K.


title: Convex pentagons that admit $i$-block transitive tilings

abstract: The problem of classifying the convex pentagons that admit tilings of the plane is a long-standing unsolved problem. Previous to this article, there were 14 known distinct kinds of convex pentagons that admit tilings of the plane. Five of these types admit tile-transitive tilings (i.e. there is a single transitivity class with respect to the symmetry group of the tiling). The remaining 9 types do not admit tile-transitive tilings, but do admit either 2-block transitive tilings or 3-block transitive tilings; these are tilings comprised of clusters of 2 or 3 pentagons such that these clusters form tile-2-transitive or tile-3-transitive tilings. In this article, we present some combinatorial results concerning pentagons that admit $i$-block transitive tilings for $i \in \mathbb{N}$. These results form the basis for an automated approach to finding all pentagons that admit $i$-block transitive tilings for each $i \in \mathbb{N}$. We will present the methods of this algorithm and the results of the computer searches so far, which includes a complete classification of all pentagons admitting 1-, 2-, and 3-block transitive tilings, among which is a new 15th type of convex pentagon that admits a tile-3-transitive tiling.


title: Static behaviour of induced seismicity

abstract: The standard paradigm to describe seismicity induced by fluid injection is to apply nonlinear diffusion dynamics in a poroelastic medium. I show that the spatiotemporal behaviour and rate evolution of induced seismicity can, instead, be expressed by geometric operations on a static stress field produced by volume change at depth. I obtain laws similar in form to the ones derived from poroelasticity while requiring a lower description length. Although fluid flow is known to occur in the ground, it is not pertinent to the behaviour of induced seismicity. The proposed model is equivalent to the static stress model for tectonic foreshocks generated by the Non- Critical Precursory Accelerating Seismicity Theory. This study hence verifies the explanatory power of this theory outside of its original scope.


title: Analytic Posteriors for Pearson's Correlation Coefficient

abstract: Pearson's correlation is one of the most common measures of linear dependence. Recently, Bernardo (2015) introduced a flexible class of priors to study this measure in a Bayesian setting. For this large class of priors we show that the (marginal) posterior for Pearson's correlation coefficient and all of the posterior moments are analytic. Our results are available in the open-source software package JASP.


title: Bayesian modelling of the emission spectrum of the JET Li-BES system

abstract: A Bayesian model of the emission spectrum of the JET lithium beam has been developed to infer the intensity of the Li I (2p-2s) line radiation and associated uncertainties. The detected spectrum for each channel of the lithium beam emission spectroscopy (Li-BES) system is here modelled by a single Li line modified by an instrumental function, Bremsstrahlung background, instrumental offset, and interference filter curve. Both the instrumental function and the interference filter curve are modelled with non-parametric Gaussian processes. All free parameters of the model, the intensities of the Li line, Bremsstrahlung background, and instrumental offset, are inferred using Bayesian probability theory with a Gaussian likelihood for photon statistics and electronic background noise. The prior distributions of the free parameters are chosen as Gaussians. Given these assumptions, the intensity of the Li line and corresponding uncertainties are analytically available using a Bayesian linear inversion technique. The proposed approach makes it possible to extract the intensity of Li line without doing a separate background subtraction through modulation of the Li beam.


title: Formulation of discontinuous Galerkin methods for relativistic astrophysics

abstract: The DG algorithm is a powerful method for solving pdes, especially for evolution equations in conservation form. Since the algorithm involves integration over volume elements, it is not immediately obvious that it will generalize easily to arbitrary time-dependent curved spacetimes. We show how to formulate the algorithm in such spacetimes for applications in relativistic astrophysics. We also show how to formulate the algorithm for equations in non-conservative form, such as Einstein's field equations themselves. We find two computationally distinct formulations in both cases, one of which has seldom been used before for flat space in curvilinear coordinates but which may be more efficient. We also give a new derivation of the ALE algorithm (Arbitrary Lagrangian-Eulerian) using 4-vector methods that is much simpler than the usual derivation and explains why the method preserves the conservation form of the equations. The various formulations are explored with some simple numerical experiments that also explore the effect of the metric identities on the results.

categories: gr-qc astro-ph.HE math.NA physics.comp-ph


title: Rational Points on the Intersection of Three Quadrics

abstract: This paper proves the Hasse principle and weak approximation for varieties defined by the smooth intersection of three quadratics in at least 19 variables, over arbitrary number fields.