title: Synthetic Gaia surveys from the FIRE cosmological simulations of Milky Way-mass galaxies

abstract: With Gaia Data Release 2, the astronomical community is entering a new era of multidimensional surveys of the Milky Way. This new phase-space view of our Galaxy demands new tools for comparing observations to simulations of Milky-Way-mass galaxies in a cosmological context, to test the physics of both dark matter and galaxy formation. We present ananke, a framework for generating synthetic phase-space surveys from high-resolution baryonic simulations, and use it to generate a suite of synthetic surveys resembling Gaia DR2 in data structure, magnitude limits, and observational errors. We use three cosmological simulations of Milky-Way-mass galaxies from the Latte suite of the Feedback In Realistic Environments (FIRE) project, which feature self-consistent clustering of star formation in dense molecular clouds and thin stellar/gaseous disks in live cosmological halos with satellite dwarf galaxies and stellar halos. We select three solar viewpoints from each simulation to generate nine synthetic Gaia-like surveys. We sample synthetic stars by assuming each star particle (of mass 7070 $M_{\odot}$) represents a single stellar population. At each viewpoint, we compute dust extinction from the simulated gas metallicity distribution and apply a simple error model to produce a synthetic Gaia-like survey that includes both observational properties and a pointer to the generating star particle. We provide the complete simulation snapshot at $z = 0$ for each simulated galaxy. We describe data access points, the data model, and plans for future upgrades. These synthetic surveys provide a tool for the scientific community to test analysis methods and interpret Gaia data.


title: One-Sided Device-Independent Certification of Unbounded Random Numbers

abstract: The intrinsic non-locality of correlations in Quantum Mechanics allow us to certify the behaviour of a quantum mechanism in a device independent way. In particular, we present a new protocol that allows an unbounded amount of randomness to be certified as being legitimately the consequence of a measurement on a quantum state. By using a sequence of non-projective measurements on single state, we show a more robust method to certify unbounded randomness than the protocol of Churchod et al., by moving to a one-sided device independent scenario. This protocol also does not assume any specific behaviour of the adversary trying to fool the participants in the protocol, which is an advantage over previous steering based protocols. We present numerical results which confirm the optimal functioning of this protocol in the ideal case. Furthermore, we also study an experimental scenario to determine the feasibility of the protocol in a realistic implementation. The effect of depolarizing noise is examined, by studying a potential state produced by a networked system of ion traps.


title: Non normal amplification of stochastic quasi-cycles

abstract: Stochastic quasi-cycles for a two species model of the excitatory-inhibitory type, arranged on a triangular loop, are studied. By increasing the strength of the inter-nodes coupling, one moves the system towards the Hopf bifurcation and the amplitude of the stochastic oscillations are consequently magnified. When the system is instead constrained to evolve on specific manifolds, selected so as to return a constant rate of deterministic damping for the perturbations, the observed amplification correlates with the degree of non normal reactivity, here quantified by the numerical abscissa. The thermodynamics of the reactive loop is also investigated and the degree of inherent reactivity shown to facilitate the out-of-equilibrium exploration of the available phase space.

categories: cond-mat.stat-mech cond-mat.dis-nn nlin.PS q-bio.PE


title: Accelerating high-throughput searches for new alloys with active learning of interatomic potentials

abstract: We propose an approach to materials prediction that uses a machine-learning interatomic potential to approximate quantum-mechanical energies and an active learning algorithm for the automatic selection of an optimal training dataset. Our approach significantly reduces the amount of DFT calculations needed, resorting to DFT only to produce the training data, while structural optimization is performed using the interatomic potentials. Our approach is not limited to one (or a small number of) lattice types (as is the case for cluster expansion, for example) and can predict structures with lattice types not present in the training dataset. We demonstrate the effectiveness of our algorithm by predicting the convex hull for the following three systems: Cu-Pd, Co-Nb-V, and Al-Ni-Ti. Our method is three to four orders of magnitude faster than conventional high-throughput DFT calculations and explores a wider range of materials space. In all three systems, we found unreported stable structures compared to the AFLOW database. Because our method is much cheaper and explores much more of materials space than high-throughput methods or cluster expansion, and because our interatomic potentials have a systematically improvable accuracy compared to empirical potentials such as EAM, etc., it will have a significant impact in the discovery of new alloy phases, particularly those with three or more components.


title: Tunable spin-orbit coupling and magnetic superstripe phase in a BEC

abstract: Superstripe phases in Bose-Einstein condensates (BECs), possessing both crystalline structure and superfluidity, opens a new avenue for exploring exotic quantum matters---supersolids. However, conclusive detection and further exploration of a superstripe is still challenging in experiments because of its short period, low visibility, fragility against magnetic field fluctuation or short lifetime. Here we propose a scheme in a spin-orbit coupled BEC which overcomes these obstacles and generates a robust magnetic superstripe phase, with only spin (no total) density modulation due to the magnetic translational symmetry, ready for direct real-space observation. In the scheme, two hyperfine spin states are individually Raman coupled with a largely-detuned third state, which induce a momentum-space separation between two lower band dispersions, yielding an effective spin-1/2 system with tunable spin-orbit coupling and Zeeman fields. Without effective Zeeman fields, spin-dependent interaction dominates, yielding a magnetic superstripe phase with a long tunable period and high visibility. Our scheme provides a platform for observing and exploring exotic properties of superstripe phases as well as novel physics with tunable spin-orbit coupling.


title: Dualities and 5-brane webs for 5d rank 2 SCFTs

abstract: We consider Type IIB 5-brane configurations for 5d rank 2 superconformal theories which are classified recently by geometry in arXiv:1801.04036. We propose all the 5-brane web diagrams for these rank 2 theories and show dualities between some of different gauge theories with explicit duality map of mass parameters and Coulomb branch moduli. In particular, we explicitly construct 5-brane configurations for $G_2$ gauge theory with six flavors and its dual $Sp(2)$ and $SU(3)$ gauge theories. We also present 5-brane webs for $SU(3)$ theories of Chern-Simons level greater than 5.


title: Modeling Majorness as a Perceptual Property in Music from Listener Ratings

abstract: For the tasks of automatic music emotion recognition, genre recognition, music recommendation it is helpful to be able to extract mode from any section of a musical piece as a perceived amount of major or minor mode (majorness) inside that section, perceived as a whole (one or several melodies and any harmony present). In this paper we take a data-driven approach (modeling directly from data without giving an explicit definition or explicitly programming an algorithm) towards modeling this property. We collect annotations from musicians and show that majorness can be understood by musicians in an intuitive way. We model this property from the data using deep learning.


title: Annihilation of Relative Equilibria in the Gravitational Field of Irregular-shaped Minor Celestial Bodies

abstract: The rotational speeds of irregular shaped minor celestial bodies can be changed by the YORP effect. This variation in speed can make the numbers, positions, stabilities, and topological cases of the minor body s relative equilibrium points vary. The numbers of relative equilibrium points can be reduced through the collision and annihilation of relative equilibrium points, or increase through the creation and separation of relative equilibrium points. Here we develop a classification system of multiple annihilation behaviors of the equilibrium points for irregular shaped minor celestial bodies. Most minor bodies have five equilibrium points; there are twice the number of annihilations per equilibrium point when the number of equilibrium points is between one and five. We present the detailed annihilation classifications for equilibria of objects which have five equilibrium points. Additionally, the annihilation classification for the seven equilibria of Kleopatra shaped objects and the nine equilibria of Bennu shaped objects are also discussed. Equilibria of different objects fall into different annihilation classifications. By letting the rotational speeds vary, we studied the annihilations and creations of relative equilibria in the gravitational field of ten minor bodies.


title: On the dynamics of pebbles in protoplanetary disks with magnetically-driven winds

abstract: We present an analytical model to investigate the production of pebbles and their radial transport through a protoplanetary disk (PPD) with magnetically driven winds. While most of the previous analytical studies in this context assume that the radial turbulent coefficient is equal to the vertical dust diffusion coefficient, in the light of the results of recent numerical simulations, we relax this assumption by adopting effective parametrisations of the turbulent coefficients involved in terms of the strength of the magnetic fields driving the wind. Theoretical studies have already pointed out that even in the absence of winds, these coefficients are not necessarily equal, though its consequences regarding pebble production have not been explored. In this paper, we investigate the evolution of the pebble production line, the radial mass flux of the pebbles and their corresponding surface density as a function of the plasma parameter at the disk midplane. Our analysis explicitly demonstrates that the presence of magnetically-driven winds in a PPD leads to considerable reduction of the rate and duration of the pebble delivery. We show that when the wind is strong, the core growth in mass due to the pebble accretion is so slow that it is unlikely that a core could reach a pebble isolation mass during a PPD lifetime. When the mass of a core reaches this critical value, pebble accretion is halted due to core-driven perturbations in the gas. With decreasing wind strength, however, pebble accretion may, in a shorter time, increase the mass of a core to the pebble isolation mass.


title: Dynamical Control of Order in a Cavity-BEC System

abstract: We demonstrate dynamical control of the superradiant transition of cavity-BEC system via periodic driving of the pump laser. We show that the dominant density wave order of the superradiant state can be suppressed, and that the subdominant competing order of Bose-Einstein condensation emerges in the steady state. Furthermore, we show that additional, non-equilibrium density wave orders, which do not exist in equilibrium, can be stabilized dynamically. Finally, for strong driving, chaotic dynamics emerges.


title: This Looks Like That: Deep Learning for Interpretable Image Recognition

abstract: When we are faced with challenging image classification tasks, we often explain our reasoning by dissecting the image, and pointing out prototypical aspects of one class or another. The mounting evidence for each of the classes helps us make our final decision. In this work, we introduce a deep network architecture -- prototypical part network (ProtoPNet), that reasons in a similar way: the network dissects the image by finding prototypical parts, and combines evidence from the prototypes to make a final classification. The model thus reasons in a way that is qualitatively similar to the way ornithologists, physicians, and others would explain to people on how to solve challenging image classification tasks. The network uses only image-level labels for training without any annotations for parts of images. We demonstrate our method on the CUB-200-2011 dataset and the Stanford Cars dataset. Our experiments show that ProtoPNet can achieve comparable accuracy with its analogous non-interpretable counterpart, and when several ProtoPNets are combined into a larger network, it can achieve an accuracy that is on par with some of the best-performing deep models. Moreover, ProtoPNet provides a level of interpretability that is absent in other interpretable deep models.


title: On the Indeterministic Nature of Star Formation on the Cloud Scale

abstract: Molecular clouds are turbulent structures whose star formation efficiency (SFE) is strongly affected by internal stellar feedback processes. In this paper we determine how sensitive the SFE of molecular clouds is to randomised inputs in the star formation feedback loop, and to what extent relationships between emergent cloud properties and the SFE can be recovered. We introduce the yule suite of 26 radiative magnetohydrodynamic (RMHD) simulations of a 10,000 solar mass cloud similar to those in the solar neighbourhood. We use the same initial global properties in every simulation but vary the initial mass function (IMF) sampling and initial cloud velocity structure. The final SFE lies between 6 and 23 percent when either of these parameters are changed. We use Bayesian mixed-effects models to uncover trends in the SFE. The number of photons emitted early in the cluster's life and the length of the cloud provide are the strongest predictors of the SFE. The HII regions evolve following an analytic model of expansion into a roughly isothermal density field. The more efficient feedback is at evaporating the cloud, the less the star cluster is dispersed. We argue that this is because if the gas is evaporated slowly, the stars are dragged outwards towards surviving gas clumps due to the gravitational attraction between the stars and gas. While star formation and feedback efficiencies are dependent on nonlinear processes, statistical models describing cloud-scale processes can be constructed.


title: Observation of the decay $\overline{B_s^0} \rightarrow \chi_{c2} K^+ K^- $ in the $\varphi$ mass region

abstract: The $\overline{B_s^0} \rightarrow \chi_{c2} K^+ K^- $ decay mode is observed and its branching fraction relative to the corresponding $\chi_{c1}$ decay mode, in a $\pm 15 \textrm{MeV}/c^2$ window around the $\phi$ mass, is found to be $\frac{\mathcal{B}(\overline{B_s^0} \rightarrow \chi_{c2} K^+ K^-) }{ \mathcal{B}(\overline{B_s^0} \rightarrow \chi_{c1} K^+ K^-)} = (17.1 \pm 3.1 \pm 0.4 \pm 0.9)\%,$ where the first uncertainty is statistical, the second systematic and the third due to the knowledge of the branching fractions of radiative $\chi_c$ decays. The decay mode $\overline{B_s^0} \rightarrow \chi_{c1} K^+ K^- $ allows the $ B_s^0$ mass to be measured as $m(B_s^0) = 5366.83 \pm 0.25 \pm 0.27 \, \textrm{MeV}/c^2,$ where the first uncertainty is statistical and the second systematic. A combination of this result with other LHCb determinations of the $B_s^0$ mass is made.


title: Adult content in Social Live Streaming Services: Characterizing deviant users and relationships

abstract: Social Live Stream Services (SLSS) exploit a new level of social interaction. One of the main challenges in these services is how to detect and prevent deviant behaviors that violate community guidelines. In this work, we focus on adult content production and consumption in two widely used SLSS, namely Live.me and Loops Live, which have millions of users producing massive amounts of video content on a daily basis. We use a pre-trained deep learning model to identify broadcasters of adult content. Our results indicate that moderation systems in place are highly ineffective in suspending the accounts of such users. We create two large datasets by crawling the social graphs of these platforms, which we analyze to identify characterizing traits of adult content producers and consumers, and discover interesting patterns of relationships among them, evident in both networks.


title: Fiber product homotopy method for multiparameter eigenvalue problems

abstract: We develop a new homotopy method for solving multiparameter eigenvalue problems (MEPs) called the fiber product homotopy method. For a $k$-parameter eigenvalue problem with matrices of sizes $n_1,\dots ,n_k = O(n)$, fiber product homotopy method requires deformation of $O(1)$ linear equations, while existing homotopy methods for MEPs require $O(n)$ nonlinear equations. We show that the fiber product homotopy method theoretically finds all eigenpairs of an MEP with probability one. It is especially well-suited for dimension-deficient singular MEPs, a weakness of all other existing methods, as the fiber product homotopy method is provably convergent with probability one for such problems as well, a fact borne out by numerical experiments. More generally, our numerical experiments indicate that the fiber product homotopy method significantly outperforms the standard Delta method in terms of accuracy, with consistent backward errors on the order of $10^{-16}$, even for dimension-deficient singular problems, and without any use of extended precision. In terms of speed, it significantly outperforms previous homotopy-based methods on all problems and outperforms the Delta method on larger problems, and is also highly parallelizable. We show that the fiber product MEP that we solve in the fiber product homotopy method, although mathematically equivalent to a standard MEP, is typically a much better conditioned problem.


title: Logarithmic Gradient Transformation and Chaos Expansion of Ito Processes

abstract: Since the seminal work of Wiener, the chaos expansion has evolved to a powerful methodology for studying a broad range of stochastic differential equations. Yet its complexity for systems subject to the white noise remains significant. The issue appears due to the fact that the random increments generated by the Brownian motion, result in a growing set of random variables with respect to which the process could be measured. In order to cope with this high dimensionality, we present a novel transformation of stochastic processes driven by the white noise. In particular, we show that under suitable assumptions, the diffusion arising from white noise can be cast into a logarithmic gradient induced by the measure of the process. Through this transformation, the resulting equation describes a stochastic process whose randomness depends only upon the initial condition. Therefore the stochasticity of the transformed system lives in the initial condition and thereby it can be treated conveniently with the chaos expansion tools.


title: Constraining the Asymptotically Safe Cosmology: cosmic acceleration without dark energy

abstract: A recently proposed Asymptotically Safe cosmology provides an elegant mechanism towards understanding the nature of dark energy and its associated cosmic coincidence problem. The underlying idea is that the accelerated expansion of the universe can occur due to infrared quantum gravity modifications at intermediate astrophysical scales (galaxies or galaxy clusters) which produce local anti-gravity sources. In this cosmological model no extra unproven energy scales or fine-tuning are used. In this study the Asymptotically Safe model is confronted with the most recent observational data from low-redshift probes, namely measurements of the Hubble parameter, standard candles (Pantheon SnIa, Quasi-stellar objects), Baryonic Acoustic Oscillations (BAOs) and high redshift probes (CMB shift parameters). Performing an overall likelihood analysis we constrain the free parameters of the model and we test its performance against the concordance model (flat $\Lambda$CDM) utilizing a large family of information criteria. We find that the Asymptotically Safe model is statistically equivalent with $\Lambda$CDM, hence it can be seen as a viable and very efficient cosmological alternative.


title: Low-dissipation Carnot-like heat engines at maximum efficient power

abstract: We study the optimal performance of Carnot-like heat engines working in low dissipation regime using the product of the efficiency and the power output, also known as the efficient power, as our objective function. Efficient power function represents the best trade-off between power and efficiency of a heat engine. We find lower and upper bounds on the efficiency in case of extreme asymmetric dissipation when the ratio of dissipation coefficients at the cold and the hot contacts approaches, respectively, zero or infinity. In addition, we obtain the form of efficiency for the case of symmetric dissipation. We also discuss the universal features of efficiency at maximum efficient power and derive the bounds on the efficiency using global linear-irreversible framework introduced recently by one of the authors.


title: A functional approach to estimation of the parameters of generalized negative binomial and gamma distributions

abstract: The generalized negative binomial distribution (GNB) is a new flexible family of discrete distributions that are mixed Poisson laws with the mixing generalized gamma (GG) distributions. This family of discrete distributions is very wide and embraces Poisson distributions, negative binomial distributions, Sichel distributions, Weibull--Poisson distributions and many other types of distributions supplying descriptive statistics with many flexible models. These distributions seem to be very promising for the statistical description of many real phenomena. GG distributions are widely applied in signal and image processing and other practical problems. The statistical estimation of the parameters of GNB and GG distributions is quite complicated. To find estimates, the methods of moments or maximum likelihood can be used as well as two-stage grid EM-algorithms. The paper presents a methodology based on the search for the best distribution using the minimization of $\ell^p$-distances and $L^p$-metrics for GNB and GG distributions, respectively. This approach, first, allows to obtain parameter estimates without using grid methods and solving systems of nonlinear equations and, second, yields not point estimates as the methods of moments or maximum likelihood do, but the estimate for the density function. In other words, within this approach the set of decisions is not a Euclidean space, but a functional space.


title: On the Error in Phase Transition Computations for Compressed Sensing

abstract: Evaluating the statistical dimension is a common tool to determine the asymptotic phase transition in compressed sensing problems with Gaussian ensemble. Unfortunately, the exact evaluation of the statistical dimension is very difficult and it has become standard to replace it with an upper-bound. To ensure that this technique is suitable, [1] has introduced an upper-bound on the gap between the statistical dimension and its approximation. In this work, we first show that the error bound in [1] in some low-dimensional models such as total variation and $\ell_1$ analysis minimization becomes poorly large. Next, we develop a new error bound which significantly improves the estimation gap compared to [1]. In particular, unlike the bound in [1] that is not applicable to settings with overcomplete dictionaries, our bound exhibits a decaying behavior in such cases.


title: Implementation of a Near-Optimal Complex Root Clustering Algorithm

abstract: We describe Ccluster, a software for computing natural $\epsilon$-clusters of complex roots in a given box of the complex plane. This algorithm from Becker et al.~(2016) is near-optimal when applied to the benchmark problem of isolating all complex roots of an integer polynomial. It is one of the first implementations of a near-optimal algorithm for complex roots. We describe some low level techniques for speeding up the algorithm. Its performance is compared with the well-known MPSolve library and Maple.


title: Micromachined force scale for optical power measurement by radiation pressure sensing

abstract: We introduce a micromachined force scale for laser power measurement by means of radiation pressure sensing. With this technique, the measured laser light is not absorbed and can be utilized while being measured. We employ silicon micromachining technology to construct a miniature force scale, opening the potential to its use for fast in-line laser process monitoring. Here we describe the mechanical sensing principle and conversion to an electrical signal. We further outline an electrostatic force substitution process for nulling of the radiation pressure force on the sensor mirror. Finally, we look at the performance of a proof-of-concept device in open-loop operation (without the nulling electrostatic force) subjected to a modulated laser at 250 W and find its response time is less than 20 ms with noise floor dominated by electronics at 2.5 W/root(Hz).


title: Approximability of Discriminators Implies Diversity in GANs

abstract: While Generative Adversarial Networks (GANs) have empirically produced impressive results on learning complex real-world distributions, recent works have shown that they suffer from lack of diversity or mode collapse. The theoretical work of Arora et al. suggests a dilemma about GANs' statistical properties: powerful discriminators cause overfitting, whereas weak discriminators cannot detect mode collapse. By contrast, we show in this paper that GANs can in principle learn distributions in Wasserstein distance (or KL-divergence in many cases) with polynomial sample complexity, if the discriminator class has strong distinguishing power against the particular generator class (instead of against all possible generators). For various generator classes such as mixture of Gaussians, exponential families, and invertible and injective neural networks generators, we design corresponding discriminators (which are often neural nets of specific architectures) such that the Integral Probability Metric (IPM) induced by the discriminators can provably approximate the Wasserstein distance and/or KL-divergence. This implies that if the training is successful, then the learned distribution is close to the true distribution in Wasserstein distance or KL divergence, and thus cannot drop modes. Our preliminary experiments show that on synthetic datasets the test IPM is well correlated with KL divergence or the Wasserstein distance, indicating that the lack of diversity in GANs may be caused by the sub-optimality in optimization instead of statistical inefficiency.


title: Light Particles with Spin in Inflation

abstract: The existence of light particles with spin during inflation is prohibited by the Higuchi bound. This conclusion can be evaded if one considers states with a sizeable coupling with the inflaton foliation, since this breaks the de Sitter isometries. The action for these states can be constructed within the Effective Field Theory of Inflation, or using a CCWZ procedure. Light particles with spin have prescribed couplings with soft inflaton perturbations, which are encoded in consistency relations. We study the phenomenology of light states with spin 2. These mix with the graviton changing the tensor power spectrum and can lead to sizeable tensor non-Gaussianities. They also give rise to a scalar bispectrum and trispectrum with a characteristic angle-dependent non-Gaussianity.


title: Supercritical causal maps : geodesics and simple random walk

abstract: We study the random planar maps obtained from supercritical Galton--Watson trees by adding the horizontal connections between successive vertices at each level. These are the hyperbolic analog of the maps studied by Curien, Hutchcroft and Nachmias in arXiv:1710.03137, and a natural model of random hyperbolic geometry. We first establish metric hyperbolicity properties of these maps: we show that they admit bi-infinite geodesics and satisfy a weak version of Gromov-hyperbolicity. We also study the simple random walk on these maps: we identify their Poisson boundary and, in the case where the underlying tree has no leaf, we prove that the random walk has positive speed. Some of the methods used here are robust, and allow us to obtain more general results about planar maps containing a supercritical Galton--Watson tree.


title: An Adaptive Volumetric Flux Boundary Condition for Lattice Boltzmann Methods

abstract: This paper presents a spatially and temporally adaptive boundary condition to specify the volumetric flux for lattice Boltzmann methods. The approach differs from standard velocity boundary conditions because it allows the velocity to vary over the boundary region provided that the total flux through the boundary satisfies a prescribed constraint, which is a typical scenario for laboratory experimental studies. This condition allows the boundary pressure to adjust dynamically to yield a specified boundary flow rate as a means to avoid unphysical mismatch between the boundary velocity and the interior flow field that can arise when a constant velocity boundary condition is applied. The method is validated for simulation of one- and two-fluid flow in complex materials, with conditions determined to match typical experiments used to study flow in porous media.


title: Strong and weak field criticality of 2D liquid gas transition

abstract: Finite size scaling (FSS) analysis of the liquid gas criticality is complicated by the absence of any broken symmetry. This, in particular, does not allow a straightforward finding of the coexistence line and the critical point. The numerical flowgram (NF) method is adapted for a controlled determination of the coexistence line and the critical point, with the critical indices $\mu,\, \nu$ measured within 1-2% of the total error. The approach based on the NF for measuring the non-analytical corrections to the diameter - the mean density of the liquid and gas phases along the coexistence line - is outlined. Our analysis is a first step toward a general evaluation of isolated critical and multicritical points.


title: On operator relations for gravitational form factors of a spin-0 hadron

abstract: The gravitational form factors for a hadron, the form factors for the hadron matrix element of the QCD energy-momentum tensor, not only describe the coupling of the hadron with a graviton, but also serve as unique quantities for describing the shape inside the hadron reflecting dynamics of quarks and gluons, such as the internal shear forces acting on the quarks/gluons and their pressure distributions. We consider the quark contribution to the gravitational form factors for a (pseudo)scalar hadron, and derive and clarify the relations satisfied by them as direct consequences of the symmetries and the equations of motion in QCD, and connections to the generalized parton distributions. Our results reveal the connections between the gravitational form factors and the higher-twist quark-gluon correlation effects inside the hadrons.


title: Electron-hole tunneling revealed by quantum oscillations in the nodal-line semimetal HfSiS

abstract: We report a study of quantum oscillations in the high-field magneto-resistance of the nodal-line semimetal HfSiS. In the presence of a magnetic field up to 31 T parallel to the c-axis, we observe quantum oscillations originating both from orbits of individual electron and hole pockets, and from magnetic breakdown between these pockets. In particular, we find an oscillation associated with a breakdown orbit enclosing one electron and one hole pocket in the form of a `figure of eight'. This observation represents an experimental confirmation of the momentum space analog of Klein tunneling. When the c-axis and the magnetic field are misaligned with respect to one another, this oscillation rapidly decreases in intensity. Finally, we extract the cyclotron masses from the temperature dependence of the oscillations, and find that the mass of the 'figure of eight' orbit corresponds to the sum of the individual pockets, consistent with theoretical predictions for Klein tunneling in topological semimetals.


title: Simple Second-Order Finite Differences for Elliptic PDEs with Discontinuous Coefficients and Interfaces

abstract: In multi-phase fluid flow, fluid-structure interaction, and other applications, partial differential equations (PDEs) often arise with discontinuous coefficients and singular sources (e.g., Dirac delta functions). These complexities arise due to changes in material properties at an immersed interface or embedded boundary, which may have an irregular shape. Consequently, the solution and its gradient can be discontinuous, and numerical methods can be difficult to design. Here a new method is presented and analyzed, using a simple formulation of one-dimensional finite differences on a Cartesian grid, allowing for a relatively easy setup for one-, two-, or three-dimensional problems. The method preserves a sharp interface with discontinuous solutions, obtained from a small number of iterations (approximately five) of solving a symmetric linear system with updates to the right- hand side. Second-order accuracy is rigorously proven in one spatial dimension and demonstrated through numerical examples in two and three spatial dimensions. The method is tested here on the variable-coefficient Poisson equation, and it could be extended for use on time-dependent problems of heat transfer, fluid dynamics, or other applications.


title: Towards Exotic Matter and Discrete Non-Abelian Symmetries in F-theory

abstract: We present a prescription in F-theory for realizing matter in "exotic" representations of product gauge groups. For 6D vacua, bifundamental hypermultiplets are engineered by starting at a singular point in moduli space which includes 6D superconformal field theories coupled to gravity. A deformation in Higgs branch moduli space takes us to a weakly coupled gauge theory description. In the corresponding elliptically fibered Calabi--Yau threefold, the minimal Weierstrass model parameters $(f,g,\Delta)$ vanish at collisions of the discriminant at least to order $(4,6,12)$, but with sufficiently high order of tangency to ensure the existence of T-brane deformations to a weakly coupled gauge theory with exotic bifundamentals. We present explicit examples including bifundamental hypermultiplets of $\mathfrak{e}_7 \times \mathfrak{su}_2$ and $\mathfrak{e}_6 \times \mathfrak{su}_3$, each of which have dual heterotic orbifold descriptions. Geometrically, these matter fields are delocalized across multiple points of an F-theory geometry. Symmetry breaking with such representations can be used to produce high dimension representations of simple gauge groups such as the four-index symmetric representation of $\mathfrak{su}_2$ and the three-index symmetric representation of $\mathfrak{su}_3$, and after further higgsing can yield discrete non-abelian symmetries.


title: Thermal Marangoni-driven Dynamics of Spinning Liquid Films

abstract: Spin coating of thin films of viscous liquids on a rotating substrate is a core technological component of semiconductor microchip fabrication. The thinning dynamics is influenced by many physical processes. Specifically temperature gradients affect thin liquid films through their influence on the local fluid surface tension. We show here experimentally and numerically that adding a static temperature gradient has a significant effect on the equilibrium film thickness and height profile reached in spin coating. Our results suggest that thermal gradients can be used to control film height profile dynamics.


title: A Hubble constant measurement from superluminal motion of the jet in GW170817

abstract: The Hubble constant ($H_0$) measures the current expansion rate of the Universe, and plays a fundamental role in cosmology. Tremendous effort has been dedicated over the past decades to measure $H_0$. Notably, Planck cosmic microwave background (CMB) and the local Cepheid-supernovae distance ladder measurements determine $H_0$ with a precision of $\sim 1\%$ and $\sim 2\%$ respectively. A $3$-$\sigma$ level of discrepancy exists between the two measurements, for reasons that have yet to be understood. Gravitational wave (GW) sources accompanied by electromagnetic (EM) counterparts offer a completely independent standard siren (the GW analogue of an astronomical standard candle) measurement of $H_0$, as demonstrated following the discovery of the neutron star merger, GW170817. This measurement does not assume a cosmological model and is independent of a cosmic distance ladder. The first joint analysis of the GW signal from GW170817 and its EM localization led to a measurement of $H_0=74^{+16}_{-8}$ km/s/Mpc (median and symmetric $68\%$ credible interval). In this analysis, the degeneracy in the GW signal between the source distance and the weakly constrained viewing angle dominated the $H_0$ measurement uncertainty. Recently, Mooley et al. (2018) obtained tight constraints on the viewing angle using high angular resolution imaging of the radio counterpart of GW170817. Here we obtain a significantly improved measurement $H_0=68.9^{+4.7}_{-4.6}$ km/s/Mpc by using these new radio observations, combined with the previous GW and EM data. We estimate that 15 more localized GW170817-like events (comparable signal-to-noise ratio, favorable orientation), having radio images and light curve data, will potentially bring resolution to the tension between the Planck and Cepheid-supernova measurements, as compared to 50-100 GW events without such data.


title: Data clustering based on Langevin annealing with a self-consistent potential

abstract: This paper introduces a novel data clustering algorithm based on Langevin dynamics, where the associated potential is constructed directly from the data. To introduce a self-consistent potential, we adopt the potential model from the established Quantum Clustering method. The first step is to use a radial basis function to construct a density distribution from the data. A potential function is then constructed such that this density distribution is the ground state solution to the time-independent Schrodinger equation. The second step is to use this potential function with the Langevin dynamics at sub-critical temperature to avoid ergodicity. The Langevin equations take a classical Gibbs distribution as the invariant measure, where the peaks of the distribution coincide with minima of the potential surface. The time dynamics of individual data points lead to different metastable states, which are interpreted as cluster centers. Clustering is therefore achieved when subsets of the data aggregate - as a result of the Langevin dynamics for a moderate period of time - in the neighborhood of a particular potential minimum. While the data points are pushed towards potential minima by the potential gradient, Brownian motion allows them to effectively tunnel through local potential barriers and escape saddle points into locations of the potential surface otherwise forbidden. The algorithm's feasibility is first established based on several illustrating examples and theoretical analyses, followed by a stricter evaluation using a standard benchmark dataset.


title: Scale Prevention By Ceramic Balls

abstract: In industrial facilities and household equipments scale formation leads to reduced efficiency and damages. Therefore various devices for anti-scale have been designed for a long time. Recently one of them, an aggregation of ceramic spheres, was proposed to prevent scale formation and its efficiency shown experimentally. The purpose of this paper is to clarify the function of this device by proposing several mathematical models and by pursuing their mathematical and numerical analysis. The first model measures the effect on nucleation of calcite of the electric potential near the surface of a ceramic sphere in natural water. The second model is based on a crystal lattice model and argues that the surface tension energy of the calcite particles is reduced by the polarization energy brought by the ceramic balls. The third model is macroscopic and numerical and studies the effect of two ceramic balls arrangements in a flow of water containing calcite particles.


title: On the accuracy of spectral solvers for micromechanics based fatigue modeling

abstract: A framework based on FFT is proposed for micromechanical fatigue modeling of polycrystals as alternative to the Finite Element method (FEM). The variational FFT approach is used with a crystal plasticity model for the cyclic behavior of the grains introduced through a FEM material subroutine, in particular an Abaqus umat. The framework also includes an alternative projection operator based on discrete differentiation to improve the microfield fidelity allowing to include second phases. The accuracy and efficiency of the FFT framework for microstructure sensitive fatigue prediction are assessed by comparing with FEM. The macroscopic cyclic response of a polycrystal obtained with both methods were indistinguishable, irrespective of the number of cycles. The microscopic fields presented small differences that decrease when using the discrete projection operator, which indeed allowed simulating accurately microstructures containing very stiff particles. Finally, the maximum differences in the fatigue life estimation from the microfields respect FEM were around 15% . In summary, this framework allows predicting fatigue life with a similar accuracy than using FEM but strongly reducing the computational cost.


title: Weight distributions of all irreducible $\mu$-constacyclic codes of length $\ell^n$

abstract: Let $\mathbb{F}_q$ be a finite field of order $q$ and integer $n\ge 1$. Let $\ell$ be a prime such that $\ell^k|(q-1)$ for some integer $k\ge 1$ and $\mu$ be an element of order $\ell^k$ in $\mathbb{F}_q$. In this paper, we determine the weight distributions of all irreducible $\mu$-constacyclic codes of length $\ell^n$ over $\mathbb{F}_q$. Explicit expressions for the generator polynomials and codewords of these codes are also obtained.


title: Cosmic ray short burst observed with the Global Muon Detector Network (GMDN) on June 22, 2015

abstract: We analyze the short cosmic ray intensity increase ("cosmic ray burst": CRB) on June 22, 2015 utilizing a global network of muon detectors and derive the global anisotropy of cosmic ray intensity and the density (i.e. the omnidirectional intensity) with 10-minute time resolution. We find that the CRB was caused by a local density maximum and an enhanced anisotropy of cosmic rays both of which appeared in association with Earth's crossing of the heliospheric current sheet (HCS). This enhanced anisotropy was normal to the HCS and consistent with a diamagnetic drift arising from the spatial gradient of cosmic ray density, which indicates that cosmic rays were drifting along the HCS from the north of Earth. We also find a significant anisotropy along the HCS, lasting a few hours after the HCS crossing, indicating that cosmic rays penetrated into the inner heliosphere along the HCS. Based on the latest geomagnetic field model, we quantitatively evaluate the reduction of the geomagnetic cut-off rigidity and the variation of the asymptotic viewing direction of cosmic rays due to a major geomagnetic storm which occurred during the CRB and conclude that the CRB is not caused by the geomagnetic storm, but by a rapid change in the cosmic ray anisotropy and density outside the magnetosphere.


title: Non-Markovian feature of the classical Hall effect

abstract: The classical Hall effect resulting from the impact of external magnetic and electric fields on the non-Markovian dynamics of charge carriers is studied. The dependence of the tangent of the Hall angle on the magnetic field is derived and compared with the experimental data for Zn. The method is proposed to determine experimentally the memory time in a system.


title: Existence and uniqueness of mild solutions for BGK models for gas mixtures of polyatomic molecules

abstract: We consider two models for a two component gas mixture with translational and internal energy degrees of freedom described by a BGK approximation assuming that the number of particles of each species remains constant. The two species are allowed to have different degrees of freedom in internal energy and are modelled by a system of kinetic BGK equations featuring two interaction terms to account for momentum and energy transfer between the species. We consider two models in the literature and prove that the models satisfy a different equation of state than the ideal gas law. Moreover, we prove that for these models there exist a unique positive solution.


title: Quark contribution to the proton spin from 2+1+1-flavor lattice QCD

abstract: We present the first chiral-continuum extrapolated up, down and strange quark spin contribution to the proton spin using lattice QCD. For the connected contributions, we use eleven ensembles of 2+1+1-flavor of Highly Improved Staggered Quarks (HISQ) generated by the MILC Collaboration. They cover four lattice spacings $a \approx \{0.15,0.12,0.09,0.06\}$ fm and three pion masses, $M_\pi \approx \{315,220,135\}$ MeV, of which two are at the physical pion mass. The disconnected strange calculations are done on seven of these ensembles, covering the four lattice spacings but only one with the physical pion mass. The disconnected light quark calculation was done on six ensembles at two values of $M_\pi \approx \{315,220\}$ MeV. High-statistics estimates on each ensemble for all three quantities allow us to quantify systematic uncertainties and perform a simultaneous chiral-continuum extrapolation in the lattice spacing and the light-quark mass. Our final results are $\Delta u \equiv \langle 1 \rangle_{\Delta u^+} = 0.777(25)(30)$, $\Delta d \equiv \langle 1 \rangle_{\Delta d^+} = -0.438(18)(30)$, and $\Delta s \equiv \langle 1 \rangle_{\Delta s^+} = -0.053(8)$, adding up to a total quark contribution to proton spin of $\sum_{q=u,d,s} (\frac{1}{2} \Delta q) = 0.143(31)(36)$. The second error is the systematic uncertainty associated with the chiral-continuum extrapolation. These results are obtained without model assumptions and are in good agreement with the recent COMPASS analysis $0.13 < \frac{1}{2} \Delta \Sigma < 0.18$, and with the $\Delta q$ obtained from various global analyses of polarized beam or target data.


title: The Impact of Binaries on the Stellar Initial Mass Function

abstract: (abridged) The stellar initial mass function (IMF) can be conveniently represented as a canonical two-part power law function and is largely invariant for star formation regions evident in the Local Group of galaxies. The IMF is a hilfskonstrukt. It is a mathematical formulation of an idealised population of stars formed together in one star formation event. The nature of the IMF (is it a probability density or an optimal sampling distribution function?) is raised. Binary stars, if unresolved, have a very significant influence at low stellar masses. Especially important is to take care of the changing binary fraction as a result of stellar-dynamical evolution of the embedded clusters which spawn the field populations of galaxies, given that the binary fraction at birth is very high and independent of primary-star mass. The high multiplicity fraction amongst massive stars leads to a substantial fraction of these being ejected out of their birth clusters and to massive stars merging. This explains the top-lightness of the IMF in star clusters in M31. In close binaries also the masses of the components can be changed due to mass transfer. A large amount of evidence points to the IMF becoming top-heavy with decreasing metallicity and above a star-formation-rate density of about 0.1 Msun/(pc^3 yr) of the cluster-forming cloud core. This is also indicated through the observed supernova rates in star-bursting galaxies. At the same time, the IMF may be bottom light at low metallicity and bottom-heavy at high metallicity, possibly accounting for the results on elliptical galaxies and ultra-faint dwarf galaxies, respectively.


title: The Dark Side of the Littlest Seesaw: freeze-in, the two right-handed neutrino portal and leptogenesis-friendly fimpzillas

abstract: We propose a minimal model to simultaneously account for a realistic neutrino spectrum through a type-I seesaw mechanism and a viable dark matter relic density. The model is an extension of the Littlest Seesaw model in which the two right-handed neutrinos of the model are coupled to a $Z_2$-odd dark sector via right-handed neutrino portal couplings. In this model, a highly constrained and direct link between dark matter and neutrino physics is achieved by considering the freeze-in production mechanism of dark matter. We show that the neutrino Yukawa couplings which describe neutrino mass and mixing may also play a dominant role in the dark matter production. We investigate the allowed regions in the parameter space of the model that provide the correct neutrino masses and mixing and simultaneously give the correct dark matter relic abundance. In certain cases the right-handed neutrino mass may be arbitrarily large, for example in the range $10^{10}-10^{11}$ GeV required for vanilla leptogenesis, with a successful relic density arising from frozen-in dark matter particles with masses around this scale, which we refer to as "fimpzillas".


title: A Universal Bound on the Strong Coupling Scale of a Gravitationally Coupled Massive Spin-2 Particle

abstract: We find a model-independent upper bound on the strong coupling scale for a massive spin-2 particle coupled to Einstein gravity. Our approach is to directly construct tree-level scattering amplitudes for these degrees of freedom and use them to find the maximum scale of perturbative unitarity violation. The highest scale is $\Lambda_3=\left(m^2M_P\right)^{1/3}$, which is saturated by ghost-free bigravity. The strong coupling scale can be further raised to $M_P$ if the kinetic term for one particle has the wrong sign, which uniquely gives the amplitudes of quadratic curvature gravity. We also discuss the generalization to massive higher-spin particles coupled to gravity.


title: Cosmological implications of ultra-light axion-like fields

abstract: Cosmological observations are used to test for imprints of an ultra-light axion-like field (ULA), with a range of potentials $V(\phi)\propto[1-\cos(\phi/f)]^n$ set by the axion-field value $\phi$ and decay constant $f$. Scalar field dynamics dictate that the field is initially frozen and then begins to oscillate around its minimum when the Hubble parameter drops below some critical value. For $n\!=\!1$, once dynamical, the axion energy density dilutes as matter; for $n\!=\!2$ it dilutes as radiation and for $n\!=\!3$ it dilutes faster than radiation. Both the homogeneous evolution of the ULA and the dynamics of its linear perturbations are included, using an effective fluid approximation generalized from the usual $n=1$ case. ULA models are parameterized by the redshift $z_c$ when the field becomes dynamical, the fractional energy density $f_{z_c} \equiv \Omega_a(z_c)/\Omega_{\rm tot}(z_c)$ in the axion field at $z_c$, and the effective sound speed $c_s^2$. Using Planck, BAO and JLA data, constraints on $f_{z_c}$ are obtained. ULAs are degenerate with dark energy for all three potentials if $1+z_c \lesssim 10$. When $3\times10^4 \gtrsim 1+z_c \gtrsim 10 $, $f_{z_c}$ is constrained to be $ \lesssim 0.004 $ for $n=1$ and $f_{z_c} \lesssim 0.02 $ for the other two potentials. The constraints then relax with increasing $z_c$. These results strongly constrain ULAs as a resolution to cosmological tensions, such as discrepant measurements of the Hubble constant, or the EDGES measurement of the global 21 cm signal.


title: A nearby superluminous supernova with a long pre-maximum 'plateau' and strong CII features

abstract: Super-luminous supernovae (SLSNe) are rare events defined as being significantly more luminous than normal terminal stellar explosions. The source of the extra powering needed to achieve such luminosities is still unclear. Discoveries in the local Universe (i.e. $z<0.1$) are scarce, but afford dense multi-wavelength observations. Additional low-redshift objects are therefore extremely valuable. We present early-time observations of the type I SLSN ASASSN-18km/SN~2018bsz. These data are used to characterise the event and compare to literature SLSNe and spectral models. Host galaxy properties are also analysed. Optical and near-IR photometry and spectroscopy were analysed. Early-time ATLAS photometry was used to constrain the rising light curve. We identified a number of spectral features in optical-wavelength spectra and tracked their time evolution. Finally, we used archival host galaxy photometry together with HII region spectra to constrain the host environment. ASASSN-18km/SN~2018bsz is found to be a type I SLSN in a galaxy at a redshift of 0.0267 (111 Mpc), making it the lowest-redshift event discovered to date. Strong CII lines are identified in the spectra. Spectral models produced by exploding a Wolf-Rayet progenitor and injecting a magnetar power source are shown to be qualitatively similar to ASASSN-18km/SN~2018bsz, contrary to most SLSNe-I that display weak/non-existent CII lines. ASASSN-18km/SN~2018bsz displays a long, slowly rising, red 'plateau' of $>$26 days, before a steeper, faster rise to maximum. The host has an absolute magnitude of --19.8 mag ($r$), a mass of M$_{*}$ = 1.5$^{+0.08}_{-0.33}$ $\times$10$^{9}$ M$_{\odot}$ , and a star formation rate of = 0.50$^{+2.22}_{-0.19}$ M$_{\odot}$ yr$^{-1}$. A nearby HII region has an oxygen abundance (O3N2) of 8.31$\pm$0.01 dex.


title: Constraints on Binary Black Hole Populations from LIGO-Virgo Detections

abstract: We reanalyse the LIGO-Virgo strain data of the 10 binary black hole mergers reported to date and compute the likelihood function in terms of chirp mass, mass ratio and effective spin. We discuss the strong degeneracy between mass ratio and spin for the three lighter events. We use this likelihood and an estimate of the horizon volume as a function of intrinsic parameters to constrain the properties of the population of merging binary black holes. The data disfavour large spins. Typical spins are constrained to $\overline a \lesssim 0.4$, even if the underlying population has randomly-oriented spins. For aligned spins the constraints are tighter, with typical spins required to be around $\overline a\sim 0.1$ and have comparable dispersion. We detect no statistically significant tendency towards a positive average spin in the direction of the orbital angular momentum. We put an upper limit on the fraction of systems where the secondary could have been tidally locked prior to the formation of the black holes (corresponding to merger times shorter than $10^8$ years) $f \lesssim 0.3$. Four events are consistent with having a maximally-spinning secondary, although one only marginally. We confirm previous findings that there is a hint of a cutoff at high mass. The data favour distributions of mass ratios with an average $\overline q \gtrsim 0.7$.


title: Galactic synchrotron distribution derived from 152 HII region absorption features in the full GLEAM survey

abstract: We derive the synchrotron distribution in the Milky Way disk from HII region absorption observations over -40{\deg} < l < 40{\deg} at six frequencies of 76.2, 83.8, 91.5, 99.2, 106.9, and 114.6 MHz with the GaLactic and Extragalactic All-sky Murchison widefield array survey (GLEAM). We develop a new method of emissivity calculation by taking advantage of the Haslam et al., (1981) map and known spectral indices, which enable us to simultaneously derive the emissivity and the optical depth of HII regions at each frequency. We show our derived synchrotron emissivities based on 152 absorption features of HII regions using both the method previously adopted in the literature and our improved method. We derive the synchrotron emissivity from HII regions to the Galactic edge along the line of sight and, for the first time, derive the emissivity from HII regions to the Sun. These results provide direct information on the distribution of the Galactic magnetic field and cosmic-ray electrons for future modelling.


title: HOD modelling of high redshift galaxies using the BLUETIDES simulation

abstract: We construct halo occupation distribution (HOD) models of high redshift ($z \gtrsim 7.5$) galaxies with $M_{*}>10^8~M_{\odot}/h$ using the BlueTides hydrodynamic simulation suite, with a particular emphasis on modelling the small scale / 1-halo clustering ($0.01\lesssim r \lesssim 1~ h^{-1}\rm{Mpc}$). Similar to low redshift studies, we find that the central and satellite mean HODs ($\left<N_{\mathrm{cen}}\right>$ and $\left<N_{\mathrm{sat}}\right>$) can be modeled by a smoothed step function and a power law respectively. The number density of satellite galaxies is however significantly suppressed compared to low redshift (satellite fractions drop from $\sim 50 \%$ at $z=0$ to $\lesssim 10 \%$ at $z=7.5$). The mean number of satellites, $\left<N_{\mathrm{sat}}\right> < 1$ for halo masses below $3 \times 10^{11} M_{\odot}/h$ (a rare halo at these redshifts). For the radial number density profiles, satellites with $10^8 \lesssim M^* \lesssim 10^{9} M_{\odot}/h$ in halos with $M_H \gtrsim 3 \times10^{11} M_{\odot}/h$ are consistent with NFW (with concentrations $c_{\mathrm{sat}} \sim 10-40$). Within halos of mass $M_H\lesssim 3 \times 10^{11} M_{\odot}/h$ satellites exhibit a power law profile with slope -3. Because these halos dominate the small scale clustering, the resulting 1-halo term is steeper than predicted using standard NFW profiles. Using this power-law profile for satellites, we can successfully reproduce the small-scale clustering exhibited by BlueTides galaxies using HOD modelling. We predict the highest probability of detecting satellites at $z>7.5$ is around centrals of $M^*\sim 3 \times 10^{10} M_{\odot}/h$ (with $M^{*}\gtrsim$ a few $10^{7} M_{\odot}/h$ ). This should be achievable with the James Webb Space Telescope (JWST).


title: Host galaxies of SNe Ic-BL with and without long gamma-ray bursts

abstract: Broad-line Ic supernovae (SNe Ic-BL) are a very rare class of core-collapse supernovae exhibiting high ejecta velocities and high kinetic energies. They are the only type of SNe that accompany long gamma-ray burst (GRB) explosions. Systematic differences found in the spectra of SNe Ic-BL with and without GRBs (GRB-SNe and SNe Ic-BL, respectively) suggest that either the progenitor or/and explosion mechanism of SNe Ic-BL with and without a GRB differ, or the difference could be only due to the viewing angle of the observer with respect to the orientation of the collimated explosion. We present the systematic comparison of the host galaxies of broad-lined SNe Ic with and without a detected GRB, the latter being detected in untargeted surveys, with the aim to find out whether there are any systematic differences between the environments in which these two classes of SNe preferentially explode. We study photometric properties of the host galaxies of a sample of 8 GRB-SNe and a sample of 28 SNe Ic-BL at z < 0.2. The two galaxy samples have indistinguishable luminosity and proper size distribution. We find indications that GRB-SNe on average occur closer to the centres of their host galaxies, i.e. the samples have a different distribution of projected offsets, normalized by the galaxy sizes. In addition we compare gas-phase metallicities of the GRB-SNe and SNe Ic-BL host samples and find that a larger fraction of super-solar metallicity hosts are found among the SNe Ic-BL without a GRB. Our results are indicative of a genuine difference between the two types of explosions and suggest that the viewing angle is not the main source of difference in the spectra of the two classes. We discuss the implications our results have on our understanding of progenitors of SNe Ic-BL with and without a GRB.


title: The Wendelstein Weak Lensing (WWL) pathfinder: Accurate weak lensing masses for Planck clusters

abstract: We present results from the Wendelstein Weak Lensing (WWL) pathfinder project, in which we have observed three intermediate redshift Planck clusters of galaxies with the new 30'$\times 30$' wide field imager at the 2m Fraunhofer Telescope at Wendelstein Observatory. We investigate the presence of biases in our shear catalogues and estimate their impact on our weak lensing mass estimates. The overall calibration uncertainty depends on the cluster redshift and is below 8.1-15 per cent for $z \approx 0.27-0.77$. It will decrease with improvements on the background sample selection and the multiplicative shear bias calibration. We present the first weak lensing mass estimates for PSZ1 G109.88+27.94 and PSZ1 G139.61+24.20, two SZ-selected cluster candidates. Based on Wendelstein colors and SDSS photometry, we find that the redshift of PSZ1 G109.88+27.94 has to be corrected to $z \approx 0.77$. We investigate the influence of line-of-sight structures on the weak lensing mass estimates and find upper limits for two groups in each of the fields of PSZ1 G109.88+27.94 and PSZ1 G186.98+38.66. We compare our results to SZ and dynamical mass estimates from the literature, and in the case of PSZ1 G186.98+38.66 to previous weak lensing mass estimates. We conclude that our pathfinder project demonstrates that weak lensing cluster masses can be accurately measured with the 2m Fraunhofer Telescope.


title: An optomechanical Bell test

abstract: Over the past few decades, experimental tests of Bell-type inequalities have been at the forefront of understanding quantum mechanics and its implications. These strong bounds on specific measurements on a physical system originate from some of the most fundamental concepts of classical physics - in particular that properties of an object are well defined independent of measurements (realism) and only affected by local interactions (locality). The violation of these bounds unambiguously shows that the measured system does not behave classically, void of any assumption on the validity of quantum theory. It has also found applications in quantum technologies for certifying the suitability of devices for generating quantum randomness, distributing secret keys and for quantum computing. Here we report on the violation of a Bell inequality involving a massive, macroscopic mechanical system. We create light-matter entanglement between the vibrational motion of two silicon optomechanical oscillators, each comprising approx. $10^{10}$ atoms, and two optical modes. This state allows us to violate a Bell inequality by more than 4 standard deviations, directly confirming the non-classical behavior of our optomechanical system under the fair sampling assumption.


title: Jet Dynamics in Compact Object Mergers: GW170817 Likely had a Successful Jet

abstract: We use relativistic hydrodynamic numerical calculations to study the interaction between a jet and a homologous outflow produced dynamically during binary neutron star mergers. We quantify how the thermal energy supplied by the jet to the ejecta and the ability of a jet to escape the homologous ejecta depend on the parameters of the jet engine and the ejecta. For collimated jets initiated at early times compared to the engine duration, we show that successful breakout of the forward cocoon shock necessitates a jet that successfully escapes the ejecta. This is because the ejecta is expanding and absorbing thermal energy, so that the forward shock from a failed jet stalls before it reaches the edge of the ejecta. This conclusion can be circumvented only for very energetic wide angle jets, with parameters that are uncomfortable given short-duration GRB observations. For successful jets, we find two regimes of jet breakout from the ejecta, early breakout on timescales shorter than the engine duration, and late breakout well after the engine shuts off. A late breakout can explain the observed delay between gravitational waves and gamma rays in GW 170817. We show that for the entire parameter space of jet parameters surveyed here (covering energies $\sim 10^{48}-10^{51}$ ergs and opening angles $\theta_j \sim 0.07-0.4$) the thermal energy deposited into the ejecta by the jet propagation is less than that produced by r-process heating on second timescales by a factor of $\gtrsim 10$. Shock heating is thus energetically subdominant in setting the luminosity of thermally powered transients coincident with neutron star mergers (kilonovae). For typical short GRB jet parameters, our conclusion is stronger: there is little thermal energy in the cocoon, much less than what is needed to explain the early blue component of the kilonova in GW 170817.


title: Pulsar observations at millimetre wavelengths

abstract: Detecting and studying pulsars above a few GHz in the radio band is challenging due to the typical faintness of pulsar radio emission, their steep spectra, and the lack of observatories with sufficient sensitivity operating at high frequency ranges. Despite the difficulty, the observations of pulsars at high radio frequencies are valuable because they can help us to understand the radio emission process, complete a census of the Galactic pulsar population, and possibly discover the elusive population in the Galactic Centre, where low-frequency observations have problems due to the strong scattering. During the decades of the 1990s and 2000s, the availability of sensitive instrumentation allowed for the detection of a small sample of pulsars above 10$\,$GHz, and for the first time in the millimetre band. Recently, new attempts between 3 and 1$\,$mm ($\approx$86$-$300$\,$GHz) have resulted in the detections of a pulsar and a magnetar up to the highest radio frequencies to date, reaching 291$\,$GHz (1.03$\,$mm). The efforts continue, and the advent of new or upgraded millimetre facilities like the IRAM 30-m, NOEMA, the LMT, and ALMA, warrants a new era of high-sensitivity millimetre pulsar astronomy in the upcoming years.


title: Bringing high-spectral resolution to VLT/SPHERE with a fibre coupling to VLT/CRIRES+

abstract: Atmospheric composition provides essential markers of the most fundamental properties of giant exoplanets, such as their formation mechanism or internal structure. New-generation exoplanet imagers, like VLT/SPHERE or Gemini/GPI, have been designed to achieve very high contrast (>15 mag) at small angular separations ($<$0.5\as) for the detection of young giant planets in the near-infrared, but they only provide very low spectral resolutions ($R<100$) for their characterization. High-dispersion spectroscopy at resolutions up to $10^5$ is one of the most promising pathways for the detailed characterization of exoplanets, but it is currently out of reach for most directly imaged exoplanets because current high-dispersion spectrographs in the near-infrared lack coronagraphs to attenuate the stellar signal and the spatial resolution necessary to resolve the planet. Project HiRISE (High-Resolution Imaging and Spectroscopy of Exoplanets) ambitions to develop a demonstrator that will combine the capabilities of two flagship instruments installed on the ESO Very Large Telescope, the high-contrast exoplanet imager SPHERE and the high-resolution spectrograph CRIRES+, with the goal of answering fundamental questions on the formation, composition and evolution of young planets. In this work, we will present the project, the first set of realistic simulations and the preliminary design of the fiber injection unit that will be implemented in SPHERE.


title: Merging Cluster Collaboration: A Panchromatic Atlas of Radio Relic Mergers

abstract: Golovich et al. 2017b presents an optical imaging and spectroscopic survey of 29 radio relic merging galaxy clusters. In this paper, we study this survey to identify substructure and quantify the dynamics of the mergers. Using a combined photometric and spectroscopic approach, we identify the minimum number of substructures in each system to describe the galaxy populations and estimate the line of sight velocity difference between likely merging subclusters. We find that the line-of-sight velocity components of the mergers are typically small compared with the maximum three dimensional relative velocity (usually $<1000$ km s$^{-1}$ and often consistent with zero). This suggests that the merger axes of these systems are generally in or near the plane of the sky matching findings in magneto-hydrodynamical simulations. In 28 of the 29 systems we identify substructures in the galaxy population aligned with the radio relic(s) and presumed associated merger induced shock. From this ensemble, we identify eight systems to include in a `gold' sample that is prime for further observation, modeling, and simulation study. Additional papers will present weak lensing mass maps and dynamical modeling for each merging system, ultimately leading to new insight into a wide range of astrophysical phenomena at some of the largest scales in the universe.


title: Steady state of interacting Floquet insulators

abstract: Floquet engineering offers tantalizing opportunities for controlling the dynamics of quantum many body systems and realizing new nonequilibrium phases of matter. However, this approach faces a major challenge: generic interacting Floquet systems absorb energy from the drive, leading to uncontrolled heating which washes away the sought after behavior. How to achieve and control a non-trivial nonequilibrium steady state is therefore of crucial importance. In this work, we study the dynamics of an interacting one-dimensional periodically-driven electronic system coupled to a phonon heat bath. Using the Floquet-Boltzmann equation (FBE) we show that the electronic populations of the Floquet eigenstates can be controlled by the dissipation. We find the regime in which the steady state features an insulatorlike filling of the Floquet bands, with a low density of additional excitations. Furthermore, we develop a simple rate equation model for the steady state excitation density that captures the behavior obtained from the numerical solution of the FBE over a wide range of parameters.

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


title: Interacting Floquet polaritons

abstract: Ordinarily, photons do not interact with one another. However, atoms can be used to mediate photonic interactions, raising the prospect of forming synthetic materials and quantum information systems from photons. One promising approach uses electromagnetically-induced transparency with highly-excited Rydberg atoms to generate strong photonic interactions. Adding an optical cavity shapes the available modes and forms strongly-interacting polaritons with enhanced light-matter coupling. However, since every atom of the same species is identical, the atomic transitions available are only those prescribed by nature. This inflexibility severely limits their utility for mediating the formation of photonic materials in cavities, as the resonator mode spectrum is typically poorly matched to the atomic spectrum. Here we use Floquet engineering to redesign the spectrum of Rubidium and make it compatible with the spectrum of a cavity, in order to explore strongly interacting polaritons in a customized space. We show that periodically modulating the energy of an atomic level redistributes its spectral weight into lifetime-limited bands separated by multiples of the modulation frequency. Simultaneously generating bands resonant with two chosen spatial modes of an optical cavity supports "Floquet polaritons" in both modes. In the presence of Rydberg dressing, we find that these polaritons interact strongly. Floquet polaritons thus provide a promising new path to quantum information technologies such as multimode photon-by-photon switching, as well as to ordered states of strongly-correlated photons, including crystals and topological fluids.


title: Joint resummation of two angularities at next-to-next-to-leading logarithmic order

abstract: Multivariate analyses are emerging as important tools to understand properties of hadronic jets which play a key role in the LHC experimental program. We take a first step towards precise and differential theory predictions by calculating the cross section for $e^+ e^- \to$ 2 jets differential in the angularities $e_\alpha$ and $e_\beta$. The logarithms of $e_\alpha$ and $e_\beta$ in the cross section are jointly resummed to next-to-next-to-leading logarithmic accuracy, using the SCET+ framework we developed, and are matched to the next-to-leading order cross section. We perform analytic one-loop calculations that serve as input for our numerical analysis, provide controlled theory uncertainties, and compare our results to Pythia. We also obtain predictions for the cross section differential in the ratio $e_\alpha/e_\beta$, which cannot be determined from a fixed-order calculation. The effect of nonperturbative corrections is also investigated. Using Event2, we validate the logarithmic structure of the single angularity cross section predicted by factorization theorems at ${\mathcal O}(\alpha_s^2)$, demonstrating that for specific angularities recoil must be taken into account when using the thrust axis, while it can be ignored if these are measured with respect to the winner-take-all axis.


title: iPTF 16hgs: A double-peaked Ca-rich gap transient in a metal poor, star forming dwarf galaxy

abstract: Calcium rich gap transients represent an intriguing new class of faint and fast evolving supernovae that exhibit strong [Ca II] emission in their nebular phase spectra. In this paper, we present the discovery and follow-up observations of iPTF 16hgs -- an intermediate luminosity and fast evolving transient that exhibited a double peaked light curve. Exhibiting a typical Type Ib spectrum in the photospheric phase and an early transition to a [Ca II] dominated nebular phase, we show that iPTF 16hgs shows properties consistent with the class of Ca-rich gap transients, with two interesting exceptions. First, while the second peak of the light curve is similar to other Ca-rich gap transients (suggesting $M_{ej}$ of 0.4 M$_\odot$ and peak luminosity of $3 \times 10^{41}$ ergs s$^{-1}$), we show that the first blue and fast declining (over $2$ days) peak is unique to this source. Second, with Integral Field Unit observations of the host galaxy, we find that iPTF 16hgs occurred in the outskirts (projected offset of $6$ kpc $ = 1.9 R_{eff}$) of a low metallicity (0.4 Z$_\odot$), star forming, dwarf spiral galaxy. Using deep late-time VLA and uGMRT observations, we place stringent limits on the local environment of the source, ruling out a large parameter space of circumstellar densities and mass loss environments of the progenitor. If iPTF 16hgs shares explosion physics with the class of Ca-rich gap transients, the presence of the first peak can be explained by enhanced mixing of 0.01 M$_\odot$ of $^{56}$Ni into the outer layers the ejecta, reminiscent of some models of He-shell detonations on WDs. On the other hand, if iPTF 16hgs is physically unrelated to the class, the first peak is consistent with shock cooling emission (of an envelope with a mass of 0.08 M$_\odot$ and radius of 13 R$_\odot$) associated with a core-collapse explosion of a highly stripped massive star in a close binary system.


title: A luminous blue kilonova and an off-axis jet from a compact binary merger at z=0.1341

abstract: The recent discovery of a faint gamma-ray burst (GRB) coincident with the gravitational wave (GW) event GW 170817 revealed the existence of a population of low-luminosity short duration gamma-ray transients produced by neutron star mergers in the nearby Universe. These events could be routinely detected by existing gamma-ray monitors, yet previous observations failed to identify them without the aid of GW triggers. Here we show that GRB150101B was an analogue of GRB170817A located at a cosmological distance. GRB 150101B was a faint short duration GRB characterized by a bright optical counterpart and a long-lived X-ray afterglow. These properties are unusual for standard short GRBs and are instead consistent with an explosion viewed off-axis: the optical light is produced by a luminous kilonova component, while the observed X-rays trace the GRB afterglow viewed at an angle of ~13 degrees. Our findings suggest that these properties could be common among future electromagnetic counterparts of GW sources.


title: Color gradients reflect an inside-out growth in early-type galaxies of the cluster MACS J1206.2-0847

abstract: We use color gradients to explore the evolution of early-type galaxies in the core of the massive galaxy cluster MACS J1206.2-0847 at z = 0.44. We used data from the CLASH and CLASH-VLT surveys to perform multiwavelength optimized model fitting using Galapagos-2 from the MegaMorph project to measure their photometric parameters. We derive color gradients for $g_{475} - I_{814}$, $r_{625} - Y_{105}$, $I_{814} - H_{160}$ , and $Y_{105} - H_{160}$ at radii ranging between 0.1 - 2 $r_e$ for 79 early-type cluster galaxies. From synthetic spectral models that use simple star formation recipes, we inferred ages and metallicities of the stellar population at different locations within each galaxy and characterized their influence on the radial color trends. We measure that galaxy sizes are $\sim$ 25% smaller in the red $H_{160}$ filter than in the blue $r_{625}$ filter but maintain a constant (within 3$\sigma$) S\'{e}rsic index $n$ with wavelength. We find negative color gradients in all colors with slopes ranging between -0.07 and -0.17 mag dex$^{-1}$ and with no obvious dependence on total magnitude, stellar mass, or location inside the cluster core. We explain the observed radial trends of color gradients as a result of the ages and metallicities of the respective stellar populations. Red galaxy cores are typically $\sim$ 3 Gyr older and more enriched in metals than the galaxy outskirts, which are of solar metallicity.


title: Sublinear-Time Quadratic Minimization via Spectral Decomposition of Matrices

abstract: We design a sublinear-time approximation algorithm for quadratic function minimization problems with a better error bound than the previous algorithm by Hayashi and Yoshida (NIPS'16). Our approximation algorithm can be modified to handle the case where the minimization is done over a sphere. The analysis of our algorithms is obtained by combining results from graph limit theory, along with a novel spectral decomposition of matrices. Specifically, we prove that a matrix $A$ can be decomposed into a structured part and a pseudorandom part, where the structured part is a block matrix with a polylogarithmic number of blocks, such that in each block all the entries are the same, and the pseudorandom part has a small spectral norm, achieving better error bound than the existing decomposition theorem of Frieze and Kannan (FOCS'96). As an additional application of the decomposition theorem, we give a sublinear-time approximation algorithm for computing the top singular values of a matrix.


title: Exact traveling wave solutions of one-dimensional models of cancer invasion

abstract: In this paper we consider continuous mathematical models of tumour growth and invasion based on the model introduced by Chaplain and Lolas \cite{Chaplain&Lolas2006}, for the case of one space dimension. The models consist of a system of three coupled nonlinear reaction-diffusion-taxis partial differential equations describing the interactions between cancer cells, the matrix degrading enzyme and the tissue. For these models under certain conditions on the model parameters we obtain exact analytical solutions in terms of traveling wave variables. These solutions are smooth positive definite functions for some of which whose profiles agree with those obtained from numerical computations \cite{Chaplain&Lolas2006} for not very large time intervals.


title: Method of model reduction and multifidelity models for solute transport in random layered porous media

abstract: This work presents a hierarchical model for solute transport in bounded layered porous media with random permeability. The model generalizes the Taylor-Aris dispersion theory to stochastic transport in random layered porous media with a known velocity covariance function. In the hierarchical model, we represent (random) concentration in terms of its cross-sectional average and a variation function. We derive a one-dimensional stochastic advection-dispersion-type equation for the average concentration and a stochastic Poisson equation for the variation function, as well as expressions for the effective velocity and dispersion coefficient. We observe that velocity fluctuations enhance dispersion in a non-monotonic fashion: the dispersion initially increases with correlation length {\lambda}, reaches a maximum, and decreases to zero at infinity. Maximum enhancement can be obtained at the correlation length about 0.25 the size of the porous media perpendicular to flow.


title: Sterile Neutrinos with Secret Interactions -- Cosmological Discord?

abstract: Several long-standing anomalies from short-baseline neutrino oscillation experiments -- most recently corroborated by new data from MiniBooNE -- have led to the hypothesis that extra, 'sterile', neutrino species might exist. Models of this type face severe cosmological constraints, and several ideas have been proposed to avoid these constraints. Among the most widely discussed ones are models with so-called 'secret interactions' in the neutrino sector. In these models, sterile neutrinos are hypothesized to couple to a new interaction, which dynamically suppresses their production in the early Universe through finite-temperature effects. Recently, it has been argued that the original calculations demonstrating the viability of this scenario need to be refined. Here, we update our earlier results from arXiv:1310.6337 [JCAP 1510 (2015) no.10, 011] accordingly. We confirm that much of the previously open parameter space for secret interactions is in fact ruled out by cosmological constraints on the sum of neutrino masses and on free-streaming of active neutrinos. We then discuss possible modifications of the vanilla scenario that would reconcile sterile neutrinos with cosmology.


title: Non-parametric decompositions of disk galaxies in S${^4}$G using DiskFit

abstract: We present photometric models of 532 disk galaxies in 3.6{\mu}m images from the Spitzer Survey of Stellar Structure in Galaxies (S$^4$G) using the non-parametric DiskFit algorithm. We first test DiskFit's performance on 400 synthetic S$^4$G-like galaxy images. DiskFit is unreliable in the bulge region, but accurately disentangles exponential disks from Ferrers bars farther out as long as their position angles differ by more than 5${^\circ}$. We then proceed to model the S$^4$G galaxies, successfully fitting 489 of them using an automated approach for initializing DiskFit, optimizing the model and deriving uncertainties using a bootstrap-resampling technique. The resulting component geometries and surface brightness profiles are compared to those derived by Salo et al. (2015) using the parametric model galfit. We find generally good agreement between the models, but discrepancies between best-fitting values for individual systems are often significant: the choice of algorithm clearly impacts the inferred disk and bar structure. In particular, we find that DiskFit typically assigns more light to the bar and less light to the disk relative to the Ferrers and exponential profiles presented by Salo et al. (2015) in the bar region. Given DiskFit's reliability at disentangling these components in our synthetic images, we conclude that the surface brightness distributions of barred S$^4$G galaxies are not well-represented by these functional forms. The results presented here underscore the importance of validating photometric decomposition algorithms before applying them to real data and the utility of DiskFit's non-parametric approach at measuring the structure of disks and bars in nearby galaxies.


title: Fractional Dp-Branes with Transverse and Tangential Dynamics in the Presence of Background Fluxes

abstract: We shall construct two boundary states which are corresponding to a dynamical fractional D$p$-brane in the presence of the fluxes of the Kalb-Ramond field and a $U(1)$ gauge potential, in the partially orbifold spacetime $\mathbb{R}^{1, 5}\times\mathbb{C}^{2}/\mathbb{Z}_{2}$. These states accurately describe the D$p$-brane in the twisted and untwisted sectors under the orbifold projection. We use them to compute the interaction of two parallel fractional D$p$-branes with the transverse velocities, tangential rotations and tangential linear motions. Various properties of the interaction, such as its long-range force, will be discussed.


title: Cable-Driven Actuation for Highly Dynamic Robotic Systems

abstract: This paper presents design and experimental evaluations of an articulated robotic limb called Capler-Leg. The key element of Capler-Leg is its single-stage cable-pulley transmission combined with a high-gap radius motor. Our cable-pulley system is designed to be as light-weight as possible and to additionally serve as the primary cooling element, thus significantly increasing the power density and efficiency of the overall system. The total weight of active elements on the leg, i.e. the stators and the rotors, contribute more than 60% of the total leg weight, which is an order of magnitude higher than most existing robots. The resulting robotic leg has low inertia, high torque transparency, low manufacturing cost, no backlash, and a low number of parts. Capler-Leg system itself, serves as an experimental setup for evaluating the proposed cable- pulley design in terms of robustness and efficiency. A continuous jump experiment shows a remarkable 96.5 % recuperation rate, measured at the battery output. This means that almost all the mechanical energy output used during push-off returned back to the battery during touch-down.


title: Low-energy antiferromagnetic spin fluctuations limit the coherent superconducting gap in cuprates

abstract: Motivated by recent attention to a potential antiferromagnetic quantum critical point at $x_c\sim 0.19$, we have used inelastic neutron scattering to investigate the low-energy spin excitations in crystals of La$_{2-x}$Sr$_x$CuO$_4$ bracketing $x_c$. We observe a peak in the normal-state spin-fluctuation weight at $\sim20$~meV for both $x=0.21$ and 0.17, inconsistent with quantum critical behavior. The presence of the peak raises the question of whether low-energy spin fluctuations limit the onset of superconducting order. Empirically evaluating the spin gap $\Delta_{\rm spin}$ in the superconducting state, we find that $\Delta_{\rm spin}$ is equal to the coherent superconducting gap $\Delta_c$ determined by electronic spectroscopies. To test whether this is a general result for other cuprate families, we have checked through the literature and find that $\Delta_c\le\Delta_{\rm spin}$ for cuprates with uniform $d$-wave superconductivity. We discuss the implications of this result.


title: The Swarm Telescope Concept

abstract: As telescope facilities become increasingly more capable they also become increasingly complex and require additional resources to operate. This is particularly true for the current and future generations of "software defined telescopes" that can support a variety of observing programs simultaneously, either through commensal observations or through support for multiple pointing centers as in the case of dipole arrays or dishes equipped with phased array feeds. At the same time, many current and future facilities are also distributed over large geographic areas, making monitoring and maintenance more difficult and costly. For these reasons we have developed a new paradigm for telescope operations called the "swarm telescope" that breaks large, single facilities into smaller groups of independent systems that can collaboratively work together to function as a single facility but with much less operational overhead. In this paper we outline the swarm telescope concept and an example of its implementation at the Long Wavelength Array. We also discuss potential advantages of using this approach for other facilities, in particular the Next Generation Very Large Array.


title: Strongly coupled electron fluids in the Poiseuille regime

abstract: In the context of describing electrons in solids as a fluid in the hydrodynamic regime, we consider a flow of electrons in a channel of finite width, i.e.~a Poiseuille flow. The electrons are accelerated by a constant electric field. We develop the appropriate relativistic hydrodynamic formalism in 2+1 dimensions and show that the fluid has a finite dc conductivity due to boundary-induced momentum relaxation, even in the absence of impurities. We use methods involving the AdS/CFT correspondence to examine the system in the strong-coupling regime. We calculate and study velocity profiles across the channel, from which we obtain the differential resistance $dV/dI$. We find that $dV/dI$ decreases with increasing current $I$ as expected for a Poiseuille flow, also at strong coupling and in the relativistic velocity regime. Moreover, we vary the coupling strength by varying $\eta/s$, the ratio of shear viscosity over entropy density. We find that $dV/dI$ decreases when the coupling is increased. We also find that strongly coupled fluids are more likely to become ultra-relativistic and turbulent. These conclusions are insensitive to the presence of impurities. In particular, we predict that in channels which are clearly in the hydrodynamic regime already at small currents, the DC channel resistance strongly depends on $\eta/s$.


title: Quasiparticle Screening near a Bosonic Superconductor-Insulator Transition Revealed by Magnetic Impurity Doping

abstract: Experiments show that the Cooper pair transport in the insulator phase that forms at thin film superconductor to insulator transitions (SIT) is simply activated. This activated behavior depends on the microscopic factors that drive the localization of the Cooper pairs. To test proposed models, we investigated how a perturbation that weakens Cooper pair binding, magnetic impurity doping, affects the characteristic activation energy, $T_0$. The data show that $T_0$ decreases monotonically with doping in films tuned farther from the SIT and increases and peaks in films that are closer to the SIT critical point. These observations provide strong evidence that the bosonic SIT in thin films is a Mott transition driven by Coulomb interactions that are screened by virtual quasi-particle excitations. This dependence on underlying fermionic degrees of freedom distinguishes these SITs from those in micro-fabricated Josephson Junction Arrays, cold atom systems, and likely in high temperature superconductors with nodes in their quasiparticle density of states.


title: Particle Camera MX-10 in Physics Education

abstract: We present several applications of the particle camera MX-10 in different radiation environments, leading to strikingly different observed patterns and consequently particle composition as well as recorded dose. We describe the measurements of the radiation background, cosmic muons direction analysis and the detection of alpha particles in a natural cave. We focus on the analysis of tracks of highly ionizing particles aboard an airplane and compute the dose rate with altitude. We also comment briefly on a test in beams of particles at an accelerator. We discuss the problem of presenting sensitive results on radiation background levels in various environments to general public. Most of the experiments presented here were performed and analyzed by high-school students, confirming the application of the device not only for class-room demonstrations but also for science projects in physics and IT.


title: Sustainable blockchain-enabled services: Smart contracts

abstract: This chapter contributes to evolving the versatility and complexity of blockchain-enabled services through extending the functionality of blockchain-enforced smart contracts. The contributions include: (i) a method for automated management of contracts with hierarchical conditionality structures through an hierarchy of intelligent agents and the use of hierarchical cryptographic key-pairs; (ii) a method for efficient and secure matching and transfer of smart- contract underlyings (entities) among disparate smart contracts/subcontracts; (iii) a method for producing an hierarchy of common secrets to facilitate hierarchical communication channels of increased security in the context of smart contracts/subcontracts/underlyings; and (iv) a method for building secure and optimized repositories through distributed hash tables in the context of contracts/ subcontracts/underlyings. These methods help providing services that allow both narrower and worldwide reach and distribution of resources. The longevity of the blockchain technology is achieved through continuous innovation. Blockchain-enabled services are potentially an efficient, secure, automated, and cost-effective alternative or complement to current service infrastructures in a range of domains (legal, medical, financial, government, IoT).


title: An Introduction to Animal Movement Modeling with Hidden Markov Models using Stan for Bayesian Inference

abstract: Hidden Markov models (HMMs) are popular time series model in many fields including ecology, economics and genetics. HMMs can be defined over discrete or continuous time, though here we only cover the former. In the field of movement ecology in particular, HMMs have become a popular tool for the analysis of movement data because of their ability to connect observed movement data to an underlying latent process, generally interpreted as the animal's unobserved behavior. Further, we model the tendency to persist in a given behavior over time. Notation presented here will generally follow the format of Zucchini et al. (2016) and cover HMMs applied in an unsupervised case to animal movement data, specifically positional data. We provide Stan code to analyze movement data of the wild haggis as presented first in Michelot et al. (2016).


title: A Diagrammatic Representation of Phase Portraits and Bifurcation Diagrams of Two-Dimensional Dynamical Systems

abstract: We treat the problem of characterizing in a systematic way the qualitative features of two-dimensional dynamical systems. To that end, we construct a representation of the topological features of phase portraits by means of diagrams that discard their quantitative information. All codimension 1 bifurcations are naturally embodied in the possible ways of transitioning smoothly between diagrams. We introduce a representation of bifurcation curves in parameter space that guides the proposition of bifurcation diagrams compatible with partial information about the system.


title: Hyperuniformity and anti-hyperuniformity in one-dimensional substitution tilings

abstract: We consider the scaling properties characterizing the hyperuniformity (or anti-hyperuniformity) of long wavelength fluctuations in a broad class of one-dimensional substitution tilings. We present a simple argument that predicts the exponent $\alpha$ governing the scaling of Fourier intensities at small wavenumbers, tilings with $\alpha>0$ being hyperuniform, and confirm with numerical computations that the predictions are accurate for quasiperiodic tilings, tilings with singular continuous spectra, and limit-periodic tilings. Tilings with quasiperiodic or singular continuous spectra can be constructed with $\alpha$ arbitrarily close to any given value between $-1$ and $3$. Limit-periodic tilings can be constructed with $\alpha$ between $-1$ and $1$ or with Fourier intensities that approach zero faster than any power law.

categories: cond-mat.stat-mech cond-mat.mtrl-sci math-ph math.DS math.MP


title: PIDS - A Behavioral Framework for Analysis and Detection of Network Printer Attacks

abstract: Nowadays, every organization might be attacked through its network printers. The malicious exploitation of printing protocols is a dangerous and underestimated threat against every printer today, as highlighted by recent published researches. This article presents PIDS (Printers' IDS), an intrusion detection system for detecting attacks on printing protocols. PIDS continuously captures various features and events obtained from traffic produced by printing protocols in order to detect attacks. As part of this research we conducted thousands of automatic and manual printing protocol attacks on various printers and recorded thousands of the printers' benign network sessions. Then we applied various supervised machine learning (ML) algorithms to classify the collected data as normal (benign) or abnormal (malicious). We evaluated several detection algorithms, feature selection methods, and the features needed in order to obtain the best detection results for protocol traffic of printers. Our empirical results suggest that the proposed framework is effective in detecting printing protocol attacks, providing an accuracy of 99.9 with negligible fall-positive rate.


title: Compact interaction potential for van der Waals nanorods

abstract: We studied the van der Waals interactions of two finite, solid, cylindrical rods at arbitrary angle and position with respect to each other. An analytic interpolative formula for the interaction potential energy is constructed, based on various asymptotic cases. The potential can be readily used for numerical and analytic description of multi-wall carbon nanotubes, metallic nanorods, rod-shaped colloids, or any other similar objects with significant van der Waals interactions.