title: An ultrahot gas-giant exoplanet with a stratosphere

abstract: Infrared radiation emitted from a planet contains information about the chemical composition and vertical temperature profile of its atmosphere. If upper layers are cooler than lower layers, molecular gases will produce absorption features in the planetary thermal spectrum. Conversely, if there is a stratosphere - where temperature increases with altitude - these molecular features will be observed in emission. It has been suggested that stratospheres could form in highly irradiated exoplanets, but the extent to which this occurs is unresolved both theoretically and observationally. A previous claim for the presence of a stratosphere remains open to question, owing to the challenges posed by the highly variable host star and the low spectral resolution of the measurements. Here we report a near-infrared thermal spectrum for the ultrahot gas giant WASP-121b, which has an equilibrium temperature of approximately 2,500 kelvin. Water is resolved in emission, providing a detection of an exoplanet stratosphere at 5-sigma confidence. These observations imply that a substantial fraction of incident stellar radiation is retained at high altitudes in the atmosphere, possibly by absorbing chemical species such as gaseous vanadium oxide and titanium oxide.


title: Complex rotation numbers: bubbles and their intersections

abstract: The construction of complex rotation numbers, due to V.Arnold, gives rise to a fractal-like set "bubbles" related to a circle diffeomorphism. "Bubbles" is a complex analogue to Arnold tongues. This article contains a survey of the known properties of bubbles, as well as a variety of open questions. In particular, we show that bubbles can intersect and self-intersect, and provide approximate pictures of bubbles for perturbations of M\"obius circle diffeomorphisms.


title: The Spatial Shape of Avalanches

abstract: In disordered elastic systems, driven by displacing a parabolic confining potential adiabatically slowly, all advance of the system is in bursts, termed avalanches. Avalanches have a finite extension in time, which is much smaller than the waiting-time between them. Avalanches also have a finite extension $\ell$ in space, i.e. only a part of the interface of size $\ell$ moves during an avalanche. Here we study their spatial shape $\left< S(x)\right>_{\ell}$ given $\ell$, as well as its fluctuations encoded in the second cumulant $\left< S^{2}(x)\right>_{\ell}^{{\rm c}}$. We establish scaling relations governing the behavior close to the boundary. We then give analytic results for the Brownian force model, in which the microscopic disorder for each degree of freedom is a random walk. Finally, we confirm these results with numerical simulations. To do this properly we elucidate the influence of discretization effects, which also confirms the assumptions entering into the scaling ansatz. This allows us to reach the scaling limit already for avalanches of moderate size. We find excellent agreement for the universal shape, its fluctuations, including all amplitudes.


title: The Gram-Schmidt Walk: A Cure for the Banaszczyk Blues

abstract: An important result in discrepancy due to Banaszczyk states that for any set of $n$ vectors in $\mathbb{R}^m$ of $\ell_2$ norm at most $1$ and any convex body $K$ in $\mathbb{R}^m$ of Gaussian measure at least half, there exists a $\pm 1$ combination of these vectors which lies in $5K$. This result implies the best known bounds for several problems in discrepancy. Banaszczyk's proof of this result is non-constructive and a major open problem has been to give an efficient algorithm to find such a $\pm 1$ combination of the vectors. In this paper, we resolve this question and give an efficient randomized algorithm to find a $\pm 1$ combination of the vectors which lies in $cK$ for $c>0$ an absolute constant. This leads to new efficient algorithms for several problems in discrepancy theory.


title: Satellitesimal Formation via Collisional Dust Growth in Steady Circumplanetary Disks

abstract: The icy satellites around Jupiter are considered to have formed in a circumplanetary disk. While previous models focused on the formation of satellites starting from satellitesimals, the question of how satellitesimals form from smaller dust particles has not been addressed so far. In this work, we study the possibility that satellitesimals form in situ in a circumplanetary disk. We calculate the radial distribution of the surface density and representative size of icy dust particles that grow by colliding with each other and drift toward the central planet in a steady circumplanetary disk with a continuous supply of gas and dust from the parent protoplanetary disk. The radial drift barrier is overcome if the ratio of the dust to gas accretion rates onto the circumplanetary disk, $\dot{M}_{\mathrm{d}}/\dot{M}_{\mathrm{g}}$, is high and the strength of turbulence, $\alpha$, is not too low. The collision velocity is lower than the critical velocity of fragmentation when $\alpha$ is low. Taken together, we find that the conditions for satellitesimal formation via dust coagulation are given by $\dot{M}_{\mathrm{d}}/\dot{M}_{\mathrm{g}}\ge1$ and $10^{-4}\le\alpha<10^{-2}$. The former condition is generally difficult to achieve, suggesting that the in-situ satellitesimal formation via particle sticking is viable only under an extreme condition. We also show that neither satellitesimal formation via the collisional growth of porous aggregates nor via streaming instability is viable as long as $\dot{M}_{\mathrm{d}}/\dot{M}_{\mathrm{g}}$ is low.


title: Lower bounds for the measurable chromatic number of the hyperbolic plane

abstract: Consider the graph $\mathbb{H}(d)$ whose vertex set is the hyperbolic plane, where two points are connected with an edge when their distance is equal to some $d>0$. Asking for the chromatic number of this graph is the hyperbolic analogue to the famous Hadwiger-Nelson problem about colouring the points of the Euclidean plane so that points at distance $1$ receive different colours. As in the Euclidean case, one can lower bound the chromatic number of $\mathbb{H}(d)$ by $4$ for all $d$. Using spectral methods, we prove that if the colour classes are measurable, then at least $6$ colours are are needed to properly colour $\mathbb{H}(d)$ when $d$ is sufficiently large.


title: Complete analysis of ensemble inequivalence in the Blume-Emery-Griffiths model

abstract: We study inequivalence of canonical and microcanonical ensembles in the mean-field Blume-Emery-Griffiths model. This generalizes previous results obtained for the Blume-Capel model. The phase diagram strongly depends on the value of the biquadratic exchange interaction K, the additional feature present in the Blume-Emery-Griffiths model. At small values of K, as for the Blume-Capel model, lines of first and second order phase transitions between a ferromagnetic and a paramagnetic phase are present, separated by a tricritical point whose location is different in the two ensembles. At higher values of K the phase diagram changes substantially, with the appearance of a triple point in the canonical ensemble which does not find any correspondence in the microcanonical ensemble. Moreover, one of the first order lines that starts from the triple point ends in a critical point, whose position in the phase diagram is different in the two ensembles. This line separates two paramagnetic phases characterized by a different value of the quadrupole moment. These features were not previously studied for other models and substantially enrich the landscape of ensemble inequivalence, identifying new aspects that had been discussed in a classification of phase transitions based on singularity theory. Finally, we discuss ergodicity breaking, which is highlighted by the presence of gaps in the accessible values of magnetization at low energies: it also displays new interesting patterns that are not present in the Blume-Capel model.


title: Target of Opportunity Observations of Blazars with H.E.S.S

abstract: The very-high-energy (VHE, E>100 GeV) extragalactic sky is dominated by blazars, a class of active galactic nuclei which show rapid variability at all wavelengths. Target of Opportunity (ToO) observations triggered by flaring activity detected at longer wavelengths are, thus, an important part of the blazar observing strategy of H.E.S.S., an array of five imaging atmospheric Cherenkov telescopes sensitive to VHE photons. In this contribution we detail the H.E.S.S. extragalactic ToO program, describing the specific procedures currently in place to follow up on multi-wavelength alerts. The program is illustrated by discussing a few recent noteworthy targets observed with the H.E.S.S. phase II array over the last two years of blazar ToO observations.


title: Square function estimates for the Bochner-Riesz means

abstract: We consider the square function (known as Stein's square function) estimate associated with the Bochner-Riesz means. The previously known range of sharp estimate is improved. Our results are based on vector valued extensions of Bennett-Carbery-Tao's multilinear (adjoint) restriction estimate and adaption of induction argument due to Bourgain-Guth. Unlike the previous work by Bourgain-Guth on $L^p$ boundedness of Bochner-Riesz means in which oscillatory operators associated to the kernel had been studied, we take more direct approach by working on the Fourier transform side. This enables us to obtain the correct order of smoothing which is essential for obtaining the sharp estimate for the square function.


title: Lorenz curves interpretations of the Bruss-Duerinckx theorem for resource dependent branching processes

abstract: The Bruss and Duerinckx theorem for resource dependent branching processes states that the survival of any society form is nested in an envelope formed by two extreme policies. The objective of this paper is to give a novel interpretation of this theorem through the use of Lorenz curves. This representation helps us visualize how the parameters interplay. Besides, as we will show, it clarifies the impact of inequality in consumption.


title: Nodal-line semimetals from Weyl superlattices

abstract: The existence and topological classification of lower-dimensional Fermi surfaces is often tied to the crystal symmetries of the underlying lattice systems. Artificially engineered lattices, such as heterostructures and other superlattices, provide promising avenues to realize desired crystal symmetries that protect lower-dimensional Fermi surface, such as nodal lines. In this work, we investigate a Weyl semimetal subjected to spatially periodic onsite potential, giving rise to several phases, including a nodal-line semimetal phase. In contrast to proposals that purely focus on lattice symmetries, the emergence of the nodal line in this setup does not require small spin-orbit coupling, but rather relies on its presence. We show that the stability of the nodal line is understood from reflection symmetry and a combination of a fractional lattice translation and charge-conjugation symmetry. Depending on the choice of parameters, this model exhibits drumhead surface states that are exponentially localized at the surface, or weakly localized surface states that decay into the bulk at all energies.


title: AWAKE readiness for the study of the seeded self-modulation of a 400\,GeV proton bunch

abstract: AWAKE is a proton-driven plasma wakefield acceleration experiment. % We show that the experimental setup briefly described here is ready for systematic study of the seeded self-modulation of the 400\,GeV proton bunch in the 10\,m-long rubidium plasma with density adjustable from 1 to 10$\times10^{14}$\,cm$^{-3}$. % We show that the short laser pulse used for ionization of the rubidium vapor propagates all the way along the column, suggesting full ionization of the vapor. % We show that ionization occurs along the proton bunch, at the laser time and that the plasma that follows affects the proton bunch. %


title: GRB Observations with H.E.S.S. II

abstract: The High Energy Stereoscopic System (H.E.S.S.) has been searching for counterparts of Gamma Ray Bursts (GRBs) for many years. In 2012 the system was upgraded with a fifth $28$ m diameter telescope (CT5) which is equipped with faster motors for rapid repointing, marking the start of the second phase of H.E.S.S. operation (H.E.S.S. II). CT5s large light collection area of $600\,{\rm m}^{2}$ improves the sensitivity to low-energy gamma-rays and even extends the energy range below $100$ GeV. The search for counterparts continues now in the energy range of tens of GeV to tens of TeV. A detection in this energy range would open a new window to the part of the spectrum of these highly energetic explosions which Fermi-LAT has only successfully detected in a reduced subset of events, with rather limited statistics. In the past years, H.E.S.S. has performed followup observations based on GRB detections by Swift-BAT and Fermi-GBM/-LAT. This Target of Opportunity observation program was carried out with a generalised Target of Opportunity Alert system. This contribution will highlight key features of the Target of Opportunity Alert system, present follow-up statistics of GRBs as well as detailed results of promising follow-up observations.


title: Using the SLEUTH urban growth model to simulate the impacts of future policy scenarios on urban land use in the Tehran metropolitan area in Iran

abstract: The SLEUTH model, based on the Cellular Automata (CA), can be applied to city development simulation in metropolitan areas. In this study the SLEUTH model was used to model the urban expansion and predict the future possible behavior of the urban growth in Tehran. The fundamental data were five Landsat TM and ETM images of 1988, 1992, 1998, 2001 and 2010. Three scenarios were designed to simulate the spatial pattern. The first scenario assumed historical urbanization mode would persist and the only limitations for development were height and slope. The second one was a compact scenario which makes the growth mostly internal and limited the expansion of suburban areas. The last scenario proposed a polycentric urban structure which let the little patches grow without any limitation and would not consider the areas beyond the specific buffer zone from the larger patches for development. Results showed that the urban growth rate was greater in the first scenario in comparison with the other two scenarios. Also it was shown that the third scenario was more suitable for Tehran since it could avoid undesirable effects such as congestion and pollution and was more in accordance with the conditions of Tehran city.


title: Topological entropy of closed sets in $[0,1]^2$

abstract: We generalize the definition of topological entropy due to Adler, Konheim, and McAndrew \cite{AKM} to set-valued functions from a closed subset $A$ of the interval to closed subsets of the interval. We view these set-valued functions, via their graphs, as closed subsets of $[0,1]^{2}$. We show that many of the topological entropy properties of continuous functions of a compact topological space to itself hold in our new setting, but not all. We also compute the topological entropy of some examples, relate the entropy to other dynamical and topological properties of the examples, and we give an example of a closed subset $G$ of $[0,1]^2$ that has $0$ entropy but $G\cup \{ (p,q) \},$ where $(p,q) \in [0,1]^2\setminus G,$ has infinite entropy.


title: Multi-conjugated adaptive optics imaging of distant galaxies -- A comparison of Gemini/GSAOI and VLT/HAWK-I data

abstract: Multi-conjugated adaptive optics (MCAO) yield nearly diffraction-limited images at 2$\mu$m wavelengths. Currently, GeMS/GSAOI at Gemini South is the only MCAO facility instrument at an 8m telescope. Using real data and for the first time, we investigate the gain in depth and S/N when MCAO is employed for $K_{\rm s}$-band observations of distant galaxies. Our analysis is based on the Frontier Fields cluster MACS J0416.1-2403, observed with GeMS/GSAOI (near diffraction-limited) and compared against VLT/HAWK-I (natural seeing) data. Using galaxy number counts, we show that the substantially increased thermal background and lower optical throughput of the MCAO unit are fully compensated for by the wavefront correction, because the galaxy images can be measured in smaller apertures with less sky noise. We also performed a direct comparison of the signal-to-noise ratios (S/N) of sources detected in both data sets. For objects with intrinsic angular sizes corresponding to half the HAWK-I image seeing, the gain in S/N is 40 per cent. Even smaller objects experience a boost in S/N by a up to a factor of 2.5 despite our suboptimal natural guide star configuration. The depth of the near diffraction limited images is more difficult to quantify than that of seeing limited images, due to a strong dependence on the intrinsic source profiles. Our results emphasize the importance of cooled MCAO systems for $K_{\rm s}$-band observations with future, extremely large telescopes.


title: Chaotic Properties of Single Element Nonlinear Chimney Model: Effect of Directionality

abstract: We generalize the chimney model by introducing nonlinear restoring and gravitational forces for the purpose of modeling swaying of trees at high wind speeds. We have derived general equations governing the system using Lagrangian formulation. We have studied the simplest case of a single element in more detail. The governing equation we arrive at for this case has not been studied so far. We study the chaotic properties of this simple building block and also the effect of directionality in the wind on the chaotic properties. We also consider the special case of two elements.


title: N\'emethi's division algorithm for zeta-functions of plumbed 3-manifolds

abstract: A polynomial counterpart of the Seiberg-Witten invariant associated with a negative definite plumbed 3-manifold has been proposed by earlier work of the authors. It is provided by a special decomposition of the zeta-function defined by the combinatorics of the manifold. In this article we give an algorithm, based on multivariable Euclidean division of the zeta-function, for the explicit calculation of the polynomial, in particular for the Seiberg--Witten invariant.


title: Equidistribution theorems on strongly pseudoconvex domains

abstract: This work consists of two parts. In the first part, we consider a compact connected strongly pseudoconvex CR manifold $X$ with a transversal CR $S^{1}$ action. We establish an equidistribution theorem on zeros of CR functions. The main techniques involve a uniform estimate of Szeg\H{o} kernel on $X$. In the second part, we consider a general complex manifold $M$ with a strongly pseudoconvex boundary $X$. By using classical result of Boutet de Monvel-Sj\"ostrand about Bergman kernel asymptotics, we establish an equidistribution theorem on zeros of holomorphic functions on $\overline M$.


title: On regular induced subgraphs of generalized polygons

abstract: The cage problem asks for the smallest number $c(k,g)$ of vertices in a $k$-regular graph of girth $g$ and graphs meeting this bound are known as cages. While cages are known to exist for all integers $k \ge 2$ and $g \ge 3$, the exact value of $c(k, g)$ is known only for some small values of $k, g$ and three infinite families where $g \in \{6, 8, 12\}$ and $k - 1$ is a prime power. These infinite families come from the incidence graphs of generalized polygons. Some of the best known upper bounds on $c(k,g)$ for $g \in \{6, 8, 12\}$ have been obtained by constructing small regular induced subgraphs of these cages. In this paper, we first use the Expander Mixing Lemma to give a general lower bound on the size of an induced $k$-regular subgraph of a regular bipartite graph in terms of the second largest eigenvalue of the host graph. We use this bound to show that the known construction of $(k,6)$-graphs using Baer subplanes of the Desarguesian projective plane is the best possible. For generalized quadrangles and hexagons, our bounds are new. In particular, we improve the known lower bound on the size of a $q$-regular induced subgraphs of the classical generalized quadrangle $\mathsf{Q}(4,q)$ and show that the known constructions are asymptotically sharp. For prime powers $q$, we also improve the known upper bounds on $c(q,8)$ and $c(q,12)$ by giving new geometric constructions of $q$-regular induced subgraphs in the symplectic generalized quadrangle $\mathsf{W}(3,q)$ and the split Cayley hexagon $\mathsf{H}(q)$, respectively. Our constructions show that \[c(q,8) \le 2(q^3 - q\sqrt{q} - q)\] for $q$ an even power of a prime, and \[c(q, 12) \le 2(q^5 - 3q^3)\] for all prime powers $q$. For $q \in \{3,4,5\}$ we also give a computer classification of all $q$-regular induced subgraphs of the classical generalized quadrangles of order $q$.


title: From sudden quench to adiabatic dynamics in the attractive Hubbard model

abstract: We study the crossover between the sudden quench limit and the adiabatic dynamics of superconducting states in the attractive Hubbard model. We focus on the dynamics induced by the change of the attractive interaction during a finite ramp time which is varied in order to track the evolution of the dynamical phase diagram from the sudden quench to the equilibrium limit. Two different dynamical regimes are realized for quenches towards weak and strong coupling interactions. At weak coupling the dynamics depends only on the energy injected into the system, whereas a dynamics retaining memory of the initial state takes place at strong coupling. We show that this is related to a sharp transition between a weak and a strong coupling quench dynamical regime, which defines the boundaries beyond which a dynamics independent from the initial state is recovered. Comparing the dynamics in the superconducting and non-superconducting phases we argue that this is due to the lack of an adiabatic connection to the equilibrium ground state for non-equilibrium superconducting states in the strong coupling quench regime.


title: Absolute Magnitude Calibration for W UMa-Type Systems Based on Gaia Data

abstract: Tycho-Gaia Astrometric Solution (TGAS) parallax data are used to determine absolute magnitudes $M_V$ for 318 W~UMa-type (EW) contact binary stars. A very steep (slope $\simeq -9$), single-parameter ($\log{P}$), linear calibration can be used to predict $M_V$ to about 0.1 -- 0.3 mag over the whole range of accessible orbital period, $0.22\!<\!P\!<\!0.88$ days. A similar calibration for the most common systems with $0.275\!<\!P\!<\!0.575$ days predicts $M_V$ values to about 0.06 -- 0.16 mag. For orbital period values both shorter and longer than the central range, the period dependence is respectively steeper and shallower, i.e.\ the binaries are fainter in $M_V$ than predicted by the whole-range, linear law. The steepness of the relation for short-period systems implies important consequences for the detectability of the faintest binaries defining the short-period cut-off of the period distribution. Although the scatter around the linear $\log{P}$-fit is fairly large (0.2 -- 0.4 mag), the current data do not support the inclusion of a $B\!-\!V$ color term in the calibration.


title: Parameter optimization for surface flux transport models

abstract: Accurate prediction of solar activity calls for precise calibration of solar cycle models. Consequently we aim to find optimal parameters for models which describe the physical processes on the solar surface, which in turn act as proxies for what occurs in the interior and provide source terms for coronal models. We use a genetic algorithm to optimize surface flux transport models using National Solar Observatory (NSO) magnetogram data for Solar Cycle 23. This is applied to both a 1D model that inserts new magnetic flux in the form of idealized bipolar magnetic regions, and also to a 2D model that assimilates specific shapes of real active regions. The genetic algorithm searches for parameter sets (meridional flow speed and profile, supergranular diffusivity, initial magnetic field, and radial decay time) that produce the best fit between observed and simulated butterfly diagrams, weighted by a latitude-dependent error structure which reflects uncertainty in observations. Due to the easily adaptable nature of the 2D model, the optimization process is repeated for Cycles 21, 22, and 24 in order to analyse cycle-to-cycle variation of the optimal solution. We find that the ranges and optimal solutions for the various regimes are in reasonable agreement with results from the literature, both theoretical and observational. The optimal meridional flow profiles for each regime are almost entirely within observational bounds determined by magnetic feature tracking, with the 2D model being able to accommodate the mean observed profile more successfully. Differences between models appear to be important in deciding values for the diffusive and decay terms. In like fashion, differences in the behaviours of different solar cycles lead to contrasts in parameters defining the meridional flow and initial field strength.


title: Generation of Low Frequency Plasma Waves After Wave-Breaking

abstract: Spatio-temporal evolution of a non-relativistic electrostatic waves in a cold plasma has been studied in the wave-breaking regime using a 1D particle-in-cell simulation. It is found that plasma gets heated after the wave-breaking but a fraction of initial energy always remains with the remnant wave in the form of BGK mode in warm plasma. An interesting finding of this work is that the frequency of the resultant BGK wave is found be below electron plasma frequency which decreases with increasing initial amplitude. Moreover, the acceleration mechanism after the wave-breaking is also found to be different from the previous work. In order to explain the results observed in the numerical experiments, a simplified theoretical model is constructed which exhibits a good agreement with the simulation. These investigations have direct relevance in wake-field acceleration experiments.


title: Spatial Adiabatic Passage of Massive Quantum Particles

abstract: By adiabatically manipulating tunneling amplitudes of cold atoms in a periodic potential with a multiple sublattice structure, we are able to coherently transfer atoms from a sublattice to another without populating the intermediate sublattice, which can be regarded as a spatial analogue of stimulated Raman adiabatic passage. A key is the existence of dark eigenstates forming a flat band in a Lieb-type optical lattice. We also successfully observe a matter-wave analogue of Autler-Townes doublet using the same setup. This work shed light on a novel kind of coherent control of cold atoms in optical potentials.


title: Learning Feature Pyramids for Human Pose Estimation

abstract: Articulated human pose estimation is a fundamental yet challenging task in computer vision. The difficulty is particularly pronounced in scale variations of human body parts when camera view changes or severe foreshortening happens. Although pyramid methods are widely used to handle scale changes at inference time, learning feature pyramids in deep convolutional neural networks (DCNNs) is still not well explored. In this work, we design a Pyramid Residual Module (PRMs) to enhance the invariance in scales of DCNNs. Given input features, the PRMs learn convolutional filters on various scales of input features, which are obtained with different subsampling ratios in a multi-branch network. Moreover, we observe that it is inappropriate to adopt existing methods to initialize the weights of multi-branch networks, which achieve superior performance than plain networks in many tasks recently. Therefore, we provide theoretic derivation to extend the current weight initialization scheme to multi-branch network structures. We investigate our method on two standard benchmarks for human pose estimation. Our approach obtains state-of-the-art results on both benchmarks. Code is available at https://github.com/bearpaw/PyraNet.


title: Launching and Controlling Gaussian Beams from Point Sources via Planar Transformation Media

abstract: Based on operations prescribed under the paradigm of Complex Transformation Optics (CTO) [1-5], it was recently shown in [5] that a complex source point (CSP) can be mimicked by a parity-time ($\mathcal{PT}$) transformation media. Such coordinate transformation has a mirror symmetry for the imaginary part, and results in a balanced loss/gain metamaterial slab. A CSP produces a Gaussian beam and, consequently, a point source placed at the center of such metamaterial slab produces a Gaussian beam propagating away from the slab. Here, we extend the CTO analysis to non-symmetric complex coordinate transformations as put forth in [6] and verify that, by using simply a (homogeneous) doubly anisotropic gain-media metamaterial slab, one can still mimic a CSP and produce Gaussian beam. In addition, we show that a Gaussian-like beams can be produced by point sources placed {\it outside} the slab as well [6]. By making use of the extra degrees of freedom (real and imaginary part of the coordinate transformation) provided by CTO, the near-zero requirement on the real part of the resulting constitutive parameters can be relaxed to facilitate potential realization of Gaussian-like beams. We illustrate how beam properties such as peak amplitude and waist location can be controlled by a proper choice of (complex-valued) CTO Jacobian elements. In particular, the beam waist location may be moved bidirectionally by allowing for negative entries in the Jacobian (equivalent to inducing negative refraction effects). These results are then interpreted in light of the ensuing CSP location.


title: Extended solutions for the biadjoint scalar field

abstract: Biadjoint scalar field theories are increasingly important in the study of scattering amplitudes in various string and field theories. Recently, some first exact nonperturbative solutions of biadjoint scalar theory were presented, with a pure power-like form corresponding to isolated monopole-like objects located at the origin of space. In this paper, we find a novel family of extended solutions, involving non-trivial form factors that partially screen the divergent field at the origin. All previous solutions emerge as special cases.


title: A glass-box interactive machine learning approach for solving NP-hard problems with the human-in-the-loop

abstract: The goal of Machine Learning to automatically learn from data, extract knowledge and to make decisions without any human intervention. Such automatic (aML) approaches show impressive success. Recent results even demonstrate intriguingly that deep learning applied for automatic classification of skin lesions is on par with the performance of dermatologists, yet outperforms the average. As human perception is inherently limited, such approaches can discover patterns, e.g. that two objects are similar, in arbitrarily high-dimensional spaces what no human is able to do. Humans can deal only with limited amounts of data, whilst big data is beneficial for aML; however, in health informatics, we are often confronted with a small number of data sets, where aML suffer of insufficient training samples and many problems are computationally hard. Here, interactive machine learning (iML) may be of help, where a human-in-the-loop contributes to reduce the complexity of NP-hard problems. A further motivation for iML is that standard black-box approaches lack transparency, hence do not foster trust and acceptance of ML among end-users. Rising legal and privacy aspects, e.g. with the new European General Data Protection Regulations, make black-box approaches difficult to use, because they often are not able to explain why a decision has been made. In this paper, we present some experiments to demonstrate the effectiveness of the human-in-the-loop approach, particularly in opening the black-box to a glass-box and thus enabling a human directly to interact with an learning algorithm. We selected the Ant Colony Optimization framework, and applied it on the Traveling Salesman Problem, which is a good example, due to its relevance for health informatics, e.g. for the study of protein folding. From studies of how humans extract so much from so little data, fundamental ML-research also may benefit.


title: Force-free collisionless current sheet models with non-uniform temperature and density profiles

abstract: We present a class of one-dimensional, strictly neutral, Vlasov-Maxwell equilibrium distribution functions for force-free current sheets, with magnetic fields defined in terms of Jacobian elliptic functions, extending the results of Abraham-Shrauner (Phys. Plasmas 20, 102117, 2013) to allow for non-uniform density and temperature profiles. To achieve this, we use an approach previously applied to the force-free Harris sheet by Kolotkov et al. (Phys. Plasmas 22, 112902, 2015). In one limit of the parameters, we recover the model of Kolotkov et al., while another limit gives a linear force-free field. We discuss conditions on the parameters such that the distribution functions are always positive, and give expressions for the pressure, density, temperature and bulk-flow velocities of the equilibrium, discussing differences from previous models. We also present some illustrative plots of the distribution function in velocity space.


title: Numerical properties of Koszul connections

abstract: We use the notation EX(S>M), EXF(S>M) and DL(S>M), where M is a smooth manifold and S is a geometric structure. EX(S>M) is the question whether S exists in M. EXF(S>M) is the question whether M admits S-foliations. DL(S>M) is the search of an invariant measuring how M is far from admitting S. For many major geometric structures, those questions are widly open. In this paper, we address EX(S>M), EXF(S>M) and DL(S>M) for affine structure and symplectic structure, left invariant affine structure, left invariant symplectic structure and bi-invariant riemanniann structure in Lie groups


title: Asymptotics of Green function for the linear waves equations in a domain with a non-uniform bottom

abstract: We consider the linear problem for water-waves created by sources on the bottom and the free surface in a 3-D basin having slowly varying profile $z=-D(x)$. The fluid verifies Euler-Poisson equations. These (non-linear) equations have been given a Hamiltonian form by Zakharov, involving canonical variables $(\xi(x,t),\eta(x,t))$ describing the dynamics of the free surface; variables $(\xi,\eta)$ are related by the free surface Dirichlet-to-Neumann (DtN) operator. For a single variable $x\in{\bf R}$ and constant depth, DtN operator was explicitely computed in terms of a convergent series. Here we neglect quadratic terms in Zakharov equations, and consider the linear response to a disturbance of $D(x)$ harmonic in time when the wave-lenght is small compared to the depth of the basin. We solve the Green function problem for a matrix-valued DtN operator, at the bottom and the free-surface.


title: Design, fabrication and characterization of Computer Generated Holograms for anti-counterfeiting applications using OAM beams as light decoders

abstract: In this paper, we present the design, fabrication and optical characterization of computer-generated holograms (CGH) encoding information for light beams carrying orbital angular momentum (OAM). Through the use of a numerical code, based on an iterative Fourier transform algorithm, a phase-only diffractive optical element (PH-DOE) specifically designed for OAM illumination has been computed, fabricated and tested. In order to shape the incident beam into a helicoidal phase profile and generate light carrying phase singularities, a method based on transmission through high-order spiral phase plates (SPPs) has been used. The phase pattern of the designed holographic DOEs has been fabricated using high-resolution Electron-Beam Lithography (EBL) over glass substrates coated with a positive photoresist layer (polymethylmethacrylate). To the best of our knowledge, the present study is the first attempt, in a comprehensive work, to design, fabricate and characterize computer-generated holograms encoding information for structured light carrying OAM and phase singularities. These optical devices appear promising as high-security optical elements for anti-counterfeiting applications.


title: Balance between quantum Markov semigroups

abstract: The concept of balance between two state preserving quantum Markov semigroups on von Neumann algebras is introduced and studied as an extension of conditions appearing in the theory of quantum detailed balance. This is partly motivated by the theory of joinings. Balance is defined in terms of certain correlated states (couplings), with entangled states as a specific case. Basic properties of balance are derived and the connection to correspondences in the sense of Connes is discussed. Some applications and possible applications, including to non-equilibrium statistical mechanics, are briefly explored.


title: Phenomenological formula for CKM matrix and its physical interpretation

abstract: We propose a phenomenological formula relating the Cabibbo--Kobayashi--Masukawa matrix $V_{CKM}$ and quark masses in the form $(\sqrt{m_d} \, \sqrt{m_s} \, \sqrt{m_b}) \propto (\sqrt{m_u} \, \sqrt{m_c} \, \sqrt{m_t})V_{CKM}$. The results of the proposed formula are in agreement with experimental data. Under the constraint of the formula, we show that the invariant amplitude of the charged current weak interactions is maximized.


title: Overview of the CSM concept and its tests

abstract: We recall the motivations for the Compositeness Standard Model(CSM) concept, its precise description, the procedures for its applications and the particular constraints that it requires. We present its most spectacular predictions for typical processes observable at present and future colliders. To the previous results we add the treatment of the inclusive processes $e^+e^- \to H,Z,W~or~t+ anything$.


title: Arc Transitive Maps with underlying Rose Window Graphs

abstract: Let ${\cal M}$ be a map with the underlying graph $\Gamma$. The automorphism group $Aut({\cal M})$ induces a natural action on the set of all vertex-edge-face incident triples, called {\em flags} of ${\cal M}$. The map ${\cal M}$ is said to be a {\em $k$-orbit} map if $Aut({\cal M})$ has $k$ orbits on the set of all flags of ${\cal M}$. It is known that there are seven different classes of $2$-orbit maps, with only four of them corresponding to arc-transitive maps, that is maps for which $Aut{\cal M}$ acts arc-transitively on the underlying graph $\Gamma$. The Petrie dual operator links these four classes in two pairs, one of which corresponds to the chiral maps and their Petrie duals. In this paper we focus on the other pair of classes of $2$-orbit arc-transitive maps. We investigate the connection of these maps to consistent cycles of the underlying graph with special emphasis on such maps of smallest possible valence, namely $4$. We then give a complete classification of such maps whose underlying graphs are arc-transitive Rose Window graphs.


title: Heden's bound on the tail of a vector space partition

abstract: A vector space partition of $\mathbb{F}_q^v$ is a collection of subspaces such that every non-zero vector is contained in a unique element. We improve a lower bound of Heden, in a subcase, on the number of elements of the smallest occurring dimension in a vector space partition. To this end, we introduce the notion of $q^r$-divisible sets of $k$-subspaces in $\mathbb{F}_q^v$. By geometric arguments we obtain non-existence results for these objects, which then imply the improved result of Heden.


title: Bubbles dynamics in microchannels: inertial and capillary migration forces

abstract: This work focuses on the dynamics of a train of unconfined bubbles flowing in microchan- nels. We investigate the transverse position of a train of bubbles, its velocity and the associated pressure drop when flowing in a microchannel depending on the internal forces due to viscosity, inertia and capillarity. Despite the small scales of the system, inertia, referred to as inertial migration force, play a crucial role in determining the transverse equilibrium position of the bubbles. Beside inertia and viscosity, other effects may also affect the transverse migration of bubbles such as the Marangoni surface stresses and the surface deformability. We look at the influence of surfactants in the limit of infinite Marangoni effect which yields rigid bubble interface. The resulting migration force may balance external body forces if present such as buoyancy, Dean or magnetic ones. This balance not only determines the transverse position of the bubbles but, consequently, the surrounding flow structure, which can be determinant for any mass/heat transfer process involved. Finally, we look at the influence of the bubble deformation on the equilibrium position and compare it to the inertial migration force at the centred position, explaining the stable or unstable character of this position accordingly. A systematic study of the influence of the parameters - such as the bubble size, uniform body force, Reynolds and capillary numbers - has been carried out using numerical simulations based on the Finite Element Method, solving the full steady Navier-Stokes equations and its asymptotic counterpart for the limits of small Reynolds and/or capillary numbers.


title: Three-body correlations and conditional forces in suspensions of active hard disks

abstract: Self-propelled Brownian particles show rich out-of-equilibrium physics, for instance, the motility-induced phase separation (MIPS). While decades of studying the structure of liquids have established a deep understanding of passive systems, not much is known about correlations in active suspensions. In this work we derive an approximate analytic theory for three-body correlations and forces in systems of active Brownian disks starting from the many-body Smoluchowski equation. We use our theory to predict the conditional forces that act on a tagged particle and their dependence on the propulsion speed of self-propelled disks. We identify preferred directions of these forces in relation to the direction of propulsion and the positions of the surrounding particles. We further relate our theory to the effective swimming speed of the active disks, which is relevant for the physics of MIPS. To test and validate our theory, we additionally run particle-resolved computer simulations, for which we explicitly calculate the three-body forces. In this context, we discuss the modeling of active Brownian swimmers with nearly hard interaction potentials. We find very good agreement between our simulations and numerical solutions of our theory, especially for the nonequilibrium pair-distribution function. For our analytical results, we carefully discuss their range of validity in the context of the different levels of approximation we applied. This discussion allows us to study the individual contribution of particles to three-body forces and to the emerging structure. Thus, our work sheds light on the collective behavior, provides the basis for further studies of correlations in active suspensions, sheds new light onto the collective behavior, and makes a step towards an emerging liquid state theory.


title: Cross-symmetry breaking of two-component discrete dipolar matter-wave solitons

abstract: We study the spontaneous symmetry breaking of dipolar Bose--Einstein condensates trapped in stacks of two-well systems, which may be effectively built as one-dimensional trapping lattices sliced by a repelling laser sheet. If the potential wells are sufficiently deep, the system is modeled by coupled discrete Gross--Pitaevskii equations with nonlocal self- and cross-interaction terms representing dipole--dipole interactions. When the dipoles are not polarized perpendicular or parallel to the lattice, the cross-interaction is asymmetric, replacing the familiar symmetric two-component solitons with a new species of cross-symmetric or -asymmetric ones. The orientation of the dipole moments and the interwell hopping rate strongly affect the shapes of the discrete two-component solitons as well as the characteristics of the cross-symmetry breaking and the associated phase transition. The sub- and super-critical types of cross-symmetry breaking can be controlled by either the hopping rate between the components or the total norm of the solitons. The effect of the interplay between the contact nonlinearity and the dipole angle on the cross-symmetry breaking is also discussed.


title: 2014 European School of High-Energy Physics

abstract: The European School of High-Energy Physics is intended to give young physicists an introduction to the the- oretical aspects of recent advances in elementary particle physics. These proceedings contain lecture notes on the theory of quantum chromodynamics, Higgs physics, Flavour physics and CP violation, and Supersymmetry.


title: CAS - CERN Accelerator School: Plasma Wake Acceleration

abstract: These proceedings are of interest to staff and students in accelerator laboratories, university departments and companies working in or having an interest in the field of new acceleration techniques. Following introductory lectures on plasma and laser physics, the course covers the different components of a plasma wake accelerator and plasma beam systems. Topical seminars and an overview about alternative new techniques like dielectric accelerators are included. Lectures on the experimental studies and their latest results, on diagnostic tools and state of the art wake acceleration facilities, both present and planned, complement the theoretical part.


title: Towards Spectroscopically Detecting the Global Latitudinal Temperature Variation on the Solar Surface

abstract: A very slight rotation-induced latitudinal temperature variation (presumably on the order of several Kelvin) on the solar surface is theoretically expected. While recent high-precision solar brightness observations reported its detection, confirmation by an alternative approach using the strengths of spectral lines is desirable, for which reducing the noise due to random fluctuation caused by atmospheric inhomogeneity is critical. Towards this difficult task, we carried out a pilot study of spectroscopically investigating the relative variation of temperature (T) at a number of points in the solar circumference region near to the limb (where latitude dependence should show up if any exists) based on the equivalent widths (W) of 28 selected lines in the 5367--5393A and 6075--6100A regions. Our special attention was paid to i) clarifying which kind of lines should be employed and ii) how much precision would be accomplished in practice. We found that lines with strong T-sensitivity (|log W/log T|) should be used and that very weak lines should be avoided because they inevitably suffer large relative fluctuations (Delta_W/W). Our analysis revealed that a precision of Delta_T/T ~ 0.003 (corresponding to ~15K) can be achieved at best by a spectral line with comparatively large |log W/log T|, though this can possibly be further improved if other more suitable lines are used. Accordingly, if many such favorable lines could be measured with sub-% precision of Delta_W/W and by averaging the resulting Delta_T/T from each line, the random noise would eventually be reduced down to <~1K and detection of very subtle amount of global T-gradient might be possible.


title: Extracting spectroscopic molecular parameters from short pulse photo-electron angular distributions

abstract: Using a quantum wave packet simulation including the nuclear and electronic degrees of freedom, we investigate the femtosecond and picosecond energy- and angle-resolved photoelectron spectra of the E($^1\Sigma_g^+$) electronic state of Li$_2$. We find that the angular distributions of the emitted photoelectrons depend strongly on the pulse duration in the regime of ultrashort laser pulses. This effect is illustrated by the extraction of a time-dependent asymmetry parameter whose variation with pulse duration can be explained by an incoherent average over different ion rotational quantum numbers. We then derive for the variation of the asymmetry parameter a simple analytical formula, which can be used to extract the asymptotic CW asymmetry parameters of individual transitions from measurements performed with ultra-short pulses.


title: Asymptotic behaviour of randomised fractional volatility models

abstract: We study the asymptotic behaviour of a class of small-noise diffusions driven by fractional Brownian motion, with random starting points. Different scalings allow for different asymptotic properties of the process (small-time and tail behaviours in particular). In order to do so, we extend some results on sample path large deviations for such diffusions. As an application, we show how these results characterise the small-time and tail estimates of the implied volatility for rough volatility models, recently proposed in mathematical finance.


title: Solving and Sampling with Many Solutions: Satisfiability and Other Hard Problems

abstract: We investigate parameterizing hard combinatorial problems by the size of the solution set compared to all solution candidates. Our main result is a uniform sampling algorithm for satisfying assignments of 2-CNF formulas that runs in expected time $O^*(\varepsilon^{-0.617})$ where $\varepsilon$ is the fraction of assignments that are satisfying. This improves significantly over the trivial sampling bound of expected $\Theta^*(\varepsilon^{-1})$, and on all previous algorithms whenever $\varepsilon = \Omega(0.708^n)$. We also consider algorithms for 3-SAT with an $\varepsilon$ fraction of satisfying assignments, and prove that it can be solved in $O^*(\varepsilon^{-2.27})$ deterministic time, and in $O^*(\varepsilon^{-0.936})$ randomized time. Finally, to further demonstrate the applicability of this framework, we also explore how similar techniques can be used for vertex cover problems.


title: Analytic approximation for eigenvalues of a class of $\mathcal{PT}$ symmetric Hamiltonians

abstract: An analytical approximation for the eigenvalues of $\mathcal{PT}$ symmetric Hamiltonian $\mathsf{H} = -d^{2}/dx^{2} - (\mathrm{i}x)^{\epsilon+2}$, $\epsilon > -1$ is developed via simple basis sets of harmonic-oscillator wave functions with variable frequencies and equilibrium positions. We demonstrate that our approximation provides high accuracy for any given energy level for all values of $\epsilon > -1$.


title: Thermoelectric efficiency of nanoscale devices in the linear regime

abstract: We study quantum transport through two-terminal nanoscale devices in contact with two particle reservoirs at different temperatures and chemical potentials. We discuss the general expressions controlling the electric charge current, heat currents and the efficiency of energy transmutation in steady conditions in the linear regime. With focus in the parameter domain where the electron system acts as a power-generator, we elaborate workable expressions for optimal efficiency and thermoelectric pameters of nanoscale devices. The general concepts are set at work in the paradigmatic cases of Lorentzian resonances and antiresonances, and the encompassing Fano transmission function: the treatments are fully analytic, in terms of the trigamma functions and Bernoulli numbers. From the general curves here reported describing transport through the above model transmission functions, useful guidelines for optimal efficiency and thermopower can be inferred for engineering nanoscale devices in energy regions where they show similar transmission functions.


title: A Unified View-Graph Selection Framework for Structure from Motion

abstract: View-graph is an essential input to large-scale structure from motion (SfM) pipelines. Accuracy and efficiency of large-scale SfM is crucially dependent on the input view-graph. Inconsistent or inaccurate edges can lead to inferior or wrong reconstruction. Most SfM methods remove `undesirable' images and pairs using several, fixed heuristic criteria, and propose tailor-made solutions to achieve specific reconstruction objectives such as efficiency, accuracy, or disambiguation. In contrast to these disparate solutions, we propose a single optimization framework that can be used to achieve these different reconstruction objectives with task-specific cost modeling. We also construct a very efficient network-flow based formulation for its approximate solution. The abstraction brought on by this selection mechanism separates the challenges specific to datasets and reconstruction objectives from the standard SfM pipeline and improves its generalization. This paper demonstrates the application of the proposed view-graph framework with standard SfM pipeline for two particular use-cases, (i) accurate and ghost-free reconstructions of highly ambiguous datasets using costs based on disambiguation priors, and (ii) accurate and efficient reconstruction of large-scale Internet datasets using costs based on commonly used priors.


title: Evolution of the interfacial perpendicular magnetic anisotropy constant of the Co$_2$FeAl/MgO interface upon annealing

abstract: We investigate thickness series of films of the Heusler alloy Co$_2$FeAl in order to study the effect of annealing on the interface with a MgO layer and on the bulk magnetic properties. Our results reveal that while the perpendicular interface anisotropy constant $K^{\perp}_{\rm S}$ is zero for the as-deposited samples, its value increases with annealing up to a value of $1.14\, \pm \,0.07$~mJ/m$^2$ for the series annealed at 320$^{\rm o}$C and of $2.07\, \pm \,0.7$~mJ/m$^2$ for the 450$^{\rm o}$C annealed series owing to a strong modification of the interface during the thermal treatment. This large value ensures a stabilization of a perpendicular magnetization orientation for a thickness below 1.7~nm. The data additionally shows that the in-plane biaxial anisotropy constant has a different evolution with thickness in as-deposited and annealed systems. The Gilbert damping parameter $\alpha$ shows minima for all series for a thickness of 40~nm and an absolute minimum value of $2.8\pm0.1\cdot10^{-3}$. The thickness dependence is explained in terms of an inhomogenous magnetization state generated by the interplay between the different anisotropies of the system and by crystalline disorder.


title: Constructing the virtual fundamental class of a Kuranishi atlas

abstract: Consider a space $X$, such as a compact space of $J$-holomorphic stable maps, that is the zero set of a Kuranishi atlas. This note explains how to define the virtual fundamental class of $X$ by representing $X$ via the zero set of a map $S_M: M\to E$, where $E$ is a finite dimensional vector space and the domain $M$ is an oriented, weighted branched topological manifold. Moreover, $S_M$ is equivariant under the action of the global isotropy group $\Gamma$ on $M$ and $E$. This tuple $(M,E, \Gamma, S_M)$ together with a homeomorphism $S_M^{-1}(0)/\Gamma \to X$ forms a single finite dimensional model (or chart) for $X$. The construction assumes only that the atlas satisfies a topological version of the index condition that can be obtained from a standard, rather than a smooth, gluing theorem. However if $X$ is presented as the zero set of an sc-Fredholm operator on a strong polyfold bundle, we outline a much more direct construction of the branched manifold $M$ that uses an sc-smooth partition of unity.


title: Fourier analysis of He 4471/Mg 4481 line profiles for separating rotational velocity and axial inclination in rapidly-rotating B-type stars

abstract: While the effect of rotation on spectral lines is complicated in rapidly-rotating stars because of the appreciable gravity-darkening effect differing from line to line, it is possible to make use of this line-dependent complexity to separately determine the equatorial rotation velocity (ve) and the inclination angle (i) of rotational axis. Although line-widths of spectral lines were traditionally used for this aim, we tried in this study to apply the Fourier method, which utilizes the unambiguously determinable first-zero frequency (sigma1) in the Fourier transform of line profile. Equipped with this technique, we analyzed the profiles of HeI 4471 and MgII 4481 lines of six rapidly-rotating (vesini~150-300km/s) late B-type stars, while comparing them with the theoretical profiles simulated on a grid of models computed for various combination of (ve, i). According to our calculation, sigma1 tends to be larger than the classical value for given vesini. This excess progressively grows with an increase in ve, and is larger for the He line than the Mg line, which leads to sigma1He > sigma1Mg. It was shown that ve and i are separately determinable from the intersection of two loci (sets of solutions reproducing the observed sigma1 for each line) on the ve vs. i plane. Yet, line profiles alone are not sufficient for their unique discrimination, for which photometric information (such as colors) needs to be simultaneously employed.


title: Analytical calculation for the gluon fragmentation into spin-triplet S-wave quarkonium

abstract: Fragmentation is the dominant mechanism for hadron production with high transverse momentum. For spin-triplet S-wave heavy quarkonium production, contribution of gluon fragmenting to color-singlet channel has been numerically calculated since 1993. However, there is still no analytic expression available up to now because of its complexity. In this paper, we calculate both polarization-summed and polarized fragmentation functions of gluon fragmenting to a heavy quark-antiquark pair with quantum number $^3S_1^{[1]}$. Our calculations are performed in two different frameworks. One is the widely used nonrelativistic QCD factorization, and the other is the newly proposed soft gluon factorization. In either case, we calculate at both leading order and next-to-leading order in velocity expansion. All of our final results are presented in terms of compact analytic expressions.


title: Efficient pattern matching in degenerate strings with the Burrows-Wheeler transform

abstract: A degenerate or indeterminate string on an alphabet $\Sigma$ is a sequence of non-empty subsets of $\Sigma$. Given a degenerate string $t$ of length $n$, we present a new method based on the Burrows--Wheeler transform for searching for a degenerate pattern of length $m$ in $t$ running in $O(mn)$ time on a constant size alphabet $\Sigma$. Furthermore, it is a hybrid pattern-matching technique that works on both regular and degenerate strings. A degenerate string is said to be conservative if its number of non-solid letters is upper-bounded by a fixed positive constant $q$; in this case we show that the search complexity time is $O(qm^2)$. Experimental results show that our method performs well in practice.


title: A computer simulation of the Volga River hydrological regime: a problem of water-retaining dam optimal location

abstract: We investigate of a special dam optimal location at the Volga river in area of the Akhtuba left sleeve beginning (7 \, km to the south of the Volga Hydroelectric Power Station dam). We claim that a new water-retaining dam can resolve the key problem of the Volga-Akhtuba floodplain related to insufficient water amount during the spring flooding due to the overregulation of the Lower Volga. By using a numerical integration of Saint-Vacant equations we study the water dynamics across the northern part of the Volga-Akhtuba floodplain with taking into account its actual topography. As the result we found an amount of water $V_A$ passing to the Akhtuba during spring period for a given water flow through the Volga Hydroelectric Power Station (so-called hydrograph which characterises the water flow per unit of time). By varying the location of the water-retaining dam $ x_d, y_d $ we obtained various values of $V_A (x_d, y_d) $ as well as various flow spatial structure on the territory during the flood period. Gradient descent method provide us the dam coordinated with the maximum value of ${V_A}$. Such approach to the dam location choice let us to find the best solution, that the value $V_A$ increases by a factor of 2. Our analysis demonstrate a good potential of the numerical simulations in the field of hydraulic works.


title: Coherence evolution and transfer supplemented by state-restoring

abstract: The evolution of quantum coherences comes with a set of conservation laws provided that the Hamiltonian governing this evolution conserves the spin-excitation number. At that, coherences do not intertwist during the evolution. Using the transmission line and the receiver in the initial ground state we can transfer the coherences to the receiver without interaction between them, { although the matrix elements contributing to each particular coherence intertwist in the receiver's state. } Therefore we propose a tool based on the unitary transformation at the receiver side to { untwist these elements and thus} restore (at least partially) the structure of the sender's initial density matrix. A communication line with two-qubit sender and receiver is considered as an example of implementation of this technique.


title: Reactor Fuel Fraction Information on the Antineutrino Anomaly

abstract: We analyzed the evolution data of the Daya Bay reactor neutrino experiment in terms of short-baseline active-sterile neutrino oscillations taking into account the theoretical uncertainties of the reactor antineutrino fluxes. We found that oscillations are disfavored at $2.6\sigma$ with respect to a suppression of the $^{235}\text{U}$ reactor antineutrino flux and at $2.5\sigma$ with respect to variations of the $^{235}\text{U}$ and $^{239}\text{Pu}$ fluxes. On the other hand, the analysis of the rates of the short-baseline reactor neutrino experiments favor active-sterile neutrino oscillations and disfavor the suppression of the $^{235}\text{U}$ flux at $3.1\sigma$ and variations of the $^{235}\text{U}$ and $^{239}\text{Pu}$ fluxes at $2.8\sigma$. We also found that both the Daya Bay evolution data and the global rate data are well-fitted with composite hypotheses including variations of the $^{235}\text{U}$ or $^{239}\text{Pu}$ fluxes in addition to active-sterile neutrino oscillations. A combined analysis of the Daya Bay evolution data and the global rate data shows a slight preference for oscillations with respect to variations of the $^{235}\text{U}$ and $^{239}\text{Pu}$ fluxes. However, the best fits of the combined data are given by the composite models, with a preference for the model with an enhancement of the $^{239}\text{Pu}$ flux and relatively large oscillations.


title: Emission from accelerating jets in gamma-ray bursts: Radiation dominated flows with increasing mass outflow rates

abstract: We study the narrowest spectra expected from gamma-ray bursts. We present an analytical function for the spectrum that is emitted from the photosphere of a radiation-dominated flow that is under acceleration. This is the narrowest possible spectrum and it differs from a Planck function. We also present numerical spectra from photospheres occurring during the transition into the coasting phase of the flow. Using these spectral models, we reanalyse Fermi observations of GRB100507 and GRB101219, which both have been reported to have very narrow spectra. The bursts can be fitted by the spectral models: For GRB101219 the spectrum is consistent with the photosphere occurring below or close to the saturation radius, while for GRB100507 the photosphere position relative to the saturation radius can be determined as a function of time. In the latter case, we find that the photosphere initially occurs in the acceleration phase and thereafter transitions into the coasting phase. We also find that this transition occurs at the same time as the change in observed cooling behaviour: the temperature is close to constant before the break and decays after. We argue that such a transition can be explained by an increasing mass outflow rate. Both analysed bursts thus give strong evidence that the jets are (initially) radiation dominated.


title: Long range forces in a performance portable Molecular Dynamics framework

abstract: Molecular Dynamics (MD) codes predict the fundamental properties of matter by following the trajectories of a collection of interacting model particles. To exploit diverse modern manycore hardware, efficient codes must use all available parallelism. At the same time they need to be portable and easily extendible by the domain specialist (physicist/chemist) without detailed knowledge of this hardware. To address this challenge, we recently described a new Domain Specific Language (DSL) for the development of performance portable MD codes based on a "Separation of Concerns": a Python framework automatically generates efficient parallel code for a range of target architectures. Electrostatic interactions between charged particles are important in many physical systems and often dominate the runtime. Here we discuss the inclusion of long-range interaction algorithms in our code generation framework. These algorithms require global communications and careful consideration has to be given to any impact on parallel scalability. We implemented an Ewald summation algorithm for electrostatic forces, present scaling comparisons for different system sizes and compare to the performance of existing codes. We also report on further performance optimisations delivered with OpenMP shared memory parallelism.


title: Lecture hall partitions and the affine hyperoctahedral group

abstract: In 1997 Bousquet-M\'elou and Eriksson introduced lecture hall partitions as the inversion vectors of elements of the parabolic quotient $\widetilde{C}/C$. We provide a new view of their correspondence that allows results in one domain to be translated into the other. We determine the equivalence between combinatorial statistics in each domain and use this correspondence to translate certain generating function formulas on lecture hall partitions to new observations about $\widetilde{C}/C$.


title: Pico-charged intermediate particles rescue dark matter interpretation of 511 keV signal

abstract: Various alleged indirect dark matter search signals, such as the 511 keV line from galaxy center or the PAMELA/AMS02 signal, are often challenged by the absence of corresponding signal from dwarf galaxies and/or from the absence of an impact on CMB through delayed recombination. We propose a novel scenario that can avoid these bounds based on the decay of dark matter, X, to a pair of intermediate particles C and \bar{C} with a lifetime much greater than the age of universe. The annihilation of these intermediate particles eventually leads to a dark matter signal. The bounds from CMB can be easily avoided by the fact that at the time of recombination, not enough C particles had been accumulated. In order to keep C particles from leaving the galaxy, we assume the particles have a small electric charge so in the galactic disk, the magnetic field keeps the C particles in the vicinity of their production. However, they can escape the dwarf galaxies and the dark matter halo where the magnetic field is weak, leading to null signal from these regions. The small charge can have interesting consequences including a signal in direct dark matter searches.


title: Impact of New Physics on the EW vacuum stability in a curved spacetime background

abstract: It has been recently shown that, contrary to an intuitive decoupling argument, the presence of new physics at very large energy scales (say around the Planck scale) can have a strong impact on the electroweak vacuum lifetime. In particular, the vacuum could be totally destabilized. This study was performed in a flat spacetime background, and it is important to extend the analysis to curved spacetime since these are Planckian-physics effects. It is generally expected that under these extreme conditions gravity should totally quench the formation of true vacuum bubbles, thus washing out the destabilizing effect of new physics. In this work we extend the analysis to curved spacetime and show that, although gravity pushes toward stabilization, the destabilizing effect of new physics is still (by far) the dominating one. In order to get model independent results, high energy new physics is parametrized in two different independent ways: as higher order operators in the Higgs field, or introducing new particles with very large masses. The destabilizing effect is observed in both cases, hinting at a general mechanism that does not depend on the parametrization details for new physics, thus maintaining the results obtained from the analysis performed in flat spacetime.


title: Feedback computability on Cantor space

abstract: We introduce the notion of feedback computable functions from $2^\omega$ to $2^\omega$, extending feedback Turing computation in analogy with the standard notion of computability for functions from $2^\omega$ to $2^\omega$. We then show that the feedback computable functions are precisely the effectively Borel functions. With this as motivation we define the notion of a feedback computable function on a structure, independent of any coding of the structure as a real. We show that this notion is absolute, and as an example characterize those functions that are computable from a Gandy ordinal with some finite subset distinguished.


title: The role of symmetry in driven propulsion at low Reynolds number

abstract: We theoretically and experimentally investigate low-Reynolds-number propulsion of geometrically achiral planar objects that possess a dipole moment and that are driven by a rotating magnetic field. Symmetry considerations (involving parity, $\widehat{P}$, and charge conjugation, $\widehat{C}$) establish correspondence between propulsive states depending on orientation of the dipolar moment. Although basic symmetry arguments do not forbid individual symmetric objects to efficiently propel due to spontaneous symmetry breaking, they suggest that the average ensemble velocity vanishes. Some additional arguments show, however, that highly symmetrical ($\widehat{P}$-even) objects exhibit no net propulsion while individual less symmetrical ($\widehat{C}\widehat{P}$-even) propellers do propel. Particular magnetization orientation, rendering the shape $\widehat{C}\widehat{P}$-odd, yields unidirectional motion typically associated with chiral structures, such as helices. If instead of a structure with a permanent dipole we consider a polarizable object, some of the arguments have to be modified. For instance, we demonstrate a truly achiral ($\widehat{P}$- and $\widehat{C}\widehat{P}$-even) planar shape with an induced electric dipole that can propel by electro-rotation. We thereby show that chirality is not essential for propulsion due to rotation-translation coupling at low Reynolds number.


title: Automatic Segmentation and Disease Classification Using Cardiac Cine MR Images

abstract: Segmentation of the heart in cardiac cine MR is clinically used to quantify cardiac function. We propose a fully automatic method for segmentation and disease classification using cardiac cine MR images. A convolutional neural network (CNN) was designed to simultaneously segment the left ventricle (LV), right ventricle (RV) and myocardium in end-diastole (ED) and end-systole (ES) images. Features derived from the obtained segmentations were used in a Random Forest classifier to label patients as suffering from dilated cardiomyopathy, hypertrophic cardiomyopathy, heart failure following myocardial infarction, right ventricular abnormality, or no cardiac disease. The method was developed and evaluated using a balanced dataset containing images of 100 patients, which was provided in the MICCAI 2017 automated cardiac diagnosis challenge (ACDC). The segmentation and classification pipeline were evaluated in a four-fold stratified cross-validation. Average Dice scores between reference and automatically obtained segmentations were 0.94, 0.88 and 0.87 for the LV, RV and myocardium. The classifier assigned 91% of patients to the correct disease category. Segmentation and disease classification took 5 s per patient. The results of our study suggest that image-based diagnosis using cine MR cardiac scans can be performed automatically with high accuracy.


title: Phase-error estimation and image reconstruction from digital-holography data using a Bayesian framework

abstract: The estimation of phase errors from digital-holography data is critical for applications such as imaging or wave-front sensing. Conventional techniques require multiple i.i.d. data and perform poorly in the presence of high noise or large phase errors. In this paper we propose a method to estimate isoplanatic phase errors from a single data realization. We develop a model-based iterative reconstruction algorithm which computes the maximum a posteriori estimate of the phase and the speckle-free object reflectance. Using simulated data, we show that the algorithm is robust against high noise and strong phase errors.


title: Three-dimensional planar model estimation using multi-constraint knowledge based on k-means and RANSAC

abstract: Plane model extraction from three-dimensional point clouds is a necessary step in many different applications such as planar object reconstruction, indoor mapping and indoor localization. Different RANdom SAmple Consensus (RANSAC)-based methods have been proposed for this purpose in recent years. In this study, we propose a novel method-based on RANSAC called Multiplane Model Estimation, which can estimate multiple plane models simultaneously from a noisy point cloud using the knowledge extracted from a scene (or an object) in order to reconstruct it accurately. This method comprises two steps: first, it clusters the data into planar faces that preserve some constraints defined by knowledge related to the object (e.g., the angles between faces); and second, the models of the planes are estimated based on these data using a novel multi-constraint RANSAC. We performed experiments in the clustering and RANSAC stages, which showed that the proposed method performed better than state-of-the-art methods.


title: Direct nonlinear Fourier transform algorithms for the computation of solitonic spectra in focusing nonlinear Schr\"odinger equation

abstract: Starting from a comparison of some established numerical algorithms for the computation of the eigenvalues (discrete or solitonic spectrum) of the non-Hermitian version of the Zakharov-Shabat spectral problem, this article delivers new algorithms that combine the best features of the existing ones and thereby allays their relative weaknesses. Our algorithm is modelled within the remit of the so-called direct nonlinear Fourier transform (NFT) associated with the focusing nonlinear Schr\"odinger equation. First, we present the data for the calibration of methods comparing the relative errors associated with the computation of the continuous NF spectrum. Then each method is paired with different numerical algorithms for finding zeros of a complex-valued function to obtain the eigenvalues. Next we describe a new class of methods based on the contour integrals evaluation for the efficient search of eigenvalues. After that, we introduce a new hybrid method, one of our main results: the method combines the advances of contour integral approach and makes use of the iterative algorithms at its second stage for the refined eigenvalues search. The veracity of our new hybrid algorithm is established by estimating the convergence speed and accuracy across three independent test profiles. Along with the development of a new approach for the computation of the eigenvalues, our study also addresses the problem of computation of the so-called norming constants associated with the eigenvalues. We show that our formalism effectively amounts to accurate and fast enough computation of residues of the reflection coefficient in the upper complex half-plane of the spectral parameter.


title: Solar wind magnetic field background spectrum from fluid to kinetic scales

abstract: The solar wind is highly structured in fast and slow flows. These two dynamical regimes remarkably differ not only for the average values of magnetic field and plasma parameters but also for the type of fluctuations they transport. Fast wind is characterized by large amplitude, incompressible fluctuations, mainly Alfv\'{e}nic, slow wind is generally populated by smaller amplitude and less Alfv\'{e}nic fluctuations, mainly compressive. The typical corotating fast stream is characterized by a stream interface, a fast wind region and a slower rarefaction region formed by the trailing expansion edge of the stream. Moving {between these two regions}, from faster to slower wind, we observe the following behavior: a) the power level of magnetic fluctuations within the inertial range largely decreases, keeping the typical Kolmogorov scaling; b) at proton scales, for about one decade right beyond the high frequency break, the spectral index becomes flatter and flatter towards a value around -2.7; c) at higher frequencies, before the electron scales, the spectral index remains around -2.7 and, {based on suitable observations available for $4$ corotating streams}, the power level does not change, irrespective of the flow speed. All these spectral features, characteristic of high speed streams, suggest the existence of a sort of magnetic field background spectrum. This spectrum would be common to both faster and slower wind but, any time the observer would cross the inner part of a fluxtube channeling the faster wind into the interplanetary space, a turbulent and large amplitude Alfv\'{e}nic spectrum would be superposed to it.


title: Preselection via Classification: A Case Study on Evolutionary Multiobjective Optimization

abstract: In evolutionary algorithms, a preselection operator aims to select the promising offspring solutions from a candidate offspring set. It is usually based on the estimated or real objective values of the candidate offspring solutions. In a sense, the preselection can be treated as a classification procedure, which classifies the candidate offspring solutions into promising ones and unpromising ones. Following this idea, we propose a classification based preselection (CPS) strategy for evolutionary multiobjective optimization. When applying classification based preselection, an evolutionary algorithm maintains two external populations (training data set) that consist of some selected good and bad solutions found so far; then it trains a classifier based on the training data set in each generation. Finally it uses the classifier to filter the unpromising candidate offspring solutions and choose a promising one from the generated candidate offspring set for each parent solution. In such cases, it is not necessary to estimate or evaluate the objective values of the candidate offspring solutions. The classification based preselection is applied to three state-of-the-art multiobjective evolutionary algorithms (MOEAs) and is empirically studied on two sets of test instances. The experimental results suggest that classification based preselection can successfully improve the performance of these MOEAs.


title: Statistical Reconstruction of Microstructures Using Entropic Descriptors

abstract: We report a multiscale approach of broad applicability to stochastic reconstruction of multiphase materials, including porous ones. The approach devised uses an optimization method, such as the simulated annealing (SA) and the so-called entropic descriptors (EDs). For a binary pattern, they quantify spatial inhomogeneity or statistical complexity at discrete length-scales. The EDs extract dissimilar structural information to that given by two-point correlation functions (CFs). Within the SA, we use an appropriate cost function consisting of EDs or comprised of EDs and CFs. It was found that the stochastic reconstruction is computationally efficient when we begin with a preliminary synthetic configuration having in part desirable features. Another option is low-cost approximate reconstructing of the entire multiphase medium beyond the SA technique. The information included in the target ED-curve was utilized for this purpose. For a given volume fraction the low-cost trial microstructures can be generated in two ways. In the first one, applied to ceramics and carbonate samples, the interpenetrating spheres generate a number of trial configurations. In the second one, with phase-EDs, here used to the sandstone sample, the overlapping superspheres do it. Both methods use a radius determined from the EDs-linked two-exponent power-law. However, the supersphere deformation parameter allows controlling of the spatial inhomogeneity of prototypical microstructures. At last, even for a hypothetical ED-curve (under reasonable assumptions), the specific microstructure can be found, if it is realizable for a given volume fraction. In general, the EDs-based methods offer a compromise between computational efficiency and the accuracy of reconstructions.


title: Some practical versions of boundary variation diminishing (BVD) algorithm

abstract: This short note presents some variant schemes of boundary variation diminishing (BVD) algorithm in one dimension with the results of numerical tests for linear advection equation to facilitate practical use. In spite of being presented in 1D fashion, all the schemes are simple and easy to implement in multi-dimensions on structured and unstructured grids for nonlinear and system equations.


title: Bi-stability of magnetic skyrmions in ultrathin multilayer nanodots induced by magnetostatic interaction

abstract: We report the results of simulations of magnetic skyrmion stability in ultrathin magnetic multilayer nanodots with interfacial Dzyaloshinskii-Moriya exchange interaction (DMI). We found that in presence of the lateral confinement the magnetostatic energy significantly influences the skyrmion stability and leads to stabilization of large radius skyrmion even at low values of the DMI strength, in addition to small radius skyrmion stabilized by DMI. In particular, stabilization of the skyrmion state with two different radii (bi-stability) is found in dipolarly-coupled (Pt/Co/Ir)$_{n}$ circular nanodots with the number of repeats of the unit cell $n$ = 3 and 5. The bi-stability range is located at the DMI strength of 0.9-1.1 mJ/m$^2$ or at the total Co-layer thickness of 2.2-2.6 nm.


title: Early-universe cosmology in Einstein-scalar-Gauss-Bonnet gravity

abstract: Theoretical arguments and cosmological observations suggest that Einstein's theory of general relativity needs to be modified at high energies. One of the best motivated higher-curvature extensions of general relativity is Einstein-scalar-Gauss-Bonnet gravity, in which a scalar field is coupled to quadratic curvature invariants. This theory is inspired by an effective string-theory model and its predictions dramatically differ from Einstein's theory in high-curvature regions - such as the interior of black holes and the early universe - where it aims at resolving curvature singularities. In this work we derive cosmological solutions in Einstein-scalar-Gauss-Bonnet gravity for quadratic and for exponential coupling functions, and for any spatial curvature. We discuss already known solutions and find new nonsingular, inflationary, and bouncing solutions. We study the linear stability of these solutions and the absence of ghosts, finding that all the aforementioned solutions are unstable against tensor perturbations. We then introduce a simple, quadratic potential for the scalar field. In some cases the presence of a mass term cures the tensor instability. The proposed model is therefore a viable and attractive candidate for inflation, and one in which the scalar field is naturally provided by the gravitational sector.


title: Robust causal structure learning with some hidden variables

abstract: We introduce a new method to estimate the Markov equivalence class of a directed acyclic graph (DAG) in the presence of hidden variables, in settings where the underlying DAG among the observed variables is sparse, and there are a few hidden variables that have a direct effect on many of the observed ones. Building on the so-called low rank plus sparse framework, we suggest a two-stage approach which first removes the effect of the hidden variables, and then estimates the Markov equivalence class of the underlying DAG under the assumption that there are no remaining hidden variables. This approach is consistent in certain high-dimensional regimes and performs favourably when compared to the state of the art, both in terms of graphical structure recovery and total causal effect estimation.


title: Existence, uniqueness and ergodic properties for time-homogeneous It\^o-SDEs with locally integrable drifts and Sobolev diffusion coefficients

abstract: Using elliptic and parabolic regularity results in $L^p$-spaces and generalized Dirichlet form theory, we construct for every starting point weak solutions to SDEs in $\mathbb{R}^d$ up to their explosion times including the following conditions. For arbitrary but fixed $p>d$ the diffusion coefficient $A=(a_{ij})_{1\le i,j\le d}$ is locally uniformly strictly elliptic with functions $a_{ij}\in H^{1,p}_{loc}(\mathbb{R}^d)$ and the drift coefficient $\mathbf{G}=(g_1,\dots, g_d)$ consists of functions $g_i\in L^p_{loc}(\mathbb{R}^d)$. The solution originates by construction from a Hunt process with continuous sample paths on the one-point compactification of $\mathbb{R}^d$ and the corresponding SDE is by a known local well-posedness result pathwise unique up to an explosion time. Just under the given assumptions we show irreducibility and the strong Feller property on $L^{1}(\mathbb{R}^d,m)+L^{\infty}(\mathbb{R}^d,m)$ of its transition function, and the strong Feller property on $L^{q}(\mathbb{R}^d,m)+L^{\infty}(\mathbb{R}^d,m)$, $q=\frac{dp}{d+p}\in (d/2,p/2)$, of its resolvent, which both include the classical strong Feller property. We present moment inequalities and classical-like non-explosion criteria for the solution which lead to pathwise uniqueness results up to infinity under presumably optimal general non-explosion conditions. We further present explicit conditions for recurrence and ergodicity, including existence as well as uniqueness of invariant probability measures.


title: Intrinsic spectra of H.E.S.S. blazars : what would we see without EBL absorption

abstract: The vast majority of extragalactic sources detected in the very high energy (E> 100 GeV) domaine are active galactic nuclei (AGN) located at cosmological distances. During their travel towards Earth, the emitted gamma-rays suffer from absorption by the extragalactic background light (EBL). The density of the EBL is not very well constrained by direct or indirect measurement which leads to uncertainties on the intrinsic spectrum of the sources. High-quality AGN spectra obtained with the High Energy Stereoscopic System (H.E.S.S.) have been used to perform a model-independent measurement of the EBL spectral energy distribution. While the precision of this measurement remains limited, it reflects the sensitivity of H.E.S.S. alone to the EBL and offers the possibility to access the intrinsic spectra of AGNs in a consistent manner, taking into account the derived uncertainties on the EBL spectral energy distribution. In this contribution, we study the intrinsic spectra as measured by H.E.S.S. and by the Fermi Large Area Telescope (LAT) of several blazars of the H.E.S.S. sky. This provides the opportunity to have new insight into the emission processes at play in the jets of AGN. The data presented consist of monitoring data of quiescent state of blazars and also bright blazar flares, such as PKS~2155-304, Mrk~421, recorded by H.E.S.S.


title: The imprint of neutrinos on clustering in redshift-space

abstract: (abridged) We investigate the signatures left by the cosmic neutrino background on the clustering of matter, CDM+baryons and halos in redshift-space using a set of more than 1000 N-body and hydrodynamical simulations with massless and massive neutrinos. We find that the effect neutrinos induce on the clustering of CDM+baryons in redshift-space on small scales is almost entirely due to the change in $\sigma_8$. Neutrinos imprint a characteristic signature in the quadrupole of the matter (CDM+baryons+neutrinos) field on small scales, that can be used to disentangle the effect of $\sigma_8$ and $M_\nu$. We show that the effect of neutrinos on the clustering of halos is very different, on all scales, to the one induced by $\sigma_8$. We find that the effects of neutrinos of the growth rate of CDM+baryons ranges from $\sim0.3\%$ to $2\%$ on scales $k\in[0.01, 0.5]~h{\rm Mpc}^{-1}$ for neutrinos with masses $M_\nu \leqslant 0.15$ eV. We compute the bias between the momentum of halos and the momentum of CDM+baryon and find it to be 1 on large scales for all models with massless and massive neutrinos considered. This point towards a velocity bias between halos and total matter on large scales that it is important to account for in order to extract unbiased neutrino information from velocity/momentum surveys such as kSZ observations. We show that baryonic effects can affect the clustering of matter and CDM+baryons in redshift-space by up to a few percent down to $k=0.5~h{\rm Mpc}^{-1}$. We find that hydrodynamics and astrophysical processes, as implemented in our simulations, only distort the relative effect that neutrinos induce on the anisotropic clustering of matter, CDM+baryons and halos in redshift-space by less than $1\%$. Thus, the effect of neutrinos in the fully non-linear regime can be written as a transfer function with very weak dependence on astrophysics.


title: Deep MR to CT Synthesis using Unpaired Data

abstract: MR-only radiotherapy treatment planning requires accurate MR-to-CT synthesis. Current deep learning methods for MR-to-CT synthesis depend on pairwise aligned MR and CT training images of the same patient. However, misalignment between paired images could lead to errors in synthesized CT images. To overcome this, we propose to train a generative adversarial network (GAN) with unpaired MR and CT images. A GAN consisting of two synthesis convolutional neural networks (CNNs) and two discriminator CNNs was trained with cycle consistency to transform 2D brain MR image slices into 2D brain CT image slices and vice versa. Brain MR and CT images of 24 patients were analyzed. A quantitative evaluation showed that the model was able to synthesize CT images that closely approximate reference CT images, and was able to outperform a GAN model trained with paired MR and CT images.


title: Structure of the solar photosphere studied from the radiation hydrodynamics code ANTARES

abstract: The ANTARES radiation hydrodynamics code is capable of simulating the solar granulation in detail unequaled by direct observation. We introduce a state-of-the-art numerical tool to the solar physics community and demonstrate its applicability to model the solar granulation. The code is based on the weighted essentially non-oscillatory finite volume method and by its implementation of local mesh refinement is also capable of simulating turbulent fluids. While the ANTARES code already provides promising insights into small-scale dynamical processes occurring in the quiet-Sun photosphere, it will soon be capable of modeling the latter in the scope of radiation magnetohydrodynamics. In this first preliminary study we focus on the vertical photospheric stratification by examining a 3-D model photosphere with an evolution time much larger than the dynamical timescales of the solar granulation and of particular large horizontal extent corresponding to $25\!" \!\! \times \, 25\!"$ on the solar surface to smooth out horizontal spatial inhomogeneities separately for up- and downflows. The highly resolved Cartesian grid thereby covers $\sim 4~\mathrm{Mm}$ of the upper convection zone and the adjacent photosphere. Correlation analysis, both local and two-point, provides a suitable means to probe the photospheric structure and thereby to identify several layers of characteristic dynamics: The thermal convection zone is found to reach some ten kilometers above the solar surface, while convectively overshooting gas penetrates even higher into the low photosphere. An $\approx 145\,\mathrm{km}$ wide transition layer separates the convective from the oscillatory layers in the higher photosphere.


title: A conformally invariant gap theorem in Yang-Mills theory

abstract: We show a sharp conformally invariant gap theorem for Yang-Mills connections in dimension 4 by exploiting an associated Yamabe-type problem.


title: Landau theory and giant room-temperature barocaloric effect in MF$_3$ metal trifluorides

abstract: The structural phase transitions of MF$_3$ (M=Al, Cr, V, Fe, Ti, Sc) metal trifluorides are studied within a simple Landau theory consisting of tilts of rigid MF$_6$ octahedra associated with soft antiferrodistoritive optic modes that are coupled to long-wavelength strain generating acoustic phonons. We calculate the temperature and pressure dependence of several quantities such as the spontaneous distortions, volume expansion and shear strains as well as $T-P$ phase diagrams. By contrasting our model to experiments we quantify the deviations from mean-field behavior and found that the tilt fluctuations of the MF$_6$ octahedra increase with metal cation size. We apply our model to predict giant barocaloric effects in Sc substituted TiF$_3$ of up to about $15\,$JK$^{-1}$kg$^{-1}$ for modest hydrostatic compressions of $0.2\,$GPa. The effect extends over a wide temperature range of over $140\,$K (including room temperature) due to a large predicted rate $dT_c/dP = 723\,$K GPa$^{-1}$, which exceeds those of typical barocaloric materials. Our results suggest that open lattice frameworks such as the trifluorides are an attractive platform to search for giant barocaloric effects.


title: Using Graph Properties to Speed-up GPU-based Graph Traversal: A Model-driven Approach

abstract: While it is well-known and acknowledged that the performance of graph algorithms is heavily dependent on the input data, there has been surprisingly little research to quantify and predict the impact the graph structure has on performance. Parallel graph algorithms, running on many-core systems such as GPUs, are no exception: most research has focused on how to efficiently implement and tune different graph operations on a specific GPU. However, the performance impact of the input graph has only been taken into account indirectly as a result of the graphs used to benchmark the system. In this work, we present a case study investigating how to use the properties of the input graph to improve the performance of the breadth-first search (BFS) graph traversal. To do so, we first study the performance variation of 15 different BFS implementations across 248 graphs. Using this performance data, we show that significant speed-up can be achieved by combining the best implementation for each level of the traversal. To make use of this data-dependent optimization, we must correctly predict the relative performance of algorithms per graph level, and enable dynamic switching to the optimal algorithm for each level at runtime. We use the collected performance data to train a binary decision tree, to enable high-accuracy predictions and fast switching. We demonstrate empirically that our decision tree is both fast enough to allow dynamic switching between implementations, without noticeable overhead, and accurate enough in its prediction to enable significant BFS speedup. We conclude that our model-driven approach (1) enables BFS to outperform state of the art GPU algorithms, and (2) can be adapted for other BFS variants, other algorithms, or more specific datasets.


title: Critical flavour number of the Thirring model in three dimensions

abstract: The Thirring model is a four-fermion theory with a current-current interaction and $U(2N)$ chiral symmetry. It is closely related to three-dimensional QED and other models used to describe properties of graphene. In addition it serves as a toy model to study chiral symmetry breaking. In the limit of flavour number $N \to 1/2$ it is equivalent to the Gross-Neveu model, which shows a parity-breaking discrete phase transition. The model was already studied with different methods, including Dyson-Schwinger equations, functional renormalisation group methods and lattice simulations. Most studies agree that there is a phase transition from a symmetric phase to a spontaneously broken phase for a small number of fermion flavours, but no symmetry breaking for large $N$. But there is no consensus on the critical flavour number $N^\text{cr}$ above which there is no phase transition anymore and on further details of the critical behaviour. Values of $N$ found in the literature vary between $2$ and $7$. All earlier lattice studies were performed with staggered fermions. Thus it is questionable if in the continuum limit the lattice model recovers the internal symmetries of the continuum model. We present new results from lattice Monte Carlo simulations of the Thirring model with SLAC fermions which exactly implement all internal symmetries of the continuum model even at finite lattice spacing. If we reformulate the model in an irreducible representation of the Clifford algebra, we find, in contradiction to earlier results, that the behaviour for even and odd flavour numbers is very different: For even flavour numbers, chiral and parity symmetry are always unbroken. For odd flavour numbers parity symmetry is spontaneously broken below the critical flavour number $N_\text{ir}^\text{cr}=9$ while chiral symmetry is still unbroken.


title: Why we like the ECI+ algorithm

abstract: Recently a measure for Economic Complexity named ECI+ has been proposed by Albeaik et al. We like the ECI+ algorithm because it is mathematically identical to the Fitness algorithm, the measure for Economic Complexity we introduced in 2012. We demonstrate that the mathematical structure of ECI+ is strictly equivalent to that of Fitness (up to normalization and rescaling). We then show how the claims of Albeaik et al. about the ability of Fitness to describe the Economic Complexity of a country are incorrect. Finally, we hypothesize how the wrong results reported by these authors could have been obtained by not iterating the algorithm.