title: Sudden Expansion of a One-Dimensional Bose Gas from Power-Law Traps

abstract: We analyze free expansion of a trapped one-dimensional Bose gas after a sudden release from the confining trap potential. By using the stationary phase and local density approximations, we show that the long-time asymptotic density profile and the momentum distribution of the gas are determined by the initial distribution of Bethe rapidities (quasimomenta) and hence can be obtained from the solutions to the Lieb-Liniger equations in the thermodynamic limit. For expansion from a harmonic trap, and in the limits of very weak and very strong interactions, we recover the self-similar scaling solutions known from the hydrodynamic approach. For all other power-law traps and arbitrary interaction strengths, the expansion is not self-similar and shows strong dependence of the density profile evolution on the trap anharmonicity. We also characterize dynamical fermionization of the expanding cloud in terms of correlation functions describing phase and density fluctuations.


title: Essential Spectra of Induced Operators on Subspaces and Quotients

abstract: Let $X$ be a complex Banach space and let $T$ be a bounded linear operator on $X$. For any closed $T$-invariant subspace $F$ of $X$, $T$ induces operators $T_{|F}:F \longrightarrow F$ and $T/F:X/F\longrightarrow X/F$. In this note, we give a simple proof of the fact that the essential spectra of $T_{|F}$ and $T/F$ are always contained in the polynomial hull of the essential spectrum of $T$.


title: Fast Estimation of Diffusion Tensors under Rician noise by the EM algorithm

abstract: This paper presents a fast computational method, the Expectation Maximization algorithm, for Maximum Likelihood (ML) estimation in diffusion tensor imaging under the Rice noise model. We further extend the ML framework to the maximum a posterior (MAP) estimation and describe the numerical similarities of both ML and MAP estimators. This novel method is implemented and applied using both synthetic and real data in a wide range of b amplitudes. The comparison with other popular methods are made in accuracy, methodology and computation.


title: Cardinal Interpolation With General Multiquadrics

abstract: This paper studies the cardinal interpolation operators associated with the general multiquadrics, $\phi_{\alpha,c}(x) = (\|x\|^2+c^2)^\alpha$, $x\in\mathbb{R}^d$. These operators take the form $$\mathscr{I}_{\alpha,c}\mathbf{y}(x) = \sum_{j\in\mathbb{Z}^d}y_jL_{\alpha,c}(x-j),\quad\mathbf{y}=(y_j)_{j\in\mathbb{Z}^d},\quad x\in\mathbb{R}^d,$$ where $L_{\alpha,c}$ is a fundamental function formed by integer translates of $\phi_{\alpha,c}$ which satisfies the interpolatory condition $L_{\alpha,c}(k) = \delta_{0,k},\; k\in\mathbb{Z}^d$. We consider recovery results for interpolation of bandlimited functions in higher dimensions by limiting the parameter $c\to\infty$. In the univariate case, we consider the norm of the operator $\mathscr{I}_{\alpha,c}$ acting on $\ell_p$ spaces as well as prove decay rates for $L_{\alpha,c}$ using a detailed analysis of the derivatives of its Fourier transform, $\widehat{L_{\alpha,c}}$.


title: Interference of short optical pulses from independent gain-switched laser diodes for quantum secure communications

abstract: Since the introduction of the decoy-state technique, phase-randomised weak coherent light pulses have been the key to increase the practicality of quantum-based communications. Their ultra-fast generation was accomplished via compact gain-switched (GS) lasers, leading to high key rates in quantum key distribution (QKD). Recently, the question arose of whether the same laser could be employed to achieve high-speed measurement-device-independent-QKD, a scheme that promises long-haul quantum communications immune to all detector attacks. For that, a challenging highvisibility interference between independent picosecond optical pulses is required. Here, we answer the above question in the affirmative by demonstrating high-visibility interference from two independent GS lasers triggered at 1GHz. The result is obtained through a careful characterization of the laser frequency chirp and time jitter. By relating these quantities to the interference visibility, we obtain a parameter-free verification of the experimental data and a numerical simulation of the achievable key rates. These findings are beneficial to other applications making use of GS lasers, including random number generation and standard QKD.


title: Predictive Cyber-security Analytics Framework: A non-homogenous Markov model for Security Quantification

abstract: Numerous security metrics have been proposed in the past for protecting computer networks. However we still lack effective techniques to accurately measure the predictive security risk of an enterprise taking into account the dynamic attributes associated with vulnerabilities that can change over time. In this paper we present a stochastic security framework for obtaining quantitative measures of security using attack graphs. Our model is novel as existing research in attack graph analysis do not consider the temporal aspects associated with the vulnerabilities, such as the availability of exploits and patches which can affect the overall network security based on how the vulnerabilities are interconnected and leveraged to compromise the system. Gaining a better understanding of the relationship between vulnerabilities and their lifecycle events can provide security practitioners a better understanding of their state of security. In order to have a more realistic representation of how the security state of the network would vary over time, a nonhomogeneous model is developed which incorporates a time dependent covariate, namely the vulnerability age. The daily transition-probability matrices are estimated using Frei's Vulnerability Lifecycle model. We also leverage the trusted CVSS metric domain to analyze how the total exploitability and impact measures evolve over a time period for a given network.


title: SmO thin films: a flexible route to correlated flat bands with nontrivial topology

abstract: Using density functional theory based calculations, we show that the correlated mixed-valent compound SmO is a 3D strongly topological semi-metal as a result of a 4$f$-5$d$ band inversion at the X point. The [001] surface Bloch spectral density reveals two weakly interacting Dirac cones that are quasi-degenerate at the M_bar-point and another single Dirac cone at the Gamma_bar-point. We also show that the topological non-triviality in SmO is very robust and prevails for a wide range of lattice parameters, making it an ideal candidate to investigate topological nontrivial correlated flat bands in thin-film form. Moreover, the electron filling is tunable by strain. In addition, we find conditions for which the inversion is of the 4f-6s type, making SmO to be a rather unique system. The similarities of the crystal symmetry and the lattice constant of SmO to the well studied ferromagnetic semiconductor EuO, makes SmO/EuO thin film interfaces an excellent contender towards realizing the quantum anomalous Hall effect in a strongly correlated electron system.


title: Interests Diffusion in Social Networks

abstract: Understanding cultural phenomena on Social Networks (SNs) and exploiting the implicit knowledge about their members is attracting the interest of different research communities both from the academic and the business side. The community of complexity science is devoting significant efforts to define laws, models, and theories, which, based on acquired knowledge, are able to predict future observations (e.g. success of a product). In the mean time, the semantic web community aims at engineering a new generation of advanced services by defining constructs, models and methods, adding a semantic layer to SNs. In this context, a leapfrog is expected to come from a hybrid approach merging the disciplines above. Along this line, this work focuses on the propagation of individual interests in social networks. The proposed framework consists of the following main components: a method to gather information about the members of the social networks; methods to perform some semantic analysis of the Domain of Interest; a procedure to infer members' interests; and an interests evolution theory to predict how the interests propagate in the network. As a result, one achieves an analytic tool to measure individual features, such as members' susceptibilities and authorities. Although the approach applies to any type of social network, here it is has been tested against the computer science research community. The DBLP (Digital Bibliography and Library Project) database has been elected as test-case since it provides the most comprehensive list of scientific production in this field.


title: Some qualitative problems in network dynamics

abstract: This paper addresses analytical aspects of deterministic, continuous-time dynamical systems defined on networks. The goal is to model and analyze certain phenomena which must be framed beyond the context of networked dynamical systems, understood as a set of interdependent dynamical systems defined on the nodes of a (possibly evolving) graph. In order to advocate for a more flexible approach to the study of network dynamics, we tackle some qualitative problems which do not fall in this working scenario. First, we address a stability problem on a network of heterogeneous agents, some of which are of dynamic nature while others just impose restrictions on the system behavior. Our second context assumes that the network is clustered, and we address a two-level stability problem involving the dynamics of both individual agents and groups. The aforementioned problems exhibit lines of non-isolated equilibria, and the analysis implicitly assumes a positiveness condition on the edge weights; the removal of this restriction complicates matters, and our third problem concerns a graph-theoretic characterization of the equilibrium set in the dynamics of certain networks with positive and negative weights, the results applying in particular to signed graphs. Our approach combines graph theory and dynamical systems theory, but also uses specific tools coming from linear algebra, including e.g. Gersgorin discs or Maxwell-type determinantal expansions. The results are of application in social, economic, and flow networks, among others, and are also aimed at motivating further research.


title: SHM method for damage localization based on substructuring and VARX models

abstract: A novel damage localization method is proposed, which is based on a substructuring approach and makes use of Vector Auto-Regressive with eXogenous input (VARX) models. The substructuring approach aims to divide the monitored structure into several multi-DOF isolated substructures. Later, each individual substructure is modeled by a VARX model, and the health of each substructure is determined analyzing the variation of the VARX model. The method allows to detect whether the isolated substructure is damaged, and besides allows to locate the damage within the substructure. Only measured displacement data is required to estimate the isolated substructure's VARX model. Moreover, it is not necessary to have a priori knowledge of the structural model. The proposed method is validated by simulations of an eight-storey shear building.


title: What is NP? - Interpretation of a Chinese paradox "white horse is not horse"

abstract: The notion of nondeterminism has disappeared from the current definition of NP, which has led to ambiguities in understanding NP, and caused fundamental difficulties in studying the relation P versus NP. In this paper, we question the equivalence of the two definitions of NP, the one defining NP as the class of problems solvable by a nondeterministic Turing machine in polynomial time, and the other defining NP as the class of problems verifiable by a deterministic Turing machine in polynomial time, and reveal cognitive biases in this equivalence. Inspired from a famous Chinese paradox white horse is not horse, we further analyze these cognitive biases. The work shows that these cognitive biases arise from the confusion between different levels of nondeterminism and determinism, due to the lack of understanding about the essence of nondeterminism. Therefore, we argue that fundamental difficulties in understanding P versus NP lie firstly at cognition level, then logic level.


title: The Force Field for Amino Acid Based Ionic Liquids: Polar Residues

abstract: Ionic liquids (ILs) constitute one of the most active fields of research nowadays. Many organic and inorganic molecules can be converted into ions via relatively simple procedures. These ions can be combined into ILs. Amino acid based ILs (AAILs) represent a specific interest due to solubilization of biological species, participation in enzymatic catalysis, and capturing toxic gases. We develop a new force field (FF) for the seven selected AAILs comprising 1-ethyl-3-methylimidazolium cation and amino acid anions with polar residues. The anions were obtained via deprotonation of carboxyl group. We account for peculiar interactions between the anion and the cation by fitting electrostatic potential for an ion pair, in contrast to isolated ions. The van der Waals interactions were transferred from the CHARMM36 FF with minor modifications, as suggested by hybrid density functional theory. Compatibility between our parameters and CHARMM36 parameters is preserved. The developed interaction model fosters computational investigation of ionic liquids.


title: The Compass for Statistical Researchers

abstract: We have hiked many miles alongside several professors as we traversed our statistical path -- a regime switching trail which changed direction following a class on the foundations of our discipline. As we play the game of research in that limbo between student and academic, one thing among Prof. Bernardi's teachings has never been more clear: to draw a route in the research map you not only need to know your destination, but you must also understand where you are and how you arrived there.


title: Z-score-based modularity for community detection in networks

abstract: Identifying community structure in networks is an issue of particular interest in network science. The modularity introduced by Newman and Girvan [Phys. Rev. E 69, 026113 (2004)] is the most popular quality function for community detection in networks. In this study, we identify a problem in the concept of modularity and suggest a solution to overcome this problem. Specifically, we obtain a new quality function for community detection. We refer to the function as Z-modularity because it measures the Z-score of a given division with respect to the fraction of the number of edges within communities. Our theoretical analysis shows that Z-modularity mitigates the resolution limit of the original modularity in certain cases. Computational experiments using both artificial networks and well-known real-world networks demonstrate the validity and reliability of the proposed quality function.


title: What is Cook's theorem?

abstract: In this paper, we make a preliminary interpretation of Cook's theorem presented in [1]. This interpretation reveals cognitive biases in the proof of Cook's theorem that arise from the attempt of constructing a formula in CNF to represent a computation of a nondeterministic Turing machine. Such cognitive biases are due to the lack of understanding about the essence of nondeterminism, and lead to the confusion between different levels of nondeterminism and determinism, thus cause the loss of nondeterminism from the NP-completeness theory. The work shows that Cook's theorem is the origin of the loss of nondeterminism in terms of the equivalence of the two definitions of NP, the one defining NP as the class of problems solvable by a nondeterministic Turing machine in polynomial time, and the other defining NP as the class of problems verifiable by a deterministic Turing machine in polynomial time. Therefore, we argue that fundamental difficulties in understanding P versus NP lie firstly at cognition level, then logic level.


title: Silicene Nanomesh

abstract: Similar to graphene, zero band gap limits the application of silicene in nanoelectronics despite of its high carrier mobility. By using first-principles calculations, we reveal that a band gap is opened in silicene nanomesh (SNM) when the width W of the wall between the neighboring holes is even. The size of the band gap increases with the reduced W and has a simple relation with the ratio of the removed Si atom and the total Si atom numbers of silicene. Quantum transport simulation reveals that the sub-10 nm single-gated SNM field effect transistors show excellent performance at zero temperature but such a performance is greatly degraded at room temperature.


title: New remarks on DFR noncommutative phase-space

abstract: The so-called canonical noncommutativity is based on a constant noncommutative parameter ($\theta$). However, this formalism breaks Lorentz invariance and one way to recover it is to define the NC parameter as a variable, an extra coordinate of the system. One approach that uses the variable $\theta$ was developed by Doplicher, Fredenhagen and Roberts (DFR) and hence, their phase-space is formed by $(x,p,\theta)$ with extra-dimensions. In this work we have demonstrated precisely that this phase-space is incomplete because the variable $\theta$ requires an associated momentum and the so-called DFR phase-space is in fact formed by $(x, p, \theta, \pi)$, where $\pi$ is an useful object. One of the models used here to demonstrate this fact brought other interesting results. We have used this complete phase-space to explain some undefined results in the $\theta$-variable literature. Finally, we have shown the importance of this DFR-momentum since with it we could fill the gap that exist in $\theta$-variable results. In other words, we have computed the field commutation relations of a QFT in this DFR phase-space. The results obtained here match exactly with the postulated (not demonstrated) values that dwell in the DFR literature.


title: Tur\'an Number of Generalized Triangles

abstract: The family $\Sigma_r$ consists of all $r$-graphs with three edges $D_1,D_2,D_3$ such that $|D_1\cap D_2|=r-1$ and $D_1 \triangle D_2 \subseteq D_3$. A generalized triangle, $\mathcal{T}_r \in \Sigma_r$ is an $r$-graph on $\{1,2,\ldots,2r-1\}$ with three edges $D_1, D_2, D_3$, such that $D_1=\{1,2,\dots,r-1, r\}, D_2= \{1, 2, \dots, r-1, r+1 \}$ and $D_3 = \{r, r+1, \dots, 2r-1\}.$ Frankl and F\"{u}redi conjectured that for all $r\geq 4$, $ex(n,\Sigma_r) = ex(n,\mathcal{T}_r )$ for all sufficiently large $n$ and they also proved it for $r=3$. Later, Pikhurko showed that the conjecture holds for $r=4$. In this paper we determine $ex(n,\mathcal{T}_5)$ and $ex(n,\mathcal{T}_6)$ for sufficiently large $n$, proving the conjecture for $r=5,6$.


title: Some New Results on Binary Relations

abstract: It is well known that if a function from set A to set B has a right inverse then the function is a surjection and the right inverse is an injection. For finite sets, the number of functions, injections, and surjections can also be counted. Relations generalize functions: do similar results exist for relations? This paper proves several new results concerning binary relations. For finite sets, we derive formulas for the number of right total, right unique, left total, and left unique relations. We also provide formulas that count the number of relations that are both right unique and left unique; right unique and right total; and left unique and left total. We conclude by discussing the probability that a relation selected at random is right unique or right total.


title: On the Vanishing Topology of Isolated Cohen-Macaulay Codimension 2 Singularities

abstract: Isolated Cohen-Macaulay codimension 2 singularities share many common features with isolated complete intersection singularities, but they also exhibit some striking new behaviour. One such instance was recently observed by Damon and Pike in their study of the vanishing topology and Euler characteristic, where they took this class of singularities as examples. In this article, we explore their findings further by determining the Betti numbers explicitly and explain the new phenomena. An important tool here is the Tjurina modification relating a Cohen-Macaulay codimension 2 singularity to a finite number of complete intersection singularities.


title: Dynamics of entanglement in systems of identical fermions undergoing decoherence

abstract: Information that is stored in quantum-mechanical systems can be easily lost because of the interaction with the environment in a process known as decoherence. Possible physical implementations of many processes in quantum information theory involve systems of identical particles, whence comprehension of the dynamics of entanglement induced by decoherence processes in identical-particle open systems becomes relevant. Here we study the effects, and concomitant entanglement evolution, arising from the interaction between a system of two identical fermions and the environment, for two paradigmatic quantum channels. New entanglement measures are introduced to quantify the entanglement between the different parties, and a study of the dynamics of entanglement for some particular examples is carried out. Our analysis, which includes also the evolution of an entanglement indicator based on an entropic criteria, offers new insights into the dynamics of entanglement in open systems of identical particles, involving the emergence of multipartite genuine entanglement. The results improve our understanding of the phenomenon of decoherence and will provide new strategies to control it.


title: On the optimal constants in Korn's and geometric rigidity estimates, in bounded and unbounded domains, under Neumann boundary conditions

abstract: We are concerned with the optimal constants: in the Korn inequality under tangential boundary conditions on bounded sets $\Omega \subset \mathbb{R}^n$, and in the geometric rigidity estimate on the whole $\mathbb{R}^2$. We prove that the latter constant equals $\sqrt{2}$, and we discuss the relation of the former constants with the optimal Korn's constants under Dirichlet boundary conditions, and in the whole $\mathbb{R}^n$, which are well known to equal $\sqrt{2}$. We also discuss the attainability of these constants and the structure of deformations/displacement fields in the optimal sets.


title: Ehrenfeucht's lemma in set theory

abstract: Ehrenfeucht's lemma (1973) asserts that whenever one element of a model of Peano arithmetic is definable from another, then they satisfy different types. We consider here the analogue of Ehrenfeucht's lemma for models of set theory. The original argument applies directly to the ordinal-definable elements of any model of set theory, and in particular, Ehrenfeucht's lemma holds fully for models of set theory satisfying $V=HOD$. We show that the lemma can fail, however, in models of set theory with $V\neq HOD$, and it necessarily fails in the forcing extension to add a generic Cohen real. We go on to formulate a scheme of natural parametric generalizations of Ehrenfeucht's lemma, namely, the principles of the form $EL(A,P,Q)$, which asserts that whenever an object $b$ is definable from some $a\in A$ using parameters in $P$, with $b\neq a$, then the types of $a$ and $b$ over $Q$ are different. We also consider various analogues of Ehrenfeucht's lemma obtained by using algebraicity in place of definability, where a set $b$ is algebraic in $a$ if it is a member of a finite set definable from $a$ (as in Hamkins, Leahy arXiv:1305.5953). Ehrenfeucht's lemma holds for the ordinal-algebraic sets, we prove, if and only if the ordinal-algebraic and ordinal-definable sets coincide. Using similar analysis, we answer two open questions posed by Hamkins and Leahy, by showing that (i) algebraicity and definability need not coincide in models of set theory and (ii) the internal and external notions of being ordinal algebraic need not coincide.


title: Paradoxes of cosmological physics in the beginning of the 21-st century

abstract: In the history of cosmology physical paradoxes played important role for development of contemporary world models. Within the modern standard cosmological model there are both observational and conceptual cosmological paradoxes which stimulate to search their solution. Confrontation of theoretical predictions of the standard cosmological model with the latest astrophysical observational data is considered. A review of conceptual problems of the Friedmann space expending models, which are in the bases of modern cosmological model, is discussed. The main paradoxes, which are discussed in modern literature, are the Newtonian character of the exact Friedmann equation, the violation of the energy conservation within any comoving local volume, violation of the limiting recession velocity of galaxies for the observed high redshift objects. Possible observational tests of the nature of the cosmological redshift are discussed


title: A Theoretical Investigation of Decay and Energy Transfer Rates and Efficiencies Near Gold Nanospheres

abstract: We consider the effect of Au nanospheres of subwavelength sizes on the decay and energy transfer rates of quantum systems placed in the proximity of these nanospheres. We find that, for the sphere sizes considered in this contribution, the radiative decay rate is barely affected by the presence of the nanosphere, whereas the non-radiative decay rate is greatly enhanced due to energy transfer from the quantum system to the nanosphere, leading to a strong quenching of the emission of the quantum system. The emission wavelength of the quantum emitter and its intrinsic quantum yield play an important role and the impact of both has to be considered together when investigating their effect on the non-radiative decay rate. The energy transfer process from the emitter to the nanosphere has a complicated distance dependence, with a 1/r^6 regime, characteristic of the F\"orster energy transfer mechanism, but also exhibiting other distance dependence regimes. In the case of a donor-acceptor pair in the presence of a Au nanosphere, the donor couples strongly to the nanosphere, acting as an enhanced dipole; the donor-acceptor energy transfer rate then follows a F\"orster trend, with an increased F\"orster radius. The angular dependence of the energy transfer efficiency between donor and acceptor has a strong dipole-dipole trend for small spheres and deviates from it for larger spheres, especially when the donor and acceptor are on opposite sides of the sphere. The spectral overlap of the donor emission, acceptor absorption and gold nanosphere extinction shows an interesting trend in that the largest F\"orster radius is obtained when the donor emission and acceptor absorption maxima are somewhat red-shifted from the localized surface plasmon peak in the extinction spectrum of the Au nanosphere, being located between it and the near-field scattering maximum.


title: Remote optical sensing on the nanometer scale with a bowtie aperture nano-antenna on a SNOM fiber tip

abstract: Plasmonic nano-antennas have proven the outstanding ability of sensing chemical and physical processes down to the nano-meter scale. Sensing is usually achieved within the highly confined optical fields generated resonantly by the nano-antennas, i.e. in contact to the nano-structures. In these paper, We demonstrate the sensing capability of nano-antennas to their larger scale environment, well beyond their plasmonic confinement volume, leading to the concept of 'remote' (non contact) sensing on the nano-meter scale. On the basis of a bowtie-aperture nano-antenna (BNA) integrated at the apex of a SNOM fiber tip, we introduce an ultra-compact, move-able and background-free optical nano-sensor for the remote sensing of a silicon surface (up to distance of 300 nm). Sensitivity of the BNA to its large scale environment is high enough to expect the monitoring and control of the spacing between the nano-antenna and a silicon surface with sub-nano-meter accuracy. This work paves the way towards a new class of nano-positionning technique, based on nano-antenna resonance monitoring, that are alternative to nano-mechanical and optical interference-based devices.


title: On the parallel sum of positive operators, forms, and functionals

abstract: The parallel sum $A:B$ of two bounded positive linear operators $A,B$ on a Hilbert space $H$ is defined to be the positive operator having the quadratic form \begin{equation*} \inf\{(A(x-y)\,|\,x-y)+(By\,|\,y)\,|\,y\in H\} \end{equation*} for fixed $x\in H$. The purpose of this paper is to provide a factorization of the parallel sum of the form $J_APJ_A^*$ where $J_A$ is the embedding operator of an auxiliary Hilbert space associated with $A$ and $B$, and $P$ is an orthogonal projection onto a certain linear subspace of that Hilbert space. We give similar factorizations of the parallel sum of nonnegative Hermitian forms, positive operators of a complex Banach space $E$ into its topological anti-dual $\bar{E}'$, and of representable positive functionals on a $^*$-algebra.


title: On the Piecewise Parabolic Method for Compressible Flow with Stellar Equations of State

abstract: The piecewise parabolic method and related schemes are widely used to model stellar flows. Several different methods for extending the validity of these methods to a general equation of state have been proposed over time, but direct comparisons amongst one-another and exact solutions with stellar equations of state are not widely available. We introduce some simple test problems with exact solutions run with a popular stellar equation of state and test how two existing codes with different approaches to incorporating general gases perform. The source code for generating the exact solutions is made available.


title: Less is More: Building Selective Anomaly Ensembles

abstract: Ensemble techniques for classification and clustering have long proven effective, yet anomaly ensembles have been barely studied. In this work, we tap into this gap and propose a new ensemble approach for anomaly mining, with application to event detection in temporal graphs. Our method aims to combine results from heterogeneous detectors with varying outputs, and leverage the evidence from multiple sources to yield better performance. However, trusting all the results may deteriorate the overall ensemble accuracy, as some detectors may fall short and provide inaccurate results depending on the nature of the data in hand. This suggests that being selective in which results to combine is vital in building effective ensembles---hence "less is more". In this paper we propose SELECT; an ensemble approach for anomaly mining that employs novel techniques to automatically and systematically select the results to assemble in a fully unsupervised fashion. We apply our method to event detection in temporal graphs, where SELECT successfully utilizes five base detectors and seven consensus methods under a unified ensemble framework. We provide extensive quantitative evaluation of our approach on five real-world datasets (four with ground truth), including Enron email communications, New York Times news corpus, and World Cup 2014 Twitter news feed. Thanks to its selection mechanism, SELECT yields superior performance compared to individual detectors alone, the full ensemble (naively combining all results), and an existing diversity-based ensemble.


title: Higher Lie and Leibniz algebras

abstract: Higher structures - infinity algebras and other objects up to homotopy, categorified algebras, `oidified' concepts, operads, higher categories, higher Lie theory, higher gauge theory... - are currently intensively investigated in Mathematics and Physics. The present thesis deals with abstractions and flexibilizations of algebraic structures, and more precisely, with algebraic operads, homotopification and categorification.


title: Solid solutions with coexisting antiferroelectric and ferroelectric phases as materials with negative refractive index

abstract: A possibility of use of the controlled decomposition of solid solutions of oxides with perovskite structure in the state of coexisting domains of the antiferroelectric (AFE) and ferroelectric (FE) phases for manufacturing of materials with negative refractive index is demonstrated. The lead zirconate titanate based solid solutions are considered as an example of substances suitable for creation of such materials. Manufactured composites constitute a dielectric matrix with a structure of conducting interphase boundaries separating domains of the FE and AFE phases. The electric conductivity of interphase boundaries occurs as a result of the local decomposition of solid solutions in the vicinity of these boundaries. The process of local decomposition and consequently the interphase boundaries conductivity can be controlled by means of external influences.


title: Momentum entanglement in relativistic quantum mechanics

abstract: I present a new group-theoretical approach to the interaction mechanism of elementary particle physics. Within an irreducible unitary two-particle representation of the Poincare group, the commutation relations of the Poincare group require that the two-particle states be momentum entangled. As in gauge theories, momentum entanglement defines a correlation between two particles that can be described as an interaction provided by the exchange of virtual (gauge) quanta. The coupling constant of this interaction is uniquely determined by the structure of the irreducible two-particle state space. For two massive spin one-half particles, the coupling constant matches the empirical value of the electromagnetic coupling constant.


title: Geometric Phase Atom Optics and Interferometry

abstract: We illustrate how geometric gauge forces and topological phase effects emerge in quantum systems without employing assumptions that rely on adiabaticity. We show how geometric magnetism may be harnessed to engineer novel quantum devices including a velocity sieve, a component in mass spectrometers, for neutral atoms or neutrons. We outline a possible experimental setup that demonstrates topological interferometry for neutral spin 1/2 systems.


title: Navigating the Space of Symmetric CMC Surfaces

abstract: In this paper we introduce a flow on the spectral data for symmetric CMC surfaces in the $3$-sphere. The flow is designed in such a way that it changes the topology but fixes the intrinsic (metric) and certain extrinsic (periods) closing conditions of the CMC surfaces. For rational times we obtain closed (possibly branched) connected CMC surfaces of higher genus. We prove the short time existence of this flow near the spectral data of (a class of) CMC tori. In particular we prove that flowing the spectral data for the Clifford torus is equivalent to the flow of Plateau solutions by varying the angle of the fundamental piece in Lawson's construction for the minimal surfaces $\xi_{g,1}.$


title: Design, Analysis, and Simulation of a Pipe-Welding Robot with Fixed Plinth

abstract: Industrial requirements concerning the increased efficiency and high rate of manufacturing result in the development of manufacturer robots, and a vast group of these types of robots is used for welding. This study presented the design, analysis, and simulation of a pipe-welding robot with fixed plinth for a constant circular welding around the pipes. Design of a welding robot capable of keeping the electrode orientation, welding speed, and distance between electrode and pipe surface constant can improve the quality of welding; thus, a five-linked articulated robot was designed for this purpose. Solving of direct and diverse kinematics and dynamics equations of the robot was done by means of Matlab software. The robot was also simulated using a program written in Matlab and the diagrams of angles, velocities, and accelerations of all the arms, and the applied force and torque of each arm required for drive the mechanism were obtained.


title: Focusing of Rayleigh waves generated by high-speed trains under the condition of ground vibration boom

abstract: In the present paper, the effects of focusing of Rayleigh waves generated by high speed trains in the supporting ground under the condition of ground vibration boom are considered theoretically. These effects are similar to the effects of focusing of sound waves radiated by aircraft under the condition of sonic boom. In particular, if a railway track has a bend to provide the possibility of changing direction of train movement, the Rayleigh surface waves generated by high-speed trains under the condition of ground vibration boom may become focused. This results in concentration of their energy along a simple caustic line at one side of the track and in the corresponding increase in ground vibration amplitudes. The effect of focusing of Rayleigh waves may occur also if a train moves along a straight line with acceleration and its current speed is higher than Rayleigh wave velocity in the ground. The obtained results are illustrated by numerical calculations.


title: Long-Distance Entanglement of Soliton Spin Qubits in Gated Nanowires

abstract: We investigate numerically charge, spin, and entanglement dynamics of two electrons confined in a gated semiconductor nanowire. The electrostatic coupling between electrons in the nanowire and the charges in the metal gates leads to a self-trapping of the electrons which results in soliton-like properties. We show that the interplay of an all-electrically controlled coherent transport of the electron solitons and of the exchange interaction can be used to realize ultrafast SWAP and entangling $\sqrt{\text{SWAP}}$ gates for distant spin qubits. We demonstrate that the latter gate can be used to generate a maximally entangled spin state of spatially separated electrons. The results are obtained by quantum mechanical time-dependent calculations with exact inclusion of electron-electron correlations.


title: Nonintrusive coupling of 3D and 2D laminated composite models based on finite element 3D recovery

abstract: In order to simulate the mechanical behavior of large structures assembled from thin composite panels, we propose a coupling technique which substitutes local 3D models for the global plate model in the critical zones where plate modeling is inadequate. The transition from 3D to 2D is based on stress and displacement distributions associated with Saint-Venant problems which are precalculated automatically for a simple 3D cell. The hybrid plate/3D model is obtained after convergence of a series of iterations between a global plate model of the structure and localized 3D models of the critical zones. This technique is nonintrusive because the global calculations can be carried out using commercial software. Evaluation tests show that convergence is fast and that the resulting hybrid model is very close to a full 3D model.


title: From Self-Assembly to Controlled-Assembly, From Optical Manipulation to AFM Manipulation

abstract: Moving nanoparticles/atoms to study the nearfield interaction between them is one of the many approaches to explore the optical and electrical properties of these assemblies. Traditional approach included the self assembly by spinning or drying nanoparticles in aqua on the substrate is well practiced. Lithography technique is another popular approach to deposit limited nano/micro patterns on substrates. Later optical and mechanical manipulations were used to have more control over moving individual elements of nano and microstructures and even atoms. Optical tweezers, optical trapping and AFM manipulation are examples of these precise approaches.


title: Incommensurate magnetic structure, Fe/Cu chemical disorder and magnetic interactions in the high-temperature multiferroic YBaCuFeO5

abstract: Motivated by the recent observations of incommensurate magnetic order and electric polarization in YBaCuFeO5 up to temperatures TN2 as high as 230K [1,2] we report here for the first time a model for the incommensurate magnetic structure of this material that we complement with ab-initio calculations of the magnetic exchange parameters. Using neutron powder diffraction we show that the appearance of polarization below TN2 is accompanied by the replacement of the high temperature collinear magnetic order by a circular inclined spiral with propagation vector ki=(1/2, 1/2, 1/2+-q). Moreover, we find that the polarization approximately scales with the modulus of the magnetic modulation vector q down to the lowest temperature investigated (T=3K). Further, we observe occupational Fe/Cu disorder in the FeO5-CuO5 bipyramids, although a preferential occupation of such units by Fe-Cu pairs is supported by the observed magnetic order and by density functional calculations. We calculate exchange coupling constants for different Fe/Cu distributions and show that, for those containing Fe-Cu dimers, the resulting magnetic order is compatible with the experimentally observed collinear magnetic structure (kc=(1/2, 1/2, 1/2), TN2 > T > TN1 = 440K). Based on these results we discuss possible origins for the incommensurate modulation and its coupling with ferroelectricity.


title: Quasi-power law ensembles

abstract: Quasi-power law ensembles are discussed from the perspective of nonextensive Tsallis distributions characterized by a nonextensive parameter $q$. A number of possible sources of such distributions are presented in more detail. It is further demonstrated that data suggest that nonextensive parameters deduced from Tsallis distributions functions $f\left(p_T\right)$, $q_1$, and from multiplicity distributions (connected with Tsallis entropy), $q_2$, are not identical and that they are connected via $q_1 + q_2 = 2$. It is also shown that Tsallis distributions can be obtained directly from Shannon information entropy, provided some special constraints are imposed. They are connected with the type of dynamical processes under consideration (additive or multiplicative). Finally, it is shown how a Tsallis distribution can accommodate the log-oscillating behavior apparently seen in some multiparticle data.


title: Calibrating an ice sheet model using high-dimensional binary spatial data

abstract: Rapid retreat of ice in the Amundsen Sea sector of West Antarctica may cause drastic sea level rise, posing significant risks to populations in low-lying coastal regions. Calibration of computer models representing the behavior of the West Antarctic Ice Sheet is key for informative projections of future sea level rise. However, both the relevant observations and the model output are high-dimensional binary spatial data; existing computer model calibration methods are unable to handle such data. Here we present a novel calibration method for computer models whose output is in the form of binary spatial data. To mitigate the computational and inferential challenges posed by our approach, we apply a generalized principal component based dimension reduction method. To demonstrate the utility of our method, we calibrate the PSU3D-ICE model by comparing the output from a 499-member perturbed-parameter ensemble with observations from the Amundsen Sea sector of the ice sheet. Our methods help rigorously characterize the parameter uncertainty even in the presence of systematic data-model discrepancies and dependence in the errors. Our method also helps inform environmental risk analyses by contributing to improved projections of sea level rise from the ice sheets.


title: New physics in $\Delta \Gamma_d$

abstract: We analyze the possibility of having new physics effects in the decay rate difference, $\Delta \Gamma_d$, of neutral $B_d$ mesons. Three different sources of enhancement are considered, CKM unitarity violations, beyond standard model effects in the tree-level dimension six operators $(\bar{d}p)(\bar{p}'b)$ with $p,p'=u,c$; and large enhancements of the almost unconstrained operators $(\bar{d}b)(\bar{\tau}\tau)$. We find that deviations of several hundred per cent from the standard model prediction of $\Delta \Gamma_d$ are not excluded by current experimental data.


title: Where Graph Topology Matters: The Robust Subgraph Problem

abstract: Robustness is a critical measure of the resilience of large networked systems, such as transportation and communication networks. Most prior works focus on the global robustness of a given graph at large, e.g., by measuring its overall vulnerability to external attacks or random failures. In this paper, we turn attention to local robustness and pose a novel problem in the lines of subgraph mining: given a large graph, how can we find its most robust local subgraph (RLS)? We define a robust subgraph as a subset of nodes with high communicability among them, and formulate the RLS-PROBLEM of finding a subgraph of given size with maximum robustness in the host graph. Our formulation is related to the recently proposed general framework for the densest subgraph problem, however differs from it substantially in that besides the number of edges in the subgraph, robustness also concerns with the placement of edges, i.e., the subgraph topology. We show that the RLS-PROBLEM is NP-hard and propose two heuristic algorithms based on top-down and bottom-up search strategies. Further, we present modifications of our algorithms to handle three practical variants of the RLS-PROBLEM. Experiments on synthetic and real-world graphs demonstrate that we find subgraphs with larger robustness than the densest subgraphs even at lower densities, suggesting that the existing approaches are not suitable for the new problem setting.


title: What flows in the chiral magnetic effect? -- Simulating the particle production with CP-breaking backgrounds

abstract: To address a question of whether the chiral magnetic current is a static polarization or a genuine flow of charged particles, we elucidate the numerical formulation to simulate the net production of right-handed particles and anomalous currents with CP-breaking background fields which cause an imbalance between particles and anti-particles. For a concrete demonstration we numerically impose pulsed electric and magnetic fields to confirm our answer to the question that the produced net particles flow in the dynamical chiral magnetic effect. The rate for the particle production and the chiral magnetic current generation is quantitatively consistent with the axial anomaly, while they appear with a finite response time. We emphasize the importance to quantify the response time that would suppress observable effects of the anomalous current.


title: Bloofi: Multidimensional Bloom Filters

abstract: Bloom filters are probabilistic data structures commonly used for approximate membership problems in many areas of Computer Science (networking, distributed systems, databases, etc.). With the increase in data size and distribution of data, problems arise where a large number of Bloom filters are available, and all them need to be searched for potential matches. As an example, in a federated cloud environment, each cloud provider could encode the information using Bloom filters and share the Bloom filters with a central coordinator. The problem of interest is not only whether a given element is in any of the sets represented by the Bloom filters, but which of the existing sets contain the given element. This problem cannot be solved by just constructing a Bloom filter on the union of all the sets. Instead, we effectively have a multidimensional Bloom filter problem: given an element, we wish to receive a list of candidate sets where the element might be. To solve this problem, we consider 3 alternatives. Firstly, we can naively check many Bloom filters. Secondly, we propose to organize the Bloom filters in a hierarchical index structure akin to a B+ tree, that we call Bloofi. Finally, we propose another data structure that packs the Bloom filters in such a way as to exploit bit-level parallelism, which we call Flat-Bloofi. Our theoretical and experimental results show that Bloofi and Flat-Bloofi provide scalable and efficient solutions alternatives to search through a large number of Bloom filters.


title: A regularized representation of the fractional Laplacian in n dimensions and its relation to Weierstrass-Mandelbrot type fractal functions

abstract: We demonstrate that the fractional Laplacian (FL) is the principal characteristic operator of harmonic systems with {\it self-similar} interparticle interactions. We show that the FL represents the "{\it fractional continuum limit}" of a discrete "self-similar Laplacian" which is obtained by Hamilton's variational principle from a discrete spring model. We deduce from generalized self-similar elastic potentials regular representations for the FL which involve convolutions of symmetric finite difference operators of even orders extending the standard representation of the FL. Further we deduce a regularized representation for the FL $-(-\Delta)^{\frac{\alpha}{2}}$ holding for $\alpha\in \R \geq 0$. We give an explicit proof that the regularized representation of the FL gives for integer powers $\frac{\alpha}{2} \in \N\_0$ a distributional representation of the standard Laplacian operator $\Delta$ including the trivial unity operator for $\alpha\rightarrow 0$. We demonstrate that self-similar {\it harmonic} systems are {\it all} governed in a distributional sense by this {\it regularized representation of the FL} which therefore can be conceived as characteristic footprint of self-similarity.


title: Stability enhancement by joint phase measurements in a single cold atomic fountain

abstract: We propose a method of joint interrogation in a single atom interferometer which overcomes the dead time between consecutive measurements in standard cold atomic fountains. The joint operation enables for a faster averaging of the Dick effect associated with the local oscillator noise in clocks and with vibration noise in cold atom inertial sensors. Such an operation allows achieving the lowest stability limit due to atom shot noise. We demonstrate a multiple joint operation in which up to five clouds of atoms are interrogated simultaneously in a single setup. The essential feature of multiple joint operation, demonstrated here for a micro-wave Ramsey interrogation, can be generalized to go beyond the current stability limit associated with dead times in present-day cold atom interferometer inertial sensors.


title: Separating subadditive Euclidean functionals

abstract: If we are given $n$ random points in the hypercube $[0,1]^d$, then the minimum length of a Traveling Salesperson Tour through the points, the minimum length of a spanning tree, and the minimum length of a matching, etc., are known to be asymptotically $\beta n^{\frac{d-1}{d}}$ a.s., where $\beta$ is an absolute constant in each case. We prove separation results for these constants. In particular, concerning the constants $\beta_{\mathrm{TSP}}^d$, $\beta_{\mathrm{MST}}^d$, $\beta_{\mathrm{MM}}^d$, and $\beta_{\mathrm{TF}}^d$ from the asymptotic formulas for the minimum length TSP, spanning tree, matching, and 2-factor, respectively, we prove that $\beta_{\mathrm{MST}}^d<\beta_{\mathrm{TSP}}^d$, $2\beta_{\mathrm{MM}}^d<\beta_{\mathrm{TSP}}^d$, and $\beta_{\mathrm{TF}}^d<\beta_{\mathrm{TSP}}^d$ for all $d\geq 2$. We also asymptotically separate the TSP from its linear programming relaxation in this setting. Our results have some computational relevance, showing that a certain natural class of simple algorithms cannot solve the random Euclidean TSP efficiently.


title: Determination of time dependent factors of coefficients in fractional diffusion equations

abstract: We consider fractional diffusion equations and study the stability of the inverse problem of determining the time-dependent parameter in a source term or a coefficient of zero-th order term from observations of the solution at one point in a bounded domain.


title: Multi-Beam RF Aperture Using Multiplierless FFT Approximation

abstract: Multiple independent radio frequency (RF) beams find applications in communications, radio astronomy, radar, and microwave imaging. An $N$-point FFT applied spatially across an array of receiver antennas provides $N$-independent RF beams at $\frac{N}{2}\log_2N$ multiplier complexity. Here, a low-complexity multiplierless approximation for the 8-point FFT is presented for RF beamforming, using only 26 additions. The algorithm provides eight beams that closely resemble the antenna array patterns of the traditional FFT-based beamformer albeit without using multipliers. The proposed FFT-like algorithm is useful for low-power RF multi-beam receivers; being synthesized in 45 nm CMOS technology at 1.1 V supply, and verified on-chip using a Xilinx Virtex-6 Lx240T FPGA device. The CMOS simulation and FPGA implementation indicate bandwidths of 588 MHz and 369 MHz, respectively, for each of the independent receive-mode RF beams.


title: The dual fibration in elementary terms

abstract: We give an elementary construction of the dual fibration of a fibration. It does not use the non-elementary notion of (pseudo-) functor into the category of categories.


title: Boundary between Stable and Unstable Regimes of Accretion. Ordered and Chaotic Unstable Regimes

abstract: We present a study of the Rayleigh-Taylor unstable regime of accretion onto rotating magnetized stars in a set of high grid resolution three-dimensional (3D) magnetohydrodynamic (MHD) simulations performed in low-viscosity discs. We find that the boundary between the stable and unstable regimes is determined almost entirely by the fastness parameter omega_s=Omega_star/Omega_K(r_m), where Omega_star is the angular velocity of the star and Omega_K(r_m) is the angular velocity of the Keplerian disc at the disc-magnetosphere boundary r=r_m. We found that accretion is unstable if omega_s < 0.6. Accretion through instabilities is present in stars with different magnetospheric sizes. However, only in stars with relatively small magnetospheres, r_m/R_star < 7, do the unstable tongues produce chaotic hot spots on the stellar surface and irregular light-curves. At even smaller values of the fastness parameter, omega_s < 0.45, multiple irregular tongues merge, forming one or two ordered unstable tongues that rotate with the angular frequency of the inner disc. This transition occurs in stars with even smaller magnetospheres, r_m/R_star < 4.2. Most of our simulations were performed at a small tilt of the dipole magnetosphere, Theta=5 degrees, and a small viscosity parameter alpha=0.02. Test simulations at higher alpha values show that many more cases become unstable, and the light-curves become even more irregular. Test simulations at larger tilts of the dipole Theta show that instability is present, however, accretion in two funnel streams dominates if Theta > 15 degrees. The results of these simulations can be applied to accreting magnetized stars with relatively small magnetospheres: Classical T Tauri stars, accreting millisecond X-ray pulsars, and cataclysmics variables.


title: $P$-Paracompact and $P$-Metrizable Spaces

abstract: Let $P$ be a directed set and $X$ a space. A collection $\mathcal{C}$ of subsets of $X$ is \emph{$P$-locally finite} if $\mathcal{C}=\bigcup \{ \mathcal{C}_p : p \in P\}$ where (i) if $p \le p'$ then $\mathcal{C}_p \subseteq \mathcal{C}_{p'}$ and (ii) each $\mathcal{C}_p$ is locally finite. Then $X$ is \emph{$P$-paracompact} if every open cover has a $P$-locally finite open refinement. Further, $X$ is \emph{$P$-metrizable} if it has a $(P \times \mathbb{N})$-locally finite base. This work provides the first detailed study of $P$-paracompact and $P$-metrizable spaces, particularly in the case when $P$ is a $\mathcal{K}(M)$, the set of all compact subsets of a separable metrizable space $M$ ordered by set inclusion.


title: Robust estimation of precision matrices under cellwise contamination

abstract: There is a great need for robust techniques in data mining and machine learning contexts where many standard techniques such as principal component analysis and linear discriminant analysis are inherently susceptible to outliers. Furthermore, standard robust procedures assume that less than half the observation rows of a data matrix are contaminated, which may not be a realistic assumption when the number of observed features is large. This work looks at the problem of estimating covariance and precision matrices under cellwise contamination. We consider using a robust pairwise covariance matrix as an input to various regularisation routines, such as the graphical lasso, QUIC and CLIME. To ensure the input covariance matrix is positive semidefinite, we use a method that transforms a symmetric matrix of pairwise covariances to the nearest covariance matrix. The result is a potentially sparse precision matrix that is resilient to moderate levels of cellwise contamination. Since this procedure is not based on subsampling it scales well as the number of variables increases.


title: Non-conservative Lagrangian method for half-dark solitons in spinor non-equilibrium Polariton condensates

abstract: In this work I introduce a powerful analytical method to analyze states of non-equilibrium polariton Bose-Einstein condensates (BEC). It is shown that the condensate wave functions carrying dark solitons and half-dark solitons can be expressed in terms of the natural parameters of the non-equilibrium system, i.e. the growth and decay rates and the interaction parameters. The derivation is based on a non-conservative Lagrangian formalism, which provides ordinary differential equations for the parameters of the dark soliton solutions in their dynamic environment. Explicit expressions for the stationary solution are stated, which significantly differ from those observed in atomic BEC. Subsequently the method is extended to spinor condensates, which yields ordinary differential equations for the parameters of half-dark solitons. Finally a stationary case with explicit expressions for half-dark solitons is presented as well.


title: Complete All-Optical Quantum Control of Electron Spins in InAs/GaAs Quantum Dot Molecule

abstract: The spin states of electrons and holes confined in InAs quantum dot molecules have recently come to fore as a promising system for the storage or manipulation of quantum information. We describe here a feasible scheme for complete quantum optical control of two electron spin qubits in two vertically-stacked singly-charged InAs quantum dots coupled by coherent electron tunneling. With an applied magnetic field transverse to the growth direction, we construct a universal set of gates that corresponds to the possible Raman transitions between the spin states. We detail the procedure to decompose a given two-qubit unitary operation, so as to realize it with a successive application of up to 8 of these gates. We give the pulse shapes for the laser pulses used to implement this universal set of gates and demonstrate the realization of the two-qubit quantum Fourier transform with fidelity of 0.881 and duration of 414 ps. Our proposal therefore offers an accessible path to universal computation in quantum dot molecules and points to the advantages of using pulse shaping incoherent manipulation of optically active quantum dots to mitigate the negative effects of unintended dynamics and spontaneous emission.


title: Low frequency gray-body factors and infrared divergences: rigorous results

abstract: Formal solutions to the mode equations for both spherically symmetric black holes and Bose-Einstein condensate acoustic black holes are obtained by writing the spatial part of the mode equation as a linear Volterra integral equation of the second kind. The solutions work for a massless minimally coupled scalar field in the s-wave or zero angular momentum sector for a spherically symmetric black hole and in the longitudinal sector of a 1D Bose-Einstein condensate acoustic black hole. These solutions are used to obtain in a rigorous way analytic expressions for the scattering coefficients and gray-body factors in the zero frequency limit. They are also used to study the infrared behaviors of the symmetric two-point function and two functions derived from it: the point-split stress-energy tensor for the massless minimally coupled scalar field in Schwarzschild-de Sitter spacetime and the density-density correlation function for a Bose-Einstein condensate acoustic black hole.


title: On financial applications of the two-parameter Poisson-Dirichlet distribution

abstract: Capital distribution curve is defined as log-log plot of normalized stock capitalizations ranked in descending order. The curve displays remarkable stability over periods of time. Theory of exchangeable distributions on set partitions, developed for purposes of mathematical genetics and recently applied in non-parametric Bayesian statistics, provides probabilistic-combinatorial approach for analysis and modeling of the capital distribution curve. Framework of the two-parameter Poisson-Dirichlet distribution contains rich set of methods and tools, including infinite-dimensional diffusion process. The purpose of this note is to introduce framework of exchangeable distributions on partitions in the financial context. In particular, it is shown that averaged samples from the Poisson-Dirichlet distribution provide approximation to the capital distribution curves in equity markets. This suggests that the two-parameter model can be employed for modelling evolution of market weights and prices fluctuating in stochastic equilibrium.


title: A Simple Construction of Recursion Operators for Multidimensional Dispersionless Integrable Systems

abstract: We present a simple novel construction of recursion operators for integrable multidimensional dispersionless systems that admit a Lax pair whose operators are linear in the spectral parameter and do not involve the derivatives with respect to the latter. New examples of recursion operators obtained using our technique include {\em inter alia} those for the general heavenly equation, which describes a class of anti-self-dual solutions of the vacuum Einstein equations, and a six-dimensional equation resulting from a system of Ferapontov and Khusnutdinova.


title: On a new integrable discretization of the derivative nonlinear Schr\"odinger (Chen-Lee-Liu) equation

abstract: We propose a general integrable lattice system involving some free parameters, which contains known integrable lattice systems such as the Ablowitz-Ladik discretization of the nonlinear Schr\"odinger (NLS) equation as special cases. With a suitable choice of the parameters, it provides a new integrable space-discretization of the derivative NLS equation known as the Chen-Lee-Liu equation. Analogously to the continuous case, the space-discrete Chen-Lee-Liu system possesses a Lax pair and admits a complex conjugation reduction between the two dependent variables. Thus, we obtain a proper space-discretization of the Chen-Lee-Liu equation defined on the three lattice sites $n-1$, $n$, $n+1$ for the first time. Considering a negative flow of the discrete Chen-Lee-Liu hierarchy, we obtain a proper discretization of the massive Thirring model in light-cone coordinates. Multicomponent generalizations of the obtained discrete equations are straightforward because the performed computations are valid for the general case where the dependent variables are vector- or matrix-valued.


title: Finite-SNR Bounds on the Sum-Rate Capacity of Rayleigh Block-Fading Multiple-Access Channels with no a Priori CSI

abstract: We provide nonasymptotic upper and lower bounds on the sum-rate capacity of Rayleigh block-fading multiple-access channels for the setup where a priori channel state information is not available. The upper bound relies on a dual formula for channel capacity and on the assumption that the users can cooperate perfectly. The lower bound is derived assuming a noncooperative scenario, where each user employs unitary space-time modulation (independently from the other users). Numerical results show that the gap between the upper and the lower bound is small already at moderate SNR values. This suggests that the sum-rate capacity gains obtainable through user cooperation are minimal.


title: Nonsingular Bouncing Model in Closed and Open universe

abstract: We investigate the cosmology of a class of model with noncanonical scalar field and matter both in FRW closed and open background. Writing the Einstein Equations in terms of dimensionless dynamical variables suitable for studying bouncing solutions, a fixed point analysis is carried out. Cosmological solutions satisfying the stability and bouncing conditions are obtained.


title: An adaptively refined phase-space element method for cosmological simulations and collisionless dynamics

abstract: N-body simulations are essential for understanding the formation and evolution of structure in the Universe. However, the discrete nature of these simulations affects their accuracy when modelling collisionless systems. We introduce a new approach to simulate the gravitational evolution of cold collisionless fluids by solving the Vlasov-Poisson equations in terms of adaptively refineable "Lagrangian phase space elements". These geometrical elements are piecewise smooth maps between Lagrangian space and Eulerian phase space and approximate the continuum structure of the distribution function. They allow for dynamical adaptive splitting to accurately follow the evolution even in regions of very strong mixing. We discuss in detail various one-, two- and three-dimensional test problems to demonstrate the performance of our method. Its advantages compared to N-body algorithms are: i) explicit tracking of the fine-grained distribution function, ii) natural representation of caustics, iii) intrinsically smooth gravitational potential fields, thus iv) eliminating the need for any type of ad-hoc force softening. We show the potential of our method by simulating structure formation in a warm dark matter scenario. We discuss how spurious collisionality and large-scale discreteness noise of N-body methods are both strongly suppressed, which eliminates the artificial fragmentation of filaments. Therefore, we argue that our new approach improves on the N-body method when simulating self-gravitating cold and collisionless fluids, and is the first method that allows to explicitly follow the fine-grained evolution in six-dimensional phase space.


title: How to bend galaxy disc profiles: the role of halo spin

abstract: The radial density profiles of stellar galaxy discs can be well approximated as an exponential. Compared to this canonical form, however, the profiles in the majority of disc galaxies show downward or upward breaks at large radii. Currently, there is no coherent explanation in a galaxy formation context of the radial profile per se, along with the two types of profile breaks. Using a set of controlled hydrodynamic simulations of disc galaxy formation, we find a correlation between the host halo's initial angular momentum and the resulting radial profile of the stellar disc: galaxies that live in haloes with a low spin parameter {\lambda} <~ 0.03 show an up-bending break in their disc density profiles, while galaxies in haloes of higher angular momentum show a down-bending break. We find that the case of pure exponential profiles ({\lambda} ~ 0.035) coincides with the peak of the spin parameter distribution from cosmological simulations. Our simulations not only imply an explanation of the observed behaviours, but also suggest that the physical origin of this effect is related to the amount of radial redistribution of stellar mass, which is anti-correlated with {\lambda}.


title: Neutrino-driven supernova of a low-mass iron-core progenitor boosted by three-dimensional turbulent convection

abstract: We present the first successful simulation of a neutrino-driven supernova explosion in three dimensions (3D), using the Prometheus-Vertex code with an axis-free Yin-Yang grid and a sophisticated treatment of three-flavor, energy-dependent neutrino transport. The progenitor is a nonrotating, zero-metallicity 9.6 Msun star with an iron core. While in spherical symmetry outward shock acceleration sets in later than 300 ms after bounce, a successful explosion starts at ~130 ms postbounce in two dimensions (2D). The 3D model explodes at about the same time but with faster shock expansion than in 2D and a more quickly increasing and roughly 10 percent higher explosion energy of >10^50 erg. The more favorable explosion conditions in 3D are explained by lower temperatures and thus reduced neutrino emission in the cooling layer below the gain radius. This moves the gain radius inward and leads to a bigger mass in the gain layer, whose larger recombination energy boosts the explosion energy in 3D. These differences are caused by less coherent, less massive, and less rapid convective downdrafts associated with postshock convection in 3D. The less violent impact of these accretion downflows in the cooling layer produces less shock heating and therefore diminishes energy losses by neutrino emission. We thus have, for the first time, identified a reduced mass accretion rate, lower infall velocities, and a smaller surface filling factor of convective downdrafts as consequences of 3D postshock turbulence that facilitate neutrino-driven explosions and strengthen them compared to the 2D case.


title: Natural Supersymmetry without Light Higgsinos

abstract: We present a mechanism that allows a large Higgsino mass without large fine-tuning. The Higgs is a pseudo Nambu-Goldstone boson (PNGB) of the global symmetry breaking pattern $SO(5) \to SO(4)$. Because of the PNGB nature of the light Higgs, the $SO(5)$ invariant Higgsino mass does not directly contribute to the Higgs mass. Large couplings in the Higgs sector that spontaneously breaks $SO(5)$ minimize the tuning, and are also motivated by the requirements of generating a sufficiently large Higgs quartic coupling and of maintaining a natural approximate global $SO(5)$ symmetry. When these conditions are imposed, theories of this type predict heavy Higgsinos. This construction differs from composite Higgs models in that no new particles are introduced to form complete $SO(5)$ multiplets involving the top quark---the stop is the only top partner. Compatibility with Higgs coupling measurements requires cancelations among contributions to the Higgs mass squared parameter at the 10% level. An important implication of this construction is that the compressed region of stop and sbottom searches can still be natural.


title: Green Valley Galaxies

abstract: The "green valley" is a wide region separating the blue and the red peaks in the ultraviolet-optical color magnitude diagram, first revealed using GALEX UV photometry. The term was coined by Christopher Martin in 2005. Green valley highlights the discriminating power of UV to very low relative levels of ongoing star formation, to which the optical colors, including u-r, are insensitive. It corresponds to massive galaxies below the star-forming "main" sequence, and therefore represents a critical tool for the study of the quenching of star formation and its possible resurgence in otherwise quiescent galaxies. This article reviews the results pertaining to morphology, structure, environment, dust content and gas properties of green valley galaxies in the local universe. Their relationship to AGN is also discussed. Attention is given to biases emerging from defining the "green valley" using optical colors. We review various evolutionary scenarios and we present evidence for a new, quasi-static view of the green valley, in which the majority of galaxies currently in the green valley were only partially quenched in the distant past and now participate in a slow cosmic decline of star formation, which also drives down the activity on the main sequence, presumably as a result of the dwindling accretion/cooling onto galaxy disks.


title: New Spatially Resolved Observations of the T Cha Transition Disk and Constraints on the Previously Claimed Substellar Companion

abstract: We present multi-epoch non-redundant masking observations of the T Cha transition disk, taken at the VLT and Magellan in H, Ks, and L' bands. T Cha is one of a small number of transition disks that host companion candidates discovered by high-resolution imaging techniques, with a putative companion at a position angle of 78 degrees, separation of 62 mas, and contrast at L' of 5.1 mag. We find comparable binary parameters in our re-reduction of the initial detection images, and similar parameters in the 2011 L', 2013 NaCo L', and 2013 NaCo Ks data sets. We find a close-in companion signal in the 2012 NaCo L' dataset that cannot be explained by orbital motion, and a non-detection in the 2013 MagAO/Clio2 L' data. However, Monte-carlo simulations show that the best fits to the 2012 NaCo and 2013 MagAO/Clio2 followup data may be consistent with noise. There is also a significant probability of false non-detections in both of these data sets. We discuss physical scenarios that could cause the best fits, and argue that previous companion and scattering explanations are inconsistent with the results of the much larger dataset presented here.


title: The First Calculation of Fractional Jets

abstract: In collider physics, jet algorithms are a ubiquitous tool for clustering particles into discrete jet objects. Event shapes offer an alternative way to characterize jets, and one can define a jet multiplicity event shape, which can take on fractional values, using the framework of "jets without jets". In this paper, we perform the first analytic studies of fractional jet multiplicity $\tilde{N}_{\rm jet}$ in the context of $e^+e^-$ collisions. We use fixed-order QCD to understand the $\tilde{N}_{\rm jet}$ cross section at order $\alpha_s^2$, and we introduce a candidate factorization theorem to capture certain higher-order effects. The resulting distributions have a hybrid jet algorithm/event shape behavior which agrees with parton shower Monte Carlo generators. The $\tilde{N}_{\rm jet}$ observable does not satisfy ordinary soft-collinear factorization, and the $\tilde{N}_{\rm jet}$ cross section exhibits a number of unique features, including the absence of collinear logarithms and the presence of soft logarithms that are purely non-global. Additionally, we find novel divergences connected to the energy sharing between emissions, which are reminiscent of rapidity divergences encountered in other applications. Given these interesting properties of fractional jet multiplicity, we advocate for future measurements and calculations of $\tilde{N}_{\rm jet}$ at hadron colliders like the LHC.


title: The Stellar Mass - Halo Mass Relation for Low Mass X-ray Groups at 0.5<z<1 in the CDFS with CSI

abstract: Since z~1, the stellar mass density locked in low mass groups and clusters has grown by a factor of ~8. Here we make the first statistical measurements of the stellar mass content of low mass X-ray groups at 0.5<z<1, enabling the calibration of stellar-to-halo mass scales for wide-field optical and infrared surveys. Groups are selected from combined Chandra and XMM-Newton X-ray observations in the Chandra Deep Field South (CDFS). These ultra-deep observations allow us to identify bona fide low mass groups at high redshift and enable measurements of their total halo masses. We compute aggregate stellar masses for these halos using galaxies from the Carnegie-Spitzer-IMACS (CSI) spectroscopic redshift survey. Stars comprise ~3-4% of the total mass of group halos with masses 10^{12.8}<M200/Msun<10^{13.5} (about the mass of Fornax and 1/50th the mass of Virgo). Complementing our sample with higher mass halos at these redshifts, we find that the stellar-to-halo mass ratio decreases toward higher halo masses, consistent with other work in the local and high redshift universe. The observed scatter about the stellar-halo mass relation is ~0.25 dex, which is relatively small and suggests that total group stellar mass can serve as a rough proxy for halo mass. We find no evidence for any significant evolution in the stellar-halo mass relation since z<1. Quantifying the stellar content in groups since this epoch is critical given that hierarchical assembly leads to such halos growing in number density and hosting increasing shares of quiescent galaxies.


title: The First Billion Years Project: The escape fraction of ionizing photons in the epoch of reionization

abstract: Proto-galaxies forming in low-mass dark matter haloes are thought to provide the majority of ionizing photons needed to reionize the Universe, due to their high escape fractions of ionizing photons. We study how the escape fraction in high-redshift galaxies relates to the physical properties of the halo in which the galaxies form, by computing escape fractions in more than 75000 haloes between redshifts 27 and 6 that were extracted from the First Billion Years project, high-resolution cosmological hydrodynamics simulations of galaxy formation. We find that the main constraint on the escape fraction is the gas column density in a radius of 10 pc around the stellar populations, causing a strong mass dependence of the escape fraction. The lower potential well in haloes with virial mass below 1e8 solar mass results in low column densities that can be penetrated by radiation from young stars (age < 5 Myr). In haloes with virial mass above 1e8 solar mass supernova feedback is important, but only 30% of the haloes in this mass range have an escape fraction higher than 1%. We find a large range of escape fractions in haloes with similar properties, caused by different distributions of the dense gas in the halo. This makes it very hard to predict the escape fraction on the basis of halo properties and results in a highly anisotropic escape fraction. The strong mass dependence, the large spread and the large anisotropy of the escape fraction may strongly affect the topology of reionization and is something current models of cosmic reionization should strive to take into account.


title: The origin of the spurious iron spread in the globular cluster NGC 3201

abstract: NGC 3201 is a globular cluster suspected to have an intrinsic spread in the iron content. We re-analysed a sample of 21 cluster stars observed with UVES-FLAMES at the Very Large Telescope and for which Simmerer et al. found a 0.4 dex wide [Fe/H] distribution with a metal-poor tail. We confirmed that when spectroscopic gravities are adopted, the derived [Fe/H] distribution spans ~0.4 dex. On the other hand, when photometric gravities are used, the metallicity distribution from Fe I lines remains large, while that derived from Fe II lines is narrow and compatible with no iron spread. We demonstrate that the metal-poor component claimed by Simmerer et al. is composed by asymptotic giant branch stars that could be affected by non local thermodynamical equilibrium effects driven by iron overionization. This leads to a decrease of the Fe I abundance, while leaving the Fe II abundance unaltered. A similar finding has been already found in asymptotic giant branch stars of the globular clusters M5 and 47 Tucanae. We conclude that NGC 3201 is a normal cluster, with no evidence of intrinsic iron spread.


title: The stellar kinematics of co-rotating spiral arms in Gaia mock observations

abstract: We have observed an N-body/Smoothed Particle Hydrodynamics simulation of a Milky Way like barred spiral galaxy. We present a simple method that samples N-body model particles into mock Gaia stellar observations and takes into account stellar populations, dust extinction and Gaia's science performance estimates. We examine the kinematics around a nearby spiral arm at a similar position to the Perseus arm at three lines of sight in the disc plane; (l,b)=(90,0), (120,0) and (150,0) degrees. We find that the structure of the peculiar kinematics around the co-rotating spiral arm, which is found in Kawata et al. (2014b), is still visible in the observational data expected to be produced by Gaia despite the dust extinction and expected observational errors of Gaia. These observable kinematic signatures will enable testing whether the Perseus arm of the Milky Way is similar to the co-rotating spiral arms commonly seen in N-body simulations.


title: 21-cm signatures of residual HI inside cosmic HII regions during reionization

abstract: We investigate the impact of sinks of ionizing radiation on the reionization-era 21-cm signal, focusing on 1-point statistics. We consider sinks in both the intergalactic medium and inside galaxies. At a fixed filling factor of HII regions, sinks will have two main effects on the 21-cm morphology: (i) as inhomogeneous absorbers of ionizing photons they result in smaller and more widespread cosmic HII patches; and (ii) as reservoirs of neutral gas they contribute a non-zero 21-cm signal in otherwise ionized regions. Both effects damp the contrast between neutral and ionized patches during reionization, making detection of the epoch of reionization with 21-cm interferometry more challenging. Here we systematically investigate these effects using the latest semi-numerical simulations. We find that sinks dramatically suppress the peak in the redshift evolution of the variance, corresponding to the midpoint of reionization. As previously predicted, skewness changes sign at midpoint, but the fluctuations in the residual HI suppress a late-time rise. Furthermore, large levels of residual HI dramatically alter the evolution of the variance, skewness and power spectrum from that seen at lower levels. In general, the evolution of the large-scale modes provides a better, cleaner, higher signal-to-noise probe of reionization.


title: A semi-classical limit for the many-body localization transition

abstract: We introduce a semi-classical limit for many-body localization in the absence of global symmetries. Microscopically, this limit is realized by disordered Floquet circuits composed of Clifford gates. In $d=1$, the resulting dynamics are always many-body localized with a complete set of strictly local integrals of motion. In $d\geq 2$, the system realizes both localized and delocalized phases separated by a continuous transition in which ergodic puddles percolate. We argue that the phases are stable to deformations away from the semi-classical limit and estimate the resulting phase boundary. The Clifford circuit model is a distinct tractable limit from that of free fermions and suggests bounds on the critical exponents for the generic transition.


title: Satellite Dwarf Galaxies in a Hierarchical Universe: Infall Histories, Group Preprocessing, and Reionization

abstract: In the Local Group (LG), almost all satellite dwarf galaxies that are within the virial radius of the Milky Way (MW) and Andromeda (M31) exhibit strong environmental influence. The orbital histories of these satellites provide the key to understanding the role of the MW/M31 halo, lower-mass groups, and cosmic reionization on the evolution of dwarf galaxies. We examine the virial-infall histories of satellites with M_star = 10^{3-9} M_sun using the ELVIS suite of cosmological zoom-in dissipationless simulations of 48 MW/M31-like halos. Satellites at z = 0 fell into the MW/M31 halos typically 5 - 8 Gyr ago at z = 0.5 - 1. However, they first fell into any host halo typically 7 - 10 Gyr ago at z = 0.7 - 1.5. This difference arises because many satellites experienced "group preprocessing" in another host halo, typically of M_vir ~ 10^{10-12} M_sun, before falling into the MW/M31 halos. Satellites with lower mass and/or those closer to the MW/M31 fell in earlier and are more likely to have experienced group preprocessing; half of all satellites with M_star < 10^6 M_sun were preprocessed in a group. Infalling groups also drive most satellite-satellite mergers within the MW/M31 halos. Finally, none of the surviving satellites at z = 0 were within the virial radius of their MW/M31 halo during reionization (z > 6), and only < 4% were satellites of any other host halo during reionization. Thus, effects of cosmic reionization versus host-halo environment on the formation histories of surviving dwarf galaxies in the LG occurred at distinct epochs, separated typically by 2 - 4 Gyr, so they are separable theoretically and, in principle, observationally.


title: WIMP and SIMP Dark Matter from the Spontaneous Breaking of a Global Group

abstract: We propose and study a scalar extension of the Standard Model which respects a $\mathbb{Z}_3$ symmetry remnant of the spontaneous breaking of a global $U(1)_\text{DM}$ symmetry. Consequently, this model has a natural dark matter candidate and a Goldstone boson in the physical spectrum. In addition, the Higgs boson properties are changed with respect to the Standard Model due to the mixing with a new particle. We explore regions in the parameter space taking into account bounds from the measured Higgs properties, dark matter direct detection as well as measurements of the effective number of neutrino species before recombination. The dark matter relic density is determined by three classes of processes: the usual self-annihilation, semi-annihilation and purely dark matter $3 \to 2$ processes. The latter has been subject of recent interest leading to the so-called `Strongly Interacting Massive Particle' (SIMP) scenario. We show under which conditions our model can lead to a concrete realization of such scenario and study the possibility that the dark matter self-interactions could address the small scale structure problems. In particular, we find that in order for the SIMP scenario to work, the dark matter mass must be in the range $7-115$ MeV, with the global symmetry energy breaking scale in the TeV range.


title: An emerging population of BL Lacs with extreme properties: towards a class of EBL and cosmic magnetic field probes?

abstract: High energy observations of extreme BL Lac objects, such as 1ES 0229+200 or 1ES 0347-121, recently focused interest both for blazar and jet physics and for the implication on the extragalactic background light and intergalactic magnetic field estimate. However, the number of these extreme highly peaked BL Lac objects (EHBL) is still rather small. Aiming at increase their number, we selected a group of EHBL candidates starting from the BL Lac sample of Plotkin et al. (2011), considering those undetected (or only barely detected) by the Large Area Telescope onboard Fermi and characterized by a high X-ray vs. radio flux ratio. We assembled the multi-wavelength spectral energy distribution of the resulting 9 sources, profiting of publicly available archival observations performed by the Swift, Galex and Fermi satellites, confirming their nature. Through a simple one-zone synchrotron self-Compton model we estimate the expected VHE flux, finding that in the majority of cases it is within the reach of present generation of Cherenkov arrays or of the forthcoming Cherenkov Telescope Array (CTA).


title: From Veneziano to Riemann: A String Theory Statement of the Riemann Hypothesis

abstract: We discuss a precise relation between the Veneziano amplitude of string theory, rewritten in terms of ratios of the Riemann zeta function, and two elementary criteria for the Riemann hypothesis formulated in terms of integrals of the logarithm and the argument of the zeta function. We also discuss how the integral criterion based on the argument of the Riemann zeta function relates to the Li criterion for the Riemann hypothesis. We provide a new generalization of this integral criterion. Finally, we comment on the physical interpretation of our recasting of the Riemann hypothesis in terms of the Veneziano amplitude.


title: The Efficiency of Stellar Reionization: Effects of Rotation, Metallicity, and Initial Mass Function

abstract: We compute the production rate of photons in the ionizing Lyman continua (LyC) of H I (lambda < 912 A), He I (lambda < 504 A), and He II (lambda < 228 A) using recent stellar evolutionary tracks coupled to a grid of non-LTE, line-blanketed (WM-basic) model atmospheres. The median LyC production efficiency is Q_LyC = (6+/-2)x10^60 LyC photons per Msun of star formation (range [3.1-9.4]x10^60) corresponding to a revised calibration of 10^{53.3+/-0.2} photons/s per Msun/yr. Efficiencies in the helium continua are Q_HeI ~ 10^60 photons/Msun and Q_HeII ~ 10^56 photons/Msun at solar metallicity and larger at low metallicity. The critical star formation rate needed to maintain reionization against recombinations at z = 7 is rho_SFR = (0.012 Msun/yr/Mpc^3) [(1+z)/8]^3 [(C_H /3) (0.2/ f_esc) for fiducial values of IGM clumping factor C_H = 3 and LyC escape fraction f_esc = 0.2. The boost in LyC production efficiency is an important ingredient, together with metallicity, C_H, and f_esc, in assessing whether IGM reionization was complete by z ~ 7. Monte-Carlo sampled spectra of coeval starbursts during the first 5 Myr have intrinsic flux ratios of F(1500)/F(900) = 0.4-0.5 and F(912^-)/F(912^+) = 0.4-0.7 in the far-UV (1500A), the LyC (900A), and at the Lyman edge (912A). These ratios can be used to calibrate the LyC escape fractions in starbursts.


title: The mass accretion rate of galaxy clusters: a measurable quantity

abstract: We explore the possibility of measuring the mass accretion rate (MAR) of galaxy clusters from their mass profiles beyond the virial radius $R_{200}$. We derive the accretion rate from the mass of a spherical shell whose inner radius is $2R_{200}$, whose thickness changes with redshift, and whose infall velocity is assumed to be equal to the mean infall velocity of the spherical shells of dark matter halos extracted from $N$-body simulations. This approximation is rather crude in hierarchical clustering scenarios where both smooth accretion and aggregation of smaller dark matter halos contribute to the mass accretion of clusters.Nevertheless, in the redshift range $z=[0,2]$, our prescription returns an average MAR within $20-40 \%$ of the average rate derived from the merger trees of dark matter halos extracted from $N$-body simulations. The MAR of galaxy clusters has been the topic of numerous detailed numerical and theoretical investigations, but so far it has remained inaccessible to measurements in the real universe. Since the measurement of the mass profile of clusters beyond their virial radius can be performed with the caustic technique applied to dense redshift surveys of the cluster outer regions, our result suggests that measuring the mean MAR of a sample of galaxy clusters is actually feasible. We thus provide a new potential observational test of the cosmological and structure formation models.


title: Swift reveals a ~5.7 day super-orbital period in the M31 globular cluster X-ray binary XB158

abstract: The M31 globular cluster X-ray binary XB158 (a.k.a. Bo 158) exhibits intensity dips on a 2.78 hr period in some observations, but not others. The short period suggests a low mass ratio, and an asymmetric, precessing disk due to additional tidal torques from the donor star since the disk crosses the 3:1 resonance. Previous theoretical 3D smoothed particle hydrodynamical modeling suggested a super-orbital disk precession period 29$\pm$1 times the orbital period, i.e. $\sim$81$\pm$3 hr. We conducted a Swift monitoring campaign of 30 observations over ~1 month in order to search for evidence of such a super-orbital period. Fitting the 0.3--10 keV Swift XRT luminosity lightcurve with a sinusoid yielded a period of 5.65+/-0.05 days, and a >5$\sigma$ improvement in $\chi^2$ over the best fit constant intensity model. A Lomb-Scargle periodogram revealed that periods 5.4--5.8 days were detected at a >3$\sigma$ level, with a peak at 5.6 days. We consider this strong evidence for a 5.65 day super-orbital period, ~70\% longer than the predicted period. The 0.3--10 keV luminosity varied by a factor ~5, consistent with variations seen in long-term monitoring from Chandra. We conclude that other X-ray binaries exhibiting similar long-term behaviour are likely to also be X-ray binaries with low mass ratios and super-orbital periods.


title: Gap formation and stability in non-isothermal protoplanetary discs

abstract: Several observations of transition discs show lopsided dust-distributions. A potential explanation is the formation of a large-scale vortex acting as a dust-trap at the edge of a gap opened by a giant planet. Numerical models of gap-edge vortices have thus far employed locally isothermal discs, but the theory of this vortex-forming or `Rossby wave' instability was originally developed for adiabatic discs. We generalise the study of planetary gap stability to non-isothermal discs using customised numerical simulations of disc-planet systems where the planet opens an unstable gap. We include in the energy equation a simple cooling function with cooling timescale $t_c=\beta\Omega_k^{-1}$, where $\Omega_k$ is the Keplerian frequency, and examine the effect of $\beta$ on the stability of gap edges and vortex lifetimes. We find increasing $\beta$ lowers the growth rate of non-axisymmetric perturbations, and the dominant azimuthal wavenumber $m$ decreases. We find a quasi-steady state consisting of one large-scale, over-dense vortex circulating the outer gap edge, typically lasting $O(10^3)$ orbits. Vortex lifetimes were found to generally increase with cooling times up to an optimal value, beyond which vortex lifetimes decrease. This non-monotonic dependence is qualitatively consistent with recent studies using strictly isothermal discs that vary the disc aspect ratio.


title: Spatially adaptive radiation-hydrodynamical simulations of galaxy formation during cosmological reionization

abstract: We present a suite of cosmological radiation-hydrodynamical simulations of the assembly of galaxies driving the reionization of the intergalactic medium (IGM) at z >~ 6. The simulations account for the hydrodynamical feedback from photoionization heating and the explosion of massive stars as supernovae (SNe). Our reference simulation, which was carried out in a box of size 25 comoving Mpc/h using 2 x 512^3 particles, produces a reasonable reionization history and matches the observed UV luminosity function of galaxies. Simulations with different box sizes and resolutions are used to investigate numerical convergence, and simulations in which either SNe or photoionization heating or both are turned off, are used to investigate the role of feedback from star formation. Ionizing radiation is treated using accurate radiative transfer at the high spatially adaptive resolution at which the hydrodynamics is carried out. SN feedback strongly reduces the star formation rates (SFRs) over nearly the full mass range of simulated galaxies and is required to yield SFRs in agreement with observations. Photoheating helps to suppress star formation in low-mass galaxies, but its impact on the cosmic SFR is small. Because the effect of photoheating is masked by the strong SN feedback, it does not imprint a signature on the UV galaxy luminosity function, although we note that our resolution is insufficient to model star-forming minihaloes cooling through molecular hydrogen transitions. Photoheating does provide a strong positive feedback on reionization because it smooths density fluctuations in the IGM, which lowers the IGM recombination rate substantially. Our simulations demonstrate a tight non-linear coupling of galaxy formation and reionization, motivating the need for the accurate and simultaneous inclusion of photoheating and SN feedback in models of the early Universe.


title: Electromagnetic and neutral-weak response functions of 4He and 12C

abstract: Ab initio calculations of the quasi-elastic electromagnetic and neutral-weak response functions of 4He and 12C are carried out for the first time. They are based on a realistic approach to nuclear dynamics, in which the strong interactions are described by two- and three-nucleon potentials and the electroweak interactions with external fields include one- and two-body terms. The Green's function Monte Carlo method is used to calculate directly the Laplace transforms of the response functions, and maximum-entropy techniques are employed to invert the resulting imaginary-time correlation functions with associated statistical errors. The theoretical results, confirmed by experiment in the electromagnetic case, show that two-body currents generate excess transverse strength from threshold to the quasi-elastic to the dip region and beyond. These findings challenge the conventional picture of quasi-elastic inclusive scattering as being largely dominated by single-nucleon knockout processes.


title: Cosmic X-ray Surveys of Distant Active Galaxies: The Demographics, Physics, and Ecology of Growing Supermassive Black Holes

abstract: We review results from cosmic X-ray surveys of active galactic nuclei (AGNs) over the past ~ 15 yr that have dramatically improved our understanding of growing supermassive black holes in the distant universe. First, we discuss the utility of such surveys for AGN investigations and the capabilities of the missions making these surveys, emphasizing Chandra, XMM-Newton, and NuSTAR. Second, we briefly describe the main cosmic X-ray surveys, the essential roles of complementary multiwavelength data, and how AGNs are selected from these surveys. We then review key results from these surveys on the AGN population and its evolution ("demographics"), the physical processes operating in AGNs ("physics"), and the interactions between AGNs and their environments ("ecology"). We conclude by describing some significant unresolved questions and prospects for advancing the field.


title: Neutrino-nucleus interactions and the short-range structure of nuclei

abstract: Improvements in theoretical modeling of Short Range structures and phenomena, and comparisons with data, will require sustained collaboration between nuclear theorists and neutrino experimentalists. The extensive history of studying this area of nuclear physics in electron- and hadron-scattering experiments, coupled with the transformative capabilities of LArTPCs to identify neutrinos, will provide a ripe opportunity for new discoveries that will further our understanding of the nucleus.