title: Second cluster integral from the spectrum of an infinite XXZ spin chain

abstract: First and second terms of the low-temperature cluster expansion for the free energy density of a magnetically polarized $XXZ$ spin chain is obtained within the propagator approach suggested by E. W. Montroll and J. C. Ward. All the calculations employ only one- and two-magnon infinite-chain spectrums. In the $XXX$-point the result reproduces the well known S. Katsura's formula obtained 50 years ago by finite-chain calculations.


title: EURETILE D7.3 - Dynamic DAL benchmark coding, measurements on MPI version of DPSNN-STDP (distributed plastic spiking neural net) and improvements to other DAL codes

abstract: The EURETILE project required the selection and coding of a set of dedicated benchmarks. The project is about the software and hardware architecture of future many-tile distributed fault-tolerant systems. We focus on dynamic workloads characterised by heavy numerical processing requirements. The ambition is to identify common techniques that could be applied to both the Embedded Systems and HPC domains. This document is the first public deliverable of Work Package 7: Challenging Tiled Applications.


title: A source of the C-metric with perfect translational inertial dragging

abstract: A new source of the C-metric is described using Israel's formalism. This source is a singular accelerated shell. By construction, perfect inertial dragging is realized inside the shell. The equation of state and energy conditions for the shell are discussed.


title: Violations of the second law of thermodynamics by a non completely positive dynamics

abstract: We consider a recently proposed model of driven open quantum microcircuit [F. Pellegrini et al., Phys. Rev. Lett. 107, 060401 (2011)] amenable to experimental investigations. We show that such an open quantum system provides a concrete physical instance where to prove that modeling its time-evolution with a dynamics lacking complete positivity conflicts with the second law of thermodynamics.


title: Gaussian Interference Channel with Intermittent Feedback

abstract: We investigate how to exploit intermittent feedback for interference management by studying the two-user Gaussian interference channel (IC). We approximately characterize (within a universal constant) the capacity region for the Gaussian IC with intermittent feedback. We exactly characterize the the capacity region of the linear deterministic version of the problem, which gives us insight into the Gaussian problem. We find that the characterization only depends on the forward channel parameters and the marginal probability distribution of each feedback link. The result shows that passive and unreliable feedback can be harnessed to provide multiplicative capacity gain in Gaussian interference channels. We find that when the feedback links are active with sufficiently large probabilities, the perfect feedback sum-capacity is achieved to within a constant gap. In contrast to other schemes developed for interference channel with feedback, our achievable scheme makes use of quantize-map-and-forward to relay the information obtained through feedback, performs forward decoding, and does not use structured codes. We also develop new outer bounds enabling us to obtain the (approximate) characterization of the capacity region.


title: Inertial effects in accelerating spacetimes

abstract: The motion of test particles along circular orbits in the vacuum $C$ metric is studied in the Frenet-Serret formalism. Special orbits and corresponding intrinsically defined geometrically relevant properties are selectively studied.


title: Strains and axial outflows in the field of a rotating black hole

abstract: We study the behaviour of an initially spherical bunch of accelerated particles emitted along trajectories parallel to the symmetry axis of a rotating black hole. We find that, under suitable conditions, curvature and inertial strains compete to model the shape of axial outflows of matter contributing to generate jet-like structures. This is of course a purely kinematical effect which does not account by itself for physical processes underlying the formation of jets. In our analysis a crucial role is played by a property of the electric and magnetic part of the Weyl tensor to be Lorentz-invariant boosting along the axis of symmetry in Kerr spacetime.


title: Massless field perturbations of the spinning C metric

abstract: A single master equation is given describing spin $s\le2$ test fields that are gauge- and tetrad-invariant perturbations of the {\it spinning C metric} spacetime representing a source with mass $M$, uniformly rotating with angular momentum per unit mass $a$ and uniformly accelerated with acceleration $A$. This equation can be separated into its radial and angular parts. The behavior of the radial functions near the horizons is studied and used to examine the influence of $A$ on the phenomenon of superradiance, while the angular equation leads to modified spin-weighted spheroidal harmonic solutions generalizing those of the Kerr spacetime. Finally the coupling between the spin of the perturbing field and the acceleration parameter $A$ is discussed.


title: Sound waves and modulational instabilities on continuous wave solutions in spinor Bose-Einstein condensates

abstract: We analyze sound waves (phonons, Bogoliubov excitations) propagating on continuous wave (cw) solutions of repulsive $F=1$ spinor Bose-Einstein condensates (BECs), such as $^{23}$Na (which is antiferromagnetic or polar) and $^{87}$Rb (which is ferromagnetic). Zeeman splitting by a uniform magnetic field is included. All cw solutions to ferromagnetic BECs with vanishing $M_F=0$ particle density and non-zero components in both $M_F=\pm 1$ fields are subject to modulational instability (MI). MI increases with increasing particle density. MI also increases with differences in the components' wavenumbers; this effect is larger at lower densities but becomes insignificant at higher particle densities. CW solutions to antiferromagnetic (polar) BECS with vanishing $M_F=0$ particle density and non-zero components in both $M_F=\pm 1$ fields do not suffer MI if the wavenumbers of the components are the same. If there is a wavenumber difference, MI initially increases with increasing particle density, then peaks before dropping to zero beyond a given particle density. The cw solutions with particles in both $M_F=\pm 1$ components and nonvanishing $M_F=0$ components do not have MI if the wavenumbers of the components are the same, but do exhibit MI when the wavenumbers are different. Direct numerical simulations of a cw with weak white noise confirm that weak noise grows fastest at wavenumbers with the largest MI, and shows some of the results beyond small amplitude perturbations. Phonon dispersion curves are computed numerically; we find analytic solutions for the phonon dispersion in a variety of limiting cases.


title: Spectral Slope Variation at Proton Scales from Fast to Slow Solar Wind

abstract: We investigated the behavior of the spectral slope of interplanetary magnetic field fluctuations at proton scales for selected high resolution time intervals from WIND and MESSENGER spacecraft at $1$ AU and $0.56$ AU, respectively. The analysis was performed within the profile of high speed streams, moving from fast to slow wind regions. The spectral slope showed a large variability between $-3.75$ and $-1.75$ and a robust tendency for this parameter to be steeper within the trailing edge where the speed is higher and to be flatter within the subsequent slower wind, following a gradual transition between these two states. The value of the spectral index seems to depend firmly on the power associated to the fluctuations within the inertial range, higher the power steeper the slope. Our result support previous analyses suggesting that there must be some response of the dissipation mechanism to the level of the energy transfer rate along the inertial range.


title: Charged massive particle at rest in the field of a Reissner-Nordstr\"om black hole. II. Analysis of the field lines and the electric Meissner effect

abstract: The properties of the electric field of a two-body system consisting of a Reissner-Nordstr\"om black hole and a charged massive particle at rest have recently been analyzed in the framework of first order perturbation theory following the standard approach of Regge, Wheeler and Zerilli. In the present paper we complete this analysis by numerically constructing and discussing the lines of force of the "effective" electric field of the sole particle with the subtraction of the dominant contribution of the black hole. We also give the total field due to the black hole and the particle. As the black hole becomes extreme an effect analogous to the Meissner effect arises for the electric field, with the "effective field" lines of the point charge being expelled by the outer horizon of the black hole. This effect existing at the level of test field approximation, i.e. by neglecting the backreaction on the background metric and electromagnetic field due to the particle's mass and charge, is here found also at the complete perturbative level. We point out analogies with similar considerations for magnetic fields by Bi{\v c}\'ak and Dvo{\v r}\'ak. We also explicitly show that the linearization of the recently obtained Belinski-Alekseev exact solution coincides with our solution in the Regge-Wheeler gauge. Our solution thus represents a "bridge" between the test field solution, which neglects all the feedback terms, and the exact two-body solution, which takes into account all the non-linearity of the interaction.


title: Dye's Theorem and Gleason's Theorem for AW*-algebras

abstract: We prove that any map between projection lattices of $AW^\ast$-algebras $A$ and $B$, where $A$ has no Type $I_2$ direct summand, that preserves orthocomplementation and suprema of arbitrary elements, is a restriction of a normal Jordan $\ast$-homomorphism between $A$ and $B$. This allows us to generalize Dye's Theorem from von Neumann algebras to $AW^\ast$-algebras. We show that Mackey-Gleason-Bunce-Wright Theorem can be extended to homogeneous $AW^\ast$-algebras of Type I. The interplay between Dye's Theorem and Gleason's Theorem is shown. As an application we prove that Jordan $\ast$-homomorphims are commutatively determined. Another corollary says that Jordan parts of $AW^\ast$-algebras can be reconstructed from posets of their abelian subalgebras.


title: Physical frames along circular orbits in stationary axisymmetric spacetimes

abstract: Three natural classes of orthonormal frames, namely Frenet-Serret, Fermi-Walker and parallel transported frames, exist along any timelike world line in spacetime. Their relationships are investigated for timelike circular orbits in stationary axisymmetric spacetimes, and illustrated for black hole spacetimes.


title: ls1 mardyn: The massively parallel molecular dynamics code for large systems

abstract: The molecular dynamics simulation code ls1 mardyn is presented. It is a highly scalable code, optimized for massively parallel execution on supercomputing architectures, and currently holds the world record for the largest molecular simulation with over four trillion particles. It enables the application of pair potentials to length and time scales which were previously out of scope for molecular dynamics simulation. With an efficient dynamic load balancing scheme, it delivers high scalability even for challenging heterogeneous configurations. Presently, multi-center rigid potential models based on Lennard-Jones sites, point charges and higher-order polarities are supported. Due to its modular design, ls1 mardyn can be extended to new physical models, methods, and algorithms, allowing future users to tailor it to suit their respective needs. Possible applications include scenarios with complex geometries, e.g. for fluids at interfaces, as well as non-equilibrium molecular dynamics simulation of heat and mass transfer.


title: Carbon Memory Assessment

abstract: The geometrical and performance scaling of silicon CMOS integrated circuit technology over the past 50 years has enabled many affordable new products for business and consumer applications. Recognizing that Flash is approaching its ultimate physical scaling limits within the next 10 years or so, the global electronics research community has begun an intense search for a new paradigm and technology for extending the functional scaling of memory technologies. Several promising nonvolatile memory concepts have emerged, based on different switching and retention mechanisms from those of Flash memory, e.g., STTRAM, RRAM, PCM and more recently, resistive memories based on carbon, which are the topic of this paper. This paper will introduce into the diverse field of carbon materials by recollecting some effects in carbon that can be used to produce a multiple time switchable, non-volatile unipolar resistive memory with potential high scalability down to atomic dimensions. Carbon-based memory is a non-volatile resistive memory, therefore, the same architectures, circuits, select transistor or diodes like in ReRAM or PCRAM can be considered as implementation. The big advantage of carbon memory might be the high temperature retention of 250 C, which makes it attractive for automotive and harsh conditions. This is a white paper for the ITRS meeting on emerging research devices (ERD) in Albuquerque, New Mexico, on August 25-26, 2014.


title: The Kerr-Schild ansatz revised

abstract: Kerr-Schild metrics have been introduced as a linear superposition of the flat spacetime metric and a squared null vector field, say $\boldsymbol{k}$, multiplied by some scalar function, say $H$. The basic assumption which led to Kerr solution was that $\boldsymbol{k}$ be both geodesic and shearfree. This condition is relaxed here and Kerr-Schild ansatz is revised by treating Kerr-Schild metrics as {\it exact linear perturbations} of Minkowski spacetime. The scalar function $H$ is taken as the perturbing function, so that Einstein's field equations are solved order by order in powers of $H$. It turns out that the congruence must be geodesic and shearfree as a consequence of third and second order equations, leading to an alternative derivation of Kerr solution.


title: The effective geometry of the $n=1$ uniformly rotating self-gravitating polytrope

abstract: The \lq\lq effective geometry" formalism is used to study the perturbations of a perfect barotropic Newtonian self-gravitating rotating and compressible fluid coupled with gravitational backreaction. The case of a uniformly rotating polytrope with index $n=1$ is investigated, due to its analytical tractability. Special attention is devoted to the geometrical properties of the underlying background acoustic metric, focusing in particular on null geodesics as well as on the analog light cone structure.


title: Measuring growth index in a universe with sterile neutrinos

abstract: Consistency tests for the general relativity (GR) can be performed by constraining the growth index $\gamma$ using the measurements of redshift-space distortions (RSD) in conjunction with other observations. In previous studies, deviations from the GR expected value of $\gamma\approx 0.55$ at the 2--3$\sigma$ level were found. In this work, we reconsider the measurement of $\gamma$ in a universe with sterile neutrinos. We constrain the sterile neutrino cosmological model using the RSD measurements combined with the cosmic microwave background data (Planck temperature data plus WMAP 9-yr polarization data), the baryon acoustic oscillation data, the Hubble constant direct measurement, the Planck Sunyaev-Zeldovich cluster counts data, and the galaxy shear data. We obtain the constraint result of the growth index, $\gamma=0.584^{+0.047}_{-0.048}$, well consistent with the GR expected value (the consistency is at the 0.6$\sigma$ level). For the parameters of sterile neutrino, we obtain $N_{\rm{eff}}=3.62^{+0.26}_{-0.42}$ and $m_{\nu,{\rm{sterile}}}^{\rm{eff}}=0.48^{+0.11}_{-0.14}$ eV. We also consider the BICEP2 data and perform an analysis on the model with tensor modes. Similar fit results are obtained, showing that once light sterile neutrino is considered in the universe, GR will become well consistent with the current observations.


title: Morphological and Kinematic Evolution of Three Interacting Coronal Mass Ejections of 2011 February 13-15

abstract: During 2011 February 13 to 15, three Earth-directed CMEs launched in successively were recorded as limb CMEs by coronagraphs (COR) of STEREO. These CMEs provided an opportunity to study their geometrical and kinematic evolution from multiple vantage points. In this paper, we examine the differences in geometrical evolution of slow and fast speed CMEs during their propagation in the heliosphere. We also study their interaction and collision using STEREO/SECCHI COR and Heliospheric Imager (HI) observations. We have found evidence of interaction and collision between the CMEs of February 15 and 14 in COR2 and HI1 FOV, respectively, while the CME of February 14 caught the CME of February 13 in HI2 FOV. By estimating the true mass of these CMEs and using their pre and post-collision dynamics, the momentum and energy exchange between them during collision phase are studied. We classify the nature of observed collision between CME of February 14 and 15 as inelastic, reaching close to elastic regime. Relating imaging observations with the in situ measurements, we find that the CMEs move adjacent to each other after their collision in the heliosphere and are recognized as distinct structures in in situ observations by WIND spacecraft at L1. Our results highlight the significance of HI observations in studying CME-CME collision for the purpose of improved space weather forecasting.


title: Constrained Dirac gluino mediation

abstract: We perform a comparison study of the Constrained Minimal Supersymmetric Standard Model and Constrained General Gauge Mediation with and without a heavy Dirac gluino. These extremely simple models have very few free parameters and exhibit the characteristic features of supersoftness and supersafeness. We determine the characteristic low energy spectra, the production cross sections of key processes at the Large Hadron Collider and the degree of fine tuning for a representative range of parameters for each model.


title: On a nonlinear model for tumor growth in a cellular medium

abstract: We investigate the dynamics of a nonlinear model for tumor growth within a cellular medium. In this setting the "tumor" is viewed as a multiphase flow consisting of cancerous cells in either proliferating phase or quiescent phase and a collection of cells accounting for the "waste" and/or dead cells in the presence of a nutrient. Here, the tumor is thought of as a growing continuum $\Omega$ with boundary $\partial \Omega$ both of which evolve in time. The key characteristic of the present model is that the total density of cancerous cells is allowed to vary, which is often the case within cellular media. We refer the reader to the articles \cite{Enault-2010}, \cite{LiLowengrub-2013} where compressible type tumor growth models are investigated. Global-in-time weak solutions are obtained using an approach based on penalization of the boundary behavior, diffusion, viscosity and pressure in the weak formulation, as well as convergence and compactness arguments in the spirit of Lions \cite{Lions-1998} (see also \cite{Feireisl-book, DT-MixedModel-2013}).


title: Josephson junction ratchet: effects of finite capacitances

abstract: We study transport in an asymmetric SQUID which is composed of a loop with three capacitively and resistively shunted Josephson junctions: two in series in one arm and the remaining one in the other arm. The loop is threaded by an external magnetic flux and the system is subjected to both a time-periodic and a constant current. We formulate the deterministic and, as well, the stochastic dynamics of the SQUID in terms of the Stewart-McCumber model and derive an equation for the phase difference across one arm, in which an effective periodic potential is of the ratchet type, i.e. its reflection symmetry is broken. In doing so, we extend and generalize earlier study by Zapata et al. [Phys. Rev. Lett. 77, 2292 (1996)] and analyze directed transport in wide parameter regimes: covering the over-damped to moderate damping regime up to its fully under-damped regime. As a result we detect the intriguing features of a negative (differential) conductance, repeated voltage reversals, noise induced voltage reversals and solely thermal noise-induced ratchet currents. We identify a set of parameters for which the ratchet effect is most pronounced and show how the direction of transport can be controlled by tailoring the external magnetic flux.


title: Electromagnetic field generation in the downstream of electrostatic shocks due to electron trapping

abstract: A new magnetic field generation mechanism in electrostatic shocks is found, which can produce fields with magnetic energy density as high as 0.01 of the kinetic energy density of the flows on time scales $ \tilde \, 10^4 \, {\omega}_{pe}^{-1}$. Electron trapping during the shock formation process creates a strong temperature anisotropy in the distribution function, giving rise to the pure Weibel instability. The generated magnetic field is well-confined to the downstream region of the electrostatic shock. The shock formation process is not modified and the features of the shock front responsible for ion acceleration, which are currently probed in laser-plasma laboratory experiments, are maintained. However, such a strong magnetic field determines the particle trajectories downstream and has the potential to modify the signatures of the collisionless shock.


title: Spatial low-discrepancy sequences, spherical cone discrepancy, and applications in financial modeling

abstract: In this paper we introduce a reproducing kernel Hilbert space defined on $\mathbb{R}^{d+1}$ as the tensor product of a reproducing kernel defined on the unit sphere $\mathbb{S}^{d}$ in $\mathbb{R}^{d+1}$ and a reproducing kernel defined on $[0,\infty)$. We extend Stolarsky's invariance principle to this case and prove upper and lower bounds for numerical integration in the corresponding reproducing kernel Hilbert space. The idea of separating the direction from the distance from the origin can also be applied to the construction of quadrature methods. An extension of the area-preserving Lambert transform is used to generate points on $\mathbb{S}^{d-1}$ via lifting Sobol' points in $[0,1)^{d}$ to the sphere. The $d$-th component of each Sobol' point, suitably transformed, provides the distance information so that the resulting point set is normally distributed in $\mathbb{R}^{d}$. Numerical tests provide evidence of the usefulness of constructing Quasi-Monte Carlo type methods for integration in such spaces. We also test this method on examples from financial applications (option pricing problems) and compare the results with traditional methods for numerical integration in $\mathbb{R}^{d}$.


title: Concentration polarization, surface currents, and bulk advection in a microchannel

abstract: We present a comprehensive analysis of salt transport and overlimiting currents in a microchannel during concentration polarization. We have carried out full numerical simulations of the coupled Poisson-Nernst-Planck-Stokes problem governing the transport and rationalized the behaviour of the system. A remarkable outcome of the investigations is the discovery of strong couplings between bulk advection and the surface current; without a surface current, bulk advection is strongly suppressed. The numerical simulations are supplemented by analytical models valid in the long channel limit as well as in the limit of negligible surface charge. By including the effects of diffusion and advection in the diffuse part of the electric double layers, we extend a recently published analytical model of overlimiting current due to surface conduction.


title: Clouds of strings in third-order Lovelock gravity

abstract: Lovelock theory is a natural extension of the Einstein theory of general relativity to higher dimensions in which the first and second orders correspond, respectively, to general relativity and Einstein-Gauss-Bonnet gravity. We present exact black hole solutions of $D\geq 4$-dimensional spacetime for first-, second-, and third-order Lovelock gravities in a string cloud background. Further, we compute the mass, temperature, and entropy of black hole solutions for the higher-dimensional general relativity and Einstein-Gauss-Bonnet theories and also perform thermodynamic stability of black holes. It turns out that the presence of the Gauss-Bonnet term and/or background string cloud completely changes the black hole thermodynamics. Interestingly, the entropy of a black hole is unaffected due to a background string cloud. We rediscover several known spherically symmetric black hole solutions in the appropriate limits.


title: Kriging based Surrogate Modeling for Fractional Order Control of Microgrids

abstract: This paper investigates the use of fractional order (FO) controllers for a microgrid. The microgrid employs various autonomous generation systems like wind turbine generator (WTG), solar photovoltaic (PV), diesel energy generator (DEG) and fuel-cells (FC). Other storage devices like the battery energy storage system (BESS) and the flywheel energy storage system (FESS) are also present in the power network. An FO control strategy is employed and the FO-PID controller parameters are tuned with a global optimization algorithm to meet system performance specifications. A kriging based surrogate modeling technique is employed to alleviate the issue of expensive objective function evaluation for the optimization based controller tuning. Numerical simulations are reported to prove the validity of the proposed methods. The results for both the FO and the integer order (IO) controllers are compared with standard evolutionary optimization techniques and the relative merits and demerits of the kriging based surrogate modeling are discussed. This kind of optimization technique is not only limited to this specific case of microgrid control but also can be ported to other computationally expensive power system optimization problems.


title: Bifurcations for a Coupled Schr\"odinger System with Multiple Components

abstract: In this paper, we study local bifurcations of an indefinite elliptic system with multiple components: \begin{equation*} \left\{\begin{array}{ll} -\Delta u_j + au_j = \mu_ju_j^3+\beta\sum_{k\ne j}u_k^2u_j, u_j>0\ \ \hbox{in}\ \Omega, u_j=0 \ \ \hbox{on}\ \partial\Omega,\ j=1,\dots,n. \end{array} \right. \end{equation*} Here $\Omega\subset{\mathbb{R}}^N$ is a smooth and bounded domain, $n\ge3$, $a<-\Lambda_1$ where $\Lambda_1$ is the principal eigenvalue of $(-\Delta, H_0^1(\Omega))$; $\mu_j$ and $\beta$ are real constants. Using the positive and non-degenerate solution of the scalar equation $-\Delta\omega-\omega=-\omega^3$, $\omega\in H_0^1(\Omega)$, we construct a synchronized solution branch $\mathcal{T}_\omega$. Then we find a sequence of local bifurcations with respect to $\mathcal{T}_\omega$, and we find global bifurcation branches of partially synchronized solutions.


title: A Generalized Damage-Based Estimate of Crash Severity

abstract: In this paper, principles of the CRASH3 damage-based collision algorithm are generalized for 3-dimensional motor vehicle crush profiles.


title: Photoproduction of doubly heavy baryon at the ILC

abstract: In the present paper, we make a detailed study on the doubly heavy baryon photoproduction in the future $e^+e^-$ International Linear Collider (ILC). The baryons $\Xi_{cc}$, $\Xi_{bc}$, and $\Xi_{bb}$ are produced via the channel $\gamma \gamma \to \Xi_{QQ'} +\bar{Q'} +\bar{Q}$, where $Q$ and $Q'$ stand for heavy $c$ or $b$ quark, respectively. As for the $\Xi_{QQ'}$-baryon production, it shall first generate a $(QQ')[n]$-diquark and then form the final baryon via fragmentation, where $[n]$ stands for the color- and spin- configurations for the $(QQ')$-diquark states. According to the non-relativistic QCD theory, four diquark configurations shall provide sizable contributions to the baryon production, e.g., $[n]$ equals $[^3S_1]_{\bar{\textbf{3}}}$, $[^1S_0]_{\textbf{6}}$, $[^3S_1]_{\textbf{6}}$, or $[^1S_0]_{\bar{\textbf{3}}}$, respectively. We adopt the improved helicity amplitude approach for the hard scattering amplitude to improve the calculation efficiency. Total and differential cross sections of those channels, as well as the theoretical uncertainties, are presented. We show that sizable amounts of baryon events can be generated at the ILC, i.e., about $2.0\times 10^{6}$ $\Xi_{cc}$, $2.2\times 10^{5}$ $\Xi_{bc}$, as well as $3.0\times 10^{3}\;\Xi_{bb}$ events are to be generated in one operation year for $\sqrt{S}=500$ GeV and ${\cal L}\simeq 10^{36}$cm$^{-2}$s$^{-1}$.


title: Exploring quantum phases by driven dissipation

abstract: Ever since the insight spreaded that tailored dissipation can be employed to control quantum systems and drive them towards pure states, the field of non-equilibrium quantum mechanics gained remarkable momentum. So far research focussed on emergent phenomena caused by the interplay and competition of unitary Hamiltonian and dissipative Markovian dynamics. In this manuscript we zero in on a so far rather understudied aspect of open quantum systems and non-equilibrium physics, namely the utilization of purely dissipative couplings to explore pure quantum phases and non-equilibrium phase transitions. To illustrate this concept, we introduce and scrutinize purely dissipative counterparts of (1) the paradigmatic transverse field Ising model and (2) the considerably more complex $\mathbb{Z}_2$ lattice gauge theory with coupled matter field. We show that, in mean field approximation, the non-equilibrium phase diagrams parallel the (thermal) phase diagrams of the Hamiltonian "blue print" theories qualitatively.


title: An all-optical nanomechanical heat engine

abstract: We propose and theoretically investigate a nanomechanical heat engine. We show how a levitated nanoparticle in a harmonic optical trap inside an optical cavity can be used to realize a Stirling cycle in the underdamped regime. The all-optical approach enables fast and exible control of all the thermodynamical parameters and the effcient optimization of the performance of the engine. We develop a systematic optimization procedure to determine optimal driving protocols. We further perform numerical simulations with realistic parameters and evaluate the maximum power and the corresponding effciency.


title: Diversification and Endogenous Financial Networks

abstract: We test the hypothesis that interconnections across financial institutions can be explained by a diversification motive. This idea stems from the empirical evidence of the existence of long-term exposures that cannot be explained by a liquidity motive (maturity or currency mismatch). We model endogenous interconnections of heterogenous financial institutions facing regulatory constraints using a maximization of their expected utility. Both theoretical and simulation-based results are compared to a stylized genuine financial network. The diversification motive appears to plausibly explain interconnections among key players. Using our model, the impact of regulation on interconnections between banks -currently discussed at the Basel Committee on Banking Supervision- is analyzed.


title: Invariant space under H\'enon renormalization : Intrinsic geometry of Cantor attractor

abstract: Three dimensional H\'non-like map $$ F(x,y,z) = (f(x) - \epsilon (x,y,z),\ x,\ \delta (x,y,z)) $$ is defined on the cubic box $ B $. An invariant space under renormalization would appear only in higher dimension. Consider renormalizable maps each of which satisfies the condition $$ \partial_y \delta \circ F(x,y,z) + \partial_z \delta \circ F(x,y,z) \cdot \partial_x \delta (x,y,z) \equiv 0 $$ for $ (x,y,z) \in B $. Denote the set of maps satisfying the above condition be $ \mathcal N $. Then the set $ \mathcal N \cap \mathcal I(\bar \epsilon) $ is invariant under the renormalization operator where $ \mathcal I(\bar \epsilon) $ is the set of infinitely renormalizable maps. H\'enon like diffeomorphism in $ \mathcal N \cap \mathcal I(\bar \epsilon) $ has universal numbers, $ b_2 \asymp | \partial_z \delta | $ and $ b_1 = b_F /b_2 $ where $ b_F $ is the average Jacobian of $ F $. The Cantor attractor of $ F \in \mathcal N \cap \mathcal I(\bar \epsilon) $, $ \mathcal O_F $ has {\em unbounded geometry} almost everywhere in the parameter space of $ b_1 $. If two maps in $ \mathcal N $ has different universal numbers $ b_1 $ and $ \widetilde b_1 $, then the homeomorphism between two Cantor attractor is at most H\"older continuous, which is called {\em non rigidity}.


title: Magnetization of fluid phonons and large-scale curvature perturbations

abstract: The quasinormal mode of a gravitating and magnetized fluid in a spatially flat, isotropic and homogeneous cosmological background is derived in the presence of the fluid sources of anisotropic stress and of the entropic fluctuations of the plasma. The obtained gauge-invariant description involves a system of two coupled differential equations whose physical content is analyzed in all the most relevant situations. The Cauchy problem of large-scale curvature perturbations during the radiation dominated stage of expansion can be neatly formulated and its general solution is shown to depend on five initial data assigned when the relevant physical wavelengths are larger than the particle horizon. The consequences of this approach are explored.


title: The detection of upwardly propagating waves channelling energy from the chromosphere to the low corona

abstract: There have been ubiquitous observations of wave-like motions in the solar atmosphere for decades. Recent improvements to space- and ground-based observatories have allowed the focus to shift to smaller magnetic structures on the solar surface. In this paper, high-resolution ground-based data taken using the Swedish 1 m Solar Telescope is combined with co-spatial and co-temporal data from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO) satellite to analyze running penumbral waves (RPWs). RPWs have always been thought to be radial wave propagation that occurs within sunspots. Recent research has suggested that they are in fact upwardly propagating field-aligned waves (UPWs). Here, RPWs within a solar pore are observed for the first time and are interpreted as UPWs due to the lack of a penumbra that is required to support RPWs. These UPWs are also observed co-spatially and co-temporally within several SDO/AIA elemental lines that sample the transition region and low corona. The observed UPWs are traveling at a horizontal velocity of around 17 +- 0.5 km s-1 and a minimum vertical velocity of 42 +- 21 km s-1. The estimated energy of the waves is around 150 W m-2, which is on the lower bound required to heat the quiet-Sun corona. This is a new, yet unconsidered source of wave energy within the solar chromosphere and low corona.


title: A new integral loss function for Bayesian optimization

abstract: We consider the problem of maximizing a real-valued continuous function $f$ using a Bayesian approach. Since the early work of Jonas Mockus and Antanas \v{Z}ilinskas in the 70's, the problem of optimization is usually formulated by considering the loss function $\max f - M_n$ (where $M_n$ denotes the best function value observed after $n$ evaluations of $f$). This loss function puts emphasis on the value of the maximum, at the expense of the location of the maximizer. In the special case of a one-step Bayes-optimal strategy, it leads to the classical Expected Improvement (EI) sampling criterion. This is a special case of a Stepwise Uncertainty Reduction (SUR) strategy, where the risk associated to a certain uncertainty measure (here, the expected loss) on the quantity of interest is minimized at each step of the algorithm. In this article, assuming that $f$ is defined over a measure space $(\mathbb{X}, \lambda)$, we propose to consider instead the integral loss function $\int_{\mathbb{X}} (f - M_n)_{+}\, d\lambda$, and we show that this leads, in the case of a Gaussian process prior, to a new numerically tractable sampling criterion that we call $\rm EI^2$ (for Expected Integrated Expected Improvement). A numerical experiment illustrates that a SUR strategy based on this new sampling criterion reduces the error on both the value and the location of the maximizer faster than the EI-based strategy.


title: Dust emission in star-forming dwarf galaxies: General properties and the nature of the sub-mm excess

abstract: We studied the global characteristics of dust emission in a large sample of emission-line star-forming galaxies. The sample consists of two subsamples. One subsample (SDSS sample) includes ~4000 compact star-forming galaxies from the SDSS, which were also detected in all four bands at 3.4, 4.6, 12, and 22 mum of the WISE all-sky survey. The second subsample (Herschel sample) is a sample of 28 compact star-forming galaxies observed with Herschel in the FIR range. Data of the Herschel sample were supplemented by the photometric data from the Spitzer observations, GALEX, SDSS, WISE, 2MASS, NVSS, and FIRST surveys, as well as optical and Spitzer spectra and data in sub-mm and radio ranges. It is found that warm dust luminosities of galaxies from the SDSS sample and cold and warm dust luminosities of galaxies from the Herschel sample are strongly correlated with Hbeta luminosities, which implies that one of the main sources of dust heating in star-forming galaxies is ionising UV radiation of young stars. Using the relation between warm and cold dust masses for estimating the total dust mass in star-forming galaxies with an accuracy better than ~0.5 dex is proposed. On the other hand, it is shown for both samples that dust temperatures do not depend on the metallicities. The dust-to-neutral gas mass ratio strongly declines with decreasing metallicity, similar to that found in other studies of local emission-line galaxies, high-redshift GRB hosts, and DLAs. On the other hand, the dust-to-ionised gas mass ratio is about one hundred times as high implying that most of dust is located in the neutral gas. It is found that thermal free-free emission of ionised gas in compact star-forming galaxies might be responsible for the sub-mm emission excess. This effect is stronger in galaxies with lower metallicities and is also positively affected by an increased star-formation rate.


title: Monte Carlo simulations of polymers with nearest- and next nearest-neighbor interactions on square and cubic lattices

abstract: We study a generalized interacting self-avoiding walk (ISAW) model with nearest- and next nearest-neighbor (NN and NNN) interactions on the square and cubic lattices. In both dimensions, the phase diagrams show coil and globule phases separated by continuous transition lines. Along these lines, we calculate the metric $\nu_t$, crossover $\phi_t$ and entropic $\gamma_t$ exponents, all of them in good agreement with the exact values of the $\Theta$ universality class. Therefore, the introduction of NNN interactions does not change the class of the ISAW model, which still exists even for repulsive forces. The growth parameters $\mu_t$ are shown to change monotonically with temperature along the $\Theta$-lines. In the square lattice, the $\Theta$-line has an almost linear behavior, which was not found in the cubic one. Although the region of repulsive NNN interactions, with attractive NN ones, leads to stiff polymers, no evidence of a transition to a crystalline phase was found.


title: Replicated Spectrographs in Astronomy

abstract: As telescope apertures increase, the challenge of scaling spectrographic astronomical instruments becomes acute. The next generation of extremely large telescopes (ELTs) strain the availability of glass blanks for optics and engineering to provide sufficient mechanical stability. While breaking the relationship between telescope diameter and instrument pupil size by adaptive optics is a clear path for small fields of view, survey instruments exploiting multiplex advantages will be pressed to find cost-effective solutions. In this review we argue that exploiting the full potential of ELTs will require the barrier of the cost and engineering difficulty of monolithic instruments to be broken by the use of large-scale replication of spectrographs. The first steps in this direction have already been taken with the soon to be commissioned MUSE and VIRUS instruments for the Very Large Telescope and the Hobby-Eberly Telescope, respectively. MUSE employs 24 spectrograph channels, while VIRUS has 150 channels. We compare the information gathering power of these replicated instruments with the present state of the art in more traditional spectrographs, and with instruments under development for ELTs. Design principles for replication are explored along with lessons learned, and we look forward to future technologies that could make massively-replicated instruments even more compelling.


title: Stress-Minimizing Orthogonal Layout of Data Flow Diagrams with Ports

abstract: We present a fundamentally different approach to orthogonal layout of data flow diagrams with ports. This is based on extending constrained stress majorization to cater for ports and flow layout. Because we are minimizing stress we are able to better display global structure, as measured by several criteria such as stress, edge-length variance, and aspect ratio. Compared to the layered approach, our layouts tend to exhibit symmetries, and eliminate inter-layer whitespace, making the diagrams more compact.


title: Self-organised criticality in dilute lattice

abstract: The variation of z in BTW model in presence of holes (dissipative sites) has been studied. The value of z decreases as the fraction of number of holes increases. Interstingly, it is observed that the variation of the rate of change of z with the fraction of number of holes is different for the two different types of distribution of holes over the lattice. When the holes are randomly distributed over the lattice then the dissipation is more compared to that of the case when the holes are present in the form of a single compact cluster with same fraction. The value of z is less in the first case than that observed in the second case.


title: Incremental Cardinality Constraints for MaxSAT

abstract: Maximum Satisfiability (MaxSAT) is an optimization variant of the Boolean Satisfiability (SAT) problem. In general, MaxSAT algorithms perform a succession of SAT solver calls to reach an optimum solution making extensive use of cardinality constraints. Many of these algorithms are non-incremental in nature, i.e. at each iteration the formula is rebuilt and no knowledge is reused from one iteration to another. In this paper, we exploit the knowledge acquired across iterations using novel schemes to use cardinality constraints in an incremental fashion. We integrate these schemes with several MaxSAT algorithms. Our experimental results show a significant performance boost for these algo- rithms as compared to their non-incremental counterparts. These results suggest that incremental cardinality constraints could be beneficial for other constraint solving domains.


title: Ring polymers in crowded environment: conformational properties

abstract: We analyze the universal size characteristics of flexible ring polymers in solutions in presence of structural obstacles (impurities) in d dimensions. One encounters such situations when considering polymers in gels, colloidal solutions, intra- and extracellular environments. A special case of extended impurities correlated on large distances r according to a power law \sim r^{-a} is considered. Applying the direct polymer renormalization scheme, we evaluate the estimates for averaged gyration radius $\langle R_{g\,{\rm ring}} \rangle$ and spanning radius $\langle R_{1/2\,{\rm ring}} \rangle$ of typical ring polymer conformation up to the first order of double \varepsilon=4-d, \delta=4-a expansion. Our results quantitatively reveal an extent of the effective size and anisotropy of closed ring macromolecules in disordered environment. In particular, the size ratio of ring and open (linear) polymers of the same molecular weight grows when increasing the strength of disorder according to $\langle R^2_{g\,{\rm ring}} \rangle / \langle R^2_{g\,{\rm chain}} \rangle =\frac{1}{2} \left(1+\frac{13}{48}\delta \right)$.


title: A non-commutative analogue of the Odlyzko bounds and bounds on performance for space-time lattice codes

abstract: This paper considers space-time coding over several independently Rayleigh faded blocks. In particular we will concentrate on giving upper bounds for the coding gain of lattice space-time codes as the number of blocks grow. This problem was previously considered in the single antenna case by Bayer et al. in 2006. Crucial to their work was Odlyzko's bound on the discriminant of an algebraic number field, as this provides an upper bound for the normalized coding gain of number field codes. In the MIMO context natural codes are constructed from division algebras defined over number fields and the coding gain is measured by the discriminant of the corresponding (non-commutative) algebra. In this paper we will develop analogues of the Odlyzko bounds in this context and show how these bounds limit the normalized coding gain of a very general family of division algebra based space-time codes. These bounds can also be used as benchmarks in practical code design and as tools to analyze asymptotic bounds of performance as the number of independently faded blocks increases.


title: Reducing quadratic forms by kneading sequences

abstract: We introduce an invertible operation on finite sequences of positive integers and call it "kneading". Kneading preserves three invariants of sequences -- the parity of the length, the sum of the entries, and one we call the "alternant". We provide a bijection between the set of sequences with alternant $a$ and parity $s$ and the set of Zagier-reduced indefinite binary quadratic forms with discriminant $a^2 + (-1)^s \cdot 4$, and show that kneading corresponds to Zagier reduction of the corresponding forms. It follows that the sum of a sequence is a class invariant of the corresponding form. We conclude with some observations and conjectures concerning this new invariant.


title: Membrane Instantons from Toda Field Theory

abstract: Four-dimensional quaternion-Kaehler metrics with an isometry are determined by solutions to the SU($\infty$) Toda equation. We set up a perturbation theory to construct approximate solutions to the latter which can be interpreted as membrane instanton corrections to the moduli space metric of the universal hypermultiplet. We compute one such solution exactly up to the five-instanton level, including all perturbative fluctuations about the instantons. The result shows a pattern that allows us to determine the asymptotic behaviour of all higher instanton corrections within this solution. We find that the generating function for the leading terms of the latter is given by the Lambert W-function.


title: Level-Rank Duality in Chern-Simons Theory from a Non-Supersymmetric Brane Configuration

abstract: We derive level-rank duality in pure Chern-Simons gauge theories from a non-supersymmetric Seiberg duality by using a non-supersymmetric brane configuration in type IIB string theory. The brane configuration consists of fivebranes, N D3 antibranes and an O3 plane. By swapping the fivebranes we derive a 3d non-supersymmetric Seiberg duality. After level shifts from loop effects, this identifies the IR of Sp(2N)_{2k-2N+2} and Sp(2k-2N+2)_{-2N} pure Chern-Simons theories, which is a level-rank pair. We also derive level-rank duality in a Chern-Simons theory based on a unitary group.


title: Further results on $B$-tensors with application to the location of real eigenvalues

abstract: In this paper, we give a further study on $B$-tensors and introduce doubly $B$-tensors that contain $B$-tensors. We show that they have similar properties, including their decompositions and strong relationship with strictly (doubly) diagonally dominated tensors. As an application, the properties of $B$-tensors are used to localize real eigenvalues of some tensors, which would be very useful in verifying the positive semi-definiteness of a tensor.


title: Quantum shells in a quantum space-time

abstract: We study the quantum motion of null shells in the quantum space-time of a black hole in loop quantum gravity. We treat the shells as test fields and use an effective dynamics for the propagation equations. The shells propagate through the region where the singularity was present in the classical black hole space-time, but is absent in the quantum space-time, eventually emerging through a white hole to a new asymptotic region of the quantum space-time. The profiles of the shells get distorted due to the quantum fluctuations in the Planckian region that replaces the singularity. The evolution of the shells is unitary throughout the whole process.


title: Do we always need a filter?

abstract: Since the groundbreaking work of the Kalman filter in the 1960s, considerable effort has been devoted to various discrete time filters for dynamic state estimation, especially including dozens of different types of suboptimal implementations of the Bayes filters. This has been accompanied by the rapid development of simulation/approximation theories and technologies. While admitting the success of filters in many cases, this study investigates the failure cases when they are in fact ineffective for state estimation. Several classic models have shown that the straightforward observation-only (O2) inference that does not need system modeling can perform better (in terms of both accuracy and computing speed) for estimation than filters. Special attention has been paid to quantitatively analyze when and why a filter will not outperform the O2 inference from the information fusion perspective. Thanks to the rapid development of advanced sensors, the O2 inference is not only engineering friendly and computationally fast but can also be very accurate and reliable by fusing the information received from multiple sensors. The statistical attributes of the multi-sensor O2 inference are analyzed and demonstrated through simulations. In the situation with limited sensors, the O2 approach can work jointly with existing clutter filtering and data association algorithms for multi-target tracking in clutter environments. Given an adequate number of sensors, the O2 approach can employ the multi-sensor data fusion to deal with clutter and can handle the very general multi-target tracking scenario with no background information.


title: Symmetric isostatic frameworks with $\ell^1$ or $\ell^\infty$ distance constraints

abstract: Combinatorial characterisations of minimal rigidity are obtained for symmetric 2-dimensional bar-joint frameworks with either $\ell^1$ or $\ell^\infty$ distance constraints. The characterisations are expressed in terms of symmetric tree packings and the number of edges fixed by the symmetry operations. The proof uses new Henneberg-type inductive construction schemes.


title: New algorithms for modular inversion and representation by binary quadratic forms arising from structure in the Euclidean algorithm

abstract: We observe structure in the sequences of quotients and remainders of the Euclidean algorithm with two families of inputs. Analyzing the remainders, we obtain new algorithms for computing modular inverses and representating prime numbers by the binary quadratic form $x^2 + 3xy + y^2$. The Euclidean algorithm is commenced with inputs from one of the families, and the first remainder less than a predetermined size produces the modular inverse or representation.


title: Integration of acoustical sensors into the KM3NeT Optical Modules

abstract: The next generation multi-cubic-kilometre water Cherenkov neutrino telescope will be build in the Mediterranean Sea. This telescope, called KM3NeT, is currently entering a first construction phase. The KM3NeT research infrastructure will comprise 690 so-called Detection Units in its final design which will be anchored to the sea bed and held upright by submerged floats. The positions of these Detection Units, several hundred metres in length, and their attached Optical Modules for the detection of Cherenkov light have to be monitored continously to provide the telescope with its desired pointing precision. A standard way to do this is the utilisation of an acoustic positioning system using emitters at fixed positions and receivers distributed along the Detection Units. The KM3NeT neutrino telescope comprises a custom-made acoustic positioning system with newly designed emitters attached to the anchors of the Detection Units and custom-designed receivers attached to the Detection Units. This article describes an approach for a receiver and its performance. The proposed Opto-Acoustical Modules combine the optical sensors for the telescope with the acoustical sensors necessary for the positioning of the module itself. This combination leads to a compact design suited for an easy deployment of the numerous Detection Units. Furthermore, the instrumented volume can be used for scientific analyses such as marine science and acoustic particle detection.


title: Bose-Einstein Condensation in an electro-pneumatically transformed quadrupole-Ioffe magnetic trap

abstract: We report a novel approach for preparing a Bose-Einstein condensate (BEC) of $^{87}$Rb atoms using electro-pneumatically driven transfer of atoms into a Quadrupole-Ioffe magnetic trap (QUIC Trap). More than 5$\times$$10^{8}$ atoms from a Magneto-optical trap are loaded into a spherical quadrupole trap and then these atoms are transferred into an Ioffe trap by moving the Ioffe coil towards the center of the quadrupole coil, thereby, changing the distance between quadrupole trap center and the Ioffe coil. The transfer efficiency is more than 80 \%. This approach is different from a conventional approach of loading the atoms into a QUIC trap wherein the spherical quadrupole trap is transformed into a QUIC trap by changing the currents in the quadrupole and the Ioffe coils. The phase space density is then increased by forced rf evaporative cooling to achieve the Bose-Einstein condensation having more than $10^{5}$ atoms.


title: The predual and John-Nirenberg inequalities on generalized BMO martingale spaces

abstract: In this paper we introduce the generalized BMO martingale spaces by stopping time sequences, which enable us to characterize the dual spaces of martingale Hardy-Lorentz spaces $H_{p,q}^s$ for $0<p\leq1, 1<q<\infty$. Moreover, by duality we obtain a John-Nirenberg theorem for the generalized BMO martingale spaces when the stochastic basis is regular. We also extend the boundedness of fractional integrals to martingale Hardy-Lorentz spaces.


title: Classification of reflection symmetry protected topological semimetals and nodal superconductors

abstract: While the topological classification of insulators, semimetals, and superconductors in terms of nonspatial symmetries is well understood, less is known about topological states protected by crystalline symmetries, such as mirror reflections and rotations. In this work, we systematically classify topological semimetals and nodal superconductors that are protected, not only by nonspatial (i.e., global) symmetries, but also by a crystal reflection symmetry. We find that the classification crucially depends on (i) the codimension of the Fermi surface (nodal line or point) of the semimetal (superconductor), (ii) whether the mirror symmetry commutes or anticommutes with the nonspatial symmetries and (iii) how the Fermi surfaces (nodal lines or points) transform under the mirror reflection and nonspatial symmetries. The classification is derived by examining all possible symmetry-allowed mass terms that can be added to the Bloch or Bogoliubov-de Gennes Hamiltonian in a given symmetry class and by explicitly deriving topological invariants. We discuss several examples of reflection symmetry protected topological semimetals and nodal superconductors, including topological crystalline semimetals with mirror $\mathbb{Z}_2$ numbers and topological crystalline nodal superconductors with mirror winding numbers.


title: Asymptotics and Concentration Bounds for Bilinear Forms of Spectral Projectors of Sample Covariance

abstract: Let $X,X_1,\dots, X_n$ be i.i.d. Gaussian random variables with zero mean and covariance operator $\Sigma={\mathbb E}(X\otimes X)$ taking values in a separable Hilbert space ${\mathbb H}.$ Let $$ {\bf r}(\Sigma):=\frac{{\rm tr}(\Sigma)}{\|\Sigma\|_{\infty}} $$ be the effective rank of $\Sigma,$ ${\rm tr}(\Sigma)$ being the trace of $\Sigma$ and $\|\Sigma\|_{\infty}$ being its operator norm. Let $$\hat \Sigma_n:=n^{-1}\sum_{j=1}^n (X_j\otimes X_j)$$ be the sample (empirical) covariance operator based on $(X_1,\dots, X_n).$ The paper deals with a problem of estimation of spectral projectors of the covariance operator $\Sigma$ by their empirical counterparts, the spectral projectors of $\hat \Sigma_n$ (empirical spectral projectors). The focus is on the problems where both the sample size $n$ and the effective rank ${\bf r}(\Sigma)$ are large. This framework includes and generalizes well known high-dimensional spiked covariance models. Given a spectral projector $P_r$ corresponding to an eigenvalue $\mu_r$ of covariance operator $\Sigma$ and its empirical counterpart $\hat P_r,$ we derive sharp concentration bounds for bilinear forms of empirical spectral projector $\hat P_r$ in terms of sample size $n$ and effective dimension ${\bf r}(\Sigma).$ Building upon these concentration bounds, we prove the asymptotic normality of bilinear forms of random operators $\hat P_r -{\mathbb E}\hat P_r$ under the assumptions that $n\to \infty$ and ${\bf r}(\Sigma)=o(n).$ In a special case of eigenvalues of multiplicity one, these results are rephrased as concentration bounds and asymptotic normality for linear forms of empirical eigenvectors. Other results include bounds on the bias ${\mathbb E}\hat P_r-P_r$ and a method of bias reduction as well as a discussion of possible applications to statistical inference in high-dimensional principal component analysis.


title: Developer Belief vs. Reality: The Case of the Commit Size Distribution

abstract: The design of software development tools follows from what the developers of such tools believe is true about software development. A key aspect of such beliefs is the size of code contributions (commits) to a software project. In this paper, we show that what tool developers think is true about the size of code contributions is different by more than an order of magnitude from reality. We present this reality, called the commit size distribution, for a large sample of open source and selected closed source projects. We suggest that these new empirical insights will help improve software development tools by aligning underlying design assumptions closer with reality.


title: Diagnosing hyperuniformity in two-dimensional disordered jammed-packings of soft spheres

abstract: Hyperuniformity characterizes a state of matter for which density fluctuations diminish towards zero at the largest length scales. However, the task of determining whether or not an experimental system is hyperuniform is experimentally challenging due to finite-resolution, noise and sample-size effects that influence characterization measurements. Here we explore these issues, employing video optical microscopy to study hyperuniformity phenomena in disordered two-dimensional jammed packings of soft spheres. Using a combination of experiment and simulation we characterize the detrimental effects of particle polydispersity, image noise, and finite-size effects on the assignment of hyperuniformity, and we develop a methodology that permits improved diagnosis of hyperuniformity from real-space measurements. The key to this improvement is a simple packing reconstruction algorithm that incorporates particle polydispersity to minimize free volume. In addition, simulations show that hyperuniformity can be ascertained more accurately in direct space than in reciprocal space as a result of finite sample-size. Finally, experimental colloidal packings of soft polymeric spheres are shown to be hyperuniform.


title: Exponential decay of eigenfunctions in a continuous multi-particle Anderson model with sub-exponentially decaying interaction

abstract: This short note is a complement to our recent paper [2] where we established strong dynamical localization for a class of multi-particle Anderson models in a Euclidean space with an alloy-type random potential and a sub-exponentially decaying interaction of infinite range. We show that the localized eigenfunctions at low energies actually decay exponentially fast. This improves the results by Fauser and Warzel who established sub-exponential decay of eigenfunctions in presence of a sub-exponentially decaying interaction.


title: Two electrons in the Dirac-bubble potential well

abstract: Within the Dirac-bubble potential model an electronic structure of the double-charged negative C60 ion has been studied by a variational method.


title: The redshift-space galaxy two-point correlation function and baryon acoustic oscillations

abstract: Future galaxy surveys will measure baryon acoustic oscillations (BAOs) with high significance, and a complete understanding of the anisotropies of BAOs in redshift space will be important to exploit the cosmological information in BAOs. Here we describe the anisotropies that arise in the redshift-space galaxy two-point correlation function (2PCF) and elucidate the origin of features that arise in the dependence of the BAOs on the angle between the orientation of the galaxy pair and the line of sight. We do so with a derivation of the configuration-space 2PCF using streaming model. We find that, contrary to common belief, the locations of BAO peaks in the redshift-space 2PCF are anisotropic even in the linear theory. Anisotropies in BAO depend strongly on the method of extracting the peak, showing maximum 3 % angular variation. We also find that extracting the BAO peak of $r^2\xi(r,\mu)$ significantly reduces the anisotropy to sub-percent level angular variation. When subtracting the tilt due to the broadband behavior of the 2PCF, the BAO bump is enhanced along the line of sight because of local infall velocities toward the BAO bump. Precise measurement of the angular dependence of the redshift-space 2PCF will allow new geometrical tests of dark energy beyond the BAO.


title: Ranking the spreading influence in complex networks

abstract: Identifying the node spreading influence in networks is an important task to optimally use the network structure and ensure the more efficient spreading in information. In this paper, by taking into account the shortest distance between a target node and the node set with the highest $k$-core value, we present an improved method to generate the ranking list to evaluate the node spreading influence. Comparing with the epidemic process results for four real networks and the Barab\'{a}si-Albert network, the parameterless method could identify the node spreading influence more accurately than the ones generated by the degree $k$, closeness centrality, $k$-shell and mixed degree decomposition methods. This work would be helpful for deeply understanding the node importance of a network.


title: First observation of flux avalanches in a-MoSi superconducting thin films

abstract: We have observed the occurrence of dendritic flux avalanches in an amorphous film of Mo$_{84}$Si$_{16}$. These events are understood to have a thermomagnetic origin and involve the abrupt penetration of bursts of magnetic flux taking place within a limited window of temperatures and magnetic fields. While dc-magnetometry allows one to determine the threshold fields for the occurrence of the thermomagnetic instabilities, magneto-optical imaging reveals the spatial distribution of magnetic flux throughout the sample. Conducting appropriate experiments, typical for this goal, avalanches were confirmed to be a characteristic of this material, ruling out the otherwise admissible possibility of an experimental artifact or a feature related to defects in the film. After the present observation, a-MoSi can be included in the gallery of superconducting materials exhibiting flux avalanches when in the form of thin films, a characteristic that must be carefully taken into consideration when one plans to employ films of those materials in applications.


title: A critical analysis of shock models for chondrule formation

abstract: In recent years many models of chondrule formation have been proposed. One of those models is the processing of dust in shock waves in protoplanetary disks. In this model, the dust and the chondrule precursors are overrun by shock waves, which heat them up by frictional heating and thermal exchange with the gas. In this paper we reanalyze the nebular shock model of chondrule formation and focus on the downstream boundary condition. We show that for large-scale plane-parallel chondrule-melting shocks the postshock equilibrium temperature is too high to avoid volatile loss. Even if we include radiative cooling in lateral directions out of the disk plane into our model (thereby breaking strict plane-parallel geometry) we find that for a realistic vertical extent of the solar nebula disk the temperature decline is not fast enough. On the other hand, if we assume that the shock is entirely optically thin so that particles can radiate freely, the cooling rates are too high to produce the observed chondrules textures. Global nebular shocks are therefore problematic as the primary sources of chondrules.


title: Ground state of Yang-Mills theory in 2+1 dimensions

abstract: Yang-Mills theory in 2+1 dimensions showed to be a research area yielding firm results in theoretical physics when compared to lattice computations. Recent analysis displayed astonishing agreement for the value of the string tension and excellent comparison for the spectrum. This successful approach can be put at test with a different theoretical framework that we devised in our preceding work for the scalar field theory in the strong coupling limit. The confirmations we get are really striking supporting it in full. As a by-product we are also able to show how AdS/CFT approach, with a description using flux tubes, is supported exactly as expected in the Isgur-Paton model with a dimensionless correction factor for the ground state of the theory already determined in lattice computations.


title: Computing convex hulls and counting integer points with polymake

abstract: The main purpose of this paper is to report on the state of the art of computing integer hulls and their facets as well as counting lattice points in convex polytopes. Using the polymake system we explore various algorithms and implementations. Our experience in this area is summarized in ten "rules of thumb".


title: On the Brezis-Lieb Lemma without pointwise convergence

abstract: Brezis-Lieb lemma is a refinement of Fatou lemma providing an evaluation of the gap between the integral for a sequence and the integral for its pointwise limit. This note studies the question if such gap can be evaluated when there is no a.e. convergence. In particular, it gives the same lower bound for the gap in L^p as the gap in the Brezis-Lieb lemma (including the case vector-valued functions) provided that p is greater or equal than 3 and the sequence converges both weakly and weakly in the sense of a duality map. It also shows that the statement is false if p<3. An application is given in form of a Brezis-Lieb lemma for gradients.


title: Evaluation of hyperelastic models for unidirectional short fibre reinforced materials using a representative volume element with refined boundary conditions

abstract: The simulation of a short fibre reinforced structure by means of the FEM requires the knowledge of the material behaviour at every Gauss point. In order to obtain such information, a representative volume element (RVE) containing unidirectional short fibres is analysed in the presented work. The findings are used to assess the applicability of several hyperelastic models describing transversal isotropic materials under consideration of large deformations. As the RVE's average response represents the homogenised behaviour at a macroscopic material point, the material models' parameters can be identified by fitting them to stress-strain curves obtained from simulations with the RVE. The application of periodic boundary conditions to the RVE in tensorial form enables a simple access to consider tilted fibres and catch the anisotropy induced by the fibres. The comparison of the calibrated material model with the RVE allows the assessment of the material model's applicability and quality. Both the modelling and the calculations are carried out with the commercial FEM software ABAQUS.


title: Non-linear power spectra in the synchronous gauge

abstract: We study the non-linear corrections to the matter and velocity power spectra in the synchronous gauge (SG). We consider the perturbations up to third order in a zero-pressure fluid in flat cosmological background, which is relevant for the non-linear growth of cosmic structure. As a result, we point out that the SG is an inappropriate coordinate choice when handling the non-linear growth of the large-scale structure. Although the equations in the SG happen to coincide with those in the comoving gauge (CG) to linear order, they differ from second order. In particular, the second order hydrodynamic equations in the the SG are apparently in the Lagrangian form, whereas those in the CG are in the Eulerian form. Thus, the non-linear power spectra naively presented in the original SG show strange behavior quite different from the result of the Newtonian theory even on sub-horizon scales. The power spectra in the SG show regularized behaviors only after we introduce convective terms in the second order so that the equations in two gauges coincide to the second order.


title: The LHC p+Pb run from the nuclear PDF perspective

abstract: The p+Pb and Pb+Pb runs at the LHC have opened a possibility to investigate the validity of collinear factorization in a clearly higher center-of-mass energy scale than earlier in nuclear collisions. Indeed, some processes that have been measured routinely in p+p($\overline {\rm p}$) collisions and utilized for years in free proton PDF fits, can now finally be reached also in the nuclear case. Such new data are expected to provide conclusive answers concerning the universality of the nuclear PDFs. In this talk, I will contrast some of the first p+Pb and Pb+Pb measurements to the predictions based on the nuclear PDFs.


title: Energy and Laplacian on Hanoi-type fractal quantum graphs

abstract: This article studies potential theory and spectral analysis on compact metric spaces, which we refer to as fractal quantum graphs. These spaces can be represented as a (possibly infinite) union of 1-dimensional intervals and a totally disconnected (possibly uncountable) compact set, which roughly speaking represents the set of junction points. Classical quantum graphs and fractal spaces such as the Hanoi attractor are included among them. We begin with proving the existence of a resistance form on the Hanoi attractor, and go on to establish heat kernel estimates and upper and lower bounds on the eigenvalue counting function of Laplacians corresponding to weakly self-similar measures on the Hanoi attractor. These estimates and bounds rely heavily on the relation between the length and volume scaling factors of the fractal. We then state and prove a necessary and sufficient condition for the existence of a resistance form on a general fractal quantum graph. Finally, we extend our spectral results to a large class of weakly self-similar fractal quantum graphs.

categories: math-ph math.FA math.MG math.MP math.PR math.SP


title: LHC data challenges the contemporary parton-to-hadron fragmentation functions

abstract: We discuss the inclusive high-pT charged-particle production in proton-proton collisions at the LHC. The experimental data are compared to the NLO perturbative QCD calculations employing various sets of parton-to-hadron fragmentation functions. Most of the theoretical predictions are found to disastrously overpredict the measured cross sections, even if the scale variations and PDF errors are accounted for. The problem appears to arise from the presently too hard gluon-to-hadron fragmentation functions.


title: Joint Hierarchical Gaussian Process Model with Application to Forecast in Medical Monitoring

abstract: A novel extrapolation method is proposed for longitudinal forecasting. A hierarchical Gaussian process model is used to combine nonlinear population change and individual memory of the past to make prediction. The prediction error is minimized through the hierarchical design. The method is further extended to joint modeling of continuous measurements and survival events. The baseline hazard, covariate and joint effects are conveniently modeled in this hierarchical structure. The estimation and inference are implemented in fully Bayesian framework using the objective and shrinkage priors. In simulation studies, this model shows robustness in latent estimation, correlation detection and high accuracy in forecasting. The model is illustrated with medical monitoring data from cystic fibrosis (CF) patients. Estimation and forecasts are obtained in the measurement of lung function and records of acute respiratory events. Keyword: Extrapolation, Joint Model, Longitudinal Model, Hierarchical Gaussian Process, Cystic Fibrosis, Medical Monitoring


title: Maser and Infrared Studies of Oxygen-Rich Late/Post-AGB Stars and Water Fountains: Development of a New Identification Method

abstract: We explored an efficient method to identify evolved stars with oxygen-rich envelopes in the late AGB or post-AGB phase of stellar evolution, which include a rare class of objects - the "water fountains". Our method considers the OH and H2O maser spectra, the near infrared Q-parameters (these are colour indices accounting for the effect of extinction), and far-infrared AKARI colours. Here we first present the results of a new survey on OH and H2O masers. There were 108 colour-selected objects: 53 of them were observed in the three OH maser lines (1612, 1665, and 1667 MHz), with 24 detections (16 new for 1612 MHz); and 106 of them were observed in the H2O maser line (22 GHz) with 24 detections (12 new). We identify a new potential water fountain source, IRAS19356+0754, with large velocity coverages of both OH and H2O maser emission. In addition, several objects with high velocity OH maser emission are reported for the first time. The Q-parameters as well as the infrared [09]-[18] and [18]-[65] AKARI colours of the surveyed objects are then calculated. We suggest that these infrared properties are effective in isolating aspherical from spherical objects, but the morphology may not necessarily be related to the evolutionary status. Nonetheless, by considering altogether the maser and infrared properties, the efficiency of identifying oxygen-rich late/post-AGB stars could be improved.


title: Light Stop Decays: Implications for LHC Searches

abstract: We investigate the flavour-changing neutral current decay of the lightest stop into a charm quark and the lightest neutralino and its four-body decay into the lightest neutralino, a down-type quark and a fermion pair. These are the relevant stop search channels in the low-mass region. The SUSY-QCD corrections to the two-body decay have been calculated for the first time and turn out to be sizeable. In the four-body decay both the contributions from diagrams with flavour-changing neutral current (FCNC) couplings and the mass effects of final state bottom quarks and $\tau$ leptons have been taken into account, which are not available in the literature so far. The resulting branching ratios are investigated in detail. We find that in either of the decay channels the branching ratios can deviate significantly from one in large parts of the allowed parameter range. Taking this into account, the experimental exclusion limits on the stop, which are based on the assumption of branching ratios equal to one, are considerably weakened. This should be taken into account in future searches for light stops at the next run of the LHC, where the probed low stop mass region will be extended.


title: Exploiting Multi-Core Architectures for Reduced-Variance Estimation with Intractable Likelihoods

abstract: Many popular statistical models for complex phenomena are intractable, in the sense that the likelihood function cannot easily be evaluated. Bayesian estimation in this setting remains challenging, with a lack of computational methodology to fully exploit modern processing capabilities. In this paper we introduce novel control variates for intractable likelihoods that can dramatically reduce the Monte Carlo variance of Bayesian estimators. We prove that our control variates are well-defined and provide a positive variance reduction. Furthermore we show how to optimise these control variates for variance reduction. The methodology is highly parallel and offers a route to exploit multi-core processing architectures that complements recent research in this direction. Indeed, our work shows that it may not be necessary to parallelise the sampling process itself in order to harness the potential of massively multi-core architectures. Simulation results presented on the Ising model, exponential random graph models and non-linear stochastic differential equation models support our theoretical findings.


title: On interpreting Patterson--Sullivan measures of geometrically finite groups as Hausdorff and packing measures

abstract: We provide a new proof of a theorem whose proof was sketched by Sullivan ('82), namely that if the Poincar\'e exponent of a geometrically finite Kleinian group $G$ is strictly between its minimal and maximal cusp ranks, then the Patterson--Sullivan measure of $G$ is not proportional to the Hausdorff or packing measure of any gauge function. This disproves a conjecture of Stratmann ('97, '06).


title: In-situ hydrogen charging of zirconium powder to study isothermal percipitation of hydrides and determination of Zr-hydride crystal structure

abstract: Zirconium alloys are widely used in the nuclear industry because of their high strength, good corrosion resistance and low neutron absorption cross-section. However, zirconium has strong affinity for hydrogen which leads to hydrogen concentration build-up over time. It is well known that the formation of hydrides will degrade the material and leads to, for example, delayed hydride cracking during high burn up. Even though zirconium hydrides have been studied for several decades, there still remain some controversies regarding the formation mechanisms, exact crystal structure, and stability of various hydride phases. This study uses high resolution synchrotron radiation as a probing tool to observe the precipitation and dissolution of hydrides in highly pure zirconium powder during in-situ hydrogen charging. The experiment enabled the direct observation of the hydride formation and phase transformations. It, also, provided high quality data for crystal structure determination.


title: Pattern Recognition in a Ring of Delayed Phase Oscillators

abstract: We show that a ring of phase oscillators coupled with transmission delays can be used as a pattern recognition system. The introduced model encodes patterns as stable periodic orbits. We present a detailed analysis of the underlying dynamics. In particular, we show that the system possesses a multitude of periodic solutions, prove stability results and present a bifurcation analysis. Furthermore, we show successful recognition results using artificial patterns and speech recordings.


title: Cavity solitons in vertical-cavity surface-emitting lasers

abstract: We investigate a control of the motion of localized structures of light by means of delay feedback in the transverse section of a broad area nonlinear optical system. The delayed feedback is found to induce a spontaneous motion of a solitary localized structure that is stationary and stable in the absence of feedback. We focus our analysis on an experimentally relevant system namely the Vertical-Cavity Surface-Emitting Laser (VCSEL). In the absence of the delay feedback we present experimental evidence of stationary localized structures in a 80 $\mu$m aperture VCSEL. The spontaneous formation of localized structures takes place above the lasing threshold and under optical injection. Then, we consider the effect of the time-delayed optical feedback and investigate analytically the role of the phase of the feedback and the carrier lifetime on the self-mobility properties of the localized structures. We show that these two parameters affect strongly the space time dynamics of two-dimensional localized structures. We derive an analytical formula for the threshold associated with drift instability of localized structures and a normal form equation describing the slow time evolution of the speed of the moving structure.


title: Interferometry from Space: A Great Dream

abstract: During some thirty years, 1980-2010, technical studies of optical interferometry from instruments in space were pursued as promising for higher spatial resolution and for higher astrometric accuracy. Nulling interferometry was studied for both high spatial resolution and high contrast. These studies were great dreams deserving further historical attention. ESA's interest in interferometry began in the early 1980s. The studies of optical interferometry for the global astrometry mission GAIA began in 1993 and ended in 1998 when interferometry was dropped as unsuited for the purpose, and the Gaia mission to be launched in 2013 is not based on interferometry. \c{opyright} Anita Publications. All rights reserved.


title: Time-distance helioseismology: A new averaging scheme for measuring flow vorticity

abstract: Time-distance helioseismology provides information about vector flows in the near-surface layers of the Sun by measuring wave travel times between points on the solar surface. Specific spatial averages of travel times have been proposed for distinguishing between flows in the east-west and north-south directions and measuring the horizontal divergence of the flows. No specific measurement technique has, however, been developed to measure flow vorticity. Here we propose a new measurement technique tailored to measuring the vertical component of vorticity. Fluid vorticity is a fundamental property of solar convection zone dynamics and of rotating turbulent convection in particular. The method consists of measuring the travel time of waves along a closed contour on the solar surface in order to approximate the circulation of the flow along this contour. Vertical vorticity is related to the difference between clockwise and counter-clockwise travel times. We applied the method to characterize the vortical motions of solar convection using helioseismic data from the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory (SDO/HMI) and from the Michelson Doppler Imager onboard the Solar and Heliospheric Observatory (SOHO/MDI). Away from the equator, a clear correlation between vertical vorticity and horizontal divergence is detected. Horizontal outflows are associated with negative vorticity in the northern hemisphere and positive vorticity in the southern hemisphere. The signal is much stronger for HMI than for MDI observations. We characterize the spatial power spectrum of the signal by comparison with a noise model. Vertical vorticity at horizontal wavenumbers below 250/R_Sun can be probed with this helioseismic technique.


title: Constraints on ultralight scalar dark matter from pulsar timing

abstract: We performed a Bayesian analysis of pulsar timing residuals from the NANOGrav pulsar timing array to search for a specific form of stochastic narrow-band signal produced by oscillating gravitational potential in the Galactic halo. Such oscillations arise in models of warm dark matter composed of an ultralight massive scalar field (m\sim 10^{-23} eV), recently considered by Khmelnitsky and Rubakov (2014). In the monochromatic approximation, the stringent upper limit (95% c.l.) on the variable gravitational potential amplitude is found to be Psi_c<1.14 10^{-15}, corresponding to the characteristic strain h_c=2\sqrt{3}\Psi_c < 4 10^{-15} at frequency f=1.75 10^{-8} Hz. In the narrow-band approximation, the upper limit of this background energy density is \Omega_{GPB}<1.27 10^{-9} at frequency f=6.2 10^{-9} Hz. These limits are an order of magnitude higher than the expected signal amplitude if the galactic dark matter consists of the ultralight scalar field. The applied analysis of pulsar timing residuals can be used to search for any narrow-band stochastic signals with different correlation properties. As a by-product, parameters of the red noise present in four NANOGrav pulsars (J1713+0747, J2145-0750, B1855+09, J1744-1134) have been evaluated.


title: Indications of Water Clouds in the Coldest Known Brown Dwarf

abstract: We present a deep near-infrared image of the newly discovered brown dwarf WISE J085510.83-071442.5 (W0855) using the FourStar imager at Las Campanas Observatory. Our detection of J3=24.8+0.33 -0.53 (J_MKO=25.0+0.33-0.53) at 2.6sigma -- or equivalently an upper limit of J3 > 23.8 (J_MKO > 24.0) at 5sigma makes W0855 the reddest brown dwarf ever categorized (J_MKO - W2 = 10.984+0.33 - 0.53 at 2.6sigma -- or equivalently an upper limit of J_MKO - W2 > 9.984 at 5sigma) and refines its position on color magnitude diagrams. Comparing the new photometry with chemical equilibrium model atmosphere predictions, we demonstrate that W0855 is 4.5sigma from models using a cloudless atmosphere and well reproduced by partly cloudy models (50%) containing sulfide and water ice clouds. Non-equilibrium chemistry or non-solar metallicity may change predictions, however using currently available model approaches, this is the first candidate outside our own solar system to have direct evidence for water clouds.


title: Detecting two-site spin-entanglement in many-body systems with local particle-number fluctuations

abstract: We derive experimentally measurable lower bounds for the two-site entanglement of the spin-degrees of freedom of many-body systems with local particle-number fluctuations. Our method aims at enabling the spatially resolved detection of spin-entanglement in Hubbard systems using high-resolution imaging in optical lattices. A possible application is the observation of entanglement generation and spreading during spin impurity dynamics, for which we provide numerical simulations. More generally, the scheme can simplify the entanglement detection in ion chains, Rydberg atoms, or similar atomic systems.


title: AFP Algorithm and a Canonical Normal Form for Horn Formulas

abstract: AFP Algorithm is a learning algorithm for Horn formulas. We show that it does not improve the complexity of AFP Algorithm, if after each negative counterexample more that just one refinements are performed. Moreover, a canonical normal form for Horn formulas is presented, and it is proved that the output formula of AFP Algorithm is in this normal form.


title: Measurement of the production and lepton charge asymmetry of $\textit{W}$ bosons in Pb+Pb collisions at $\sqrt{s_{\mathrm{\mathbf{NN}}}}=$ 2.76 TeV with the ATLAS detector

abstract: A measurement of $\textit{W}$ boson production in lead-lead collisions at $\sqrt{s_{\mathrm{NN}}}=$2.76 TeV is presented. It is based on the analysis of data collected with the ATLAS detector at the LHC in 2011 corresponding to an integrated luminosity of 0.14 $\mathrm{nb}^{-1}$ and 0.15 $\mathrm{nb}^{-1}$ in the muon and electron decay channels, respectively. The differential production yields and lepton charge asymmetry are each measured as a function of the average number of participating nucleons $< N_{\mathrm{part}} >$ and absolute pseudorapidity of the charged lepton. The results are compared to predictions based on next-to-leading-order QCD calculations. These measurements are, in principle, sensitive to possible nuclear modifications to the parton distribution functions and also provide information on scaling of $\textit{W}$ boson production in multi-nucleon systems.


title: Full-wave nonlinear ultrasound simulation on distributed clusters with applications in high-intensity focused ultrasound

abstract: Model-based treatment planning and exposimetry for high-intensity focused ultrasound (HIFU) requires the numerical simulation of nonlinear ultrasound propagation through heterogeneous and absorbing media. This is a computationally demanding problem due to the large distances travelled by the ultrasound waves relative to the wavelength of the highest frequency harmonic. Here, the k-space pseudospectral method is used to solve a set of coupled partial differential equations equivalent to a generalised Westervelt equation. The model is implemented in C++ and parallelised using the message passing interface (MPI) for solving large-scale problems on distributed clusters. The domain is partitioned using a 1D slab decomposition, and global communication is performed using a sparse communication pattern. Operations in the spatial frequency domain are performed in transposed space to reduce the communication burden imposed by the 3D fast Fourier transform. The performance of the model is evaluated using grid sizes up to 4096 x 2048 x 2048 grid points distributed over a cluster using up to 1024 compute cores. Given the global nature of the gradient calculation, the model shows good strong scaling behaviour, with a speed-up of 1.7x whenever the number of cores is doubled. This means large-scale simulations can be distributed across high numbers of cores on a cluster to minimise execution times with a relatively small computational overhead. The efficacy of the model is demonstrated by simulating the ultrasound beam pattern for a HIFU sonication of the kidney.


title: The Swift UVOT Stars Survey: I. Methods and Test Clusters

abstract: We describe the motivations and background of a large survey of nearby stel- lar populations using the Ultraviolet Optical Telescope (UVOT) aboard the Swift Gamma-Ray Burst Mission. UVOT, with its wide field, NUV sensitivity, and 2.3 spatial resolution, is uniquely suited to studying nearby stellar populations and providing insight into the NUV properties of hot stars and the contribution of those stars to the integrated light of more distant stellar populations. We review the state of UV stellar photometry, outline the survey, and address problems spe- cific to wide- and crowded-field UVOT photometry. We present color-magnitude diagrams of the nearby open clusters M 67, NGC 188, and NGC 2539, and the globular cluster M 79. We demonstrate that UVOT can easily discern the young- and intermediate-age main sequences, blue stragglers, and hot white dwarfs, pro- ducing results consistent with previous studies. We also find that it characterizes the blue horizontal branch of M 79 and easily identifies a known post-asymptotic giant branch star.


title: Plant clonal morphologies and spatial patterns as self-organized responses to resource-limited environments

abstract: We propose here to interpret and model peculiar plant morphologies (cushions, tussocks) observed in the Andean altiplano as localized structures. Such structures resulting in a patchy, aperiodic aspect of the vegetation cover are hypothesized to self-organize thanks to the interplay between facilitation and competition processes occurring at the scale of basic plant components biologically referred to as 'ramets'. (Ramets are often of clonal origin.) To verify this interpretation, we applied a simple, fairly generic model (one integro-differential equation) emphasizing via Gaussian kernels non-local facilitative and competitive feedbacks of the vegetation biomass density on its own dynamics. We show that under realistic assumptions and parameter values relating to ramet scale, the model can reproduce some macroscopic features of the observed systems of patches and predict values for the inter-patch distance that match the distances encountered in the reference area (Sajama National Park in Bolivia). Prediction of the model can be confronted in the future to data on vegetation patterns along environmental gradients as to anticipate the possible effect of global change on those vegetation systems experiencing constraining environmental conditions.


title: A Computer-Assisted Stability Proof for a Stationary Solution of Reaction-Diffusion Equations

abstract: The main subject of this paper is a computer assisted stability proof for a stationary solution of reaction diffusion equations in one dimensional space. We use Nakao's numerical verification method to enclose a stationary solution of reaction-diffusion equations. Considering the linearized stability of the solution, a method of excluding eigenvalues in a half plane is adopted. We first focus on the eigenvalues for an operator linearized at an approximate solution. An excluding theorem is presented such that we know under some condition, and there is no eigenvalue in some disks. Some computable criteria are constructed to apply the theorem in a computer. And also the invertibility of some operator is proved theoretically in the paper. However, we need the information of the eigenvalues for the operator linearized at the exact solution. This can be obtained by combining with the verification results of the solution. Then we judge the stability of the solution from the domain where the eigenvalues are located. At last there are some verification results.


title: Evolution of Spacetime arises due to the departure from Holographic Equipartition in all Lanczos-Lovelock Theories of Gravity

abstract: In the case of general relativity one can interpret the Noether charge in any bulk region as the heat content $TS$ of its boundary surface. Further, the time evolution of spacetime metric in Einstein's theory arises due to the difference $(N_{sur}-N_{bulk})$ of suitably defined surface and bulk degrees of freedom. We show that this thermodynamic interpretation generalizes in a natural fashion to all Lanczos-Lovelock models of gravity. The Noether charge, related to time evolution vector field, in a bulk region of space is equal to the heat content $TS$ of the boundary surface with the temperature $T$ defined using local Rindler observers and $S$ being the Wald entropy. Using the Wald entropy to define the surface degrees of freedom $N_{sur}$ and Komar energy density to define the bulk degrees of freedom $N_{bulk}$, we can also show that the time evolution of the geometry is sourced by $(N_{sur}-N_{bulk})$. When it is possible to choose the foliation of spacetime such that metric is independent of time, the above dynamical equation yields the holographic equipartition for Lanczos-Lovelock gravity with $N_{sur}=N_{bulk}$. The implications are discussed.


title: Matter wave Fourier optics with a strongly interacting two-dimensional Fermi gas

abstract: We demonstrate and characterize an experimental technique to directly image the momentum distribution of a strongly interacting two-dimensional quantum gas with high momentum resolution. We apply the principles of Fourier optics to investigate three main operations on the expanding gas: focusing, collimation and magnification. We focus the gas in the radial plane using a harmonic confining potential and thus gain access to the momentum distribution. We pulse a different harmonic potential to stop the rapid axial expansion which allows us to image the momentum distribution with high resolution. Additionally, we propose a method to magnify the mapped momentum distribution to access interesting momentum scales. All these techniques can be applied to a wide range of experiments and in particular to study many-body phases of quantum gases.


title: Non-predetermined Model Theory

abstract: This article introduce a new model theory call non-predetermined model theory where functions and relations need not to be determined already and they are determined through time.


title: Effects of grain growth on the interstellar polarization curve

abstract: We apply the time evolution of grain size distributions by accretion and coagulation found in our previous work to the modelling of the wavelength dependence of interstellar linear polarization. We especially focus on the parameters of the Serkowski curve $K$ and $\lambda_{\max}$ characterizing the width and the maximum wavelength of this curve, respectively. We use aligned silicate and non-aligned carbonaceous spheroidal particles with different aspect ratios $a/b$. The imperfect alignment of grains with sizes larger than a cut-off size $r_{V,\rm cut}$ is considered. We find that the evolutionary effects on the polarization curve are negligible in the original model with commonly used material parameters (hydrogen number density $n_\mathrm{H}=10^3$ cm$^{-3}$, gas temperature $T_\mathrm{gas}=10$~K, and the sticking probability for accretion $S_\mathrm{acc}=0.3$). Therefore, we apply the tuned model where the coagulation threshold of silicate is removed. In this model, $\lambda_{\max}$ displaces to the longer wavelengths and the polarization curve becomes wider ($K$ reduces) on time-scales $\sim (30 - 50) (n_\mathrm{H}/10^3 \mathrm{cm}^{-3})^{-1}$ Myr. The tuned models at $T < 30 (n_\mathrm{H}/10^3 \mathrm{cm}^{-3})^{-1} $ Myr and different values of the parameters $r_{V,\rm cut}$ can also explain the observed trend between $K$ and $\lambda_{\max}$. It is significant that the evolutionary effect appears in the perpendicular direction to the effect of $r_{V,\rm cut}$ on the $K$ - $\lambda_{\max}$ diagram. Very narrow polarization curves can be reproduced if we change the type of particles (prolate/oblate) and/or vary $a/b$.