title: BayesSummaryStatLM: An R package for Bayesian Linear Models for Big Data and Data Science

abstract: Recent developments in data science and big data research have produced an abundance of large data sets that are too big to be analyzed in their entirety, due to limits on either computer memory or storage capacity. Here, we introduce our R package 'BayesSummaryStatLM' for Bayesian linear regression models with Markov chain Monte Carlo implementation that overcomes these limitations. Our Bayesian models use only summary statistics of data as input; these summary statistics can be calculated from subsets of big data and combined over subsets. Thus, complete data sets do not need to be read into memory in full, which removes any physical memory limitations of a user. Our package incorporates the R package 'ff' and its functions for reading in big data sets in chunks while simultaneously calculating summary statistics. We describe our Bayesian linear regression models, including several choices of prior distributions for unknown model parameters, and illustrate capabilities and features of our R package using both simulated and real data sets.


title: Millimeter polarisation of the protoplanetary nebula OH 231.8+4.2: A follow-up study with CARMA

abstract: In order to investigate the characteristics and influence of the magnetic field in evolved stars, we performed a follow-up investigation of our previous submillimeter analysis of the proto-planetary nebula (PPN) OH 231.8+4.2 (Sabin et al. 2014), this time at 1.3mm with the CARMA facility in polarisation mode for the purpose of a multi-scale analysis. OH 231.8+4.2 was observed at ~2.5" resolution and we detected polarised emission above the 3-sigma threshold (with a mean polarisation fraction of 3.5 %). The polarisation map indicates an overall organised magnetic field within the nebula. The main finding in this paper is the presence of a structure mostly compatible with an ordered toroidal component that is aligned with the PPN's dark lane. We also present some alternative magnetic field configuration to explain the structure observed. These data complete our previous SMA submillimeter data for a better investigation and understanding of the magnetic field structure in OH 231.8+4.2.


title: Parametric solutions of Pell equations

abstract: This short paper is concerned with polynomial Pell equations \[P^2-DQ^2=1,\] with $P,Q,D\in\Bbb C[X]$ and ${deg}(D)=2$. The main result shows that the polynomials $P$ and $Q$ are closely related to Chebyshev polynomials. We then investigate the existence of such polynomials in $\Bbb Z[X]$ specializing to fixed solutions of ordinary Pell equations over the integers.


title: A compression scheme for radio data in high performance computing

abstract: We present a procedure for efficiently compressing astronomical radio data for high performance applications. Integrated, post-correlation data are first passed through a nearly lossless rounding step which compares the precision of the data to a generalized and calibration-independent form of the radiometer equation. This allows the precision of the data to be reduced in a way that has an insignificant impact on the data. The newly developed Bitshuffle lossless compression algorithm is subsequently applied. When the algorithm is used in conjunction with the HDF5 library and data format, data produced by the CHIME Pathfinder telescope is compressed to 28% of its original size and decompression throughputs in excess of 1 GB/s are obtained on a single core.


title: Scattering of electromagnetic waves by small impedance particles of an arbitrary shape

abstract: An explicit formula is derived for the electromagnetic (EM) field scattered by one small impedance particle $D$ of an arbitrary shape. If $a$ is the characteristic size of the particle, $\lambda$ is the wavelength, $a<<\lambda$ and $\zeta$ is the boundary impedance of $D$, $[N,[E,N]]=\zeta [N,H]$ on $S$, where $S$ is the surface of the particle, $N$ is the unit outer normal to $S$, and $E$, $H$ is the EM field, then the scattered field is $E_{sc}=[\nabla g(x,x_1), Q]$. Here $g(x,y)=\frac{e^{ik|x-y|}}{4\pi |x-y|}$, $k$ is the wave number, $x_1\in D$ is an arbitrary point, and $Q=-\frac{\zeta |S|}{i\omega \mu}\tau \nabla \times E_0$, where $E_0$ is the incident field, $|S|$ is the area of $S$, $\omega$ is the frequency, $\mu$ is the magnetic permeability of the space exterior to $D$, and $\tau$ is a tensor which is calculated explicitly. The scattered field is $O(|\zeta| a^2)>> O(a^3)$ as $a\to 0$ when $\lambda$ is fixed and $\zeta$ does not depend on $a$. Thus, $|E_{sc}|$ is much larger than the classical value $O(a^3)$ for the field scattered by a small particle. It is proved that the effective field in the medium, in which many small particles are embedded, has a limit as $a\to 0$ and the number $M=M(a)$ of the particles tends to $\infty$ at a suitable rate. Thislimit solves a linear integral equation. The refraction coefficient of the limiting medium is calculated analytically. This yields a recipe for creating materials with a desired refraction coefficient.


title: A route to room temperature ferromagnetic ultrathin SrRuO$_3$ films

abstract: Experimental efforts to stabilize ferromagnetism in ultrathin films of transition metal oxides have so far failed, despite expectations based on density functional theory (DFT) and DFT+U. Here, we investigate one of the most promising materials, SrRuO$_3$, and include correlation effects beyond DFT by means of dynamical mean field theory. In agreement with experiment we find an intrinsic thickness limitation for metallic ferromagnetism in SrRuO$_3$ thin films. Indeed, we demonstrate that the realization of ultrathin ferromagnetic films is out of reach of standard thin-film techniques. Proposing charge carrier doping as a new route to manipulate thin films, we predict room temperature ferromagnetism in electron-doped SrRuO$_3$ ultra thin films.


title: Characteristic analysis for integrable soliton models on two-dimensional target spaces

abstract: We investigate the evolutionary aspects of some integrable soliton models whose Lagrangians are derived from the pullback of a volume-form to a two-dimensional target space. These models are known to have infinitely many conserved quantities and support various types of exact analytic solutions with nontrivial topology. In particular, we show that, in spite of the fact that they admit nice smooth solutions, wave propagation about these solutions will always be ill-posed. This is related to the fact that the corresponding Euler-Lagrange equations are not of hyperbolic type.


title: A short proof of the existence of the solution to elliptic boundary problem

abstract: There are several methods for proving the existence of the solution to the elliptic boundary problem $Lu=f \text{\,\, in\,\,} D,\quad u|_S=0,\quad (*)$. Here $L$ is an elliptic operator of second order, $f$ is a given function, and uniqueness of the solution to problem (*) is assumed. The known methods for proving the existence of the solution to (*) include variational methods, integral equation methods, method of upper and lower solutions. In this paper a method based on functional analysis is proposed. This method is conceptually simple and technically is easy. It requires some known a priori estimates and a continuation in a parameter method.


title: Two samples test for discrete power-law distributions

abstract: Power-law distributions occur in wide variety of physical, biological, and social phenomena. In this paper, we propose a statistical hypothesis test based on the log-likelihood ratio to assess whether two samples of discrete data are drawn from the same power-law distribution.

categories: stat.ME physics.data-an physics.soc-ph stat.AP


title: Low-Thrust Transfer Between Circular Orbits Using Natural Precession

abstract: The minimum-fuel low-thrust transfer between circular orbits is formulated using the Edelbaum's averaged dynamics with the addition of the nodal precession due to the first zonal term. The extremal analysis shows that an optimal transfer is composed of three sequences in the regular case. The optimal control problem is solved by a shooting method with a costate guess derived from an approximate solution.


title: Direct correlation between aromatization of carbon- rich organic matter and its visible electronic absorption edge

abstract: The evolution of the electronic absorption edge of type I, II and III kerogen is studied by diffuse reflectance UV-Visible absorption spectroscopy. The functional form of the electronic absorption edge for all kerogens measured is in excellent agreement with the "Urbach tail" phenomenology. The Urbach decay width extracted from the exponential fit within the visible range is strongly correlated with the aliphatic/aromatic ratio in isolated kerogen, regardless of the kerogen type. No correlation is found between the decay width and the average size of aromatic clusters, which is explained in terms of a non-linear increase in optical absorption with increasing size of the aromatic clusters determined by 13C NMR. Further, absorption spectra calculated with density functional theory calculations on proxy ensemble models of kerogen are in excellent agreement with the experimental results. The correlation of the decay width with conventional maturity indicators such as vitrinite reflectance is found to be good within a particular kerogen type, but not consistent across different kerogen types, reflecting systematic variations in bulk composition for different type kerogen types with the same vitrinite reflectance. Thus, diffuse reflectance visible absorption spectroscopy is presented as a rapid, calibrated and non-destructive method to monitor both the maturity and the chemical composition of kerogen. The chemical insight of kerogen in relation to its optical absorption provided by this methodology may serve for rapid screening of kerogen for electronics and optical devices in place of functionalized produced carbon.


title: Scanning tunneling spectroscopy reveals a silicon dangling bond charge state transition

abstract: We report the study of single dangling bonds (DB) on the hydrogen terminated silicon (100) surface using a low temperature scanning tunneling microscope (LT-STM). By investigating samples prepared with different annealing temperatures, we establish the critical role of subsurface arsenic dopants on the DB electronic properties. We show that when the near surface concentration of dopants is depleted as a result of $1250{\deg}C$ flash anneals, a single DB exhibits a sharp conduction step in its I(V) spectroscopy that is not due to a density of states effect but rather corresponds to a DB charge state transition. The voltage position of this transition is perfectly correlated with bias dependent changes in STM images of the DB at different charge states. Density functional theory (DFT) calculations further highlight the role of subsurface dopants on DB properties by showing the influence of the DB-dopant distance on the DB state. We discuss possible theoretical models of electronic transport through the DB that could account for our experimental observations.


title: Deeply virtual Compton scattering at small-$x$ in future Electron - Ion Colliders

abstract: The study of exclusive processes in the future electron-ion ($eA$) colliders will be an important tool to investigate the QCD dynamics at high energies as they are in general driven by the gluon content of the target which is strongly subject to parton saturation effects. In this paper we compute the coherent and incoherent cross sections for the deeply virtual Compton scattering (DVCS) process relying on the color dipole approach and considering different models for the dipole - proton scattering amplitude. The dependencies of the cross sections with the energy, photon virtuality, nuclear mass number and squared momentum transfer are analysed in detail. We demonstrate that the ratio between the incoherent and coherent cross sections decreases at smaller values of $Q^2$ and increases at smaller values of $A$. We show that the coherent cross section dominates at small $t$ and exhibits the typical diffractive pattern, with the number of dips in the range $|t| \le 0.3$ GeV$^2$ increasing with the mass atomic number. Our results indicate that the position of the dips are independent of the model used to treat the dipole - proton interaction as well as of the center-of-mass energy.


title: Offloading on the Edge: Analysis and Optimization of Local Data Storage and Offloading in HetNets

abstract: The rapid increase in data traffic demand has overloaded existing cellular networks. Planned upgrades in the communication architecture (e.g. LTE), while helpful, are not expected to suffice to keep up with demand. As a result, extensive densification through small cells, caching content closer to or even at the device, and device-to-device (D2D) communications are seen as necessary components for future heterogeneous cellular networks to withstand the data crunch. Nevertheless, these options imply new CAPEX and OPEX costs, extensive backhaul support, contract plan incentives for D2D, and a number of interesting tradeoffs arise for the operator. In this paper, we propose an analytical model to explore how much local storage and communication through "edge" nodes could help offload traffic in various heterogeneous network (HetNet) setups and levels of user tolerance to delays. We then use this model to optimize the storage allocation and access mode of different contents as a tradeoff between user satisfaction and cost to the operator. Finally, we validate our findings through realistic simulations and show that considerable amounts of traffic can be offloaded even under moderate densification levels.


title: Representation of vector fields

abstract: A simple proof is given for the explicit formula which allows one to recover a $C^2-$smooth vector field $A=A(x)$ in $\mathbb{R}^3$, decaying at infinity, from the knowledge of its $\nabla \times A$ and $\nabla \cdot A$. The representation of $A$ as a sum of the gradient field and a divergence-free vector fields is derived from this formula. Similar results are obtained for a vector field in a bounded $C^2-$smooth domain.


title: Consistent Answers of Conjunctive Queries on Graphs

abstract: During the past decade, there has been an extensive investigation of the computational complexity of the consistent answers of Boolean conjunctive queries under primary key constraints. Much of this investigation has focused on self-join-free Boolean conjunctive queries. In this paper, we study the consistent answers of Boolean conjunctive queries involving a single binary relation, i.e., we consider arbitrary Boolean conjunctive queries on directed graphs. In the presence of a single key constraint, we show that for each such Boolean conjunctive query, either the problem of computing its consistent answers is expressible in first-order logic, or it is polynomial-time solvable, but not expressible in first-order logic.


title: An Anisotropic Landau-Lifschitz-Gilbert model of dissipation in qubits

abstract: We derive a microscopic model for dissipative dynamics in a system of mutually interacting qubits coupled to a thermal bath that generalises the dissipative model of Landau-Lifschitz-Gilbert to the case of anisotropic bath couplings. We show that the dissipation acts to bias the quantum trajectories towards a reduced phase space. This model applies to a system of superconducting flux qubits whose coupling to the environment is necessarily anisotropic. We study the model in the context of the D-Wave computing device and show that the form of environmental coupling in this case produces dynamics that are closely related to several models proposed on phenomenological grounds.


title: Integral equations and applications

abstract: The goal of this Section is to formulate some of the basic results on the theory of integral equations and mention some of its applications. The literature of this subject is very large. Proofs are not given due to the space restriction. The results are taken from the works mentioned in the references.


title: Rigidity and sharp stability estimates for hypersurfaces with constant and almost-constant nonlocal mean curvature

abstract: We prove that the boundary of a (not necessarily connected) bounded smooth set with constant nonlocal mean curvature is a sphere. More generally, and in contrast with what happens in the classical case, we show that the Lipschitz constant of the nonlocal mean curvature of such a boundary controls its $C^2$-distance from a single sphere. The corresponding stability inequality is obtained with a sharp decay rate.


title: Powerful radiative jets in super-critical accretion disks around non-spinning black holes

abstract: We describe a set of simulations of super-critical accretion onto a non-rotating supermassive BH. The accretion flow is radiation pressure dominated and takes the form of a geometrically thick disk with twin low-density funnels around the rotation axis. For accretion rates $\gtrsim 10 \dot M_{\rm Edd}$, there is sufficient gas in the funnel to make this region optically thick. Radiation from the disk first flows into the funnel, after which it accelerates the optically thick funnel gas along the axis. The resulting jet is baryon-loaded and has a terminal density-weighted velocity $\approx 0.3c$. Much of the radiative luminosity is converted into kinetic energy by the time the escaping gas becomes optically thin. For an observer viewing down the axis, the isotropic equivalent luminosity of total energy is as much as $10^{48}\,\rm erg\,s^{-1}$ for a $10^7 M_\odot$ BH accreting at $10^3$ Eddington. Therefore, energetically, the simulated jets are consistent with observations of the most powerful tidal disruption events, e.g., Swift J1644. The jet velocity is, however, too low to match the Lorentz factor $\gamma > 2$ inferred in J1644. There is no such conflict in the case of other tidal disruption events. Since favorably oriented observers see isotropic equivalent luminosities that are highly super-Eddington, the simulated models can explain observations of ultra-luminous X-ray sources, at least in terms of luminosity and energetics, without requiring intermediate mass black holes. Finally, since the simulated jets are baryon-loaded and have mildly relativistic velocities, they match well the jets observed in SS433. The latter are, however, more collimated than the simulated jets. This suggests that, even if magnetic fields are not important for acceleration, they may perhaps still play a role in confining the jet.


title: Approximation of the Scattering Amplitude using Nonsymmetric Saddle Point Matrices

abstract: In this paper we examine iterative methods for solving the forward ($A{\bf x}={\bf b}$) and adjoint ($A^{T}{\bf y}={\bf g}$) systems of linear equations used to approximate the scattering amplitude, defined by ${\bf g}^{T}{\bf x}={\bf y}^{T}{\bf b}$. Based on an idea first proposed by Gene Golub, we use a conjugate gradient-like iteration for a nonsymmetric saddle point matrix that is constructed so as to have a real positive spectrum. Numerical experiments show that this method is more consistent than known methods for computing the scattering amplitude such as GLSQR or QMR. We then demonstrate that when combined with known preconditioning techniques, the proposed method exhibits more rapid convergence than state-of-the-art iterative methods for nonsymmetric systems.


title: Space-Constrained Massive MIMO: Hitting the Wall of Favorable Propagation

abstract: The recent development of the massive multiple-input multiple-output (MIMO) paradigm, has been extensively based on the pursuit of favorable propagation: in the asymptotic limit, the channel vectors become nearly orthogonal and inter-user interference tends to zero [1]. In this context, previous studies have considered fixed inter-antenna distance, which implies an increasing array aperture as the number of elements increases. Here, we focus on a practical, space-constrained topology, where an increase in the number of antenna elements in a fixed total space imposes an inversely proportional decrease in the inter-antenna distance. Our analysis shows that, contrary to existing studies, inter-user interference does not vanish in the massive MIMO regime, thereby creating a saturation effect on the achievable rate.


title: SDSSJ14584479+3720215: A Benchmark JHK Blazar Light Curve from the 2MASS Calibration Scans

abstract: Active galactic nuclei (AGNs) are well-known to exhibit flux variability across a wide range of wavelength regimes, but the precise origin of the variability at different wavelengths remains unclear. To investigate the relatively unexplored near-IR variability of the most luminous AGNs, we conduct a search for variability using well sampled JHKs-band light curves from the 2MASS survey calibration fields. Our sample includes 27 known quasars with an average of 924 epochs of observation over three years, as well as one spectroscopically confirmed blazar (SDSSJ14584479+3720215) with 1972 epochs of data. This is the best-sampled NIR photometric blazar light curve to date, and it exhibits correlated, stochastic variability that we characterize with continuous auto-regressive moving average (CARMA) models. None of the other 26 known quasars had detectable variability in the 2MASS bands above the photometric uncertainty. A blind search of the 2MASS calibration field light curves for AGN candidates based on fitting CARMA(1,0) models (damped-random walk) uncovered only 7 candidates. All 7 were young stellar objects within the {\rho} Ophiuchus star forming region, five with previous X-ray detections. A significant {\gamma}-ray detection (5{\sigma}) for the known blazar using 4.5 years of Fermi photon data is also found. We suggest that strong NIR variability of blazars, such as seen for SDSSJ14584479+3720215, can be used as an efficient method of identifying previously-unidentified {\gamma}-ray blazars, with low contamination from other AGN.


title: More Analysis of Double Hashing for Balanced Allocations

abstract: With double hashing, for a key $x$, one generates two hash values $f(x)$ and $g(x)$, and then uses combinations $(f(x) +i g(x)) \bmod n$ for $i=0,1,2,...$ to generate multiple hash values in the range $[0,n-1]$ from the initial two. For balanced allocations, keys are hashed into a hash table where each bucket can hold multiple keys, and each key is placed in the least loaded of $d$ choices. It has been shown previously that asymptotically the performance of double hashing and fully random hashing is the same in the balanced allocation paradigm using fluid limit methods. Here we extend a coupling argument used by Lueker and Molodowitch to show that double hashing and ideal uniform hashing are asymptotically equivalent in the setting of open address hash tables to the balanced allocation setting, providing further insight into this phenomenon. We also discuss the potential for and bottlenecks limiting the use this approach for other multiple choice hashing schemes.


title: Emerging heterogeneities in Italian customs and comparison with nearby countries

abstract: In this work we apply techniques and modus operandi typical of Statistical Mechanics to a large dataset about key social quantifiers and compare the resulting behaviours of five European nations, namely France, Germany, Italy, Spain and Switzerland. The social quantifiers considered are $i.$ the evolution of the number of autochthonous marriages (i.e. between two natives) within a given territorial district and $ii.$ the evolution of the number of mixed marriages (i.e. between a native and an immigrant) within a given territorial district. Our investigations are twofold. From a theoretical perspective, we develop novel techniques, complementary to classical methods (e.g. historical series and logistic regression), in order to detect possible collective features underlying the empirical behaviours; from an experimental perspective, we evidence a clear outline for the evolution of the social quantifiers considered. The comparison between experimental results and theoretical predictions is excellent and allows speculating that France, Italy and Spain display a certain degree of {\em internal heterogeneity}, that is not found in Germany and Switzerland; such heterogeneity, quite mild in France and in Spain, is not negligible in Italy and highlights quantitative differences in the customs of Northern and Southern regions. These findings may suggest the persistence of two culturally distinct communities, long-term lasting heritages of different and well-established cultures.


title: Bell's theorem and the issue of determinism and indeterminism

abstract: The paper considers the claim that quantum theories with a deterministic dynamics of objects in ordinary space-time, such as Bohmian mechanics, contradict the assumption that the measurement settings can be freely chosen in the EPR experiment. That assumption is one of the premises of Bell's theorem. I first argue that only a premise to the effect that what determines the choice of the measurement settings is independent of what determines the past state of the measured system is needed for the derivation of Bell's theorem. Determinism as such does not undermine that independence (unless there are particular initial conditions of the universe that would amount to conspiracy). Only entanglement could do so. However, generic entanglement without collapse on the level of the universal wave function can go together with effective wave functions for subsystems of the universe, as in Bohmian mechanics. The paper argues that such effective wave functions are sufficient for the mentioned independence premise to hold.


title: Composite Photon Theory Versus Elementary Photon Theory

abstract: The purpose of this paper is to show that the composite photon theory measures up well against the Standard Model's elementary photon theory. This is done by comparing the two theories area by area. Although the predictions of quantum electrodynamics are in excellent agreement with experiment (as in the anomalous magnetic moment of the electron), there are some problems, such as the difficulty in describing the electromagnetic field with the four-component vector potential because the photon has only two polarization states. In most areas the two theories give similar results, so it is impossible to rule out the composite photon theory. Pryce's arguments in 1938 against a composite photon theory are shown to be invalid or irrelevant. Recently, it has been realized that in the composite theory the antiphoton does not interact with matter because it is formed of a neutrino and an antineutrino with the wrong helicity. This leads to experimental tests that can determine which theory is correct.


title: Generating the local oscillator "locally" in continuous-variable quantum key distribution based on coherent detection

abstract: Continuous-variable quantum key distribution (CV-QKD) protocols based on coherent detection have been studied extensively in both theory and experiment. In all the existing implementations of CV-QKD, both the quantum signal and the local oscillator (LO) are generated from the same laser and propagate through the insecure quantum channel. This arrangement may open security loopholes and also limit the potential applications of CV-QKD. In this paper, we propose and demonstrate a pilot-aided feedforward data recovery scheme which enables reliable coherent detection using a "locally" generated LO. Using two independent commercial laser sources and a spool of 25 km optical fiber, we construct a coherent communication system. The variance of the phase noise introduced by the proposed scheme is measured to be 0.04 (rad^2), which is small enough to enable secure key distribution. This technology also opens the door for other quantum communication protocols, such as the recently proposed measurement-device-independent (MDI) CV-QKD where independent light sources are employed by different users.


title: Three versions of categorical crossed-product duality

abstract: In this partly expository paper we compare three different categories of C*-algebras in which crossed-product duality can be formulated, both for actions and for coactions of locally compact groups. In these categories, the isomorphisms correspond to C*-algebra isomorphisms, imprimitivity bimodules, and outer conjugacies, respectively. In each case, a variation of the fixed-point functor that arises from classical Landstad duality is used to obtain a quasi-inverse for a crossed-product functor. To compare the various cases, we describe in a formal way our view of the fixed-point functor as an "inversion" of the process of forming a crossed product. In some cases, we obtain what we call "good" inversions, while in others we do not. For the outer-conjugacy categories, we generalize a theorem of Pedersen to obtain a fixed-point functor that is quasi-inverse to the reduced-crossed-product functor for actions, and we show that this gives a good inversion. For coactions, we prove a partial version of Pedersen's theorem that allows us to define a fixed-point functor, but the question of whether it is a quasi-inverse for the crossed-product functor remains open.


title: Constraining the Mass and the Non-Radial Drag Coefficient of a Solar Coronal Mass Ejection

abstract: Decades of observations show that CMEs can deflect from a purely radial trajectory yet no consensus exists as to the cause of these deflections. Many of theories attribute the CME deflection to magnetic forces. We developed ForeCAT (Kay et al. 2013, Kay et al. 2015), a model for CME deflections based solely on magnetic forces, neglecting any reconnection effects. Here we compare ForeCAT predictions to the observed deflection of the 2008 December 12 CME and find that ForeCAT can accurately reproduce the observations. Multiple observations show that this CME deflected nearly 30{\deg} in latitude (Byrne et al. 2010, Gui et al. 2011) and 4.4{\deg} in longitude (Gui et al. 2011). From the observations, we are able to constrain all of the ForeCAT input parameters (initial position, radial propagation speed, and expansion) except the CME mass and the drag coefficient that affects the CME motion. By minimizing the reduced chi-squared, $\chi^2_{\nu}$, between the ForeCAT results and the observations we determine an acceptable mass range between 4.5x10$^{14}$ and 1x10$^{15}$ g and the drag coefficient less than 1.4 with a best fit at 7.5x10$^{14}$ g and 0 for the mass and drag coefficient. ForeCAT is sensitive to the magnetic background and we are also able to constrain the rate at which the quiet sun magnetic field falls to be similar or to or fall slightly slower than the Potential Field Source Surface model.


title: Computing cobordism maps in link Floer homology and the reduced Khovanov TQFT

abstract: We study the maps induced on link Floer homology by elementary decorated link cobordisms. We compute these for births, deaths, stabilizations, and destabilizations, and show that saddle cobordisms can be computed in terms of maps in a decorated skein exact triangle that extends the oriented skein exact triangle in knot Floer homology. In particular, we completely determine the Alexander and Maslov grading shifts. As a corollary, we compute the maps induced by elementary cobordisms between unlinks. We show that these give rise to a $(1+1)$-dimensional TQFT that coincides with the reduced Khovanov TQFT. Hence, when applied to the cube of resolutions of a marked link diagram, it gives the complex defining the reduced Khovanov homology of the knot. Finally, we define a spectral sequence from (reduced) Khovanov homology using these cobordism maps, and we prove that it is an invariant of the (marked) link.


title: Homothety and Quasi Self Similarity in Asymptotically Flat and AdS Spacetimes

abstract: For a spherically symmetric self-gravitating scalar field we study self similar and quasi-self similar solutions in asymptotically flat and AdS spacetimes in various dimensions. Our main approach relies on reducing the Einstein-Klein-Gordon system of equations under the self similarity assumption to a dynamical system. Applying standard techniques in the analysis of dynamical system we study the corresponding phase diagrams, gaining some insight into the short and long-term evolution of the system. In the asymptotically flat case the fixed points of the dynamical system correspond to conformally flat spacetimes. We observed a new transition from a fixed point to a black hole and we compute the corresponding exponent analytically. In the case of AdS spacetime, we observed the quasi self similar scalar field scattering to infinity after interfering near the center corresponding to the dynamical system with no fixed points.


title: Minimal universal metric spaces

abstract: Let $\mathfrak{M}$ be a class of metric spaces. A metric space $Y$ is minimal $\mathfrak{M}$-universal if every $X\in\mathfrak{M}$ can be isometrically embedded in $Y$ but there are no proper subsets of $Y$ satisfying this property. We find conditions under which, for given metric space $X$, there is a class $\mathfrak{M}$ of metric spaces such that $X$ is minimal $\mathfrak{M}$-universal. We generalize the notion of minimal $\mathfrak{M}$-universal metric space to notion of minimal $\mathfrak{M}$-universal class of metric spaces and prove the uniqueness, up to an isomorphism, for these classes. The necessary and sufficient conditions under which the disjoint union of the metric spaces belonging to a class $\mathfrak{M}$ is minimal $\mathfrak{M}$-universal are found. Examples of minimal universal metric spaces are constructed for the classes of the three-point metric spaces and $n$-dimensional normed spaces. Moreover minimal universal metric spaces are found for some subclasses of the class of metric spaces $X$ which possesses the following property. Among every three distinct points of $X$ there is one point lying between the other two points.


title: The $R_{\infty}$-property for Chevalley groups of types $B_l, C_l, D_l$ over integral domains

abstract: We prove that Chevalley groups of the classical series $B_l, C_l, D_l$ over an integral domain of zero characteristic, which has torsion automorphism group, possess the $R_{\infty}$-property.


title: A Hebbian/Anti-Hebbian Neural Network for Linear Subspace Learning: A Derivation from Multidimensional Scaling of Streaming Data

abstract: Neural network models of early sensory processing typically reduce the dimensionality of streaming input data. Such networks learn the principal subspace, in the sense of principal component analysis (PCA), by adjusting synaptic weights according to activity-dependent learning rules. When derived from a principled cost function these rules are nonlocal and hence biologically implausible. At the same time, biologically plausible local rules have been postulated rather than derived from a principled cost function. Here, to bridge this gap, we derive a biologically plausible network for subspace learning on streaming data by minimizing a principled cost function. In a departure from previous work, where cost was quantified by the representation, or reconstruction, error, we adopt a multidimensional scaling (MDS) cost function for streaming data. The resulting algorithm relies only on biologically plausible Hebbian and anti-Hebbian local learning rules. In a stochastic setting, synaptic weights converge to a stationary state which projects the input data onto the principal subspace. If the data are generated by a nonstationary distribution, the network can track the principal subspace. Thus, our result makes a step towards an algorithmic theory of neural computation.


title: Valley properties of doped graphene in a magnetic field

abstract: The aim of this work is to describe the electronic properties of graphene in a constant magnetic field in the long wavelength approximation with random binary disorder, by solving the Soven equation self-consistently. Density of state contributions for different valleys in each sublattice sites are obtained for different values of magnetic field strength showing remarkable differences between K and K' valleys. A band gap is obtained by an asymmetric on-site impurity concentration and the graphene electrons acquire an anomalous magnetic moment, which is opposite in different valleys, which depend highly in the interplay between the impurity band, the band edges and the broadening of the Landau levels. In turn, magnetization as a function of B for different on-site random impurities is computed showing that by decreasing the on-site impurity energy values, maximum magnetization is shifted towards higher values of B which can be used to create and manipulate polarized valley currents. Finally, conductivity and local vertex function are obtained as a function of energy showing that scattering contributions from A and B sublattices differ significantly. Effective medium local two-irreducible vertex is computed showing that scattering from sublattices A and B do not contribute equally, which can be related to weak anti-localization. From these results, it could be possible to explore how the valley pseudospin can be used to create polarized currents by populating asymmetrically the sublattice sites, where the population can be tuned with the applied magnetic field strength.


title: Some singular minimizers in low dimensions in the calculus of variations

abstract: We construct a singular minimizing map ${\bf u}$ from $\mathbb{R}^3$ to $\mathbb{R}^2$ of a smooth uniformly convex functional of the form $\int_{B_1} F(D{\bf u})\,dx$.


title: Solitons and frequency combs in silica microring resonators: Interplay of the Raman and higher-order dispersion effects

abstract: The influence of Raman scattering and higher order dispersions on solitons and frequency comb generation in silica microring resonators is investigated. The Raman effect introduces a threshold value in the resonator quality factor above which the frequency locked solitons can not exist and, instead, a rich dynamics characterized by generation of self-frequency shift- ing solitons and dispersive waves is observed. A mechanism of broadening of the Cherenkov radiation through Hopf instability of the frequency locked solitons is also reported.


title: What Is an Emerging Technology?

abstract: There is considerable and growing interest in the emergence of novel technologies, especially from the policy-making perspective. Yet as an area of study, emerging technologies lacks key foundational elements, namely a consensus on what classifies a technology as 'emergent' and strong research designs that operationalize central theoretical concepts. The present paper aims to fill this gap by developing a definition of 'emerging technologies' and linking this conceptual effort with the development of a framework for the operationalisation of technological emergence. The definition is developed by combining a basic understanding of the term and in particular the concept of 'emergence' with a review of key innovation studies dealing with definitional issues of technological emergence. The resulting definition identifies five attributes that feature in the emergence of novel technologies. These are: (i) radical novelty, (ii) relatively fast growth, (iii) coherence, (iv) prominent impact, and (v) uncertainty and ambiguity. The framework for operationalising emerging technologies is then elaborated on the basis of the proposed attributes. To do so, we identify and review major empirical approaches (mainly in, although not limited to, the scientometric domain) for the detection and study of emerging technologies (these include indicators and trend analysis, citation analysis, co-word analysis, overlay mapping, and combinations thereof) and elaborate on how these can be used to operationalise the different attributes of emergence.


title: 5G Cellular: Key Enabling Technologies and Research Challenges

abstract: The evolving fifth generation (5G) cellular wireless networks are envisioned to provide higher data rates, enhanced end-user quality-of-experience (QoE), reduced end-to-end latency, and lower energy consumption. This article presents several emerging technologies, which will enable and define the 5G mobile communications standards. The major research problems, which these new technologies breed, as well as the measurement and test challenges for 5G systems are also highlighted.


title: The quantum field theory interpretation of quantum mechanics

abstract: It is shown that adopting the \emph{Quantum Field} ---extended entity in space-time build by dynamic appearance propagation and annihilation of virtual particles--- as the primary ontology the astonishing features of quantum mechanics can be rendered intuitive. This interpretation of quantum mechanics follows from the formalism of the most successful theory in physics: quantum field theory.


title: Inertial-Hall effect: the influence of rotation on the Hall conductivity

abstract: Inertial effects play an important role in classical mechanics but have been largely overlooked in quantum mechanics. Nevertheless, the analogy between inertial forces on mass particles and electromagnetic forces on charged particles is not new. In this paper, we consider a rotating non-interacting planar two-dimensional electron gas with a perpendicular uniform magnetic field and investigate the effects of the rotation in the Hall conductivi


title: Near-optimal quantum tomography: estimators and bounds

abstract: We give bounds on the average fidelity achievable by any quantum state estimator, which is arguably the most prominently used figure of merit in quantum state tomography. Moreover, these bounds can be computed online---that is, while the experiment is running. We show numerically that these bounds are quite tight for relevant distributions of density matrices. We also show that the Bayesian mean estimator is ideal in the sense of performing close to the bound without requiring optimization. Our results hold for all finite dimensional quantum systems.


title: Atom- photon entanglement beyond the multi-photon resonance condition

abstract: The density matrix equations of motion in near-degenerate three-level V-type closed-loop atomic system are calculated numerically in Floquet frame. The dynamical behavior of atom- photon entanglement between the dressed atom and its spontaneous emission is studied in semi classical approach beyond the two-photon resonance condition in such a system. The quantum entropy of these two subsystems is investigated by using the von Neumann entropy. It is shown that, the degree of entanglement measure (DEM) can be controlled via the intensity and the detuning of coupling optical field and quantum interference induced by spontaneous emission. Moreover in the absence of quantum interference the steady state behavior of DEM can be achieved even in beyond the two- photon resonance condition. Furthermore in the absence of quantum interference for special parameters of Rabi frequency and detuning of driving laser field disentanglement can be occurred. Also the electromagnetically induced transparency condition can be obtained when the system is disentangled.


title: Bulk photovoltaic effect enhancement via electrostatic control in layered ferroelectrics

abstract: The correlation between the shift current mechanism for the bulk photovoltaic effect (BPVE) and the structural and electronic properties of ferroelectric perovskite oxides is not well understood. Here, we study and engineer the shift current photovoltaic effect using a visible-light-absorbing ferroelectric Pb(Ni$_{x}$Ti$_{1-x}$)O$_{3-x}$ solid solution from first principles. We show that the covalent orbital character dicates the direction, magnitude, and onset energy of shift current in a predictable fashion. In particular, we find that the shift current response can be enhanced via electrostatic control in layered ferroelectrics, as bound charges face a stronger impetus to screen the electric field in a thicker material, delocalizing electron densities. This heterogeneous layered structure with alternative photocurrent generating and insulating layers is ideal for BPVE applications.


title: A Hebbian/Anti-Hebbian Network Derived from Online Non-Negative Matrix Factorization Can Cluster and Discover Sparse Features

abstract: Despite our extensive knowledge of biophysical properties of neurons, there is no commonly accepted algorithmic theory of neuronal function. Here we explore the hypothesis that single-layer neuronal networks perform online symmetric nonnegative matrix factorization (SNMF) of the similarity matrix of the streamed data. By starting with the SNMF cost function we derive an online algorithm, which can be implemented by a biologically plausible network with local learning rules. We demonstrate that such network performs soft clustering of the data as well as sparse feature discovery. The derived algorithm replicates many known aspects of sensory anatomy and biophysical properties of neurons including unipolar nature of neuronal activity and synaptic weights, local synaptic plasticity rules and the dependence of learning rate on cumulative neuronal activity. Thus, we make a step towards an algorithmic theory of neuronal function, which should facilitate large-scale neural circuit simulations and biologically inspired artificial intelligence.


title: Measurement of charged jet production cross sections and nuclear modification in p-Pb collisions at $\sqrt{s_\rm{NN}} = 5.02$ TeV

abstract: Charged jet production cross sections in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV measured with the ALICE detector at the LHC are presented. Using the anti-$k_{\rm T}$ algorithm, jets have been reconstructed in the central rapidity region from charged particles with resolution parameters $R = 0.2$ and $R = 0.4$. The reconstructed jets have been corrected for detector effects and the underlying event background. To calculate the nuclear modification factor, $R_{\rm pPb}$, of charged jets in p-Pb collisions, a pp reference was constructed by scaling previously measured charged jet spectra at $\sqrt{s} = 7$ TeV. In the transverse momentum range $20 \le p_{\rm T,ch\ jet} \le 120$ GeV/$c$, $R_{\rm pPb}$ is found to be consistent with unity, indicating the absence of strong nuclear matter effects on jet production. Major modifications to the radial jet structure are probed via the ratio of jet production cross sections reconstructed with the two different resolution parameters. This ratio is found to be similar to the measurement in pp collisions at $\sqrt{s} = 7$ TeV and to the expectations from PYTHIA pp simulations and NLO pQCD calculations at $\sqrt{s_{\rm NN}} = 5.02$ TeV.


title: Well-Posedness of the Nonlinear Unsteady Prandtl Equations with Robin Boundary Condition in Weighted Sobolev Spaces

abstract: In this paper, we study the well-posedness of classical solutions to the nonlinear unsteady Prandtl equations with Robin boundary condition in half space in weighted Sobolev spaces. We firstly investigate the monotonic shear flow with Robin boundary condition and the linearized Prandtl-type equations with Robin boundary condition in weighted Sobolev spaces. Due to the degeneracy of the Prandtl equations and the loss of regularity, we apply the Nash-Moser-Hormander iteration scheme to prove the existence of classical solutions to the nonlinear Prandtl equations with Robin boundary condition when the initial data is a small perturbation of a monotonic shear flow satisfying Robin boundary condition. The uniqueness and stability are also proved in the weighted Sobolev spaces. The nonlinear Prandtl equations with Robin boundary arise in the inviscid limit of incompressible Navier-Stokes equations with Navier-slip boundary condition for which the slip length is square root of viscosity. Our results are also valid for the Dirichlet boundary case.


title: Asymptotic state lumping in transport and diffusion problems on networks

abstract: One of the aims of systems biology is to build multiple layered and multiple scale models of living systems which can efficiently describe phenomena occurring at various level of resolution. Such models should consist of layers of various microsystems interconnected by a network of pathways, to form a macrosystem in a consistent way; that is, the observable characteristics of the macrosystem should be, at least asymptotically, derivable by aggregation of the appropriate features of the microsystems forming it, and from the properties of the network. In this paper we consider a general macromodel describing a population consisting of several interacting with each other subgroups, with the rules of interactions given by a system of ordinary differential equations, and we construct two different micromodels whose aggregated dynamics is approximately the same as that of the original macromodel. The micromodels offer a more detailed description of the original macromodel's dynamics by considering an internal structure of each subgroup. Here, each subgroup is represented by an edge of a graph with diffusion or transport occurring along it, while the interactions between the edges are described by interface conditions at the nodes joining them. We prove that with an appropriate scaling of such models, roughly speaking, with fast diffusion or transport combined with slow exchange at the nodes, the solutions of the micromodels are close to the solution to the macromodel.


title: First-Principles Calculation of the Bulk Photovoltaic Effect in KNbO$_{3}$ and (K,Ba)(Ni,Nb)O$_{3-\delta}$

abstract: The connection between noncentrosymmetric materials' structure, electronic structure, and bulk photovoltaic performance remains not well understood. In particular, it is still unclear which photovoltaic (PV) mechanism(s) are relevant for the recently demonstrated visible-light ferroelectric photovoltaic (K,Ba)(Ni,Nb)O$_{3-\delta}$. In this paper, we study the bulk photovoltaic effect (BPVE) of (K,Ba)(Ni,Nb)O$_{3-\delta}$ and KNbO$_{3}$ by calculating the shift current from first principles. The effects of structural phase, lattice distortion, oxygen vacancies, cation arrangement, composition, and strain on BPVE are systematically studied. We find that (K,Ba)(Ni,Nb)O$_{3-\delta}$ has a comparable BPVE with that of the broadly explored BiFeO$_{3}$, but for a much lower photon energy. In particular, the Glass coefficient of (K,Ba)(Ni,Nb)O$_{5}$ in a simply layered structure can be as large as 12 times that of BiFeO$_{3}$. Furthermore, the nature of the wavefunctions dictates the eventual shift current yield, which can be significantly affected and engineered by changing the O vacancy location, cation arrangement, and strain. This is not only helpful for understanding other PV mechanisms that relate to the motion of the photocurrent carriers, but also provides guidelines for the design and optimization of PV converters.


title: Magnetic field induced helical mode and topological transitions in a quasi-ballistic topological insulator nanoribbon with circumferentially quantized surface state sub-bands

abstract: A topological insulator (TI) nanowire (NW), where the core is insulating and the 2D spin-helical Dirac fermion topological surface states (TSS) are circumferentially quantized into a series of 1D sub-bands, promises novel topological physics and applications. An axial magnetic flux ({\Phi}) through the core drives periodic topological transitions in the surface sub-bands, changing from being all doubly-degenerate with a gapped Dirac point (DP) at integer (including zero) flux quanta ({\Phi}$_0$ = h/e, with h being the Planck constant and e the electron charge), to contain a topologically-protected, non-degenerate 1D spin helical mode with restored DP at half-integer flux quanta. The resulting magnetoconductance is predicted to exhibit Aharonov-Bohm oscillations (ABO) with maxima occurring alternatively at half-integer or integer flux quanta (referred to as {\pi}-ABO or 0-ABO), depending periodically on the Fermi wavevector (k$_F$, with period 2{\pi}/C, C being the NW circumference). Here, we report a clear observation of such k$_F$-periodic alternations between 0-ABO and {\pi}-ABO in Bi$_2$Te$_3$ TI nanoribbon (NR, a rectangular cross sectional NW) field effect devices, which exhibit quasi-ballistic transport over ~2 {\mu}m (as manifested in length-independent conductance, exponential decaying ABO amplitude with increasing temperature (T), and an 1/T dependence of the extracted phase coherence length). The conductances as functions of the gate voltage at half and zero flux quanta also exhibit clear, but anti-correlated oscillations periodic in k$_F$ (with period 2{\pi}/C, equivalently when C encloses an integer multiples of Fermi wavelength 2{\pi}/k$_F$), consistent with the circumferentially quantized surface sub-bands. We also extract the minimal Fermi energy and momentum for TSS to emerge out of the bulk valence band, in agreement with the known Bi$_2$Te$_3$ bandstructure.


title: Semigroup approach to diffusion and transport problems on networks

abstract: Models describing transport and diffusion processes occurring along the edges of a graph and interlinked by its vertices have been recently receiving a considerable attention. In this paper we generalize such models and consider a network of transport or diffusion operators defined on one dimensional domains and connected through boundary conditions linking the end-points of these domains in an arbitrary way (not necessarily as the edges of a graph are connected). We prove the existence of $C_0$-semigroups solving such problems and provide conditions fully characterizing when they are positive.


title: A review of mean-shift algorithms for clustering

abstract: A natural way to characterize the cluster structure of a dataset is by finding regions containing a high density of data. This can be done in a nonparametric way with a kernel density estimate, whose modes and hence clusters can be found using mean-shift algorithms. We describe the theory and practice behind clustering based on kernel density estimates and mean-shift algorithms. We discuss the blurring and non-blurring versions of mean-shift; theoretical results about mean-shift algorithms and Gaussian mixtures; relations with scale-space theory, spectral clustering and other algorithms; extensions to tracking, to manifold and graph data, and to manifold denoising; K-modes and Laplacian K-modes algorithms; acceleration strategies for large datasets; and applications to image segmentation, manifold denoising and multivalued regression.


title: Photoplethysmography-Based Heart Rate Monitoring in Physical Activities via Joint Sparse Spectrum Reconstruction

abstract: Goal: A new method for heart rate monitoring using photoplethysmography (PPG) during physical activities is proposed. Methods: It jointly estimates spectra of PPG signals and simultaneous acceleration signals, utilizing the multiple measurement vector model in sparse signal recovery. Due to a common sparsity constraint on spectral coefficients, the method can easily identify and remove spectral peaks of motion artifact (MA) in PPG spectra. Thus, it does not need any extra signal processing modular to remove MA as in some other algorithms. Furthermore, seeking spectral peaks associated with heart rate is simplified. Results: Experimental results on 12 PPG datasets sampled at 25 Hz and recorded during subjects' fast running showed that it had high performance. The average absolute estimation error was 1.28 beat per minute and the standard deviation was 2.61 beat per minute. Conclusion and Significance: These results show that the method has great potential to be used for PPG-based heart rate monitoring in wearable devices for fitness tracking and health monitoring.


title: Degenerate Hessian structures on radiant manifolds

abstract: We present a rigorous mathematical treatment of Ruppeiner geometry, by considering degenerate Hessian metrics defined on radiant manifolds. A manifold $M$ is said to be radiant if it is endowed with a symmetric, flat connection $\bar\nabla$ and a global vector field $\rho$ whose covariant derivative is the identity mapping. A degenerate Hessian metric on $M$ is a degenerate metric tensor $g$ that can locally be written as the covariant Hessian of a function, called potential. A function on $M$ is said to be extensive if its Lie derivative with respect to $\rho$ is the function itself. We show that the Hessian metrics appearing in equilibrium thermodynamics are necessarily degenerate, owing to the fact that their potentials are extensive (up to an additive constant). Manifolds having degenerate Hessian metrics always contain embedded Hessian submanifolds, which generalize the manifolds defined by constant volume in which Ruppeiner geometry is usually studied. By means of examples, we illustrate that linking scalar curvature to microscopic interactions within a thermodynamic system is inaccurate under this approach. In contrast, thermodynamic critical points seem to arise as geometric singularities.


title: A Hebbian/Anti-Hebbian Network for Online Sparse Dictionary Learning Derived from Symmetric Matrix Factorization

abstract: Olshausen and Field (OF) proposed that neural computations in the primary visual cortex (V1) can be partially modeled by sparse dictionary learning. By minimizing the regularized representation error they derived an online algorithm, which learns Gabor-filter receptive fields from a natural image ensemble in agreement with physiological experiments. Whereas the OF algorithm can be mapped onto the dynamics and synaptic plasticity in a single-layer neural network, the derived learning rule is nonlocal - the synaptic weight update depends on the activity of neurons other than just pre- and postsynaptic ones - and hence biologically implausible. Here, to overcome this problem, we derive sparse dictionary learning from a novel cost-function - a regularized error of the symmetric factorization of the input's similarity matrix. Our algorithm maps onto a neural network of the same architecture as OF but using only biologically plausible local learning rules. When trained on natural images our network learns Gabor-filter receptive fields and reproduces the correlation among synaptic weights hard-wired in the OF network. Therefore, online symmetric matrix factorization may serve as an algorithmic theory of neural computation.


title: Simplified dark matter top-quark interactions at the LHC

abstract: Stringent limits on the interactions between dark matter (DM) and the standard model can be set by studying how initial-state or final-state particles recoil against missing transverse energy (MET). In this work, we improve, extend and correct LHC constraints on the interactions between DM and top quarks that are mediated by the exchange of spin-0 s-channel resonances. A comparison of the LHC run-1 sensitivity of the two main search channels is presented, which shows that mono-jet searches are typically more restrictive than the MET + tbar t searches. We furthermore explore the reach of the 14 TeV LHC. The collider constraints are compared to the restrictions arising from direct and indirect detection as well as the DM relic abundance, and we also reflect on effective field theory interpretations of the LHC exclusions.


title: One of the closest exoplanet pairs to the 3:2 Mean Motion Resonance: K2-19b \& c

abstract: The K2 mission has recently begun to discover new and diverse planetary systems. In December 2014 Campaign 1 data from the mission was released, providing high-precision photometry for ~22000 objects over an 80 day timespan. We searched these data with the aim of detecting further important new objects. Our search through two separate pipelines led to the independent discovery of K2-19b \& c, a two-planet system of Neptune sized objects (4.2 and 7.2 $R_\oplus$), orbiting a K dwarf extremely close to the 3:2 mean motion resonance. The two planets each show transits, sometimes simultaneously due to their proximity to resonance and alignment of conjunctions. We obtain further ground based photometry of the larger planet with the NITES telescope, demonstrating the presence of large transit timing variations (TTVs), and use the observed TTVs to place mass constraints on the transiting objects under the hypothesis that the objects are near but not in resonance. We then statistically validate the planets through the \texttt{PASTIS} tool, independently of the TTV analysis.


title: Bayesian Optimization of Text Representations

abstract: When applying machine learning to problems in NLP, there are many choices to make about how to represent input texts. These choices can have a big effect on performance, but they are often uninteresting to researchers or practitioners who simply need a module that performs well. We propose an approach to optimizing over this space of choices, formulating the problem as global optimization. We apply a sequential model-based optimization technique and show that our method makes standard linear models competitive with more sophisticated, expensive state-of-the-art methods based on latent variable models or neural networks on various topic classification and sentiment analysis problems. Our approach is a first step towards black-box NLP systems that work with raw text and do not require manual tuning.


title: Consistent Probabilistic Social Choice

abstract: Two fundamental axioms in social choice theory are consistency with respect to a variable electorate and consistency with respect to components of similar alternatives. In the context of traditional non-probabilistic social choice, these axioms are incompatible with each other. We show that in the context of probabilistic social choice, these axioms uniquely characterize a function proposed by Fishburn (Rev. Econ. Stud., 51(4), 683--692, 1984). Fishburn's function returns so-called maximal lotteries, i.e., lotteries that correspond to optimal mixed strategies of the underlying plurality game. Maximal lotteries are guaranteed to exist due to von Neumann's Minimax Theorem, are almost always unique, and can be efficiently computed using linear programming.


title: Predicted Impacts of Proton Temperature Anisotropy on Solar Wind Turbulence

abstract: Particle velocity distributions measured in the weakly collisional solar wind are frequently found to be non-Maxwellian, but how these non-Maxwellian distributions impact the physics of plasma turbulence in the solar wind remains unanswered. Using numerical solutions of the linear dispersion relation for a collisionless plasma with a bi-Maxwellian proton velocity distribution, we present a unified framework for the four proton temperature anisotropy instabilities, identifying the associated stable eigenmodes, highlighting the unstable region of wavevector space, and presenting the properties of the growing eigenfunctions. Based on physical intuition gained from this framework, we address how the proton temperature anisotropy impacts the nonlinear dynamics of the \Alfvenic fluctuations underlying the dominant cascade of energy from large to small scales and how the fluctuations driven by proton temperature anisotropy instabilities interact nonlinearly with each other and with the fluctuations of the large-scale cascade. We find that the nonlinear dynamics of the large-scale cascade is insensitive to the proton temperature anisotropy, and that the instability-driven fluctuations are unlikely to cause significant nonlinear evolution of either the instability-driven fluctuations or the turbulent fluctuations of the large-scale cascade.


title: Origin of the enhancement of tunneling probability in the nearly integrable system

abstract: The enhancement of tunneling probability in the nearly integrable system is closely examined, focusing on tunneling splittings plotted as a function of the inverse of the Planck's constant. On the basis of the analysis using the absorber which efficiently suppresses the coupling creating spikes in the plot, we found that the splitting curve should be viewed as the staircase-shaped skeleton accompanied by spikes. We further introduce renormalized integrable Hamiltonians, and explore the origin of such a staircase structure by investigating the nature of eigenfunctions closely. It is found that the origin of the staircase structure could trace back to the anomalous structure in tunneling tail which manifests itself in the representation using renormalized action bases. This also explains the reason why the staircase does not appear in the completely integrable system.


title: Millimeter Wave Cellular Networks: A MAC Layer Perspective

abstract: The millimeter wave (mmWave) frequency band is seen as a key enabler of multi-gigabit wireless access in future cellular networks. In order to overcome the propagation challenges, mmWave systems use a large number of antenna elements both at the base station and at the user equipment, which lead to high directivity gains, fully-directional communications, and possible noise-limited operations. The fundamental differences between mmWave networks and traditional ones challenge the classical design constraints, objectives, and available degrees of freedom. This paper addresses the implications that highly directional communication has on the design of an efficient medium access control (MAC) layer. The paper discusses key MAC layer issues, such as synchronization, random access, handover, channelization, interference management, scheduling, and association. The paper provides an integrated view on MAC layer issues for cellular networks, identifies new challenges and tradeoffs, and provides novel insights and solution approaches.


title: Fault Analysis Using Gegenbauer Multiresolution Analysis

abstract: This paper exploits the multiresolution analysis in the fault analysis on transmission lines. Faults were simulated using the ATP (Alternative Transient Program), considering signals at 128/cycle. A nonorthogonal multiresolution analysis was provided by Gegenbauer scaling and wavelet filters. In the cases where the signal reconstruction is not required, orthogonality may be immaterial. Gegenbauer filter banks are thereby offered in this paper as a tool for analyzing fault signals on transmission lines. Results are compared to those ones derived from a 4-coefficient Daubechies filter. The main advantages in favor of Gegenbauer filters are their smaller computational effort and their constant group delay, as they are symmetric filters.


title: On the Possibility of Measuring the Single-tagged Exclusive Jets at the LHC

abstract: The feasibility studies of the measurement of the central exclusive jet production at the LHC using the proton tagging technique are presented. In order to reach the low jet-mass region, single tagged events were considered. The studies were performed at the c.m. energy of 14 TeV and the ATLAS detector, but are also applicable for the CMS-TOTEM experiments. Four data-taking scenarios were considered: AFP and ALFA detectors as forward proton taggers and $\beta^*$ = 0.55 m and $\beta^*$ = 90 m optics. After the event selection, the signal-to-background ratio ranges between 5 and $10^4$. Finally, the expected precision of the central exclusive dijet cross-section measurement for data collection period of 100 h is estimated.


title: An Efficient User-Side Nulling Calibration for Quantum Annealing Computers

abstract: (Withdrawn) This work describes an efficient user-side method of calibrating and correcting quantum annealing computers. For quantum annealing computers based on the Ising model, the method measures the residual bias of the h and J coefficients. Once measured, these biases can then be nulled in subsequent runs for any problem of interest. This method also returns a temperature for each qubit based on the measured versus the expected qubit distributions computed from a Boltzmann distribution model.


title: Better Than Earth

abstract: Do We Inhabit The Best O All Possible Worlds? German mathematician Gottfried Leibniz thought so, writing in 1710 that our planet, warts and all, must be the most optimal one imaginable. Leibniz's idea was roundly scorned as unscientific wishful thinking, most notably by French author Voltaire in his magnum opus, Candide. Yet Leibniz might find sympathy from at least one group of scientists - the astronomers who have for decades treated Earth as a golden standard as they search for worlds beyond our own solar system. Because earthlings still know of just one living world - our own - it makes some sense to use Earth as a template in the search for life elsewhere, such as in the most Earth-like regions of Mars or Jupiter's watery moon Europa. Now, however, discoveries of potentially habitable planets orbiting stars other than our sun - exoplanets, that is - are challenging that geocentric approach.


title: Inflation from Asymptotically Safe Theories

abstract: We investigate models in which inflation is driven by an ultraviolet safe and interacting scalar sector stemming from a new class of nonsupersymmetric gauge field theories. These new theories, differently from generic scalar models, are well defined to arbitrary short distances because of the existence of a controllable ultraviolet interacting fixed point. The scalar couplings at the ultraviolet fixed point and their overall running are predicted by the geometric structure of the underlying theory. We analyse the minimal and non-minimal coupling to gravity of these theories and the consequences for inflation. In the minimal coupling case the theory requires large non-perturbative quantum corrections to the quantum potential for the theory to agree with data, while in the non- minimal coupling case the perturbative regime in the couplings of the theory is preferred. Requiring the theory to reproduce the observed amplitude of density perturbations constrain the geometric data of the theory such as the number of colors and flavors for generic values of the non-minimal coupling.


title: Absorption Coefficients of the Methane-Nitrogen Binary Ice System: Implications for Pluto

abstract: The methane-nitrogen phase diagram of Prokhvatilov and Yantsevich (1983) indicates that at temperatures relevant to the surfaces of icy dwarf planets like Pluto, two phases contribute to the methane absorptions: nitrogen saturated with methane $\bf{\bar{N_{2}}}$:CH$_{4}$ and methane saturated with nitrogen $\bf{\bar{CH_{4}}}$:N$_{2}$. No optical constants are available so far for the latter component limiting construction of a proper model, in compliance with thermodynamic equilibrium considerations. New optical constants for solid solutions of methane diluted in nitrogen (N$_{2}$:CH$_{4}$) and nitrogen diluted in methane (CH$_{4}$:N$_{2}$) are presented at temperatures between 40 and 90 K, in the wavelength range 1.1-2.7 $\mu$m at different mixing ratios. These optical constants are derived from transmission measurements of crystals grown from the liquid phase in closed cells. A systematic study of the changes of methane and nitrogen solid mixtures spectral behavior with mixing ratio and temperature is presented.


title: Polynomial kernelization for removing induced claws and diamonds

abstract: A graph is called (claw,diamond)-free if it contains neither a claw (a $K_{1,3}$) nor a diamond (a $K_4$ with an edge removed) as an induced subgraph. Equivalently, (claw,diamond)-free graphs can be characterized as line graphs of triangle-free graphs, or as linear dominoes, i.e., graphs in which every vertex is in at most two maximal cliques and every edge is in exactly one maximal clique. In this paper we consider the parameterized complexity of the (claw,diamond)-free Edge Deletion problem, where given a graph $G$ and a parameter $k$, the question is whether one can remove at most $k$ edges from $G$ to obtain a (claw,diamond)-free graph. Our main result is that this problem admits a polynomial kernel. We complement this finding by proving that, even on instances with maximum degree $6$, the problem is NP-complete and cannot be solved in time $2^{o(k)}\cdot |V(G)|^{O(1)}$ unless the Exponential Time Hypothesis fail


title: Measuring the magnetic moment density in patterned ultrathin ferromagnets with submicron resolution

abstract: We present a new approach to infer the surface density of magnetic moments $I_s$ in ultrathin ferromagnetic films with perpendicular anisotropy. It relies on quantitative stray field measurements with an atomic-size magnetometer based on the nitrogen-vacancy center in diamond. The method is applied to microstructures patterned in a 1-nm-thick film of CoFeB. We report measurements of $I_s$ with a few percent uncertainty and a spatial resolution in the range of $(100$ nm)$^2$, an improvement by several orders of magnitude over existing methods. As an example of application, we measure the modifications of $I_s$ induced by local irradiation with He$^+$ ions in an ultrathin ferromagnetic wire. This method offers a new route to study variations of magnetic properties at the nanoscale.


title: On the predictability of rogue events

abstract: Using experimental data from three different rogue wave supporting systems, determinism and predictability of the underlying dynamics are evaluated with methods of nonlinear time series analysis. We included original records from the Draupner platform in the North Sea as well as time series from two optical systems in our analysis. One of the latter was measured in the infrared tail of optical fiber supercontinua, the other in the fluence profiles of multifilaments. All three data sets exhibit extreme-value statistics and exceed the significant wave height in the respective system by a factor larger than two. Nonlinear time series analysis indicates a different degree of determinism in the systems. The optical fiber scenario is found to be driven by quantum noise whereas rogue waves emerge as a consequence of turbulence in the others. With the large number of rogue events observed in the multifilament system, we can systematically explore the predictability of such events in a turbulent system. We observe that rogue events do not necessarily appear without a warning, but are often preceded by a short phase of relative order. This surprising finding sheds some new light on the fascinating phenomenon of rogue waves.


title: Adaptive Kinetic-Fluid Solvers for Heterogeneous Computing Architectures

abstract: We show feasibility and benefits of porting an adaptive multi-scale kinetic-fluid code to CPU-GPU systems. Challenges are due to the irregular data access for adaptive Cartesian mesh, vast difference of computational cost between kinetic and fluid cells, and desire to evenly load all CPUs and GPUs. Our Unified Flow Solver (UFS) combines Adaptive Mesh Refinement (AMR) with automatic cell-by-cell selection of kinetic or fluid solvers based on continuum breakdown criteria. Using GPUs enables hybrid simulations of mixed rarefied-continuum flows with a million of Boltzmann cells with 24x24x24 velocity mesh. We describe the implementation of CUDA kernels for three modules in UFS: the direct Boltzmann solver using discrete velocity method, the Direct Simulation Monte Carlo (DSMC) solver, and a mesoscopic solver based on Lattice Boltzmann Method, all using octree Cartesian mesh. Double digit speedups on single GPU and good scaling for multi-GPUs have been demonstrated.


title: Instanton-noninstanton transition in nonintegrable tunneling processes: A renormalized perturbation approach

abstract: The instanton-noninstanton (I-NI) transition in the tunneling process, which has been numerically observed in classically nonintegrable quantum maps, can be described by a perturbation theory based on an integrable Hamiltonian renormalized so as to incorporate the integrable part of the map. The renormalized perturbation theory is successfully applied to the two quantum maps, the H\'enon and standard maps. In spite of different nature of tunneling in the two systems, the I-NI transition exhibits very common characteristics. In particular, the manifestation of I-NI transition is obviously explained by a remarkable quenching of the renormalized transition matrix element. The enhancement of tunneling probability after the transition can be understood as a sudden change of the tunneling mechanism from the instanton to quite a different mechanism supported by classical flows just outside of the stable-unstable manifolds of the saddle on the top of the potential barrier.


title: Some new insights into information decomposition in complex systems based on common information

abstract: We take a closer look at the structure of bivariate dependency induced by a pair of predictor random variables $(X_1, X_2)$ trying to synergistically, redundantly or uniquely encode a target random variable $Y$. We evaluate a recently proposed measure of redundancy based on the G\'acs-K\"{o}rner common information (Griffith et al., Entropy 2014) and show that the measure, in spite of its elegance is degenerate for most non-trivial distributions. We show that Wyner's common information also fails to capture the notion of redundancy as it violates an intuitive monotonically non-increasing property. We identify a set of conditions when a conditional version of G\'acs and K\"{o}rner's common information is an ideal measure of unique information. Finally, we show how the notions of approximately sufficient statistics and conditional information bottleneck can be used to quantify unique information.


title: Cataland: Why the Fuss?

abstract: The three main objects in noncrossing Catalan combinatorics associated to a finite Coxeter system are noncrossing partitions, clusters, and sortable elements. The first two of these have known Fuss-Catalan generalizations. We provide new viewpoints for both and introduce the missing generalization of sortable elements by lifting the theory from the Coxeter system to the associated positive Artin monoid. We show how this new perspective ties together all three generalizations, providing a uniform framework for noncrossing Fuss-Catalan combinatorics. Having developed the combinatorial theory, we provide an interpretation of our generalizations in the language of the representation theory of hereditary Artin algebras.


title: Smelling out Code Clones: Clone Detection Tool Evaluation and Corresponding Challenges

abstract: Software clones have been an active area of research for the past two decades. However, although numerous clone detection tools are now available, only a small fraction of the literature has focused on tool evaluation, and this is in fact still an open problem. This is mostly due to the fact that standard information retrieval metrics such as recall and precision require a priori knowledge of clones already in the system. Detection tools also typically have a large number of parameters which are difficult to fine-tune for optimal performance on a particular software system, and different outputs produced by different tools add to the complexity of comparing one tool to another. In this review, we further explore the reasons why tool evaluation is still an open challenge, and present the current tools and frameworks targeted at mitigating these problems, focusing on the current standard benchmarks used to evaluate modern clone detection tools, and also presenting a recent method aimed at finding optimal tool configurations.


title: Emergent Spacetime and Cosmic Inflation

abstract: We present a novel background-independent framework for cosmic inflation, starting with a matrix model. In this framework, inflation is portrayed as a dynamic process responsible for the generation of both space and time. This stands in contrast to conventional inflation, which is characterized as a mere (exponential) expansion of an already existing spacetime, driven by the vacuum energy associated with an inflaton field. We observe that the cosmic inflation is triggered by the condensate of Planck energy into vacuum and responsible for the dynamical emergence of spacetime. The emergent spacetime picture admits a background-independent formulation so that the inflation is described by a conformal Hamiltonian system which requires neither an inflaton field nor an ad hoc inflation potential. This implies that the emergent spacetime may incapacitate all the rationales to introduce the multiverse hypothesis.


title: The ideal energy of classical lattice dynamics

abstract: We define, as local quantities, the least energy and momentum allowed by quantum mechanics and special relativity for physical realizations of some classical lattice dynamics. These definitions depend on local rates of finite-state change. In two example dynamics, we see that these rates evolve like classical mechanical energy and momentum.


title: On the consistency of tachyon warm inflation with viscous pressure

abstract: We obtain conditions for the existence of an attractor in the system of equations describing a tachyon warm inflationary model with bulk viscosity taken into account. When these conditions are met the evolution approaches slow-roll regime. We present the primordial power spectrum for the tachyon field by considering a dissipation coefficient depending on the scalar field and temperature.


title: Irreducible components of the space of foliations by surfaces

abstract: Let $\mathcal{F}$ be written as $ f^{*}(\mathcal{G})$, where $\mathcal{G}$ is a $1$-dimensional foliation on $ {\mathbb P^{n-1}}$ and $f:{\mathbb P^n}--->{\mathbb P^{n-1}}$ a non-linear generic rational map. We use local stability results of singular holomorphic foliations, to prove that: if $n\geq 4$, a foliation $\mathcal{F}$ by complex surfaces on $\mathbb P^n$ is globally stable under holomorphic deformations. As a consequence, we obtain irreducible components for the space of two-dimensional foliations in $\mathbb P^n$. We present also a result which characterizes holomorphic foliations on ${\mathbb P^n}, n\geq 4$ which can be obtained as a pull back of 1- foliations in ${\mathbb P^{n-1}}$ of degree $d\geq2$.


title: 3D Printing Meets Computational Astrophysics: Deciphering the Structure of Eta Carinae's Inner Colliding Winds

abstract: We present the first 3D prints of output from a supercomputer simulation of a complex astrophysical system, the colliding stellar winds in the massive (>120 M_Sun), highly eccentric (e ~ 0.9) binary star system Eta Carinae. We demonstrate the methodology used to incorporate 3D interactive figures into a PDF journal publication and the benefits of using 3D visualization and 3D printing as tools to analyze data from multidimensional numerical simulations. Using a consumer-grade 3D printer (MakerBot Replicator 2X), we successfully printed 3D smoothed particle hydrodynamics (SPH) simulations of Eta Carinae's inner (r ~ 110 au) wind-wind collision interface at multiple orbital phases. The 3D prints and visualizations reveal important, previously unknown 'finger-like' structures at orbital phases shortly after periastron (phi ~ 1.045) that protrude radially outward from the spiral wind-wind collision region. We speculate that these fingers are related to instabilities (e.g. thin-shell, Rayleigh-Taylor) that arise at the interface between the radiatively-cooled layer of dense post-shock primary-star wind and the fast (3000 km/s), adiabatic post-shock companion-star wind. The success of our work and easy identification of previously unrecognized physical features highlight the important role 3D printing and interactive graphics can play in the visualization and understanding of complex 3D time-dependent numerical simulations of astrophysical phenomena.


title: Anonymous broadcasting of classical information with a continuous-variable topological quantum code

abstract: Broadcasting information anonymously becomes more difficult as surveillance technology improves, but remarkably, quantum protocols exist that enable provably traceless broadcasting. The difficulty is making scalable entangled resource states that are robust to errors. We propose an anonymous broadcasting protocol that uses a continuous-variable surface-code state that can be produced using current technology. High squeezing enables large transmission bandwidth and strong anonymity, and the topological nature of the state enables local error mitigation.


title: Perturbative extraction of gravitational waveforms generated with Numerical Relativity

abstract: We derive an analytical expression for extracting the gravitational waveforms at null infinity using the Weyl scalar $\psi_4$ measured at a finite radius. Our expression is based on a series solution in orders of 1/r to the equations for gravitational perturbations about a spinning black hole. We compute this expression to order $1/r^2$ and include the spin parameter $a$ of the Kerr background. We test the accuracy of this extraction procedure by measuring the waveform for a merging black-hole binary at ten different extraction radii (in the range r/M=75-190) and for three different resolutions in the convergence regime. We find that the extraction formula provides a set of values for the radiated energy and momenta that at finite extraction radii converges towards the expected values with increasing resolution, which is not the case for the `raw' waveform at finite radius. We also examine the phase and amplitude errors in the waveform as a function of observer location and again observe the benefits of using our extraction formula. The leading corrections to the phase are ${\cal O}(1/r)$ and to the amplitude are ${\cal O}(1/r^2)$. This method provides a simple and practical way of estimating the waveform at infinity, and may be especially useful for scenarios such as well separated binaries, where the radiation zone is far from the sources, that would otherwise require extended simulation grids in order to extrapolate the `raw' waveform to infinity. Thus this method saves important computational resources and provides an estimate of errors.


title: Mid-J CO Shock Tracing Observations of Infrared Dark Clouds I

abstract: Infrared dark clouds (IRDCs) are dense, molecular structures in the interstellar medium that can harbour sites of high-mass star formation. IRDCs contain supersonic turbulence, which is expected to generate shocks that locally heat pockets of gas within the clouds. We present observations of the CO J = 8-7, 9-8, and 10-9 transitions, taken with the Herschel Space Observatory, towards four dense, starless clumps within IRDCs (C1 in G028.37+00.07, F1 and F2 in G034.43+0007, and G2 in G034.77-0.55). We detect the CO J = 8-7 and 9-8 transitions towards three of the clumps (C1, F1, and F2) at intensity levels greater than expected from photodissociation region (PDR) models. The average ratio of the 8-7 to 9-8 lines is also found to be between 1.6 and 2.6 in the three clumps with detections, significantly smaller than expected from PDR models. These low line ratios and large line intensities strongly suggest that the C1, F1, and F2 clumps contain a hot gas component not accounted for by standard PDR models. Such a hot gas component could be generated by turbulence dissipating in low velocity shocks.


title: A Comprehensive Archival Search for Counterparts to Ultra-Compact High Velocity Clouds: Five Local Volume Dwarf Galaxies

abstract: We report five Local Volume dwarf galaxies (two of which are presented here for the first time) uncovered during a comprehensive archival search for optical counterparts to ultra-compact high velocity clouds (UCHVCs). The UCHVC population of HI clouds are thought to be candidate gas-rich, low mass halos at the edge of the Local Group and beyond, but no comprehensive search for stellar counterparts to these systems has been presented. Careful visual inspection of all publicly available optical and ultraviolet imaging at the position of the UCHVCs revealed six blue, diffuse counterparts with a morphology consistent with a faint dwarf galaxy beyond the Local Group. Optical spectroscopy of all six candidate dwarf counterparts show that five have an H$\alpha$-derived velocity consistent with the coincident HI cloud, confirming their association, the sixth diffuse counterpart is likely a background object. The size and luminosity of the UCHVC dwarfs is consistent with other known Local Volume dwarf irregular galaxies. The gas fraction ($M_{HI}/M_{star}$) of the five dwarfs are generally consistent with that of dwarf irregular galaxies in the Local Volume, although ALFALFA-Dw1 (associated with ALFALFA UCHVC HVC274.68+74.70$-$123) has a very high $M_{HI}/M_{star}$$\sim$40. Despite the heterogenous nature of our search, we demonstrate that the current dwarf companions to UCHVCs are at the edge of detectability due to their low surface brightness, and that deeper searches are likely to find more stellar systems. If more sensitive searches do not reveal further stellar counterparts to UCHVCs, then the dearth of such systems around the Local Group may be in conflict with $\Lambda$CDM simulations.


title: The Black Ring is Unstable

abstract: We study non-axisymmetric linearised gravitational perturbations of the Emparan-Reall black ring using numerical methods. We find an unstable mode whose onset lies within the "fat" branch of the black ring and continues into the "thin" branch. Together with previous results using Penrose inequalities that fat black rings are unstable, this provides numerical evidence that the entire black ring family is unstable.


title: The Role of Bulge Formation in the Homogenization of Stellar Populations at $z\sim2$ as revealed by Internal Color Dispersion in CANDELS

abstract: We use data from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey to study how the spatial variation in the stellar populations of galaxies relate to the formation of galaxies at $1.5 < z < 3.5$. We use the Internal Color Dispersion (ICD), measured between the rest-frame UV and optical bands, which is sensitive to age (and dust attenuation) variations in stellar populations. The ICD shows a relation with the stellar masses and morphologies of the galaxies. Galaxies with the largest variation in their stellar populations as evidenced by high ICD have disk-dominated morphologies (with S\'{e}rsic indexes $< 2$) and stellar masses between $10 < \mathrm{Log~M/ M_\odot}< 11$. There is a marked decrease in the ICD as the stellar mass and/or the S\'ersic index increases. By studying the relations between the ICD and other galaxy properties including sizes, total colors, star-formation rate, and dust attenuation, we conclude that the largest variations in stellar populations occur in galaxies where the light from newly, high star-forming clumps contrasts older stellar disk populations. This phase reaches a peak for galaxies only with a specific stellar mass range, $10 < \mathrm{Log~M/ M_\odot} < 11$, and prior to the formation of a substantial bulge/spheroid. In contrast, galaxies at higher or lower stellar masses, and/or higher S\'{e}rsic index ($n > 2$) show reduced ICD values, implying a greater homogeneity of their stellar populations. This indicates that if a galaxy is to have both a quiescent bulge along with a star forming disk, typical of Hubble Sequence galaxies, this is most common for stellar masses $10 < \mathrm{Log~M/M_\odot} < 11$ and when the bulge component remains relatively small ($n<2$).


title: The Milky Way and Andromeda galaxies in a constrained hydrodynamical simulation: morphological evolution

abstract: We study the two main constituent galaxies of a constrained simulation of the Local Group as candidates for the Milky Way (MW) and Andromeda (M31). We focus on the formation of the stellar discs and its relation to the formation of the group as a rich system with two massive galaxies, and investigate the effects of mergers and accretion as drivers of morphological transformations. We use a state-of-the-art hydrodynamical code which includes star formation, feedback and chemical enrichment to carry out our study. We run two simulations, where we include or neglect the effects of radiation pressure from stars, to investigate the impact of this process on the morphologies and star formation rates of the simulated galaxies. We find that the simulated M31 and MW have different formation histories, even though both inhabit, at z=0, the same environment. These differences directly translate into and explain variations in their star formation rates, in-situ fractions and final morphologies. The M31 candidate has an active merger history, as a result of which its stellar disc is unable to survive unaffected until the present time. In contrast, the MW candidate has a smoother history with no major mergers at late times, and forms a disc that grows steadily; at z=0 the simulated MW has an extended, rotationally-supported disc which is dominant over the bulge. Our two feedback implementations predict similar evolution of the galaxies and their discs, although some variations are detected, the most important of which is the formation time of the discs: in the model with weaker/stronger feedback the discs form earlier/later. In summary, by comparing the formation histories of the two galaxies, we conclude that the particular merger/accretion history of a galaxy rather than its environment at the LG-scales is the main driver of the formation and subsequent growth or destruction of galaxy discs.


title: The M-dwarfs in Multiples (MinMs) survey - I. Stellar multiplicity among low-mass stars within 15 pc

abstract: We present a large-scale, volume-limited companion survey of 245 late-K to mid-M (K7-M6) dwarfs within 15 pc. Infrared adaptive optics (AO) data were analysed from the Very Large Telescope, Subaru Telescope, Canada-France-Hawaii Telescope, and MMT Observatory to detect close companions to the sample from $\sim$1 au to 100 au, while digitised wide-field archival plates were searched for wide companions from $\sim$100 au to 10,000 au. With sensitivity to the bottom of the main sequence over a separation range of 3 au to 10,000 au, multiple AO and wide-field epochs allow us to confirm candidates with common proper motions, minimize background contamination, and enable a measurement of comprehensive binary statistics. We detected 65 co-moving stellar companions and find a companion star fraction of $23.5 \pm 3.2$ per cent over the 3 au to 10,000 au separation range. The companion separation distribution is observed to rise to a higher frequency at smaller separations, peaking at closer separations than measured for more massive primaries. The mass ratio distribution across the $q = 0.2 - 1.0$ range is flat, similar to that of multiple systems with solar-type primaries. The characterisation of binary and multiple star frequency for low-mass field stars can provide crucial comparisons with star forming environments and hold implications for the frequency and evolutionary histories of their associated disks and planets.


title: Intrinsic random walks and sub-Laplacians in sub-Riemannian geometry

abstract: On a sub-Riemannian manifold we define two type of Laplacians. The \emph{macroscopic Laplacian} $\Delta_\omega$, as the divergence of the horizontal gradient, once a volume $\omega$ is fixed, and the \emph{microscopic Laplacian}, as the operator associated with a sequence of geodesic random walks. We consider a general class of random walks, where \emph{all} sub-Riemannian geodesics are taken in account. This operator depends only on the choice of a complement $\mathbf{c}$ to the sub-Riemannian distribution, and is denoted $L^c$. We address the problem of equivalence of the two operators. This problem is interesting since, on equiregular sub-Riemannian manifolds, there is always an intrinsic volume (e.g. Popp's one $P$) but not a canonical choice of complement. The result depends heavily on the type of structure under investigation. On contact structures, for every volume $\omega$, there exists a unique complement $c$ such that $\Delta_\omega=L^c$. On Carnot groups, if $H$ is the Haar volume, then there always exists a complement $c$ such that $\Delta_H=L^c$. However this complement is not unique in general. For quasi-contact structures, in general, $\Delta_P \neq L^c$ for any choice of $c$. In particular, $L^c$ is not symmetric w.r.t. Popp's measure. This is surprising especially in dimension 4 where, in a suitable sense, $\Delta_P$ is the unique intrinsic macroscopic Laplacian. A crucial notion that we introduce here is the N-intrinsic volume, i.e. a volume that depends only on the set of parameters of the nilpotent approximation. When the nilpotent approximation does not depend on the point, a N-intrinsic volume is unique up to a scaling by a constant and the corresponding N-intrinsic sub-Laplacian is unique. This is what happens for dimension smaller or equal than 4, and in particular in the 4-dimensional quasi-contact structure mentioned above.


title: Radial distributions of sub-populations in the globular cluster M15: a more centrally concentrated primordial population

abstract: We examine the radial distributions of stellar populations in the globular cluster (GC) M15, using HST/WFC3 photometry of red giants in the nitrogen-sensitive F343N-F555W color. Surprisingly, we find that giants with "primordial" composition (i.e., N abundances similar to those in field stars) are the most centrally concentrated within the WFC3 field. We then combine our WFC3 data with SDSS u, g photometry and find that the trend reverses for radii >1' (3 pc) where the ratio of primordial to N-enhanced giants increases outwards, as already found by Lardo et al. The ratio of primordial to enriched stars thus has a U-shaped dependency on radius with a minimum near the half-light radius. N-body simulations show that mass segregation might produce a trend resembling the observed one, but only if the N-enhanced giants are ~0.25 Mo less massive than the primordial giants, which requires extreme He enhancement (Y~0.40). However, such a large difference in Y is incompatible with the negligible optical color differences between primordial and enriched giants which suggest Delta Y < 0.03 and thus a difference in turn-off mass of Delta M < 0.04 Mo between the different populations. The radial trends in M15 are thus unlikely to be of dynamical origin and presumably reflect initial conditions, a result that challenges all current GC formation scenarios. We note that population gradients in the central regions of GCs remain poorly investigated and may show a more diverse behavior than hitherto thought.