title: Anticipatory Radio Resource Management for Mobile Video Streaming with Linear Programming

abstract: In anticipatory networking, channel prediction is used to improve communication performance. This paper describes a new approach for allocating resources to video streaming traffic while accounting for quality of service. The proposed method is based on integrating a model of the user's local play-out buffer into the radio access network. The linearity of this model allows to formulate a Linear Programming problem that optimizes the trade-off between the allocated resources and the stalling time of the media stream. Our simulation results demonstrate the full power of anticipatory optimization in a simple, yet representative, scenario. Compared to instantaneous adaptation, our anticipatory solution shows impressive gains in spectral efficiency and stalling duration at feasible computation time while being robust against prediction errors.


title: On the parity of the multiplicative order of certain products of integers related to Gauss factorials

abstract: In this note, we prove that under some conditions, certain products of integers related to Gauss factorials are always quadratic residues.


title: Sparse spanning $k$-connected subgraphs in tournaments

abstract: In 2009, Bang-Jensen asked whether there exists a function $g(k)$ such that every strongly $k$-connected $n$-vertex tournament contains a strongly $k$-connected spanning subgraph with at most $kn + g(k)$ arcs. In this paper, we answer the question by showing that every strongly $k$-connected $n$-vertex tournament contains a strongly $k$-connected spanning subgraph with at most $kn + 750k^2\log(k+1)$ arcs.


title: Self-organisation in protoplanetary disks: global, non-stratified Hall-MHD simulations

abstract: Recent observations revealed organised structures in protoplanetary disks, such as axisymmetric rings or horseshoe concen- trations evocative of large-scale vortices. These structures are often interpreted as the result of planet-disc interactions. However, these disks are also known to be unstable to the magneto-rotational instability (MRI) which is believed to be one of the dominant angular momentum transport mechanism in these objects. It is therefore natural to ask if the MRI itself could produce these structures without invoking planets. The nonlinear evolution of the MRI is strongly affected by the low ionisation fraction in protoplanetary disks. The Hall effect in particular, which is dominant in dense and weakly ionised parts of these objects, has been shown to spontaneously drive self- organising flows in shearing box simulations. Here, we investigate the behaviour of global MRI-unstable disc models dominated by the Hall effect and characterise their dynamics. We perform 3D unstratified Hall-MHD simulations of keplerian disks for a broad range of Hall, ohmic and ambipolar Elsasser numbers. We confirm the transition from a turbulent to an organised state as the intensity of the Hall effect is increased. We observe the formation of zonal flows, their number depending on the available magnetic flux and on the intensity of the Hall effect. For intermediate Hall intensity, the flow self-organises into long-lived magnetised vortices. Neither the addition of a toroidal field nor ohmic or ambipolar diffusion drastically change this picture in the range of parameters we have explored. The ability of these structures to trap dust particles in this configuration is demonstrated. We conclude that Hall-MRI driven organisation is a plausible scenario which could explain some of the structures found in recent observations.


title: Fair Scheduling for Data Collection in Mobile Sensor Networks with Energy Harvesting

abstract: We consider the problem of data collection from a continental-scale network of energy harvesting sensors, applied to tracking mobile assets in rural environments. Our application constraints favour a highly asymmetric solution, with heavily duty-cycled sensor nodes communicating with powered base stations. We study a novel scheduling optimisation problem for energy harvesting mobile sensor network, that maximises the amount of collected data under the constraints of radio link quality and energy harvesting efficiency, while ensuring a fair data reception. We show that the problem is NP-complete and propose a heuristic algorithm to approximate the optimal scheduling solution in polynomial time. Moreover, our algorithm is flexible in handling progressive energy harvesting events, such as with solar panels, or opportunistic and bursty events, such as with Wireless Power Transfer. We use empirical link quality data, solar energy, and WPT efficiency to evaluate the proposed algorithm in extensive simulations and compare its performance to state-of-the-art. We show that our algorithm achieves high data reception rates, under different fairness and node lifetime constraints.


abstract: Tunable photonic elements at the surface of an optical fiber with piezoelectric core are proposed and analyzed theoretically. These elements are based on whispering gallery modes whose propagation along the fiber is fully controlled by nanoscale variation of the effective fiber radius, which can be tuned by means of a piezoelectric actuator embedded into the core. The developed theory allows one to express the introduced effective radius variation through the shape of the actuator and the voltage applied to it. In particular, the design of a miniature tunable optical delay line and a miniature tunable dispersion compensator is presented. The potential application of the suggested model to the design of a miniature optical buffer is discussed.


title: An Introduction to Mechanized Reasoning

abstract: Mechanized reasoning uses computers to verify proofs and to help discover new theorems. Computer scientists have applied mechanized reasoning to economic problems but -- to date -- this work has not yet been properly presented in economics journals. We introduce mechanized reasoning to economists in three ways. First, we introduce mechanized reasoning in general, describing both the techniques and their successful applications. Second, we explain how mechanized reasoning has been applied to economic problems, concentrating on the two domains that have attracted the most attention: social choice theory and auction theory. Finally, we present a detailed example of mechanized reasoning in practice by means of a proof of Vickrey's familiar theorem on second-price auctions.


title: Evaluation of the performance of two state-transfer Hamiltonians in the presence of static disorder

abstract: We analyse the performance of two quantum-state-transfer Hamiltonians in the presence of diagonal and off-diagonal disorder, and in terms of different measures. The first Hamiltonian pertains to a fully-engineered chain and the second to a chain with modified boundary couplings. The task is to find which Hamiltonian is the most robust to given levels of disorder and irrespective of the input state. In this respect, it is shown that the performance of the two protocols are approximately equivalent.


title: A Hybrid Transmission Grid Architecture Enabling Efficient Optimal Power Flow

abstract: The recent rise of electricity generation based on renewable energy sources increases the demand for transmission capacity. Capacity expansion via the upgrade of transmission line capacity, e.g., by conversion to a high-voltage direct current (HVDC) line, is an attractive option. In this paper, it is shown that if the upgrade to HVDC is applied systematically to selected transmission lines across the grid, a hybrid architecture is obtained that enables an efficient and globally optimal solution of the optimal power flow (OPF) problem. More precisely, for conventional meshed AC transmission grids the OPF problem is nonconvex and in general NP-hard, rendering it hard to solve. We prove that after the upgrade to the proposed hybrid architecture, the same mesh topology facilitates an exact convex relaxation of the OPF problem, enabling its globally optimal solution with efficient polynomial time algorithms. This OPF method is then employed in simulations, which demonstrate that the hybrid architecture can increase the effective transmission capacity and substantially reduce the generation costs, even compared to the AC grid with optimal transmission switching.


title: A Software Package for Chemically Inspired Graph Transformation

abstract: Chemical reaction networks can be automatically generated from graph grammar descriptions, where rewrite rules model reaction patterns. Because a molecule graph is connected and reactions in general involve multiple molecules, the rewriting must be performed on multisets of graphs. We present a general software package for this type of graph rewriting system, which can be used for modelling chemical systems. The package contains a C++ library with algorithms for working with transformation rules in the Double Pushout formalism, e.g., composition of rules and a domain specific language for programming graph language generation. A Python interface makes these features easily accessible. The package also has extensive procedures for automatically visualising not only graphs and rewrite rules, but also Double Pushout diagrams and graph languages in form of directed hypergraphs. The software is available as an open source package, and interactive examples can be found on the accompanying webpage.


title: Resonantly Enhanced Multiphoton Ionization under XUV FEL radiation: A case study of the role of harmonics

abstract: We provide a detailed quantitative study of the possible role of a small admixture of harmonics on resonant two-photon ionization. The motivation comes from the occasional presence of 2nd and 3rd harmonics in FEL radiation. We obtain the dependence of ionic yields on the intensity of the fundamental, the percentage of 2nd harmonic and the detuning of the fundamental from resonance. Having examined the cases of one and two intermediate resonances, we arrive at results of general validity and global behavior, showing that even a small amount of harmonic may seem deceptively innocuous.


title: Translational and rotational dynamics of a large buoyant sphere in turbulence

abstract: We report experimental measurements of the translational and rotational dynamics of a large buoyant sphere in isotropic turbulence. We introduce an efficient method to simultaneously determine the position and (absolute) orientation of a spherical body from visual observation. The method employs a minimization algorithm to obtain the orientation from the 2D projection of a specific pattern drawn onto the surface of the sphere. This has the advantages that it does not require a database of reference images, is easily scalable using parallel processing, and enables accurate absolute orientation reference. Analysis of the sphere's translational dynamics reveals clear differences between the streamwise and transverse directions. The translational auto-correlations and PDFs provide evidence for periodicity in the particle's dynamics even under turbulent conditions. The angular autocorrelations show weak periodicity. The angular accelerations exhibit wide tails, however without a directional dependence.


title: The Exploratory Research of the Effect Communication Model and Effect Improving Strategy of Interactive Advertising

abstract: Interactive advertising which characterized by interactivity has become the mainstream of advertising by gradually replacing traditional one-way advertising during the new media era. This paper has obeyed the research outline below, literature review, background description, assumption proposed, and empirical analysis. Therefore, this paper proposed the communication model of interactive advertising and drawn two related important conclusions, interactivity has brought positive effect to advertising communication, different type of consumers tend to use different interactive options in different ways. Furthermore, this paper also presented three related optimization strategies to improving the communication of interactive advertising, namely, 1.changing communication model from one-way to two-way, 2.renovating new communication process and effect-generated path, 3.renovating new strategy portfolio to improving the communication effect of interactive advertising.


title: The Block Pseudo-Marginal Sampler

abstract: The pseudo-marginal (PM) approach is increasingly used for Bayesian inference in statistical models, where the likelihood is intractable but can be estimated unbiasedly. %Examples include random effect models, state-space models and data subsampling in big-data settings. Deligiannidis et al. (2016) show how the PM approach can be made much more efficient by correlating the underlying Monte Carlo (MC) random numbers used to form the estimate of the likelihood at the current and proposed values of the unknown parameters. Their approach greatly speeds up the standard PM algorithm, as it requires a much smaller number of samples or particles to form the optimal likelihood estimate. Our paper presents an alternative implementation of the correlated PM approach, called the block PM, which divides the underlying random numbers into blocks so that the likelihood estimates for the proposed and current values of the parameters only differ by the random numbers in one block. We show that this implementation of the correlated PM can be much more efficient for some specific problems than the implementation in Deligiannidis et al. (2016); for example when the likelihood is estimated by subsampling or the likelihood is a product of terms each of which is given by an integral which can be estimated unbiasedly by randomised quasi-Monte Carlo. Our article provides methodology and guidelines for efficiently implementing the block PM. A second advantage of the the block PM is that it provides a direct way to control the correlation between the logarithms of the estimates of the likelihood at the current and proposed values of the parameters than the implementation in Deligiannidis et al. (2016). We obtain methods and guidelines for selecting the optimal number of samples based on idealized but realistic assumptions.


title: Permutation Representations of the Orbits of the Automorphism Group of a Finite Module over Discrete Valuation Ring

abstract: Consider a discrete valuation ring $R$ whose residue field is finite of cardinality at least $3$. For a finite torsion module, we consider transitive subsets $O$ under the action of the automorphism group of the module. We prove that the associated permutation representation on the complex vector space $C[O]$ is multiplicity free. This is achieved by obtaining a complete description of the transitive subsets of $O\times O$ under the diagonal action of the automorphism group.


title: Consensus versus persistence of disagreement in opinion formation: the role of zealots

abstract: We consider a general class of three--state models where individuals hold one of two opposite opinions, or are neutral, and exchange opinions in generic pairwise interactions. We show that when opinions spread in a population where a fraction of individuals (zealots) maintain unshakably their opinion, one of four qualitatively distinct kinds of collective dynamics arise, depending on the specific rules governing the social interactions. Unsurprisingly, when their density is high, zealots drive the whole population to consensus on their opinion. For low densities a rich phase diagram emerges: a finite population of dissenters can survive and be the only stationary state or may need a critical mass or dissenters to be sustained; the critical mass may vanish or not as the density is reduced; the transition to the high density regime can be smooth or abrupt, and shows interesting hysteretic effects. For each choice of the interaction rules we calculate the critical density of zealots above which diverse opinions cannot survive.


title: Extracting Arabic Relations from the Web

abstract: The goal of this research is to extract a large list or table from named entities and relations in a specific domain. A small set of a handful of instance relations is required as input from the user. The system exploits summaries from Google search engine as a source text. These instances are used to extract patterns. The output is a set of new entities and their relations. The results from four experiments show that precision and recall varies according to relation type. Precision ranges from 0.61 to 0.75 while recall ranges from 0.71 to 0.83. The best result is obtained for (player, club) relationship, 0.72 and 0.83 for precision and recall respectively.


title: Modal operators for meet-complemented lattices

abstract: We investigate some modal operators of necessity and possibility in the context of meet-complemented (not necessarily distributive) lattices. We proceed in stages. We compare our operators with others.


title: Challenges in Nuclear Structure Theory

abstract: The goal of nuclear structure theory is to build a comprehensive microscopic framework in which properties of nuclei and extended nuclear matter, and nuclear reactions and decays can all be consistently described. Due to novel theoretical concepts, breakthroughs in the experimentation with rare isotopes, increased exchange of ideas across different research areas, and the progress in computer technologies and numerical algorithms, nuclear theorists have been quite successful in solving various bits and pieces of the nuclear many-body puzzle and the prospects are exciting. This article contains a brief, personal perspective on the status of the field.


title: Effective Capacity in Broadcast Channels with Arbitrary Inputs

abstract: We consider a broadcast scenario where one transmitter communicates with two receivers under quality-of-service constraints. The transmitter initially employs superposition coding strategies with arbitrarily distributed signals and sends data to both receivers. Regarding the channel state conditions, the receivers perform successive interference cancellation to decode their own data. We express the effective capacity region that provides the maximum allowable sustainable data arrival rate region at the transmitter buffer or buffers. Given an average transmission power limit, we provide a two-step approach to obtain the optimal power allocation policies that maximize the effective capacity region. Then, we characterize the optimal decoding regions at the receivers in the space spanned by the channel fading power values. We finally substantiate our results with numerical presentations.


title: Photometric evolution of the 2016 outburst of recurrent Nova LMC 1968: the first three weeks

abstract: Optical (BVRI) photometry of the first three weeks of the 2016 outburst of the recurrent Nova LMC 1968 is presented and discussed. The 2016 I-band light-curve is an exact replica, even in the most minute details, of that for the 2010 eruption. The maximum is inferred to have occurred on 2016 Jan 21.2 at I=11.5 mag, corresponding to an absolute magnitude M(I)=-7.15. A ~1 day plateau is present in all bands about six days past optical maximum, simultaneous with the emergence of super-soft X-ray emission in Swift observations, signalling the widespread ionization of the ejecta. The nova entered a much longer plateau about 9 days past maximum, governed by the brightness of the white dwarf, now directly visible and still nuclearly burning on its surface. An outburst recurrence mean period of ~955 days (2.6 yrs) would fit both the OGLE inter-season gaps and the observed intervals between previous outburts.


title: Clifford theory for glider representations

abstract: Classical Clifford theory studies the decomposition of simple $G$-modules into simple $H$-modules for some normal subgroup $H \triangleleft G$. In this paper we deal with chains of normal subgroups $1 \triangleleft G_1 \triangleleft \cdots \triangleleft G_d =G$, which allow to consider fragments and in particular glider representations. These are given by a descending chain of vector spaces over some field $K$ and relate different representations of the groups appearing in the chain. Picking some normal subgroup $H \triangleleft G$ one obtains a normal subchain and one can construct an induced fragment structure. Moreover, a notion of irreducibility of fragments is introduced, which completes the list of ingredients to perform a Clifford theory.


title: A Bayesian non-parametric method for clustering high-dimensional binary data

abstract: In many real life problems, objects are described by large number of binary features. For instance, documents are characterized by presence or absence of certain keywords; cancer patients are characterized by presence or absence of certain mutations etc. In such cases, grouping together similar objects/profiles based on such high dimensional binary features is desirable, but challenging. Here, I present a Bayesian non parametric algorithm for clustering high dimensional binary data. It uses a Dirichlet Process (DP) mixture model and simulated annealing to not only cluster binary data, but also find optimal number of clusters in the data. The performance of the algorithm was evaluated and compared with other algorithms using simulated datasets. It outperformed all other clustering methods that were tested in the simulation studies. It was also used to cluster real datasets arising from document analysis, handwritten image analysis and cancer research. It successfully divided a set of documents based on their topics, hand written images based on different styles of writing digits and identified tissue and mutation specificity of chemotherapy treatments.


title: Weak (anti)localization in tubular semiconductor nanowires with spin-orbit coupling

abstract: We compute analytically the weak (anti)localization correction to the Drude conductivity for electrons in tubular semiconductor systems of zinc blende type. We include linear Rashba and Dresselhaus spin-orbit coupling (SOC) and compare wires of standard growth directions $\langle100\rangle$, $\langle111\rangle$, and $\langle110\rangle$. The motion on the quasi-two-dimensional surface is considered diffusive in both directions: transversal as well as along the cylinder axis. It is shown that Dresselhaus and Rashba SOC similarly affect the spin relaxation rates. For the $\langle110\rangle$ growth direction, the long-lived spin states are of helical nature. We detect a crossover from weak localization to weak anti-localization depending on spin-orbit coupling strength as well as dephasing and scattering rate. The theory is fitted to experimental data of an undoped $\langle111\rangle$ InAs nanowire device which exhibits a top-gate-controlled crossover from positive to negative magnetoconductivity. Thereby, we extract transport parameters where we quantify the distinct types of SOC individually.


title: Which spectral distortions does $\Lambda$CDM actually predict?

abstract: Ever refined cosmological measurements have established the $\Lambda$CDM concordance model, with the key cosmological parameters being determined to percent-level precision today. This allows us to make explicit predictions for the spectral distortions of the cosmic microwave background (CMB) created by various processes occurring in the early Universe. Here, we summarize all guaranteed CMB distortions and assess their total uncertainty within $\Lambda$CDM. We also compare simple methods for approximating them, highlighting some of the subtle aspects when it comes to interpreting future distortion measurements. Under simplified assumptions, we briefly study how well a PIXIE-like experiment may measure the main distortion parameters (i.e., $\mu$ and $y$). Next generation CMB spectrometers are expected to detect the distortion caused by reionization and structure formation at extremely high significance. They will also be able to constrain the small-scale power spectrum through the associated $\mu$-distortion, improving limits on running of the spectral index. Distortions from the recombination era, adiabatic cooling of matter relative to the CMB and dark matter annihilation require a higher sensitivity than PIXIE in its current design. The crucial next step is an improved modeling of foregrounds and instrumental aspects, as we briefly discuss here.


title: Transit times and mean ages for nonautonomous and autonomous compartmental systems

abstract: We develop a theory for transit times and mean ages for nonautonomous compartmental systems. Using the McKendrick-von F\"orster equation, we show that the mean ages of mass in a compartmental system satisfy a linear nonautonomous ordinary differential equation that is exponentially stable. We then define a nonautonomous version of transit time as the mean age of mass leaving the compartmental system at a particular time and show that our nonautonomous theory generalises the autonomous case. We apply these results to study a nine-dimensional nonautonomous compartmental system modeling the terrestrial carbon cycle, which is a modification of the Carnegie-Ames-Stanford approach (CASA) model, and we demonstrate that the nonautonomous versions of transit time and mean age differ significantly from the autonomous quantities when calculated for that model.


title: Analysis of the strong decay $X(5568) \to B_s^0\pi^+$ with QCD sum rules

abstract: In this article, we take the $X(5568)$ to be the scalar diquark-antidiquark type tetraquark state, study the hadronic coupling constant $g_{XB_s\pi}$ with the three-point QCD sum rules by carrying out the operator product expansion up to the vacuum condensates of dimension-6 and including both the connected and disconnected Feynman diagrams, then calculate the partial decay width of the strong decay $ X(5568) \to B_s^0 \pi^+$ and obtain the value $\Gamma_X=\left(20.5\pm8.1\right)\,\rm{MeV}$, which is consistent with the experimental data $\Gamma_X = \left(21.9 \pm 6.4 {}^{+5.0}_{-2.5}\right)\,\rm{MeV}$ from the D0 collaboration.


title: A microscopic continuum model for defect dynamics in metallic glasses

abstract: Motivated by results of the topological theory of glasses accounting for geometric frustration, we develop the simplest possible continuum mechanical model of defect dynamics in metallic glasses that accounts for topological, energetic, and kinetic ideas. The model is aimed towards the development of a microscopic understanding of the plasticity of such materials. We discuss the expected predictive capabilities of the model vis-a-vis some observed physical behaviors of metallic glasses.


title: Elementary equivalences and accessible functors

abstract: We introduce the notion of $\lambda$-equivalence and $\lambda$-embeddings of objects in suitable categories. This notion specializes to $L_{\infty\lambda}$-equivalence and $L_{\infty\lambda}$-elementary embedding for categories of structures in a language of arity less than $\lambda$, and interacts well with functors and $\lambda$-directed colimits. We recover and extend results of Feferman and Eklof on "local functors" without fixing a language in advance. This is convenient for formalizing Lefschetz's principle in algebraic geometry, which was one of the main applications of the work of Eklof.


title: Mixture Proportion Estimation via Kernel Embedding of Distributions

abstract: Mixture proportion estimation (MPE) is the problem of estimating the weight of a component distribution in a mixture, given samples from the mixture and component. This problem constitutes a key part in many "weakly supervised learning" problems like learning with positive and unlabelled samples, learning with label noise, anomaly detection and crowdsourcing. While there have been several methods proposed to solve this problem, to the best of our knowledge no efficient algorithm with a proven convergence rate towards the true proportion exists for this problem. We fill this gap by constructing a provably correct algorithm for MPE, and derive convergence rates under certain assumptions on the distribution. Our method is based on embedding distributions onto an RKHS, and implementing it only requires solving a simple convex quadratic programming problem a few times. We run our algorithm on several standard classification datasets, and demonstrate that it performs comparably to or better than other algorithms on most datasets.


title: Jumping Jupiter can explain Mercury's orbit

abstract: The orbit of Mercury has large values of eccentricity and inclination that cannot be easily explained if this planet formed on a circular and coplanar orbit. Here, we study the evolution of Mercury's orbit during the instability related to the migration of the giant planets in the framework of the jumping Jupiter model. We found that some instability models are able to produce the correct values of Mercury's eccentricity and inclination, provided that relativistic effects are included in the precession of Mercury's perihelion. The orbital excitation is driven by the fast change of the normal oscillation modes of the system corresponding to the perihelion precession of Jupiter (for the eccentricity), and the nodal regression of Uranus (for the inclination).


title: Particle identification performance of the prototype Aerogel RICH counter for the Belle II experiment

abstract: We have developed a new type of particle identification device, called an Aerogel Ring Imaging Cherenkov (ARICH) counter, for the Belle II experiment. It uses silica aerogel tiles as Cherenkov radiators. For detection of Cherenkov photons, Hybrid Avalanche Photo-Detectors (HAPDs) are used. The designed HAPD has a high sensitivity to single photons under a strong magnetic field. We have confirmed that the HAPD provides high efficiency for single-photon detection even after exposure to neutron and gamma-ray radiation that exceeds the levels expected in the 10-year Belle II operation. In order to confirm the basic performance of the ARICH counter system, we carried out a beam test at the DESY using a prototype of the ARICH counter with six HAPD modules. The results are in agreement with our expectations and confirm the suitability of the ARICH counter for the Belle II experiment. Based on the in-beam performance of the device, we expect that the identification efficiency at 3.5 GeV/c is 97.4% and 4.9% for pions and kaons, respectively. This paper summarizes the development of the HAPD for the ARICH and the evaluation of the performance of the prototype ARICH counter built with the final design components.


title: The Erdos-Posa Property for Directed Graphs

abstract: A classical result by Erdos and Posa states that there is a function $f: {\mathbb N} \rightarrow {\mathbb N}$ such that for every $k$, every graph $G$ contains $k$ pairwise vertex disjoint cycles or a set $T$ of at most $f(k)$ vertices such that $G-T$ is acyclic. The generalisation of this result to directed graphs is known as Younger's conjecture and was proved by Reed, Robertson, Seymour and Thomas in 1996. This so-called Erdos-Posa-property can naturally be generalised to arbitrary graphs and digraphs. Robertson and Seymour proved that a graph $H$ has the Erdos-Posa-property if, and only if, $H$ is planar. In this paper we study the corresponding problem for digraphs. We obtain a complete characterisation of the class of strongly connected digraphs which have the Erdos-Posa-property (both for topological and butterfly minors). We also generalise this result to classes of digraphs which are not strongly connected. In particular, we study the class of vertex-cyclic digraphs (digraphs without trivial strong components). For this natural class of digraphs we obtain a nearly complete characterisation of the digraphs within this class with the Erdos-Posa-property. In particular we give positive and algorithmic examples of digraphs with the Erdos-Posa-property by using directed tree decompositions in a novel way.


title: A giant quiescent solar filament observed with high-resolution spectroscopy

abstract: A giant, quiet-Sun filament was observed with the high-resolution Echelle spectrograph at the Vacuum Tower Telescope at Observatorio del Teide on 2011 November 15. A mosaic of spectra (10 maps of 100" X 182") was recorded simultaneously in the chromospheric absorption lines H-alpha and Na I D2. Physical parameters of the filament plasma were derived using Cloud Model (CM) inversions and line core fits. The spectra were complemented with full-disk filtergrams (He I 10830 A, H-alpha, and Ca II K) of the Chromspheric Telescope (ChroTel) and full-disk magnetograms of HMI. The filament had extremely large linear dimensions (817"), which corresponds to about 658 Mm along a great circle on the solar surface. A total amount of 175119 H-alpha contrast profiles were inverted using the CM approach. The inferred mean line-of-sight (LOS) velocity, Doppler width, and source function were similar to previous works of smaller quiescent filaments. However, the derived optical thickness was larger. LOS velocity trends inferred from the H-alpha line core fits were in accord, but smaller, than the ones obtained with CM inversions. Signatures of counter-streaming flows were detected in the filament.The largest conglomerates of brightenings in the line core of Na I D2 coincided well with small-scale magnetic fields as seen by HMI. Mixed magnetic polarities were detected close to the ends of barbs. The computation of photospheric horizontal flows based on HMI magnetograms revealed flow kernels with a size of 5-8 Mm and velocities of 0.30-0.45 km/s at the ends of the filament. The physical properties of extremely large filaments are similar to their smaller counterparts, except for the optical thickness which in our sample was found to be larger. We found that a part of the filament, which erupted the day before, is in the process of reestablishing its initial configuration.


title: Joint law of the hitting time, overshoot and undershoot for a L\'evy process

abstract: Let be $(X_t, t\geq 0)$ be a L\'evy process which is the sum of a Brownian motion with drift and a compound Poisson process. We consider the first passage time $\tau_x$ at a fixed level $x>0$ by $(X_t, t\geq 0)$ and $K_x:= X_{\tau_x}-x$ the overshoot and $L_x:= x-X_{\tau_x^-}$ the undershoot. We first study the regularity of the density of the first passage time. Secondly, we calculate the joint law of $(\tau_x, K_x, L_x).$


title: Phonon instability of insulating states in optical lattices

abstract: The influence of collective phonon excitations, due to intersite atomic interactions, on the stability of optical lattices is analyzed. These phonon excitations are shown to essentially reduce the ability of atoms to be localized. The states that seem to be insulating in the absence of the phonon excitations can become delocalized when the latter are present. The delocalization effect exists for both long-range as well as local atomic interactions.


title: Jacobi's Inversion Problem for Genus Two Hyperelliptic Integral

abstract: Hinted by the elliptic parameterization of the Ising model, the addition formula of the elliptic function forms to give the integrable SU(2) group relation in the previous paper. We then expect that the addition formula of the Abelian function with any genus will form to give some integrable Lie group structure. In this paper, we study Jacobi's inversion problem for hyperelliptic integral with genus two and we expect some SU(2) structure for the addition formula of the hyperelliptic theta function with genus two.


title: Coherent states on quaternion Slices and a measurable field of Hilbert spaces

abstract: A set of reproducing kernel Hilbert spaces are obtained on Hilbert spaces over quaternion slices with the aid of coherent states. It is proved that the so obtained set forms a measurable field of Hilbert spaces and their direct integral appears again as a reproducing kernel Hilbert space for a bigger Hilbert space over the whole quaternions. Hilbert spaces over quaternion slices are identified as representation spaces for a set of irreducible unitary group representations and their direct integral is shown to be a reducible representation for the Hilbert space over the whole quaternion field.


title: Sunspot group tilt angle measurements from historical observations

abstract: Sunspot positions from various historical sets of solar drawings are analysed with respect to the tilt angles of bipolar sunspot groups. Data by Scheiner, Hevelius, Staudacher, Zucconi, Schwabe, and Spoerer deliver a series of average tilt angles spanning a period of 270 years, additional to previously found values for 20th-century data obtained by other authors. We find that the average tilt angles before the Maunder minimum were not significantly different from the modern values. However, the average tilt angles of a period 50 years after the Maunder minimum, namely for cycles 0 and 1, were much lower and near zero. The normal tilt angles before the Maunder minimum suggest that it was not abnormally low tilt angles which drove the solar cycle into a grand minimum.


title: Plasma-driven ultrashort bunch diagnostic

abstract: Ultrashort electron bunches are crucial for an increasing number of applications, however, diagnosing their longitudinal phase space remains a challenge. We propose a new method that harnesses the strong electric fields present in a laser driven plasma wakefield. By transversely displacing driver laser and witness bunch, a streaking field is applied to the bunch. This field maps the time information to a transverse momentum change and, consequently, to a change of transverse position. We illustrate our method with simulations where we achieve a time resolution in the attosecond range.


title: Population synthesis of RR Lyrae stars in the original $Kepler$ and $K2$ Fields of View

abstract: It is interesting to ask what fraction of the total available RR Lyrae (RRL) sample that falls in the $Kepler$ and $K2$ Fields of View (FoV) is known or discovered. In order assess the completeness of our sample we compared the known RRL sample in the $Kepler$ and $K2$ fields with synthetic Galactic models. The Catalina Sky Survey RRL sample was used to calibrate our method. We found that a large number of faint RRL stars is missing from $Kepler$ and $K2$ fields.


title: Tapering the sky response for angular power spectrum estimation from low-frequency radio-interferometric data

abstract: It is important to correctly subtract point sources from radio-interferometric data in order to measure the power spectrum of diffuse radiation like the Galactic synchrotron or the Epoch of Reionization 21-cm signal. It is computationally very expensive and challenging to image a very large area and accurately subtract all the point sources from the image. The problem is particularly severe at the sidelobes and the outer parts of the main lobe where the antenna response is highly frequency dependent and the calibration also differs from that of the phase center. Here we show that it is possible to overcome this problem by tapering the sky response. Using simulated 150 MHz observations, we demonstrate that it is possible to suppress the contribution due to point sources from the outer parts by using the Tapered Gridded Estimator to measure the angular power spectrum C_l of the sky signal. We also show from the simulation that this method can self-consistently compute the noise bias and accurately subtract it to provide an unbiased estimation of C_l.


title: Variational Autoencoders for Semi-supervised Text Classification

abstract: Although semi-supervised variational autoencoder (SemiVAE) works in image classification task, it fails in text classification task if using vanilla LSTM as its decoder. From a perspective of reinforcement learning, it is verified that the decoder's capability to distinguish between different categorical labels is essential. Therefore, Semi-supervised Sequential Variational Autoencoder (SSVAE) is proposed, which increases the capability by feeding label into its decoder RNN at each time-step. Two specific decoder structures are investigated and both of them are verified to be effective. Besides, in order to reduce the computational complexity in training, a novel optimization method is proposed, which estimates the gradient of the unlabeled objective function by sampling, along with two variance reduction techniques. Experimental results on Large Movie Review Dataset (IMDB) and AG's News corpus show that the proposed approach significantly improves the classification accuracy compared with pure-supervised classifiers, and achieves competitive performance against previous advanced methods. State-of-the-art results can be obtained by integrating other pretraining-based methods.


title: A systematic study of the initial state in heavy ion collisions based on the quark participant assumption

abstract: We investigate the initial state geometric quantities of heavy ion collisions based on the quark participant assumption in the Glauber multiple scattering approach. A systematic comparison to the nucleon participant assumption has been presented and confronted with the charged multiplicity measurements in various collision systems. It is found that the quark participant based assumption can be important to understand the data in multiplicity production and the initial spatial eccentricity in small systems.


title: Nested sampling, statistical physics and the Potts model

abstract: We present a systematic study of the nested sampling algorithm based on the example of the Potts model. This model, which exhibits a first order phase transition for $q>4$, exemplifies a generic numerical challenge in statistical physics: The evaluation of the partition function and thermodynamic observables, which involve high dimensional sums of sharply structured multi-modal density functions. It poses a major challenge to most standard numerical techniques, such as Markov Chain Monte Carlo. In this paper we will demonstrate that nested sampling is particularly suited for such problems and it has a couple of advantages. For calculating the partition function of the Potts model with $N$ sites: a) one run stops after $O(N)$ moves, so it takes $O(N^{2})$ operations for the run, b) only a single run is required to compute the partition function along with the assignment of confidence intervals, c) the confidence intervals of the logarithmic partition function decrease with $1/\sqrt{N}$ and d) a single run allows to compute quantities for all temperatures while the autocorrelation time is very small, irrespective of temperature. Thermodynamic expectation values of observables, which are completely determined by the bond configuration in the representation of Fortuin and Kasteleyn, like the Helmholtz free energy, the internal energy as well as the entropy and heat capacity, can be calculated in the same single run needed for the partition function along with their confidence intervals. In contrast, thermodynamic expectation values of magnetic properties like the magnetization and the magnetic susceptibility require sampling the additional spin degree of freedom. Results and performance are studied in detail and compared with those obtained with multi-canonical sampling. Eventually the implications of the findings on a parallel implementation of nested sampling are outlined.


title: Quasiparticle scattering interference in iron pnictides: A probe of the origin of nematicity

abstract: In this paper, we investigate the quasiparticle scattering interference(QPI) in the nematic phase of iron pnictides, based on the magnetic and orbital scenarios of nematicity, respectively. In the spin density wave(SDW) state, the QPI pattern exhibits a dimer structure in the energy region of the SDW gap, with its orientation along the ferromagnetic direction of the SDW order. When the energy is increased to be near the Fermi level, it exhibits two sets of dimers along the same direction. The dimer structure of the QPI patterns persists with the decrease of the SDW correlation length in the magnetic driven nematic phase, although it tends to merge together for the scattering patterns with energies close to the Fermi level. While in the orbital scenario, the QPI patterns exhibit a dimer structure in a wide energy region. It undergoes a {\pi}/2 rotation with the increasing of energy, which is associated with the inequivalent energies of the two Dirac nodes induced by the orbital order. These distinct features may be used to probe or distinguish two kinds of scenarios of the nemeticity.


title: A New Method to Visualize Deep Neural Networks

abstract: We present a method for visualising the response of a deep neural network to a specific input. For image data for instance our method will highlight areas that provide evidence in favor of, and against choosing a certain class. The method overcomes several shortcomings of previous methods and provides great additional insight into the decision making process of convolutional networks, which is important both to improve models and to accelerate the adoption of such methods in e.g. medicine. In experiments on ImageNet data, we illustrate how the method works and can be applied in different ways to understand deep neural nets.


title: Ballistic conductivity of graphene channel with p-n junction on ferroelectric domain wall

abstract: We study the impact of the ferroelectric domain wall on the ballistic conductance of the single-layer graphene channel in the heterostructure graphene / physical gap / ferroelectric film using Wentzel-Kramers-Brillouin approximation. Both self-consistent numerical modeling of the electric field and space charge dynamics in the heterostructure and approximate analytical theory show that the domain wall contact with the surface creates p-n junction in graphene channel. We calculated that the carriers' concentration induced in graphene by uncompensated ferroelectric dipoles originated from the spontaneous polarization abrupt near the surface can reach the values of 1019 m-2 order, which is in two orders higher than it can be obtained for the gate doped graphene on non-ferroelectric substrates. Therefore we predict that graphene channel with the p-n junction caused by ferroelectric domain wall would be characterized by rather high ballistic conductivity.


title: Antiferromagnetic spin chain behavior and a transition to 3D magnetic order in Cu(D,L-alanine)2: Roles of H-bonds

abstract: We study the spin chain behavior, a transition to 3D magnetic order and the magnitudes of the exchange interactions for the metal-amino acid complex Cu(D,L-alanine)2.H2O, a model compound to investigate exchange couplings supported by chemical paths characteristic of biomolecules. Thermal and magnetic data were obtained as a function of temperature (T) and magnetic field (B0). The magnetic contribution to the specific heat, measured between 0.48 and 30 K, displays above 1.8 K a 1D spin-chain behavior that can be fitted with an intrachain antiferromagnetic (AFM) exchange coupling constant 2J0 = (-2.12 $\pm$ 0.08) cm$^{-1}$, between neighbor coppers at 4.49 {\AA} along chains connected by non-covalent and H-bonds. We also observe a narrow specific heat peak at 0.89 K indicating a phase transition to a 3D magnetically ordered phase. Magnetization curves at fixed T = 2, 4 and 7 K with B0 between 0 and 9 T, and at T between 2 and 300 K with several fixed values of B0 were globally fitted by an intrachain AFM exchange coupling constant 2J0 = (-2.27 $\pm$ 0.02) cm$^{-1}$ and g = 2.091 $\pm$ 0.005. Interchain interactions J1 between coppers in neighbor chains connected through long chemical paths with total length of 9.51 {\AA} are estimated within the range 0.1 < |2J1| < 0.4 cm$^{-1}$, covering the predictions of various approximations. We analyze the magnitudes of 2J0 and 2J1 in terms of the structure of the corresponding chemical paths. The main contribution in supporting the intrachain interaction is assigned to H-bonds while the interchain interactions are supported by paths containing H-bonds and carboxylate bridges, with the role of the H-bonds being predominant. We compare the obtained intrachain coupling with studies of compounds showing similar behavior and discuss the validity of the approximations allowing to calculate the interchain interactions.


title: Tensor Product Multiplicities via Upper Cluster Algebras

abstract: For each valued quiver $Q$ of Dynkin type, we construct a valued ice quiver $\Delta_Q^2$. Let $G$ be a simple connected Lie group with Dynkin diagram the underlying valued graph of $Q$. The upper cluster algebra of $\Delta_Q^2$ is graded by the triple dominant weights $(\mu,\nu,\lambda)$ of $G$. We prove that when $G$ is simply-laced, the dimension of each graded component counts the tensor multiplicity $c_{\mu,\nu}^\lambda$. We conjecture that this is also true if $G$ is not simply-laced, and sketch a possible approach. Using this construction, we improve Berenstein-Zelevinsky's model, or in some sense generalize Knutson-Tao's hive model in type $A$.


title: Decoherence by spontaneous emission: a single-atom analog of superradiance

abstract: We show that the stationary decoherence rate of an open quantum system can be decomposed as a sum of local and nonlocal contributions, respectively related to the strength of the coupling between system and environment, and to the quality of the information about the system leaking into the environment. Both terms arise naturally in the framework of the influence functional, with the nonlocal term arising from the coupling between the backward and forward histories describing the open quantum system time evolution. While the local contribution always yields a positive decoherence rate, the nonlocal one may lead to recoherence when only partial information about the system is obtained from the disturbed environment. We illustrate these concepts in the framework of interferometry with trapped atoms, by analyzing the effect of spontaneous emission on the coherence of the center-of-mass.


title: Verification of an analytic fit for the vortex core profile in superfluid Fermi gases

abstract: A characteristic property of superfluidity and -conductivity is the presence of quantized vortices in rotating systems. To study the BEC-BCS crossover the two most common methods are the Bogoliubov-De Gennes theory and the usage of an effective field theory. In order to simplify the calculations for one vortex, it is often assumed that the hyperbolic tangent yields a good approximation for the vortex structure. The combination of a variational vortex structure, together with cylindrical symmetry yields analytic (or numerically simple) expressions. The focus of this article is to investigate to what extent this analytic fit truly reflects the vortex structure throughout the BEC-BCS crossover at finite temperatures. The vortex structure will be determined using the effective field theory presented in [Eur. Phys. Journal B 88, 122 (2015)] and compared to the variational analytic solution. By doing this it is possible to see where these two structures agree, and where they differ. This comparison results in a range of applicability where the hyperbolic tangent will be a good fit for the vortex structure.


title: Robust multigrid for high-order discontinuous Galerkin methods: A fast Poisson solver suitable for high-aspect ratio Cartesian grids

abstract: We present a polynomial multigrid method for nodal interior penalty and local discontinuous Galerkin formulations of the Poisson equation on Cartesian grids. For smoothing we propose two classes of overlapping Schwarz methods. The first class comprises element-centered and the second face-centered methods. Within both classes we identify methods that achieve superior convergence rates, prove robust with respect to the mesh spacing and the polynomial order, at least up to ${P=32}$. Consequent structure exploitation yields a computational complexity of $O(PN)$, where $N$ is the number of unknowns. Further we demonstrate the suitability of the face-centered method for element aspect ratios up to 32.


title: A minimum-change version of the Chung-Feller theorem for Dyck paths

abstract: A Dyck path with $2k$ steps and $e$ flaws is a path in the integer lattice that starts at the origin and consists of $k$ many $\nearrow$-steps and $k$ many $\searrow$-steps that change the current coordinate by $(1,1)$ or $(1,-1)$, respectively, and that has exactly $e$ many $\searrow$-steps below the line $y=0$. Denoting by $D_{2k}^e$ the set of Dyck paths with $2k$ steps and $e$ flaws, the Chung-Feller theorem asserts that the sets $D_{2k}^0,D_{2k}^1,\ldots,D_{2k}^k$ all have the same cardinality $\frac{1}{k+1}\binom{2k}{k}=C_k$, the $k$-th Catalan number. The standard combinatorial proof of this classical result establishes a bijection $f'$ between $D_{2k}^e$ and $D_{2k}^{e+1}$ that swaps certain parts of the given Dyck path $x$, with the effect that $x$ and $f'(x)$ may differ in many positions. In this paper we strengthen the Chung-Feller theorem by presenting a simple bijection $f$ between $D_{2k}^e$ and $D_{2k}^{e+1}$ which has the additional feature that $x$ and $f(x)$ differ in only two positions (the least possible number). We also present an algorithm that allows to compute a sequence of applications of $f$ in constant time per generated Dyck path. As an application, we use our minimum-change bijection $f$ to construct cycle-factors in the odd graph $O_{2k+1}$ and the middle levels graph $M_{2k+1}$ --- two intensively studied families of vertex-transitive graphs --- that consist of $C_k$ many cycles of the same length.


title: Testing fine-grained parallelism for the ADMM on a factor-graph

abstract: There is an ongoing effort to develop tools that apply distributed computational resources to tackle large problems or reduce the time to solve them. In this context, the Alternating Direction Method of Multipliers (ADMM) arises as a method that can exploit distributed resources like the dual ascent method and has the robustness and improved convergence of the augmented Lagrangian method. Traditional approaches to accelerate the ADMM using multiple cores are problem-specific and often require multi-core programming. By contrast, we propose a problem-independent scheme of accelerating the ADMM that does not require the user to write any parallel code. We show that this scheme, an interpretation of the ADMM as a message-passing algorithm on a factor-graph, can automatically exploit fine-grained parallelism both in GPUs and shared-memory multi-core computers and achieves significant speedup in such diverse application domains as combinatorial optimization, machine learning, and optimal control. Specifically, we obtain 10-18x speedup using a GPU, and 5-9x using multiple CPU cores, over a serial, optimized C-version of the ADMM, which is similar to the typical speedup reported for existing GPU-accelerated libraries, including cuFFT (19x), cuBLAS (17x), and cuRAND (8x).


title: Large deviations for the two-dimensional stochastic Navier-Stokes equation with vanishing noise correlation

abstract: We are dealing with the validity of a large deviation principle for the two-dimensional Navier-Stokes equation, with periodic boundary conditions, perturbed by a Gaussian random forcing. We are here interested in the regime where both the strength of the noise and its correlation are vanishing, on a length scale $\e$ and $\d(\e)$, respectively, with $0<\e,\ \d(\e)<<1$. Depending on the relationship between $\e$ and $\d(\e)$ we will prove the validity of the large deviation principle in different functional spaces.


title: Index-antiguiding in narrow-ridge GaN-based laser diodes investigated by measurements of the current-dependent gain and index spectra and by self-consistent simulation

abstract: The threshold current density of narrow (1.5 {\mu}m) ridge-waveguide InGaN multi-quantum-well laser diodes, as well as the shape of their lateral far-field patterns, strongly depend on the etch depth of the ridge waveguide. Both effects can be attributed to strong index-antiguiding. A value of the antiguiding factor R = 10 is experimentally determined near threshold by measurements of the current-dependent gain and refractive index spectra. The device performances are simulated self-consistently solving the Schr\"odinger-Poisson equations and the equations for charge transport and waveguiding. Assuming a carrier-induced index change which matches the experimentally determined antiguiding factor, both the measured high threshold current and the shape of the far-field pattern of lasers with shallow ridges can be reproduced theoretically.


title: Chemical and physical characterization of collapsing low-mass prestellar dense cores

abstract: The first hydrostatic core, also called the first Larson core, is one of the first steps in low-mass star formation, as predicted by theory. With recent and future high performance telescopes, details of these first phases become accessible, and observations may confirm theory and even bring new challenges for theoreticians. In this context, we study from a theoretical point of view the chemical and physical evolution of the collapse of prestellar cores until the formation of the first Larson core, in order to better characterize this early phase in the star formation process. We couple a state-of-the-art hydrodynamical model with full gas-grain chemistry, using different assumptions on the magnetic field strength and orientation. We extract the different components of each collapsing core (i.e., the central core, the outflow, the disk, the pseudodisk, and the envelope) to highlight their specific physical and chemical characteristics. Each component often presents a specific physical history, as well as a specific chemical evolution. From some species, the components can clearly be differentiated. The different core models can also be chemically differentiated. Our simulation suggests some chemical species as tracers of the different components of a collapsing prestellar dense core, and as tracers of the magnetic field characteristics of the core. From this result, we pinpoint promising key chemical species to be observed.


title: Canonical energy and hairy AdS black holes

abstract: We propose the modified version of the canonical energy which was introduced originally by Hollands and Wald. Our construction depends only on the Euler-Lagrange expression of the system and thus is independent of the ambiguity in the Lagrangian. After some comments on our construction, we briefly mention on the relevance of our construction to the boundary information metric in the context of the AdS/CFT correspondence. We also study the stability of three-dimensional hairy extremal black holes by using our construction.


title: On the pole content of coupled channels chiral approaches used for the $\bar{K}N$ system

abstract: Several theoretical groups describe the antikaon-nucleon interaction at low energies within approaches based on the chiral SU(3) dynamics and including next-to-leading order contributions. We present a comparative analysis of the pertinent models and discuss in detail their pole contents. It is demonstrated that the approaches lead to very different predictions for the $K^{-}p$ amplitude extrapolated to subthreshold energies as well as for the $K^{-}n$ amplitude. The origin of the poles generated by the models is traced to the so-called zero coupling limit, in which the inter-channel couplings are switched off. This provides new insights into the pole contents of the various approaches. In particular, different concepts of forming the $\Lambda(1405)$ resonance are revealed and constraints related to the appearance of such poles in a given approach are discussed.


title: On the inconsistency of $\ell_1$-penalised sparse precision matrix estimation

abstract: Various $\ell_1$-penalised estimation methods such as graphical lasso and CLIME are widely used for sparse precision matrix estimation. Many of these methods have been shown to be consistent under various quantitative assumptions about the underlying true covariance matrix. Intuitively, these conditions are related to situations where the penalty term will dominate the optimisation. In this paper, we explore the consistency of $\ell_1$-based methods for a class of sparse latent variable -like models, which are strongly motivated by several types of applications. We show that all $\ell_1$-based methods fail dramatically for models with nearly linear dependencies between the variables. We also study the consistency on models derived from real gene expression data and note that the assumptions needed for consistency never hold even for modest sized gene networks and $\ell_1$-based methods also become unreliable in practice for larger networks.


title: Invariant manifolds of parabolic fixed points (I). Existence and dependence on parameters

abstract: In this paper we study the existence and regularity of stable manifolds associated to fixed points of parabolic type in the differentiable and analytic cases, using the parametrization method. The parametrization method relies on a suitable approximate solution of a functional equation. In the case of parabolic points, if the manifolds have dimension two or higher, in general this approximation cannot be obtained in the ring of polynomials but as a sum of homogeneous functions and it is given in~\cite{BFM2015b}. Assuming a sufficiently good approximation is found, here we provide an "a posteriori" result which gives a true invariant manifold close to the approximated one. In the differentiable case, in some cases, there is a loss of regularity. We also consider the case of parabolic periodic orbits of periodic vector fields and the dependence of the manifolds on parameters. Examples are provided. We apply our method to prove that in several situations, namely, related to the parabolic infinity in the elliptic spatial three body problem, these invariant manifolds exist and do have polynomial approximations.


title: On Burenkov's extension operator preserving Sobolev-Morrey spaces on Lipschitz domains

abstract: We prove that Burenkov's Extension Operator preserves Sobolev spaces built on general Morrey spaces, including classical Morrey spaces. The analysis concerns bounded and unbounded open sets with Lipschitz boundaries in the n-dimensional Euclidean space.


title: Invariant manifolds of parabolic fixed points (II). Approximations by sums of homogeneous functions

abstract: We study the computation of local approximations of invariant manifolds of parabolic fixed points and parabolic periodic orbits of periodic vector fields. If the dimension of these manifolds is two or greater, in general, it is not possible to obtain polynomial approximations. Here we develop an algorithm to obtain them as sums of homogeneous functions by solving suitable cohomological equations. We deal with both the differentiable and analytic cases. We also study the dependence on parameters. In the companion paper, Invariant manifolds of parabolic fixed points (I), these approximations are used to obtain the existence of true invariant manifolds close by. Examples are provided.


title: Effective Vinogradov's mean value theorem via efficient boxing

abstract: We combine Wooley's efficient congruencing method with earlier work of Vinogradov and Hua to get effective bounds on Vinogradov's mean value theorem.


title: Estimation of the lifetime of small helium bubbles near tungsten surfaces -- a methodological study

abstract: Under low energy and high flux/fluence irradiation of helium (He) atoms, the formation and bursting of He bubbles on tungsten (W) surfaces play important roles in the morphological evolution of component surfaces and impurity production in fusion reactors. Microscopically, the bursting of He bubbles is a stochastic process, and He bubbles have statistically average lifetimes. In the present paper, a molecular dynamics-based method was developed to extract, for the first time, the lifetime of He bubbles near tungsten surfaces. It was found that He bubble bursting can be treated as an activated event. Its occurrence frequency or, equivalently, the average lifetime of bubbles follows the Arrhenius equation. For a given bubble size, the activation energy exhibits a good linear dependence with the depth, and the pre-exponential factor obeys the Meyer-Nedle rule. These results are useful for establishing a model in multi-scale simulations of impurity production in the fusion plasma and of the morphological evolution of component surfaces.


title: Consistency of multi-time Dirac equations with general interaction potentials

abstract: In 1932, Dirac proposed a formulation in terms of multi-time wave functions as candidate for relativistic many-particle quantum mechanics. A well-known consistency condition that is necessary for existence of solutions strongly restricts the possible interaction types between the particles. It was conjectured by Petrat and Tumulka that interactions described by multiplication operators are generally excluded by this condition, and they gave a proof of this claim for potentials without spin-coupling. Under smoothness assumptions of possible solutions we show that there are potentials which are admissible, give an explicit example, however, show that none of them fulfills the physically desirable Poincar\'e invariance. We conclude that in this sense Dirac's multi-time formalism does not allow to model interaction by multiplication operators, and briefly point out several promising approaches to interacting models one can instead pursue.


title: A Characteristic Dynamic Mode Decomposition

abstract: Temporal or spatial structures are readily extracted from complex data by modal decompositions like Proper Orthogonal Decomposition (POD) or Dynamic Mode Decomposition (DMD). Subspaces of such decompositions serve as reduced order models and define either spatial structures in time or temporal structures in space. On the contrary, convecting phenomena pose a major problem to those decompositions. A structure traveling with a certain group velocity will be perceived as a plethora of modes in time or space respectively. This manifests itself for example in poorly decaying singular values when using a POD. The poor decay is counter-intuitive, since a single structure is expected to be represented by a few modes. The intuition proves to be correct and we show that in a properly chosen reference frame along the characteristics defined by the group velocity, a POD or DMD reduces moving structures to a few modes, as expected. Beyond serving as a reduced model, the resulting entity can be used to define a constant or minimally changing structure in turbulent flows. This can be interpreted as an empirical counterpart to exact coherent structures. We present the method and its application to a head vortex of a compressible starting jet.


title: Single gap superconductivity in beta-Bi2Pd

abstract: beta-Bi2Pd compound has been proposed as another example of a multi-gap superconductor [Y. Imai et al., J. Phys. Soc. Jap. 81, 113708 (2012)]. Here, we report on measurements of several important physical quantities capable to show a presence of multiple energy gaps on our superconducting single crystals of beta-Bi2Pd with the critical temperature Tc close to 5 K. The calorimetric study via a sensitive ac technique shows a sharp anomaly at the superconducting transition, however only a single energy gap is detected. Also other characteristics inferred from calorimetric measurements as the field dependence of the Sommerfeld coefficient and the temperature and angular dependence of the upper critical magnetic field point unequivocally to standard single s-wave gap superconductivity. The Hall-probe magnetometry provides the same result from the analysis of the temperature dependence of the lower critical field. A single-gapped BCS density of states is detected by the scanning tunneling spectroscopy measurements. Then, the bulk as well as the surface sensitive probes evidence a standard conventional superconductivity in this system where the topologically protected surface states have been recently detected by ARPES [M. Sakano et al., Nature Comm. 6, 8595 (2015)] .


title: Bohmian Mechanics, Collapse Models and the emergence of Classicality

abstract: We discuss the emergence of classical trajectories in Bohmian Mechanics (BM), when a macroscopic object interacts with an external environment. We show that in such a case the conditional wave function of the system follows a dynamics which, under reasonable assumptions, corresponds to that of the Ghirardi-Rimini-Weber (GRW) collapse model. As a consequence, Bohmian trajectories evolve classically. Our analysis also shows how the GRW (istantaneous) collapse process can be derived by an underlying continuous interaction of a quantum system with an external agent, thus throwing a light on how collapses can emerge from a deeper level theory.


title: Invariant measures for piecewise continuous maps

abstract: We say that $f:[0,1]\to [0,1]$ is a {\it piecewise continuous interval map} if there exists a partition $0=x_0<x_1<\cdots<x_{d}<x_{d+1}=1$ of $[0,1]$ such that $f\vert_{(x_{i-1},x_i)}$ is continuous and the lateral limits $w_0^+=\lim_{x\to 0^+} f(x)$, $w_{d+1}^-=\lim_{x\to 1^-} f(x)$, \mbox{$w_i^{-}=\lim_{x\to x_i^-} f(x)$} and $w_i^{+}=\lim_{x\to x_i^+} f(x)$ exist for each $i$. We prove that every piecewise continuous interval map without connections admits an invariant Borel probability measure. We also prove that every injective piecewise continuous interval map with no connections and no periodic orbits is topologically semi-conjugate to an interval exchange transformation.


title: A novel mass generation scheme for an Abelian vector field

abstract: A novel mass generation procedure for an Abelian vector field is proposed. This procedure is based on the construction of a class of gauge theories whose free limit describes a free massless vector field and a set of massless real scalar fields by means of the antifield-BRST deformation technique. The relationship between our results and those arising from the Higgs mechanism based on the spontaneous symmetry breaking of an Abelian gauge symmetry is emphasized. Some examples with one, two, and three scalars are given.


title: Nonlinear hydrodynamical evolution of eccentric Keplerian discs in two dimensions: validation of secular theory

abstract: We perform global two-dimensional hydrodynamical simulations of Keplerian discs with free eccentricity over thousands of orbital periods. Our aim is to determine the validity of secular theory in describing the evolution of eccentric discs, and to explore their nonlinear evolution for moderate eccentricities. Linear secular theory is found to correctly predict the structure and precession rates of discs with small eccentricities. However, discs with larger eccentricities (and eccentricity gradients) are observed to precess faster (retrograde relative to the orbital motion), at a rate that depends on their eccentricities (and eccentricity gradients). We derive analytically a nonlinear secular theory for eccentric gas discs, which explains this result as a modification of the pressure forces whenever eccentric orbits in a disc nearly intersect. This effect could be particularly important for highly eccentric discs produced in tidal disruption events, or for narrow gaseous rings; it might also play a role in causing some of the variability in superhump binary systems. In two dimensions, the eccentricity of a moderately eccentric disc is long-lived and persists throughout the duration of our simulations. Eccentric modes are however weakly damped by their interaction with non-axisymmetric spiral density waves (driven by the Papaloizou-Pringle instability, which occurs in our idealised setup with solid walls), as well as numerical diffusion.

categories: astro-ph.EP astro-ph.HE astro-ph.SR physics.flu-dyn


title: Irreducible Green Functions Method applied to nanoscopic systems

abstract: The equation of motion method (EOM) for Green functions is one of the tools used in the analysis of quantum dot system coupled with metallic and superconducting leads. We investigate modified EOM, based on differentiation of double-time temperature dependent Green functions both after primary time t and secondary time t'. Our EOM approach allows us to obtain the Abrikosov-Suhl resonance both in the particle-hole symmetric case but also in the asymmetric cases. We will apply the irreducible Green functions technique to analyses the EOM for dot system. This method give a workable decoupling scheme breaking the infinite set of Green function equations. We apply this technique for calculating the density of the states and the differential conductance of single-level quantum dot with Coulomb repulsion attached to one metallic and one superconducting leads (N-QD-SC). Our results are compared with the experimental data and previous calculations.


title: Comparing the properties of the X-shaped bulges of NGC 4710 and the Milky Way with MUSE

abstract: We used the new ESO VLT instrument MUSE to obtain spectral and imaging coverage of NGC 4710. The wide area and excellent sampling of the MUSE integral field spectrograph allows us to investigate the dynamical properties of the X-shaped bulge of NGC 4710 and compare it with the properties of the Milky Way's own X-shaped bulge. We measured the radial velocities, velocity dispersion, and stellar populations using a penalized pixel full spectral fitting technique adopting simple stellar populations models, on a 1' x 1' area centred on the bulge of NGC 4710. We have constructed the velocity maps of the bulge of NGC 4710 and we investigated the presence of vertical metallicity gradients. These properties were compared to those of the Milky Way bulge and as well as to a simulated galaxy with boxy/peanut bulge. We find the line-of-sight velocity maps and 1D rotation curves of the bulge of NGC 4710 to be remarkably similar to those of the Milky Way bulge. Some specific differences that were identified are in good agreement with the expectations from variations in the bar orientation angle. The bulge of NGC 4710 has a boxy-peanut morphology with a pronounced X-shape, showing no indication of any additional spheroidally distributed bulge population, in which we measure a vertical metallicity gradient of 0.35 dex/kpc. The general properties of NGC 4710 are very similar to those observed in the MW bulge. However, it has been suggested that the MW bulge has an additional component that is comprised of the oldest, most metal-poor stars and which is not part of the boxy-peanut bulge structure. Such a population is not observed in NGC 4710, but could be hidden in the integrated light we observed.


title: Ligand-Surface Interactions and Surface Oxidation of Colloidal PbSe Quantum Dots Revealed by Thin-film Positron Annihilation Methods

abstract: Positron Two Dimensional Angular Correlation of Annihilation Radiation (2D-ACAR) measurements reveal modifications of the electronic structure and composition at the surfaces of PbSe quantum dots (QDs), deposited as thin films, produced by various ligands containing either oxygen or nitrogen atoms. In particular, the 2D-ACAR measurements on thin films of colloidal PbSe QDs capped with oleic acid ligands yield an increased intensity in the electron momentum density (EMD) at high momenta compared to PbSe quantum dots capped with oleylamine. Moreover, the EMD of PbSe QDs is strongly affected by the small ethylediamine ligands, since these molecules lead to small distances between QDs and favor neck formation between near neighbor QDs, inducing electronic coupling between neighboring QDs. The high sensitivity to the presence of oxygen atoms at the surface can be also exploited to monitor the surface oxidation of PbSe QDs upon exposure to air. Our study clearly demonstrates that positron annihilation spectroscopy applied to thin films can probe surface transformations of colloidal semiconductor QDs embedded in functional layers.


title: Challenges in Design of Kitaev Materials: Magnetic Interactions from Competing Energy Scales

abstract: In this study, we reanalyze the magnetic interactions in the Kitaev spin liquid candidate materials Na$_2$IrO$_3$, $\alpha$-RuCl$_3$, and $\alpha$-Li$_2$IrO$_3$ using nonperturbative exact diagonalization methods. These methods are more appropriate given the relatively itinerant nature of the systems suggested in previous works. We treat all interactions up to third neighbours on equal footing. The computed terms reveal significant long range coupling, bond-anisotropy, and/or off-diagonal couplings which we argue naturally explain the observed ordered phases in these systems. Given these observations, the potential for realizing the spin-liquid state in real materials is analyzed, and synthetic challenges are defined and explained.


title: Observation of Majorana fermions with spin selective Andreev reflection in the vortex of topological superconductor

abstract: Majorana fermion (MF) whose antiparticle is itself has been predicted in condensed matter systems. Signatures of the MFs have been reported as zero energy modes in various systems. More definitive evidences are highly desired to verify the existence of the MF. Very recently, theory has predicted MFs to induce spin selective Andreev reflection (SSAR), a novel magnetic property which can be used to detect the MFs. Here we report the first observation of the SSAR from MFs inside vortices in Bi2Te3/NbSe2 hetero-structure, in which topological superconductivity was previously established. By using spin-polarized scanning tunneling microscopy/spectroscopy (STM/STS), we show that the zero-bias peak of the tunneling differential conductance at the vortex center is substantially higher when the tip polarization and the external magnetic field are parallel than anti-parallel to each other. Such strong spin dependence of the tunneling is absent away from the vortex center, or in a conventional superconductor. The observed spin dependent tunneling effect is a direct evidence for the SSAR from MFs, fully consistent with theoretical analyses. Our work provides definitive evidences of MFs and will stimulate the MFs research on their novel physical properties, hence a step towards their statistics and application in quantum computing.


title: From dots to stripes to sheets - Shape control of lead sulfide nanostructures

abstract: Controlling anisotropy in nanostructures is a challenging but rewarding task since confinement in one or more dimensions influences the physical and chemical properties of the items decisively. In particular, semiconducting nanostructures can be tailored to gain optimized properties to work as transistors or absorber material in solar cells. We demonstrate that the shape of colloidal lead sulfide nanostructures can be tuned from spheres to stripes to sheets by means of the precursor concentrations, the concentration of a chloroalkane co-ligand and the synthesis temperature. All final structures still possess at least one dimension in confinement. Electrical transport measurements complement the findings.


title: Maximum Efficiency of Low-Dissipation Heat Engines at Arbitrary Power

abstract: We investigate maximum efficiency at a given power for low-dissipation heat engines. Close to maximum power, the maximum gain in efficiency scales as a square root of relative loss in power and this scaling is universal for a broad class of systems. For the low-dissipation engines, we calculate the maximum gain in efficiency for an arbitrary fixed power. We show that the engines working close to maximum power can operate at considerably larger efficiency compared to the efficiency at maximum power. Furthermore, we introduce universal bounds on maximum efficiency at a given power for low-dissipation heat engines. These bounds represent direct generalization of the bounds on efficiency at maximum power obtained by Esposito et al. Phys. Rev. Lett. 105, 150603 (2010). We derive the bounds analytically in the regime close to maximum power and for small power values. For the intermediate regime we present strong numerical evidence for the validity of the bounds.


title: The Magnetic Classification of Solar Active Regions 1992 - 2015

abstract: The purpose of this letter is to address a blind-spot in our knowledge of solar active region statistics. To the best of our knowledge there are no published results showing the variation of the Mount Wilson magnetic classifications as a function of solar cycle based on modern observations. We show statistics for all active regions reported in the daily Solar Region Summary from 1992 January 1 to 2015 December 31. We find that the $\alpha$ and $\beta$ class active regions (including all sub-groups e.g. $\beta\gamma$, $\beta\delta$) make up fractions of approximately 20% and 80% of the sample respectively. This fraction is relatively constant during high levels of activity, however, an increase in the $\alpha$ fraction to about 35% and and a decrease in the $\beta$ fraction to about 65% can be seen near each solar minimum and is statistically significant at the 2-$\sigma$ level. Over 30% of all active regions observed during the years of solar maxima were appended with the classifications $\gamma$ and/or $\delta$, while these classifications account for only a fraction of a percent during the years near the solar minima. This variation in the active region types indicates that the formation of complex active regions may be due to the pileup of frequent emergence of magnetic flux during solar maximum, rather than the emergence of complex, monolithic flux structures.


title: Inability of the entropy vector method to certify nonclassicality in linelike causal structures

abstract: Bell's theorem shows that our intuitive understanding of causation must be overturned in light of quantum correlations. Nevertheless, quantum mechanics does not permit signalling and hence a notion of cause remains. Understanding this notion is not only important at a fundamental level, but also for technological applications such as key distribution and randomness expansion. It has recently been shown that a useful way to decide which classical causal structures could give rise to a given set of correlations is to use entropy vectors. These are vectors whose components are the entropies of all subsets of the observed variables in the causal structure. The entropy vector method employs causal relationships among the variables to restrict the set of possible entropy vectors. Here, we consider whether the same approach can lead to useful certificates of non-classicality within a given causal structure. Surprisingly, we find that for a family of causal structures that include the usual bipartite Bell structure they do not. For all members of this family, no function of the entropies of the observed variables gives such a certificate, in spite of the existence of nonclassical correlations. It is therefore necessary to look beyond entropy vectors to understand cause from a quantum perspective.


title: A generally-relativistic gauge classification of the Dirac fields

abstract: We consider generally-relativistic gauge transformations for the spinorial fields finding two mutually exclusive but together exhaustive classes in which fermions are placed adding supplementary information to the results obtained by Lounesto, and identifying quantities analogous to the momentum vector and the Pauli-Lubanski axial vector we discuss how our results are similar to those obtained by Wigner; by taking into account the system of Dirac field equations we will investigate the consequences for the dynamics: and in particular we shall address the problem of getting the non-relativistic approximation in a consistent way. We are going to comment on extensions.


title: Search for s channel single top quark production in pp collisions at sqrt(s) = 7 and 8 TeV

abstract: A search is presented for single top quark production in the s channel in proton-proton collisions with the CMS detector at the CERN LHC in decay modes of the top quark containing a muon or an electron in the final state. The signal is extracted through a maximum-likelihood fit to the distribution of a multivariate discriminant defined using boosted decision trees to separate the expected signal contribution from background processes. The analysis uses data collected at centre-of-mass energies of 7 and 8 TeV and corresponding to integrated luminosities of 5.1 and 19.7 inverse femtobarns, respectively. The measured cross sections of 7.1 +/- 8.1 pb (at 7 TeV) and 13.4 +/- 7.3 pb (at 8 TeV) result in a best fit value of 2.0 +/- 0.9 for the combined ratio of the measured and expected values. The signal significance is 2.5 standard deviations, and the upper limit on the rate relative to the standard model expectation is 4.7 at 95% confidence level.


title: Local Lipschitz Stability for Inverse Robin Problems in Some Elliptic and Parabolic Systems

abstract: In this work, we shall study the nonlinear inverse problems of recovering the Robin coefficients in elliptic and parabolic systems of second order, and establish their local Lipschitz stabilities. Some local Lipschitz stability was derived for an elliptic inverse Robin problem. We shall first restructure the arguments in \cite{chou04} for the local Lipschitz stability so that the stability follows from three basic conditions for the elliptic inverse Robin problem. The new arguments are then generalized to help establish a novel local Lipschitz stability for parabolic inverse Robin problems.


title: Structure of K\"ahler potential for D-term inflationary attractor models

abstract: Minimal chaotic models of D-term inflation predicts too large primordial tensor perturbations. Although it can be made consistent with observations utilizing higher order terms in the K\"ahler potential, expansion is not controlled in the absence of symmetries. We comprehensively study the conditions of K\"ahler potential for D-term plateau-type potentials and discuss its symmetry. They include the alpha-attractor model with a massive vector supermultiplet and its generalization leading to pole inflation of arbitrary order. We extend the models so that it can describe Coulomb phase, gauge anomaly is cancelled, and fields other than inflaton are stabilized during inflation. We also point out a generic issue for large-field D-term inflation that the masses of the non-inflaton fields tend to exceed the Planck scale.


title: Multidimensional factorization through helical mapping

abstract: This paper proposes a new perspective on the problem of multidimensional spectral factorization, through helical mapping: $d$-dimensional ($d$D) data arrays are vectorized, processed by $1$D cepstral analysis and then remapped onto the original space. Partial differential equations (PDEs) are the basic framework to describe the evolution of physical phenomena. We observe that the minimum phase helical solution asymptotically converges to the $d$D semi-causal solution, and allows to decouple the two solutions arising from PDEs describing physical systems. We prove this equivalence in the theoretical framework of cepstral analysis, and we also illustrate the validity of helical factorization through a $2$D wave propagation example and a $3$D application to helioseismology.


title: Modifications of the Lifshitz-Kosevich formula in two-dimensional Dirac systems

abstract: Starting from the Luttinger-Ward functional we derive an expression for the oscillatory part of the grand potential of a two dimensional Dirac system in a magnetic field. We perform the computation for the clean and the disordered system, and we study the effect of electron-electron interactions on the oscillations. Unlike in the two dimensional electron gas (2DEG), a finite temperature and impurity scattering also affects the oscillation frequency. Furthermore, we find that in graphene, compared to the 2DEG, additional interaction induced damping effects occur: to two-loop order electron-electron interactions do lead to an additional damping factor in the amplitude of the Lifshitz-Kosevich-formula.


title: Predicting Mercury's Precession using Simple Relativistic Newtonian Dynamics

abstract: We present a new simple relativistic model for planetary motion describing accurately the anomalous precession of the perihelion of Mercury and its origin. The model is based on transforming Newton's classical equation for planetary motion from absolute to real spacetime influenced by the gravitational potential and introducing the concept of influenced direction.


title: Resonant enhancement of Raman scattering in metamaterials with hybrid electromagnetic and plasmonic resonances

abstract: A tri-layer metamaterial perfect absorber of light, consisting of (Al/ZnS/Al) films with the top aluminium layer patterned as an array of circular disk nanoantennas, is investigated for resonantly enhancing Raman scattering from C-60 fullerene molecules deposited on the metamaterial. The metamaterial is designed to have resonant bands due to plasmonic and electromagnetic resonances at the Raman pump frequency (725 nm) as well as Stokes emission bands. The Raman scattering from C60 on the metamaterial with resonantly matched bands is measured to be enhanced by an order of magnitude more than from C60 on metamaterials with off-resonant absorption bands peaked at 1090 nm. The Raman pump is significantly enhanced due to the resonance with a propagating surface plasmon band, while the highly impedance matched electromagnetic resonance is expected to couple out the Raman emission efficiently. The nature and hybridization of the plasmonic and electromagnetic resonances to form compound resonances are investigated by numerical simulations.


title: On Resolvability of a Graph Associated to a Finite Vector Space

abstract: The metric dimension of non-component graph, associated to a finite vector space, is determined. It is proved that the exchange property holds for resolving sets of the graph, except a special case. Some results are also related to an intersection graph.


title: A Jamming-resilient Algorithm for Self-triggered Network Coordination

abstract: The issue of cyber-security has become ever more prevalent in the analysis and design of cyber-physical systems. In this paper, we investigate self-triggered consensus networks in the presence of communication failures caused by Denialof- Service (DoS) attacks. A general framework is considered in which the network links can fail independent of each other. By introducing a notion of Persistency-of-Communication (PoC), we provide an explicit characterization of DoS frequency and duration under which consensus can be preserved by suitably designing time-varying control and communication policies. An explicit characterization of the effects of DoS on the consensus time is also provided. The considered notion of PoC is compared with classic average connectivity conditions that are found in pure continuous-time consensus networks. Finally, examples are given to substantiate the analysis.


title: Resilient Control under Denial-of-Service:Robust Design

abstract: In this paper, we study networked control systems in the presence of Denial-of-Service (DoS) attacks, namely attacks that prevent transmissions over the communication network. The control objective is to maximize frequency and duration of the DoS attacks under which closed-loop stability is not destroyed. Analog and digital predictor-based controllers with state resetting are proposed, which achieve the considered control objective for a general class of DoS signals. An example is given to illustrate the proposed solution approach.


title: Large Slowdown of Water Dynamics at Stacked Phospholipid Membranes for Increasing Hydration Level: All-Atoms Molecular Dynamics

abstract: Water hydrating phospholipid membranes determine their stability and function, as well as their interaction with other molecules. In this article we study, using all-atom molecular dynamics simulations, the rotational and translational dynamical properties of water molecules confined in stacked phospholipid membranes at different levels of hydration, from poorly hydrated to a completely hydrated membrane. We find that both the translational and the reorientation dynamics of water are dramatically slowed down as the hydration is reduced. Analyzing in details the structure and dynamics of the hydrogen bond at the interface, we show that both those among water molecules and those between water and lipids slow down by decreasing the hydration, however the latter are always slower than the former. By increasing hydration, water saturates all the possible hydrogen bonds with the lipids and, by further increase of hydration, the hydrogen bonds among waters becomes the majority. However, the dynamics of the water-lipids hydrogen bonds becomes up to one order of magnitude slower than that of the water-water hydrogen bonds, inducing a large slowing down of the dynamics of the entire system even at large hydration level.


title: A coordinate ascent method for solving semidefinite relaxations of non-convex quadratic integer programs

abstract: We present a coordinate ascent method for a class of semidefinite programming problems that arise in non-convex quadratic integer optimization. These semidefinite programs are characterized by a small total number of active constraints and by low-rank constraint matrices. We exploit this special structure by solving the dual problem, using a barrier method in combination with a coordinate-wise exact line search. The main ingredient of our algorithm is the computationally cheap update at each iteration and an easy computation of the exact step size. Compared to interior point methods, our approach is much faster in obtaining strong dual bounds. Moreover, no explicit separation and reoptimization is necessary even if the set of primal constraints is large, since in our dual approach this is covered by implicitly considering all primal constraints when selecting the next coordinate.


title: Production of ultracold ground state LiRb molecules by photoassociation through a resonantly coupled state

abstract: We report on a resonantly coupled 2(1) - 4 (1) photoassociation resonance in LiRb. This pair of states displays interesting decay physics that hint at interference-like effects caused by two different decay paths. We observe decay to predominantly $X \ ^1\Sigma^+ \ v=43$, with significant numbers of molecules produced in $X \ ^1\Sigma^+ \ v=2$ - 12. This photoassociation resonance also produces $\sim$300 $X \ ^1\Sigma^+ \ v=0 \ J=0$ molecules/second. Finally, we identify a stimulated Raman adiabatic passage transfer pathway from $v=43$ to $v=0$ that has the potential to produce up to $2 \times 10^5$ LiRb molecules/second in the $X \ ^1\Sigma^+ \ v=0 \ J=0$ state.


title: The effects of turbulence generated in Big Bang nucleosynthesis

abstract: The continuity equation requires that energy released in nuclear fusion flows away from the point of interaction and is not immediately thermalised into the CMB. It is seen that consequently the bulk of thermo-nuclear energy released in BBN is converted by radiation pressure to kinetic energy of gas motions and generates instability in homogeneous initial conditions. High levels of turbulence follow, triggering the Jeans mechanism and leading to the formation of galaxies and large-scale structure in the order of 10 Myrs. Friedmann cosmology is perturbed, creating a fractal Voronoi matter distribution and providing energy for reionisation. In keeping with observations, general relativity predicts that there is no bound on the size of the structures produced. The true (overall) rate of expansion is substantially below the observed (local) value of Hubble's constant, and is consistent with published estimates from the CMB for void models. The CMB is predicted to show near flatness over the observed region. It is seen from observation of fluid flows, from established physics, and from numerical solution of the Navier-Stokes equations that bisymmetric spiral vortices form where gas flows meet. These vortices necessarily have flat rotation curves at large radii and create conditions for spiral galaxy formation. It is shown by N-body simulation that enduring Sa, Sb and Sc spiral galaxies result from appropriate initial conditions. The simulations and analysis of spiral potential show that orbital velocities are greater than would be the case in an axisymmetric potential.


title: Geodesics of McVittie Spacetime with a Phantom Cosmological Background

abstract: We investigate the geodesics of a Schwarzschild spacetime embedded in an isotropic expanding cosmological background (McVittie metric). We focus on bound particle geodesics in a background including matter and phantom dark energy with constant dark energy equation of state parameter $w<-1$ involving a future Big Rip singularity at a time $t_{\ast}$. Such geodesics have been previously studied in the Newtonian approximation and found to lead to dissociation of bound systems at a time $t_{rip}<t_{\ast}$ which for fixed background $w$, depends on a single dimensionless parameter $\bar{\omega}_0$ related to the angular momentum and depending on the mass and the size of the bound system. We extend this analysis to large massive bound systems where the Newtonian approximation is not appropriate and we compare the derived dissociation time with the corresponding time in the context of the Newtonian approximation. By identifying the time when the effective potential minimum disappears due to the repulsive force of dark energy we find that the dissociation time of bound systems occurs earlier than the prediction of the Newtonian approximation. However, the effect is negligible for all existing cosmological bound systems and it would become important only in hypothetical bound extremely massive ($10^{20} M_\odot$) and large ($100Mpc$) bound systems. We verify this result by explicit solution of the geodesic equations. This result is due to an interplay between the repulsive phantom dark energy effects and the existence of the well known innermost stable orbits of Schwarzschild spacetimes.