title: Communicating Correlated Sources Over an Interference Channel

abstract: A new coding technique, based on \textit{fixed block-length} codes, is proposed for the problem of communicating a pair of correlated sources over a $2-$user interference channel. Its performance is analyzed to derive a new set of sufficient conditions. The latter is proven to be strictly less binding than the current known best, which is due to Liu and Chen [Dec, 2011]. Our findings are inspired by Dueck's example [March, 1981].


title: Impacts of nuclear-physics uncertainty in stellar temperatures on the s-process nucleosynthesis

abstract: We evaluated the uncertainty relevant to s-process nucleosynthesis using a Monte-Carlo centred approach. We are based on a realistic and general prescription of temperature dependent uncertainty for the reactions. We considered massive stars for the weak s-process and AGB stars for the main s-process. We found that the adopted uncertainty for (n,$\gamma$) rates, tens of per cent on average, affect the production of s-process nuclei along the $\beta$-stability line, while for $\beta$-decay, for which contributions from excited states enhances the uncertainty, has the strongest impact on branching points.


title: Recent results from NA48/2 and NA62 experiments at CERN

abstract: The NA48/2 and $\text{NA62}\text{-}{\text{R}_\text{K}}$ experiments at the CERN SPS collected a large sample of charged kaon decays in flight. $\text{NA62}\text{-}{\text{R}_\text{K}}$ was running in 2007-08 with a highly efficient minimum bias trigger for decays into electrons. A preliminary measurement of the electromagnetic transition form factor slope of the $\pi^0$ from $1.05\times 10^{6}$ fully reconstructed $\pi^0$ Dalitz decays is presented. The obtained value $a = (3.70 \pm 0.53_\text{stat} \pm 0.36_\text{syst})\times 10^{-2}$ represents a $5.8\sigma$ observation of a non-zero slope in the time-like region of momentum transfer. An upper limit on the rate of a lepton number violating decay $K^\pm\to\pi^\mp\mu^\pm\mu^\pm$ is reported from $\sim 1.6\times 10^{11} K^\pm$ decays at NA48/2 in 2003-04: $\mathcal{B} < 8.6\times 10^{-11}$ at $90 \%$ CL. Searches for heavy sterile neutrino $N_4$ and inflaton $\chi$ resonances in $K^\pm\to\pi\mu\mu$ decays are reported. No signal is observed and upper limits on the products $\mathcal{B}(K^\pm\to\mu^\pm N_4)\mathcal{B}(N_4\to\pi^\mp\mu^\pm)$ and $\mathcal{B}(K^\pm\to\pi^\pm \chi)\mathcal{B}(\chi\to\mu^+\mu^-)$ are set in the range $10^{-10}-10^{-9}$ for resonance lifetimes up to $100~\text{ps}$. The result of a search for dark photon with the same sample of decays is also reported. In the absence of observed signal, the limits on the mixing parameter $\varepsilon^2$ in the range $9-70~\text{MeV/c}^{2}$ are improved.


title: Tunable plasmon-enhanced birefringence in ribbon array of anisotropic 2D materials

abstract: We explore the far-field scattering properties of anisotropic 2D materials in ribbon array configuration. Our study reveals the plasmon-enhanced linear birefringence in these ultrathin metasurfaces, where linearly polarized incident light can be scattered into its orthogonal polarization or be converted into circular polarized light. We found wide modulation in both amplitude and phase of the scattered light via tuning the operating frequency or material's anisotropy and develop models to explain the observed scattering behavior.


title: Multi-Layer Generalized Linear Estimation

abstract: We consider the problem of reconstructing a signal from multi-layered (possibly) non-linear measurements. Using non-rigorous but standard methods from statistical physics we present the Multi-Layer Approximate Message Passing (ML-AMP) algorithm for computing marginal probabilities of the corresponding estimation problem and derive the associated state evolution equations to analyze its performance. We also give the expression of the asymptotic free energy and the minimal information-theoretically achievable reconstruction error. Finally, we present some applications of this measurement model for compressed sensing and perceptron learning with structured matrices/patterns, and for a simple model of estimation of latent variables in an auto-encoder.

categories: cs.IT cond-mat.stat-mech math.IT stat.ML


title: Reissner-Nordstr{\o}m-de Sitter Manifold : Photon Sphere and Maximal Analytic Extension

abstract: This paper is devoted to the study of the Reissner-Nordstr{\o}m-de Sitter black holes and their maximal analytic extensions. In particular, we study some of their properties that lays the groundwork for separate papers where we obtain decay results and construct conformal scattering theories for test fields on such spacetimes. Here, we find the necessary and sufficient conditions on the parameters of the Reissner-Nordstr{\o}m-de Sitter metric -namely, the mass, the charge, and the cosmological constant- to have three horizons. Under this conditions, we prove that there is only one photon sphere and we locate it. We then give a detailed construction of the maximal analytic extension of the Reissner-Nordstr{\o}m-de Sitter manifold in the case of three horizons.


title: Influence of external magnetic field, finite-size effects and chemical potential on the phase transition of a complex scalar field

abstract: A scalar model is built, as a quantum field theory defined on toroidal topologies, to describe phase transition in films subjected to periodic boundary conditions and influenced by an external and constant magnetic field. Criticality is studied and the relations between the critical temperature, the film thickness, the magnetic field strength and the chemical potential are investigated. Since the model describes a second-order phase transition a comparison with the Ginzburg-Landau theory is made.


title: Bielliptic curves of genus three and the Torelli problem for certain elliptic surfaces

abstract: We study the Hodge structure of elliptic surfaces which are canonically defined from bielliptic curves of genus three. We prove that the period map for the second cohomology has one dimensional fibers, and the period map for the total cohomology is of degree twelve, and moreover, by adding the information of the Hodge structure of the canonical divisor, we prove a generic Torelli theorem for these elliptic surfaces. Finally, we give explicit examples of the pair of non-isomorphic elliptic surfaces which have the same Hodge structure on themselves and the same Hodge structure on their canonical divisors.


title: Fine-Grained Parameterized Complexity Analysis of Graph Coloring Problems

abstract: The $q$-Coloring problem asks whether the vertices of a graph can be properly colored with $q$ colors. Lokshtanov et al. [SODA 2011] showed that $q$-Coloring on graphs with a feedback vertex set of size $k$ cannot be solved in time $\mathcal{O}^*((q-\varepsilon)^k)$, for any $\varepsilon > 0$, unless the Strong Exponential-Time Hypothesis (SETH) fails. In this paper we perform a fine-grained analysis of the complexity of $q$-Coloring with respect to a hierarchy of parameters. We show that even when parameterized by the vertex cover number, $q$ must appear in the base of the exponent: Unless ETH fails, there is no universal constant $\theta$ such that $q$-Coloring parameterized by vertex cover can be solved in time $\mathcal{O}^*(\theta^k)$ for all fixed $q$. We apply a method due to Jansen and Kratsch [Inform. & Comput. 2013] to prove that there are $\mathcal{O}^*((q - \varepsilon)^k)$ time algorithms where $k$ is the vertex deletion distance to several graph classes $\mathcal{F}$ for which $q$-Coloring is known to be solvable in polynomial time. We generalize earlier ad-hoc results by showing that if $\mathcal{F}$ is a class of graphs whose $(q+1)$-colorable members have bounded treedepth, then there exists some $\varepsilon > 0$ such that $q$-Coloring can be solved in time $\mathcal{O}^*((q-\varepsilon)^k)$ when parameterized by the size of a given modulator to $\mathcal{F}$. In contrast, we prove that if $\mathcal{F}$ is the class of paths - some of the simplest graphs of unbounded treedepth - then no such algorithm can exist unless SETH fails.


title: Quantum-trajectory analysis for charge transfer in solid materials induced by strong laser fields

abstract: We investigate the dependence of charge transfer on the intensity of driving laser field when SiO2 crystal is irradiated by an 800 nm laser. It is surprising that the direction of charge transfer undergoes a sudden reversal when the driving laser intensity exceeds critical values with different carrier envelope phases. By applying quantum-trajectory analysis, we find that the Bloch oscillation plays an important role in charge transfer in solid. Also, we study the interaction of strong laser with gallium nitride (GaN) that is widely used in optoelectronics. A pump-probe scheme is applied to control the quantum trajectories of the electrons in the conduction band. The signal of charge transfer is controlled successfully by means of theoretically proposed approach.


title: The configuration category of a product

abstract: A construction related to the Boardman-Vogt tensor product of operads allows us to describe the configuration category of a product manifold $M\times N$ in terms of the configuration categories of the factors $M$ and $N$.


title: Astrophysical interpretation of the anisotropies in the unresolved gamma-ray background

abstract: Recently, a new measurement of the auto- and cross-correlation angular power spectrum (APS) of the isotropic gamma-ray background was performed, based on 81 months of data of the Fermi Large-Area Telescope (LAT). Here, we fit, for the first time, the new APS data with a model describing the emission of unresolved blazars. These sources are expected to dominate the anisotropy signal. The model we employ in our analysis reproduces well the blazars resolved by Fermi LAT. When considering the APS obtained by masking the sources in the 3FGL catalogue, we find that unresolved blazars under-produce the measured APS below $\sim$1 GeV. Contrary to past results, this suggests the presence of a new contribution to the low-energy APS, with a significance of, at least, 5$\sigma$. The excess can be ascribed to a new class of faint gamma-ray emitters. If we consider the APS obtained by masking the sources in the 2FGL catalogue, there is no under-production of the APS below 1 GeV, but the new source class is still preferred over the blazars-only scenario (with a significance larger than 10$\sigma$). The properties of the new source class and the level of anisotropies induced in the isotropic gamma-ray background are the same, independent of the APS data used. In particular, the new gamma-ray emitters must have a soft energy spectrum, with a spectral index ranging, approximately, from 2.7 to 3.2. This complicates their interpretation in terms of known sources, since, normally, star-forming and radio galaxies are observed with a harder spectrum. The new source class identified here is also expected to contribute significantly to the intensity of the isotropic gamma-ray background.


title: Non-Adaptive Group Testing Framework based on Concatenation Code

abstract: We consider an efficiently decodable non-adaptive group testing (NAGT) problem that meets theoretical bounds. The problem is to find a few specific items (at most $d$) satisfying certain characteristics in a colossal number of $N$ items as quickly as possible. Those $d$ specific items are called \textit{defective items}. The idea of NAGT is to pool a group of items, which is called \textit{a test}, then run a test on them. If the test outcome is \textit{positive}, there exists at least one defective item in the test, and if it is \textit{negative}, there exists no defective items. Formally, a binary $t \times N$ measurement matrix $\mathcal{M} = (m_{ij})$ is the representation for $t$ tests where row $i$ stands for test $i$ and $m_{ij} = 1$ if and only if item $j$ belongs to test $i$. There are three main objectives in NAGT: minimize the number of tests $t$, construct matrix $\mathcal{M}$, and identify defective items as quickly as possible. In this paper, we present a strongly explicit construction of $\mathcal{M}$ for when the number of defective items is at most 2, with the number of tests $t \simeq 16 \log{N} = O(\log{N})$. In particular, we need only $K \simeq N \times 16\log{N} = O(N\log{N})$ bits to construct such matrices, which is optimal. Furthermore, given these $K$ bits, any entry in the matrix can be constructed in time $O \left(\ln{N}/ \ln{\ln{N}} \right)$. Moreover, $\mathcal{M}$ can be decoded with high probability in time $O\left( \frac{\ln^2{N}}{\ln^2{\ln{N}}} \right)$. When the number of defective items is greater than 2, we present a scheme that can identify at least $(1-\epsilon)d$ defective items with $t \simeq 32 C(\epsilon) d \log{N} = O(d \log{N})$ in time $O \left( d \frac{\ln^2{N}}{\ln^2{\ln{N}}} \right)$ for any close-to-zero $\epsilon$, where $C(\epsilon)$ is a constant that depends only on $\epsilon$.


title: Pumping of magnons in a Dzyaloshinskii-Moriya ferromagnet

abstract: We formulate a microscopic linear response theory of nonequilibrium magnonic torques and magnon pumping applicable to multiple-magnonic-band uniform ferromagnets with Dzyaloshinskii-Moriya interactions. From the linear response theory, we identify the extrinsic and intrinsic contributions where the latter is expressed via the Berry curvature of magnonic bands. We observe that in the presence of a time-dependent magnetization Dzyaloshinskii-Moriya interactions can act as fictitious electric fields acting on magnons. We study various current responses to this fictitious field and analyze the role of Berry curvature. After identifying the magnon-mediated contribution to the equilibrium Dzyaloshinskii-Moriya interaction, we also establish the Onsager reciprocity between the magnon-mediated torques and heat pumping. We apply our theory to the magnonic heat pumping and torque responses in honeycomb and kagome lattice ferromagnets.


title: Optimal Context Aware Transmission Strategy for non-Orthogonal D2D Communications

abstract: The increasing traffic demand in cellular networks has recently led to the investigation of new strategies to save precious resources like spectrum and energy. A possible solution employs direct device-to-device (D2D) communications, which is particularly promising when the two terminals involved in the communications are located in close proximity. The D2D communications should coexist with other transmissions, so they must be careful scheduled in order to avoid harmful interference impacts. In this paper, we analyze how a distributed context-awareness, obtained by observing few local channel and topology parameters, can be used to adaptively exploit D2D communications. We develop a rigorous theoretical analysis to quantify the balance between the gain offered by a D2D transmission, and its impact on the other network communications. Based on this analysis, we derive two theorems that define the optimal strategy to be employed, in terms of throughput maximization, when a single or multiple transmit power levels are available for the D2D communications. We compare this strategy to the state-of-the-art in the same network scenario, showing how context awareness can be exploited to achieve a higher sum throughput and an improved fairness.


title: Reconstructing CMB fluctuations and the mean reionization optical depth

abstract: The Thomson optical depth from reionization is a limiting factor in measuring the amplitude of primordial fluctuations, and hence in measuring physics that affects the low-redshift amplitude, such as the neutrino masses. Current constraints on the optical depth, based on directly measuring large-scale cosmic microwave background (CMB) polarization, are challenging due to foregrounds and systematic effects. Here, we consider an indirect measurement of large-scale polarization, using observed maps of small-scale polarization together with maps of fields that distort the CMB, such as CMB lensing and patchy reionization. We find that very futuristic CMB surveys will be able to reconstruct large-scale polarization, and thus the mean optical depth, using only measurements on small scales.


title: Trends in U. S. Storminess 1949--2009

abstract: We use an extensive NOAA database of hourly precipitation data from 5995 stations in the 48 contiguous United States over the period 1949--2009 to investigate possible trends in the frequency and severity of extreme weather events, defined as periods of intense precipitation. The frequency and intensity of these events are quantified by a dimensionless storminess, defined as the variance of the hourly rainfall at a site normalized by the square of the mean rainfall at that site. For 1722 stations with sufficient data, we compute the rate of change of the logarithm of the storminess at each station and set bounds on its mean (over stations) trend; use of the logarithms weights trends at calm stations equally to those at stormy stations and enhances the statistical power of the mean. These results are confirmed by reversing the order of averaging: first computing, for each year, the mean (over stations) logarithm of the storminess, and then fitting to its trend. We set $2\sigma$ upper bounds of 0.001/y (doubling or halving time scales $> 1000$ y) on any trend, increasing or decreasing.


title: $\mathbb{C}P^N$ sigma models via the $SU(2)$ coherent states approach

abstract: In this paper we present results obtained from the unification of $SU(2)$ coherent states with $\mathbb{C}P^N$ sigma models defined on the Riemann sphere having finite actions. The set of coherent states generated by a vector belonging to a carrier space of an irreducible representation of the group gives rise to a map from the sphere into the set of rank-1 Hermitian projectors in that space. The map can be identified with a particular solution of the $\mathbb{C}P^N$ sigma model, where $N+1$ is equal to the dimension of the representation space. In particular a choice of the generating vector as the highest weight vector of the representation gives rise to the map known as a Veronese immersion. Using a description of the matrix elements of these representations in terms of Jacobi polynomials, we obtain an explicit parametrization of the solutions of the $\mathbb{C}P^N$ models, which has not been previously found. We relate the analytical properties of the solutions, which are known to belong to separate classes --- holomorphic, anti-holomorphic and various types of mixed ones --- to the weight corresponding to the chosen weight vector. Some examples of the described constructions are elaborated in detail in this paper.


title: Modelling the Milky Way's globular cluster system

abstract: We construct a model for the Galactic globular cluster system based on a realistic gravitational potential and a distribution function (DF) analytic in the action integrals. The DF comprises disc and halo components whose functional forms resemble those recently used to describe the stellar discs and stellar halo. We determine the posterior distribution of our model parameters using a Bayesian approach. This gives us an understanding of how well the globular cluster data constrain our model. The favoured parameter values of the disc and halo DFs are similar to values previously obtained from fits to the stellar disc and halo, although the cluster halo system shows clearer rotation than does the stellar halo. Our model reproduces the generic features of the globular cluster system, namely the density profile, the mean rotation velocity. The fraction of disc clusters coincides with the observed fraction of metal-rich clusters. However, the data indicate either incompatibility between catalogued cluster distances and current estimates of distance to the Galactic Centre, or failure to identify clusters behind the bulge. As the data for our Galaxy's components increase in volume and precision over the next few years, it will be rewarding to revisit the present analysis.


title: A non-linear theory of infrahyperfunctions

abstract: We develop a nonlinear theory for infrahyperfunctions (also referred to as quasianalytic (ultra)distributions by L. H\"{o}rmander). In the hyperfunction case our work can be summarized as follows. We construct a differential algebra that contains the space of hyperfunctions as a linear differential subspace and in which the multiplication of real analytic functions coincides with their ordinary product. Moreover, by proving an analogue of Schwartz's impossibility result for hyperfunctions, we show that this embedding is optimal. Our results fully solve an earlier question raised by M. Oberguggenberger.


title: Notes on Spinning Operators in Fermionic CFT

abstract: The Gross-Neveu model defines a unitary CFT of interacting fermions in $2<d<4$ which has perturbative descriptions in the $1/N$ expansion and in the epsilon-expansion near two and four dimensions. In each of these descriptions, the CFT has an infinite tower of nearly conserved currents of all spins. We determine the structure of the non-conservation equations both at large $N$ and in the epsilon-expansion, and use it to find the leading order anomalous dimensions of the broken currents. Similarly, we use the fact that the CFT spectrum includes a nearly free fermion to fix the leading anomalous dimensions of a few scalar composite operators. We also compute the scaling dimensions of double-trace spinning operators in the large $N$ expansion, which correspond to interaction energies of two-particle states in the AdS dual higher-spin theory. We first derive these anomalous dimensions by a direct Feynman diagram calculation, and then show that the result can be exactly reproduced by analytic bootstrap methods, provided the sum over the tower of weakly broken higher-spin currents is suitably regularized. Finally, we apply the analytic bootstrap approach to derive the anomalous dimensions of the double-trace spinning operators in the 3d bosonic and fermion vector models coupled to Chern-Simons theory, to leading order in $1/N$ but exactly in the `t Hooft coupling.


title: Wetting in the presence of the electric field: The classical density functional theory study for a model system

abstract: We discuss the effect of an external electric field on the wetting of a solid surface by liquid. To this end, we use a model of the two-level-atom fluid for which the changes in interatomic interactions due to the presence of the field can be found using quantum-mechanical perturbation theory. Constructing the grand potential functional, we perform the standard calculations of Young's equilibrium contact angle. The switching on of the electric field $|{\bf {E}}| > 0$ may increase noticeably the contact angle $\theta$.


title: Low temperature features in the heat capacity of unary metals and intermetallics for the example of bulk aluminum and Al$_3$Sc

abstract: We explore the competition and coupling of vibrational and electronic contributions to the heat capacity of Al and Al$_3$Sc at temperatures below 50 K combining experimental calorimetry with highly converged finite temperature density functional theory calculations. We find that semilocal exchange correlation functionals accurately describe the rich feature set observed for these temperatures, including electron-phonon coupling. Using different representations of the heat capacity, we are therefore able to identify and explain deviations from the Debye behaviour in the low-temperature limit and in the temperature regime 30 - 50 K as well as the reduction of these features due to the addition of Sc.


title: Strain engineering of Schottky barriers in single- and few-layer MoS2 vertical devices

abstract: We study the effect of local strain in the electronic transport properties of vertical metal-atomically thin MoS2-metal structures. We use a conductive atomic force microscope tip to apply different load forces to monolayer and few-layer MoS2 crystals deposited onto a conductive indium tin oxide (ITO) substrate while measuring simultaneously the I-V characteristics of the vertical tip/MoS2/ITO structures. The structures show rectifying I-V characteristics, with rectification ratios strongly dependent on the applied load. To understand these results, we compare the experimental I-Vs with a double Schottky barrier model, which is in good agreement with our experimental results and allows us to extract quantitative information about the electronic properties of the tip/MoS2/ITO structures and their dependence on the applied load. Finally, we test the stability of the studied structures using them as mechanically tunable current rectifiers.


title: Shock-Cloud Interaction in the Solar Corona

abstract: Flare associated coronal shock waves sometimes interact with solar prominences leading to large amplitude prominence oscillations. Such prominence activation gives us unique opportunity to track time evolution of shock-cloud interaction in cosmic plasmas. Although the dynamics of interstellar shock-cloud interaction is extensively studied, coronal shock-solar prominence interaction is rarely studied in the context of shock-cloud interaction. Associated with X5.4 class solar flare occurred on 7 March, 2012, a globally propagated coronal shock wave interacted with a polar prominence leading to large amplitude prominence oscillation. In this paper, we studied bulk acceleration and excitation of internal flow of the shocked prominence using three-dimensional MHD simulations. We studied eight magnetohydrodynamic (MHD) simulation runs with different mass density structure of the prominence, and one hydrodynamic simulation run, and compared the result. In order to compare observed motion of activated prominence with corresponding simulation, we also studied prominence activation by injection of triangular shaped coronal shock. We found that magnetic tension force mainly accelerate (and then decelerate) the prominence. The internal flow, on the other hand, is excited during the shock front sweeps through the the prominence and damps almost exponentially. We construct phenomenological model of bulk momentum transfer from shock to the prominence, which agreed quantitatively with all the simulation results. Based on the phenomenological prominence-activation model, we diagnosed physical parameters of coronal shock wave. The estimated energy of the coronal shock is several percent of total energy released during the X5.4 flare.


title: HESS J1826$-$130: A Very Hard $\gamma$-Ray Spectrum Source in the Galactic Plane

abstract: HESS J1826$-$130 is an unidentified hard spectrum source discovered by H.E.S.S. along the Galactic plane, the spectral index being $\Gamma$ = 1.6 with an exponential cut-off at about 12 TeV. While the source does not have a clear counterpart at longer wavelengths, the very hard spectrum emission at TeV energies implies that electrons or protons accelerated up to several hundreds of TeV are responsible for the emission. In the hadronic case, the VHE emission can be produced by runaway cosmic-rays colliding with the dense molecular clouds spatially coincident with the H.E.S.S. source.


title: Black Hole Continuum Spectra as a Test of General Relativity: Quadratic Gravity

abstract: Observations of the continuum spectrum emitted by accretion disks around black holes allows us to infer their properties, including possibly whether black holes are described by the Kerr metric. Some modified gravity theories do not admit the Kerr metric as a solution, and thus, continuum spectrum observations could be used to constrain these theories. We here investigate whether current and next generation X-Ray observations of the black hole continuum spectrum can constrain such deviations from Einstein's theory, focusing on two well-motivated modified quadratic gravity theories: dynamical Chern-Simons gravity and Einstein-dilaton-Gauss-Bonnet gravity. We do so by determining whether the non-Kerr deviations in the continuum spectrum introduced by these theories are larger than the observational error intrinsic to the observations. We find that dynamical Chern-Simons gravity cannot be constrained better than current bounds with current or next generation continuum spectrum observations. Einstein-dilaton-Gauss-Bonnet gravity, however, may be constrained better than current bounds with next generation telescopes, as long as the systematic error inherent in the accretion disk modeling is decreased below the predicted observational error.


title: Tunneling of the hard-core model on finite triangular lattices

abstract: We consider the hard-core model on finite triangular lattices with Metropolis dynamics. Under suitable conditions on the triangular lattice dimensions, this interacting particle system has three maximum-occupancy configurations and we investigate its high-fugacity behavior by studying tunneling times, i.e., the first hitting times between between these maximum-occupancy configurations, and the mixing time. The proof method relies on the analysis of the corresponding state space using geometrical and combinatorial properties of the hard-core configurations on finite triangular lattices, in combination with known results for first hitting times of Metropolis Markov chains in the equivalent zero-temperature limit. In particular, we show how the order of magnitude of the expected tunneling times depends on the triangular lattice dimensions in the low-temperature regime and prove the asymptotic exponentiality of the rescaled tunneling time leveraging the intrinsic symmetry of the state space.


title: A Halo Model Approach to the $\rm{21\, cm}$ and $\rm{Ly}\alpha$ Cross-correlation

abstract: We present a halo-model-based approach to calculate the cross-correlation between $\rm{21\,cm}$ HI intensity fluctuations and $\rm{Ly}\alpha$ emitters (LAE) during the epoch of reionization (EoR). Ionizing radiation around dark matter halos are modeled as bubbles with the size and growth determined based on the reionization photon production, among other physical parameters. The cross-correlation shows a clear negative-to-positive transition, associated with transition from ionized to neutral hydrogen in the intergalactic medium during EoR. The cross-correlation is subject to several foreground contaminants, including foreground radio point sources important for $\rm{21\,cm}$ experiments and low-$z$ interloper emission lines, such as $\rm{H}\alpha$, OIII, and OII, for $\rm{Ly}\alpha$ experiments. Our calculations show that by masking out high fluxes in the $\rm{Ly}\alpha$ measurement, the correlated foreground contamination on the $\rm{21\,cm}$-$\rm{Ly}\alpha$ cross-correlation can be dramatically reduced. We forecast the detectability of $\rm{21\,cm}$-$\rm{Ly}\alpha$ cross-correlation at different redshifts and adopt a Fisher matrix approach to estimate uncertainties on the key EoR parameters that have not been well constrained by other observations of reionization. This halo-model-based approach enables us to explore the EoR parameter space rapidly for different $\rm{21\,cm}$ and $\rm{Ly}\alpha$ experiments.


title: Understanding EROS2 observations toward the spiral arms within a classical Galactic model framework

abstract: EROS has searched for microlensing toward four directions in the Galactic plane away from the Galactic center. The interpretation of the catalog optical depth is complicated by the spread of the source distance distribution. We compare the EROS microlensing observations with Galactic models, tuned to fit the EROS source catalogs, and take into account all observational data such as the microlensing optical depth, the Einstein crossing durations, and the color and magnitude distributions of the catalogued stars. We simulated EROS-like source catalogs using the Hipparcos database, the Galactic mass distribution, and an interstellar extinction table. Taking into account the EROS star detection efficiency, we were able to produce simulated color-magnitude diagrams that fit the observed diagrams. This allows us to estimate average microlensing optical depths and event durations that are directly comparable with the measured values. Both the Besancon model and our Galactic model allow us to fully understand the EROS color-magnitude data. The average optical depths and mean event durations calculated from these models are in reasonable agreement with the observations; consequently our simulation allows a better understanding of the lens and source spatial distributions in the microlensing events. Varying the Galactic structure parameters through simulation, we were also able to deduce contraints on the kinematics of the disk, the disk stellar mass function (IMF at a few kpc distance from the Sun), and the maximum contribution of a thick disk of compact objects in the Galactic plane (Mthick < 5 - 7 x 10**10 Msun at 95%, depending on the model). We also show that the microlensing data toward one of our monitored directions are significantly sensitive to the galactic bar parameters, although much larger statistics are needed to provide competitive constraints.


title: A New Wiretap Channel Model and its Strong Secrecy Capacity

abstract: In this paper, a new wiretap channel model is proposed, where the legitimate transmitter and receiver communicate over a discrete memoryless channel. The wiretapper has perfect access to a fixed-length subset of the transmitted codeword symbols of her choosing. Additionally, she observes the remainder of the transmitted symbols through a discrete memoryless channel. This new model subsumes the classical wiretap channel and wiretap channel II with noisy main channel as its special cases. The strong secrecy capacity of the proposed channel model is identified. Achievability is established by solving a dual secret key agreement problem in the source model, and converting the solution to the original channel model using probability distribution approximation arguments. In the dual problem, a source encoder and decoder, who observe random sequences independent and identically distributed according to the input and output distributions of the legitimate channel in the original problem, communicate a confidential key over a public error-free channel using a single forward transmission, in the presence of a compound wiretapping source who has perfect access to the public discussion. The security of the key is guaranteed for the exponentially many possibilities of the subset chosen at wiretapper by deriving a lemma which provides a doubly-exponential convergence rate for the probability that, for a fixed choice of the subset, the key is uniform and independent from the public discussion and the wiretapping source's observation. The converse is derived by using Sanov's theorem to upper bound the secrecy capacity of the new wiretap channel model by the secrecy capacity when the tapped subset is randomly chosen by nature.


title: Simplified dark matter models with a spin-2 mediator at the LHC

abstract: We consider simplified dark matter models where a dark matter candidate couples to the standard model (SM) particles via an $s$-channel spin-2 mediator, and study constraints on the model parameter space from the current LHC data. Our focus lies on the complementarity among different searches, in particular monojet and multijet plus missing energy searches and resonance searches. For universal couplings of the mediator to SM particles, missing-energy searches can give stronger constraints than $WW$, $ZZ$, dijet, dihiggs, $t\bar t$, $b\bar b$ resonance searches in the low-mass region and/or when the coupling of the mediator to dark matter is much larger than its couplings to SM particles. The strongest constraints however come from diphoton and dilepton resonance searches. Only if these modes are suppressed, missing-energy searches can be competitive in constraining dark matter models with a spin-2 mediator.


title: An involution on Dyck paths that preserves the rise composition and interchanges the number of returns and the position of the first double fall

abstract: Motivated by a recent paper of Adin, Bagno and Roichman, we present an involution on Dyck paths that preserves the rise composition and interchanges the number of returns and the position of the first double fall.


title: Infinite Mixtures of Infinite Factor Analysers

abstract: Factor-analytic Gaussian mixture models are often employed as a model-based approach to clustering high-dimensional data. Typically, the numbers of clusters and latent factors must be specified in advance of model fitting, and remain fixed. The pair which optimises some model selection criterion is then chosen. For computational reasons, models in which the number of latent factors differ across clusters are rarely considered. Here the infinite mixture of infinite factor analysers (IMIFA) model is introduced. IMIFA employs a Pitman-Yor process prior to facilitate automatic inference of the number of clusters using the stick-breaking construction and a slice sampler. Furthermore, IMIFA employs multiplicative gamma process shrinkage priors to allow cluster-specific numbers of factors, automatically inferred via an adaptive Gibbs sampler. IMIFA is presented as the flagship of a family of factor-analytic mixture models, providing flexible approaches to clustering high-dimensional data. Applications to a benchmark data set, metabolomic spectral data, and a manifold learning handwritten digit example illustrate the IMIFA model and its advantageous features. These include obviating the need for model selection criteria, reducing the computational burden associated with the search of the model space, improving clustering performance by allowing cluster-specific numbers of factors, and quantifying uncertainty in the numbers of clusters and cluster-specific factors.


title: Controlling false discoveries in Bayesian gene networks with lasso regression p-values

abstract: Bayesian networks can represent directed gene regulations and therefore are favored over co-expression networks. However, hardly any Bayesian network study concerns the false discovery control (FDC) of network edges, leading to low accuracies due to systematic biases from inconsistent false discovery levels in the same study. We design four empirical tests to examine the FDC of Bayesian networks from three p-value based lasso regression variable selections --- two existing and one we originate. Our method, lassopv, computes p-values for the critical regularization strength at which a predictor starts to contribute to lasso regression. Using null and Geuvadis datasets, we find that lassopv obtains optimal FDC in Bayesian gene networks, whilst existing methods have defective p-values. The FDC concept and tests extend to most network inference scenarios and will guide the design and improvement of new and existing methods. Our novel variable selection method with lasso regression also allows FDC on other datasets and questions, even beyond network inference and computational biology. Lassopv is implemented in R and freely available at https://github.com/lingfeiwang/lassopv and https://cran.r-project.org/package=lassopv


title: IR-safe and UV-safe integrands in the EFTofLSS with exact time dependence

abstract: Because large-scale structure surveys may very well be the next leading sources of cosmological information, it is important to have a precise understanding of the cosmological observables; for this reason, the Effective Field Theory of Large-Scale Structure (EFTofLSS) was developed. So far, most results in the EFTofLSS have used the so-called Einstein-de Sitter approximation, an approximation of the time dependence which is known to be accurate to better than one percent. However, in order to reach even higher accuracy, the full time dependence must be used. The computation with exact time dependence is sensitive to both infrared (IR) and ultraviolet (UV) effects in the loop integrands, and while these effects must cancel because of diffeomorphism invariance, they make numerical computation much less efficient. We provide a formulation of the one-loop, equal-time exact-time-dependence power spectrum of density perturbations which is manifestly free of these spurious IR and UV divergences at the level of the integrand. We extend our results to the total matter mode with clustering quintessence, show that IR and UV divergences cancel, and provide the associated IR- and UV-safe integrand. This also establishes that the consistency conditions are satisfied in this system. We then use our one-loop result to do an improved precision comparison of the two-loop dark-matter power spectrum with the Dark Sky N-body simulation.


title: Higher degree S-lemma and the stability of quadratic modules

abstract: In this work we will investigate a certain generalization of the so called S-lemma in higher degrees. The importance of this generalization is, that it is closely related to Hilbert's 1888 theorem about tenary quartics. In fact, if such a generalization exits, then one can state a Hilbert-like theorem, where positivity is only demanded on some semi-algebraic set. We will show that such a generalization is not possible, at least not without additional conditions. To prove this, we will use and generalize certain tools developed by Netzer ([Ne]). These new tools will allow us to conclude that this generalization of the S-lemma is not possible because of geometric reasons. Furthermore, we are able to establish a link between geometric reasons and algebraic reasons. This will be accomplished within the framework of quadratic modules.


title: Maximal Spectra of rings consisting of regulated functions

abstract: In 1969 H\"ochster proved that for every quasi-compact T1-space $X$ we can find a commutative ring $R$ such that $X$ is homeomorphic to the maximal spectrum $\mathrm{Specm}(R)$ of $R$. This result implies the existence of a commutative ring $R$ that admits a non-Hausdorff and totally disconnected maximal spectrum $\mathrm{Specm}(R)$. However, there has not been an example of such a commutative ring yet. The aim of this paper is to investigate rings that almost admit a non-Hausdorff and totally disconnected maximal spectrum.


title: Local analogues of high-redshift star-forming galaxies: integral field spectroscopy of green peas

abstract: We use integral field spectroscopy, from the SWIFT and Palm3K instruments, to perform a spatially-resolved spectroscopic analysis of four nearby highly star-forming `green pea' (GP) galaxies, that are likely analogues of star-forming systems at z~2.5-3. By studying emission-line maps in H$\alpha$, [NII]$\lambda \lambda$6548,6584 and [SII]$\lambda$$\lambda$6716,6731, we explore the kinematic morphology of these systems and constrain properties such as gas-phase metallicities, electron densities and gas-ionization mechanisms. Two of our GPs are rotationally-supported while the others are dispersion-dominated systems. The rotationally-supported galaxies both show evidence for recent or ongoing mergers. However, given that these systems have intact disks, these interactions are likely to have low mass ratios (i.e. minor mergers), suggesting that the minor-merger process may be partly responsible for the high SFRs seen in these GPs. Nevertheless, the fact that the other two GPs appear morphologically undisturbed suggests that mergers (including minor mergers) are not necessary for driving the high star formation rates in such galaxies. We show that the GPs are metal-poor systems (25-40 per cent of solar) and that the gas ionization is not driven by AGN in any of our systems, indicating that AGN activity is not co-eval with star formation in these starbursting galaxies.


title: Factors of sums involving $q$-binomial coefficients and powers of $q$-integers

abstract: We show that, for all positive integers $n_1, \ldots, n_m$, $n_{m+1}=n_1$, and any non-negative integers $j$ and $r$ with $j\leqslant m$, the expression $$ \frac{1}{[n_1]}{n_1+n_{m}\brack n_1}^{-1} \sum_{k=1}^{n_1}[2k][k]^{2r}q^{jk^2-(r+1)k}\prod_{i=1}^{m} {n_i+n_{i+1}\brack n_i+k} $$ is a Laurent polynomial in $q$ with integer cofficients, where $[n]=1+q+\cdots+q^{n-1}$ and ${n\brack k}=\prod_{i=1}^k(1-q^{n-i+1})/(1-q^i)$. This gives a $q$-analogue of a divisibility result on the Catalan triangle obtained by the first author and Zeng, and also confirms a conjecture of the first author and Zeng. We further propose several related conjectures.


title: Accelerated Magnetic Resonance Spectroscopy with Vandermonde Factorization

abstract: Multi-dimensional magnetic resonance spectroscopy is an important tool for studying molecular structures, interactions and dynamics in bio-engineering. The data acquisition time, however, is relatively long and non-uniform sampling can be applied to reduce this time. To obtain the full spectrum,a reconstruction method with Vandermonde factorization is proposed.This method explores the general signal property in magnetic resonance spectroscopy: Its time domain signal is approximated by a sum of a few exponentials. Results on synthetic and realistic data show that the new approach can achieve faithful spectrum reconstruction and outperforms state-of-the-art low rank Hankel matrix method.


title: Approximation of Generalized Ridge Functions in High Dimensions

abstract: This paper studies the approximation of generalized ridge functions, namely of functions which are constant along some submanifolds of $\mathbb{R}^N$. We introduce the notion of linear-sleeve functions, whose function values only depend on the distance to some unknown linear subspace $L$. We propose two effective algorithms to approximate linear-sleeve functions $f(x)=g(\text{dist}(x,L)^2)$, when both the linear subspace $L\subset \mathbb{R}^N$ and the function $g\in C^s[0,1]$ are unknown. We will prove error bounds for both algorithms and provide an extensive numerical comparison of both. We further propose an approach of how to apply these algorithms to capture general sleeve functions, which are constant along some lower dimensional submanifolds.


title: Imaging Spectropolarimeter for Multi-application Solar Telescope at Udaipur Solar Observatory: Characterization of polarimeter and preliminary observations

abstract: Multi-Application Solar Telescope (MAST) is a 50 cm off-axis Gregorian telescope and started operational recently at Udaipur Solar Observatory (USO). For understanding the evolution and dynamics of solar magnetic and velocity fields, an imaging spectropolarimeter is being developed as one of the back-end instruments of MAST. This system consists of a narrow-band filter and a polarimeter. Polarimeter includes a linear polarizer and two sets of Liquid Crystal Variable Retarders (LCVRs). The instrument is intended for the simultaneous observations in the spectral lines 617.3 nm and 854.2 nm, which are formed in photosphere and chromosphere, respectively. In this paper, we present results from the characterization of the LCVRs for the spectral lines of interest and response matrix of the polarimeter. We also present preliminary observations of an active region obtained using the spectropolarimeter. For verification, we compare the Stokes observations of the active region obtained from Helioseismic Magnetic Imager (HMI) onboard Solar Dynamics Observatory (SDO) with that of MAST observations in the spectral line 617.3 nm. We found good agreement between both the observations, considering the fact that MAST observations are seeing limited.


title: Self-adjoint Matrices are Equivariant

abstract: In this short note we prove that a matrix $A\in\mathbb{R}^{n,n}$ is self-adjoint if and only if it is equivariant with respect to the action of a group $\Gamma\subset {\bf O}(n)$ which is isomorphic to $\otimes_{k=1}^n\mathbf{Z}_2$. Moreover we discuss potential applications of this result, and we use it in particular for the approximation of higher order derivatives for smooth real valued functions of several variables.


title: Short Term Power Demand Prediction Using Stochastic Gradient Boosting

abstract: Power prediction demand is vital in power system and delivery engineering fields. By efficiently predicting the power demand, we can forecast the total energy to be consumed in a certain city or district. Thus, exact resources required to produce the demand power can be allocated. In this paper, a Stochastic Gradient Boosting (aka Treeboost) model is used to predict the short term power demand for the Emirate of Sharjah in the United Arab Emirates (UAE). Results show that the proposed model gives promising results in comparison to the model used by Sharjah Electricity and Water Authority (SEWA).


title: Experimentally Validated Bond Graph Model of a Brazed-Plate Heat Exchanger (BPHE)

abstract: The paper deals with the Bond Graph (BG) modeling and the model validation of a brazed-plate heat exchanger. This device is an important part of a mechanically pumped cooling loop. A thermo hydraulic BG model is developed and compared with experimental data. Optimization is performed to determine the best value of the convection heat exchange coefficients to be fixed in the model.


title: On the monotonicity of weighted power means of matrices

abstract: Let $\mu_p(A,B,t)=(tA^p+(1-t)B^p)^{1/p}$ denote the weighted power mean between positive operators $A$ and $B$. We show that the function $t\to \|A-\mu_p(A,B,t)\|_2$ is monotonically decreasing whenever $1/2 \leq p \leq 1$. Hence showing that the weighted power means satisfy Audenaert's "in-betweenness" property for positive operators for power satisfying $1/2 \leq p \leq 1$. We also show that when $p>2$ there exist operators for which the weighted power mean does not satisfy this "in-betweenness" property with respect to the Euclidean metric.


title: A maximal Boolean sublattice that is not the range of a Banaschewski function

abstract: We construct a countable bounded sublattice of the lattice of all subspaces of a vector space with two non-isomorphic maximal Boolean sublattice. We represent one of them as the range of a Banschewski function and we prove that this is not the case of the other. Hereby we solve a problem of F. Wehrung.


title: Unveiling the AGN in IC 883: discovery of a parsec-scale radio jet

abstract: IC883 is a luminous infrared galaxy (LIRG) classified as a starburst-active galactic nucleus (AGN) composite. In a previous study we detected a low-luminosity AGN (LLAGN) radio candidate. Here we report on our radio follow-up at three frequencies which provides direct and unequivocal evidence of the AGN activity in IC883. Our analysis of archival X-ray data, together with the detection of a transient radio source with luminosity typical of bright supernovae, give further evidence of the ongoing star formation activity, which dominates the energetics of the system. At sub-parsec scales, the radio nucleus has a core-jet morphology with the jet being a newly ejected component showing a subluminal proper motion of 0.6c-1c. The AGN contributes less than two per cent of the total IR luminosity of the system. The corresponding Eddington factor is ~1E-3, suggesting this is a low-accretion rate engine, as often found in LLAGNs. However, its high bolometric luminosity (~10E44erg/s) agrees better with a normal AGN. This apparent discrepancy may just be an indication of the transition nature of the nucleus from a system dominated by star-formation, to an AGN-dominated system. The nucleus has a strongly inverted spectrum and a turnover at ~4.4GHz, thus qualifying as a candidate for the least luminous (L_5.0GHz ~ 6.3E28erg/s/Hz) and one of the youngest (~3000yr) gigahertz-peaked spectrum (GPS) sources. If the GPS origin for the IC883 nucleus is confirmed, then advanced mergers in the LIRG category are potentially key environments to unveil the evolution of GPS sources into more powerful radio galaxies.


title: Analysing the 21cm signal from the Epoch of Reionization with artificial neural networks

abstract: The 21 cm signal from the Epoch of Reionization should be observed within the next decade. While a simple statistical detection is expected with SKA pathfinders, the SKA will hopefully produce a full 3D mapping of the signal. To extract from the observed data constraints on the parameters describing the underlying astrophysical processes, inversion methods must be developed. For example, the Markov Chain Monte Carlo method has been successfully applied. Here we test another possible inversion method: artificial neural networks (ANN). We produce a training set which consists of 70 individual sample. Each sample is made of the 21 cm power spectrum at different redshifts produced with the 21cmFast code plus the value of three parameters used in the semi-numerical simulations that describe astrophysical processes. Using this set we train the network to minimize the error between the parameter values it produces as an output and the true values. We explore the impact of the architecture of the network on the quality of the training. Then we test the trained network on the new set of 54 test samples with different values of the parameters. We find that the quality of the parameter reconstruction depends on the sensitivity of the power spectrum to the different parameters at a given redshift, that including thermal noise and sample variance decreases the quality of the reconstruction and that using the power spectrum at several redshifts as an input to the ANN improves the quality of the reconstruction. We conclude that ANNs are a viable inversion method whose main strength is that they require a sparse exploration of the parameter space and thus should be usable with full numerical simulations.


title: Analytic Closures for M1 Neutrino Transport

abstract: Carefully accounting for neutrino transport is an essential component of many astrophysical studies. Solving the full transport equation is too expensive for most realistic applications, especially those involving multiple spatial dimensions. For such cases, resorting to approximations is often the only viable option for obtaining solutions. One such approximation, which recently became popular, is the M1 method. It utilizes the system of the lowest two moments of the transport equation and closes the system with an ad hoc closure relation. The accuracy of the M1 solution depends on the quality of the closure. Several closures have been proposed in the literature and have been used in various studies. We carry out an extensive study of these closures by comparing the results of M1 calculations with precise Monte Carlo calculations of the radiation field around spherically-symmetric protoneutron star models. We find that no closure performs consistently better or worse than others in all cases. The level of accuracy a given closure yields depends on the matter configuration, neutrino type, and neutrino energy. Given this limitation, the maximum entropy closure by Minerbo (1978) on average yields relatively accurate results in the broadest set of cases considered in this work.


title: Magnetic field formation in the Milky Way-like disk galaxies of the Auriga project

abstract: The magnetic fields observed in the Milky~Way and nearby galaxies appear to be in equipartition with the turbulent, thermal, and cosmic ray energy densities, and hence are expected to be dynamically important. However, the origin of these strong magnetic fields is still unclear, and most previous attempts to simulate galaxy formation from cosmological initial conditions have ignored them altogether. Here, we analyse the magnetic fields predicted by the simulations of the Auriga Project, a set of 30 high-resolution cosmological zoom simulations of Milky~Way-like galaxies, carried out with a moving-mesh magneto-hydrodynamics code and a detailed galaxy formation physics model. We find that the magnetic fields grow exponentially at early times owing to a small-scale dynamo with an e-folding time of roughly $100\,\rm{Myr}$ in the center of halos until saturation occurs around $z=2-3$, when the magnetic energy density reaches about $10\%$ of the turbulent energy density with a typical strength of $10-50\,\rm{\mu G}$. In the galactic centers the ratio between magnetic and turbulent energy remains nearly constant until $z=0$. At larger radii, differential rotation in the disks leads to linear amplification that typically saturates around $z=0.5$ to $z=0$. The final radial and vertical variations of the magnetic field strength can be well described by two joint exponential profiles, and are in good agreement with observational constraints. Overall, the magnetic fields have only little effect on the global evolution of the galaxies as it takes too long to reach equipartition. We also demonstrate that our results are well converged with numerical resolution.


title: Are pulsars spun up or down by SASI spiral modes?

abstract: Pulsars may either be spun up or down by hydrodynamic instabilities during the supernova explosion of massive stars. Besides rapidly-rotating cases related to bipolar explosions, stellar rotation may affect the explosion of massive stars in the more common situations where the centrifugal force is minor. Using 2D simulations of a simplified setup in cylindrical geometry, we examine the impact of rotation on the Standing Accretion Shock Instability (SASI) and the corotation instability, also known as low-T/|W|. The influence of rotation on the saturation amplitude of these instabilities depends on the specific angular momentum in the accretion flow and the ratio of the shock to the neutron star radii. The spiral mode of SASI becomes more vigorous with faster rotation only if this ratio is large enough. A corotation instability develops at large rotation rates and impacts the dynamics more dramatically, leading to a strong one-armed spiral wave. Non-axisymmetric instabilities are able to redistribute angular momentum radially and affect the pulsar spin at birth. A systematic study of the relationship between the core rotation period of the progenitor and the initial pulsar spin is performed. Stellar rotation rates for which pulsars are spun up or down by SASI are estimated. Rapidly spinning progenitors are modestly spun down by spiral modes, less than $\sim 30\%$, when a corotation instability develops. Given the observational constraints on pulsar spin periods at birth, this suggests that rapid rotation might not play a significant hydrodynamic role in most core-collapse supernovae.


title: Unifying the micro and macro properties of AGN feeding and feedback

abstract: We unify the feeding and feedback of supermassive black holes with the global properties of galaxies, groups, and clusters, by linking for the first time the physical mechanical efficiency at the horizon and Mpc scale. The macro hot halo is tightly constrained by the absence of overheating and overcooling as probed by X-ray data and hydrodynamic simulations ($\varepsilon_{\rm BH} \simeq$ 10$^{-3}\, T_{\rm x,7.4}$). The micro flow is shaped by general relativistic effects tracked by state-of-the-art GR-RMHD simulations ($\varepsilon_\bullet \simeq$ 0.03). The SMBH properties are tied to the X-ray halo temperature $T_{\rm x}$, or related cosmic scaling relation (as $L_{\rm x}$). The model is minimally based on first principles, as conservation of energy and mass recycling. The inflow occurs via chaotic cold accretion (CCA), the rain of cold clouds condensing out of the quenched cooling flow and recurrently funneled via inelastic collisions. Within 100s gravitational radii, the accretion energy is transformed into ultrafast 10$^4$ km s$^{-1}$ outflows (UFOs) ejecting most of the inflowing mass. At larger radii the energy-driven outflow entrains progressively more mass: at roughly kpc scale, the velocities of the hot/warm/cold outflows are a few 10$^3$, 1000, 500 km s$^{-1}$, with median mass rates ~10, 100, several 100 M$_\odot$ yr$^{-1}$, respectively. The unified CCA model is consistent with the observations of nuclear UFOs, and ionized, neutral, and molecular macro outflows. We provide step-by-step implementation for subgrid simulations, (semi)analytic works, or observational interpretations which require self-regulated AGN feedback at coarse scales, avoiding the a-posteriori fine-tuning of efficiencies.

categories: astro-ph.HE astro-ph.GA physics.comp-ph physics.flu-dyn


title: Lyman-Werner Escape Fractions from the First Galaxies

abstract: Direct collapse black holes forming in pristine, atomically-cooling haloes at $z \approx 10-20$ may act as the seeds of supermassive black holes (BH) at high redshifts. In order to create a massive BH seed, the host halo needs to be prevented from forming stars. H$_2$ therefore needs to be irradiated by a large flux of Lyman-Werner (LW) UV photons in order to suppress H$_2$ cooling. A key uncertainty in this scenario is the escape fraction of LW radiation from first galaxies, the dominant source of UV photons at this epoch. To better constrain this escape fraction, we have performed radiation-hydrodynamical simulations of the growth of HII regions and their associated photodissociation regions in the first galaxies using the ZEUS-MP code. We find that the LW escape fraction crucially depends on the propagation of the ionisation front (I-front). For an R-type I-front overrunning the halo, the LW escape fraction is always larger than 95%. If the halo recombines later from the outside--in, due to a softened and weakened spectrum, the LW escape fraction in the rest-frame of the halo (the near-field) drops to zero. A detailed and careful analysis is required to analyse slowly moving, D-type I-fronts, where the escape fraction depends on the microphysics and can be as small as 3% in the near-field and 61% in the far-field or as large as 100% in both the near-field and the far-field.


title: Delay-time distribution of core-collapse supernovae with late events resulting from binary interaction

abstract: Most massive stars, the progenitors of core-collapse supernovae, are in close binary systems and may interact with their companion through mass transfer or merging. We undertake a population synthesis study to compute the delay-time distribution of core-collapse supernovae, that is, the supernova rate versus time following a starburst, taking into account binary interactions. We test the systematic robustness of our results by running various simulations to account for the uncertainties in our standard assumptions. We find that a significant fraction, $15^{+9}_{-8}$%, of core-collapse supernovae are `late', that is, they occur 50-200 Myrs after birth, when all massive single stars have already exploded. These late events originate predominantly from binary systems with at least one, or, in most cases, with both stars initially being of intermediate mass ($4-8M_{\odot}$). The main evolutionary channels that contribute often involve either the merging of the initially more massive primary star with its companion or the engulfment of the remaining core of the primary by the expanding secondary that has accreted mass at an earlier evolutionary stage. Also, the total number of core-collapse supernovae increases by $14^{+15}_{-14}$% because of binarity for the same initial stellar mass. The high rate implies that we should have already observed such late core-collapse supernovae, but have not recognized them as such. We argue that $\phi$ Persei is a likely progenitor and that eccentric neutron star - white dwarf systems are likely descendants. Late events can help explain the discrepancy in the delay-time distributions derived from supernova remnants in the Magellanic Clouds and extragalactic type Ia events, lowering the contribution of prompt Ia events. We discuss ways to test these predictions and speculate on the implications for supernova feedback in simulations of galaxy evolution.


title: Kinetic and radiative power from optically thin accretion flows

abstract: We perform a set of general relativistic, radiative, magneto-hydrodynamical simulations (GR-RMHD) to study the transition from radiatively inefficient to efficient state of accretion on a non-rotating black hole. We study ion to electron temperature ratios ranging from $T_{\rm i}/T_{\rm e}=10$ to $100$, and simulate flows corresponding to accretion rates as low as $10^{-6}\dot M_{\rm Edd}$, and as high as $10^{-2}\dot M_{\rm Edd}$. We have found that the radiative output of accretion flows increases with accretion rate, and that the transition occurs earlier for hotter electrons (lower $T_{\rm i}/T_{\rm e}$ ratio). At the same time, the mechanical efficiency hardly changes and accounts to ${\approx}\,3\%$ of the accreted rest mass energy flux, even at the highest simulated accretion rates. This is particularly important for the mechanical AGN feedback regulating massive galaxies, groups, and clusters. Comparison with recent observations of radiative and mechanical AGN luminosities suggests that the ion to electron temperature ratio in the inner, collisionless accretion flow should fall within $10<T_{\rm i}/T_{\rm e}<30$, i.e., the electron temperature should be several percent of the ion temperature.

categories: astro-ph.HE astro-ph.CO astro-ph.GA physics.comp-ph physics.plasm-ph


title: Constraining the Milky Way assembly history with Galactic Archaeology. Ludwig Biermann Award Lecture 2015

abstract: The aim of Galactic Archaeology is to recover the evolutionary history of the Milky Way from its present day kinematical and chemical state. Because stars move away from their birth sites, the current dynamical information alone is not sufficient for this task. The chemical composition of stellar atmospheres, on the other hand, is largely preserved over the stellar lifetime and, together with accurate ages, can be used to recover the birthplaces of stars currently found at the same Galactic radius. In addition to the availability of large stellar samples with accurate 6D kinematics and chemical abundance measurements, this requires detailed modeling with both dynamical and chemical evolution taken into account. An important first step is to understand the variety of dynamical processes that can take place in the Milky Way, including the perturbative effects of both internal (bar and spiral structure) and external (infalling satellites) agents. We discuss here (1) how to constrain the Galactic bar, spiral structure, and merging satellites by their effect on the local and global disc phase-space, (2) the effect of multiple patterns on the disc dynamics, and (3) the importance of radial migration and merger perturbations for the formation of the Galactic thick disc. Finally, we discuss the construction of Milky Way chemo-dynamical models and relate to observations.


title: Variational optimization algorithms for uniform matrix product states

abstract: We combine the Density Matrix Renormalization Group (DMRG) with Matrix Product State tangent space concepts to construct a variational algorithm for finding ground states of one dimensional quantum lattices in the thermodynamic limit. A careful comparison of this variational uniform Matrix Product State algorithm (VUMPS) with infinite Density Matrix Renormalization Group (IDMRG) and with infinite Time Evolving Block Decimation (ITEBD) reveals substantial gains in convergence speed and precision. We also demonstrate that VUMPS works very efficiently for Hamiltonians with long range interactions and also for the simulation of two dimensional models on infinite cylinders. The new algorithm can be conveniently implemented as an extension of an already existing DMRG implementation.


title: HI absorption in nearby compact radio galaxies

abstract: HI absorption studies yield information on both AGN feeding and feedback processes. This AGN activity interacts with the neutral gas in compact radio sources, which are believed to represent the young or recently re-triggered AGN population. We present the results of a survey for HI absorption in a sample of 66 compact radio sources at 0.040 < z < 0.096 with the Australia Telescope Compact Array. In total, we obtained seven detections, five of which are new, with a large range of peak optical depths (3% to 87%). Of the detections, 71% exhibit asymmetric, broad ($\Delta{v}_{\mathrm{FWHM}}$ > 100 km s$^{-1}$) features, indicative of disturbed gas kinematics. Such broad, shallow and offset features are also found within low-excitation radio galaxies which is attributed to disturbed circumnuclear gas, consistent with early-type galaxies typically devoid of a gas-rich disk. Comparing mid-infrared colours of our galaxies with HI detections indicates that narrow and deep absorption features are preferentially found in late-type and high-excitation radio galaxies in our sample. These features are attributed to gas in galactic disks. By combining XMM-Newton archival data with 21-cm data, we find support that absorbed X-ray sources may be good tracers of HI content within the host galaxy. This sample extends previous HI surveys in compact radio galaxies to lower radio luminosities and provides a basis for future work exploring the higher redshift universe.


title: Quasar emission lines as probes of orientation: implications for disc wind geometries and unification

abstract: The incidence of broad absorption lines (BALs) in quasar samples is often interpreted in the context of a geometric unification model consisting of an accretion disc and an associated outflow. We use the the Sloan Digital Sky Survey (SDSS) quasar sample to test this model by examining the equivalent widths (EWs) of CIV 1550\AA, Mg II 2800\AA, [OIII] 5007\AA\ and C III] 1909\AA. We find that the emission line EW distributions in BAL and non-BAL quasars are remarkably similar -- a property that is inconsistent with scenarios in which a BAL outflow rises equatorially from a geometrically thin, optically thick accretion disc. We construct simple models to predict the distributions from various geometries; these models confirm the above finding and disfavour equatorial geometries. We show that obscuration, line anisotropy and general relativistic effects on the disc continuum are unlikely to hide an EW inclination dependence. We carefully examine the radio and polarisation properties of BAL quasars. Both suggest that they are most likely viewed (on average) from intermediate inclinations, between type 1 and type 2 AGN. We also find that the low-ionization BAL quasars in our sample are not confined to one region of `Eigenvector I' parameter space. Overall, our work leads to one of the following conclusions, or some combination thereof: (i) the continuum does not emit like a geometrically thin, optically thick disc; (ii) BAL quasars are viewed from similar angles to non-BAL quasars, i.e. low inclinations; (iii) geometric unification does not explain the fraction of BALs in quasar samples.


title: Can we observe neutrino flares in coincidence with explosive transients?

abstract: The new generation of powerful instruments is reaching sensitivities and temporal resolutions that will allow multi-messenger astronomy of explosive transient phenomena, with high-energy neutrinos as a central figure. We derive general criteria for the detectability of neutrinos from powerful transient sources, for given instrument sensitivities. In practice, we provide the minimal photon flux necessary for neutrino detection based on two main observables, the bolometric luminosity and the time variability of the emission. This limit can be compared to the observations in specified wavelengths in order to target the most promising sources for follow-ups. Our criteria can also help discriminating false associations of neutrino events with a flaring source. We find that relativistic transient sources such as High- and Low-Luminosity GRBs, Blazar flares, Tidal Disruption Events and magnetar flares could be observed with IceCube, as they have a good chance to occur within a detectable distance. Among non-relativistic transient sources, only luminous supernovae appear as promising candidates. We caution that our criterion should not be directly applied to low-luminosity GRBs and Type Ibc supernovae, as these objects could have hosted a choked GRB, leading to neutrino emission without a relevant counterpart radiation. We treat a set of concrete examples and show that several transients, some of which are being monitored by IceCube, are far from meeting the criterion for detectability (e.g., Crab flares or Swift J1644+57).


title: Power Spectra Based Planck Constraints on Compensated Isocurvature, and Forecasts for LiteBIRD and CORE Space Missions

abstract: Compensated isocurvature perturbations (CIP), where the baryon and cold dark matter perturbations cancel, do not cause total matter isocurvature perturbation. Consequently, at the linear order in the baryon density contrast $\Delta$, CIP is not detectable by the CMB power spectra. At the second order CIP smoothes the power spectra in a similar manner as lensing, causing a degeneracy between the CIP variance $\Delta^2_{rms}=<\Delta^2>$ and lensing parameter $A_L$. We show that the CMB lensing data breaks this degeneracy. Nested sampling of the LCDM+CIP(+$A_L$) model, the Planck 2015 temperature, polarization, and lensing data give $\Delta^2_{rms}=0.0069\pm0.0030$ at 68% CL. A non-zero value is favored at 2.3$\sigma$. CIP with $\Delta^2_{rms}=0.007$ improves the bestfit $\chi^2$ by 3.6 compared to the adiabatic LCDM model. In contrast, although the temperature data favor $A_L=1.22$, allowing $A_L\ne1$ does not improve the joint fit, since the lensing data disfavor $A_L\ne1$. Indeed, CIP provides a rare example of a simple model, which can reduce the Planck lensing anomaly by fitting well simultaneously the high multipole temperature and lensing data, as well as the polarization data. Finally, we derive forecasts for future satellite missions (LiteBIRD proposal to JAXA and Exploring Cosmic Origins with CORE proposal to ESA's M5 call). Due to its coarse angular resolution, LiteBIRD is not able to improve the constraints on CIP or $A_L$, but CORE-M5 approaches the cosmic variance limit and improves the CIP constraint to $\Delta^2_{rms}<0.0006\ (0.0014)$ at 68% (95%) CL, which is 9 times better than the current trispectrum based upper bound and 6 times better than obtained from the simulated Planck data. In addition, CORE-M5 will exquisitely distinguish between CIP and $A_L$. No matter whether CIP is allowed for or not, the uncertainty of the lensing parameter will be $\sigma(A_L)=0.012$.


title: Simultaneous, full characterization of a single-photon state using semiconductor quantum-dot light

abstract: As single-photon sources become more mature and are used more often in quantum information, communications and measurement applications, their characterization becomes more important. Single-photon-like light is often characterized by its brightness, and two quantum properties: the single-photon composition and the photon indistinguishability. While it is desirable to obtain these quantities from a single measurement, currently two or more measurements are required. Here, we simultaneously determine the brightness, the single photon purity, the indistinguishability, and the statistical distribution of Fock states to third order for a quantum light source. The measurement uses a pair of two-photon (n = 2) number-resolving detectors. n > 2 number-resolving detectors provide no additional advantage in the single-photon characterization. The new method extracts more information per experimental trial than a conventional measurement for all input states, and is particularly more e cient for statistical mixtures of photon states. Thus, using this n=2, number- resolving detector scheme will provide advantages in a variety of quantum optics measurements and systems.


title: The Structure of the Distant Kuiper Belt in a Nice Model Scenario

abstract: This work explores the orbital distribution of minor bodies in the outer Solar System emplaced as a result of a Nice model migration from the simulations of Brasser & Morbidelli (2013). This planetary migration scatters a planetesimal disk from between 29-34 AU and emplaces a population of objects into the Kuiper belt region. From the 2:1 Neptune resonance and outward, the test particles analyzed populate the outer resonances with orbital distributions consistent with trans-Neptunian object (TNO) detections in semi-major axis, inclination, and eccentricity, while capture into the closest resonances is too efficient. The relative populations of the simulated scattering objects and resonant objects in the 3:1 and 4:1 resonances are also consistent with observed populations based on debiased TNO surveys, but the 5:1 resonance is severely underpopulated compared to population estimates from survey results. Scattering emplacement results in the expected orbital distribution for the majority of the TNO populations, however the origin of the large observed population in the 5:1 resonance remains unexplained.


title: Riesz Bases of Exponentials on Unbounded Multi-tiles

abstract: We prove the existence of Riesz bases of exponentials of L^2(Omega), provided that Omega in R^d is a measurable set of finite and positive measure, not necessarily bounded, that satisfies a multi-tiling condition and an arithmetic property that we call admissibility. This property is satisfied for any bounded domain, so our results extend the known case of bounded multi-tiles. We also extend known results for submulti-tiles and frames of exponentials to the unbounded case.


title: Quartified Leptonic Color, Bound States, and Future Electron-Positron Collider

abstract: The $[SU(3)]^4$ quartification model of Babu, Ma, and Willenbrock (BMW), proposed in 2003, predicts a confining leptonic color $SU(2)$ gauge symmetry, which becomes strong at the keV scale. It also predicts the existence of three families of half-charged leptons (hemions) below the TeV scale. These hemions are confined to form bound states which are not so easy to discover at the Large Hadron Collider (LHC). However, just as $J/\psi$ and $\Upsilon$ appeared as sharp resonances in $e^-e^+$ colliders of the 20th century, the corresponding 'hemionium' states are expected at a future $e^-e^+$ collider of the 21st century.


title: Experimentally feasible protocol for semiquantum key distribution

abstract: Quantum key distribution (QKD) protocols make it possible for two quantum parties to generate a secret shared key. Semiquantum key distribution (SQKD) protocols, such as "QKD with classical Bob" and "QKD with classical Alice" (that have both been proven robust), achieve this goal even if one of the parties is classical. However, existing SQKD protocols are not experimentally feasible with current technology. Here we suggest a new protocol, "Classical Alice with a controllable mirror", that can be experimentally implemented with current technology (using 4-level systems instead of qubits), and we prove it to be robust.


title: Polynomial Wolff axioms and Kakeya-type estimates in $\mathbb{R}^4$

abstract: We establish new linear and trilinear bounds for collections of tubes in $\mathbb{R}^4$ that satisfy the polynomial Wolff axioms. In brief, a collection of $\delta$-tubes satisfies the Wolff axioms if not too many tubes can be contained in the $\delta$-neighborhood of a plane. A collection of tubes satisfies the polynomial Wolff axioms if not too many tubes can be contained in the $\delta$-neighborhood of a low degree algebraic variety. First, we prove that if a set of $\delta^{-3}$ tubes in $\mathbb{R}^4$ satisfies the polynomial Wolff axioms, then the union of the tubes must have volume at least $\delta^{1-1/40}$. We also prove a more technical statement which is analogous to a maximal function estimate at dimension $3+1/40$. Second, we prove that if a collection of $\delta^{-3}$ tubes in $\mathbb{R}^4$ satisfies the polynomial Wolff axioms, and if most triples of intersecting tubes point in three linearly independent directions, then the union of the tubes must have volume at least $\delta^{3/4}$. Again, we also prove a slightly more technical statement which is analogous to a maximal function estimate at dimension $3+1/4$. We conjecture that every Kakeya set satisfies the polynomial Wolff axioms, but we are unable to prove this. If our conjecture is correct, it implies a Kakeya maximal function estimate at dimension $3+1/40$, and in particular this implies that every Kakeya set in $\mathbb{R}^4$ must have Hausdorff dimension at least $3+1/40$. This would be an improvement over the current best bound of 3, which was established by Wolff in 1995.


title: Large Scale Novel Object Discovery in 3D

abstract: We present a method for discovering never-seen-before objects in 3D point clouds obtained from sensors like Microsoft Kinect. We generate supervoxels directly from the point cloud data and use them with a Siamese network, built on a recently proposed 3D convolutional neural network architecture. We use known objects to train a non-linear embedding of supervoxels, by optimizing the criteria that supervoxels which fall on the same object should be closer than those which fall on different objects, in the embedding space. We test on unknown objects, which were not seen during training, and perform clustering in the learned embedding space of supervoxels to effectively perform novel object discovery. We validate the method with extensive experiments, quantitatively showing that it can discover numerous unseen objects while being trained on only a few dense 3D models. We also show very good qualitative results of object discovery in point cloud data when the test objects, either specific instances or even categories, were never seen during training.


title: NLO electroweak radiative corrections for four-fermionic process at Belle II

abstract: We discuss the next to the leading order (NLO) electroweak radiative corrections to the $e^- e^+ \rightarrow f^- f^+ (\gamma)$ cross section asymmetry, for polarized and unpolarized beam scenario. The left-right and forward-backward amplitudes, with and without radiative corrections, are evaluated and compared for various kinematics. The hard bremsstrahlung is included for arbitrary energy cuts. The radiative corrections are shown to be significant and having a non-trivial dependency on the kinematic conditions. The calculations are relevant for the ultra-precise low-energy experiment Belle~II planned at SuperKEKB.


title: A Connection Between Orthogonal Polynomials and Shear Instabilities in the Quasi-geostrophic Shallow Water Equations

abstract: In this paper we demonstrate a connection between the roots of a certain sequence of orthogonal polynomials on the real line and the linear instability of a $x$-directionally homogeneous background velocity profile $u^b(x,y) = \cos(y)$ in the quasi-geostrophic shallow water (QG) equation in a domain with periodic boundaries in the $y$-direction. Using the relationship we establish, we then prove that there exists a unique unstable mode for each horizontal wave number $0<k<1$ and provide mathematically rigorous estimates of the associated growth rate.


title: Spin mediated enhanced negative magnetoresistance in Ni80Fe20 and p-silicon bilayer

abstract: In this work, we present an experimental study of spin mediated enhanced negative magnetoresistance in Ni80Fe20 (50 nm)/p-Si (350 nm) bilayer. The resistance measurement shows a reduction of ~2.5% for the bilayer specimen as compared to 1.3% for Ni80Fe20 (50 nm) on oxide specimen for an out-of-plane applied magnetic field of 3T. In the Ni80Fe20-only film, the negative magnetoresistance behavior is attributed to anisotropic magnetoresistance. We propose that spin polarization due to spin-Hall effect is the underlying cause of the enhanced negative magnetoresistance observed in the bilayer. Silicon has weak spin orbit coupling so spin Hall magnetoresistance measurement is not feasible. We use V2{\omega} and V3{\omega} measurement as a function of magnetic field and angular rotation of magnetic field in direction normal to electric current to elucidate the spin-Hall effect. The angular rotation of magnetic field shows a sinusoidal behavior for both V2{\omega} and V3{\omega}, which is attributed to the spin phonon interactions resulting from the spin-Hall effect mediated spin polarization. We propose that the spin polarization leads to a decrease in hole-phonon scattering resulting in enhanced negative magnetoresistance.


title: Exogeneity tests, incomplete models, weak identification and non-Gaussian distributions: invariance and finite-sample distributional theory

abstract: We study the distribution of Durbin-Wu-Hausman (DWH) and Revankar-Hartley (RH) tests for exogeneity from a finite-sample viewpoint, under the null and alternative hypotheses. We consider linear structural models with possibly non-Gaussian errors, where structural parameters may not be identified and where reduced forms can be incompletely specified (or nonparametric). On level control, we characterize the null distributions of all the test statistics. Through conditioning and invariance arguments, we show that these distributions do not involve nuisance parameters. In particular, this applies to several test statistics for which no finite-sample distributional theory is yet available, such as the standard statistic proposed by Hausman (1978). The distributions of the test statistics may be non-standard -- so corrections to usual asymptotic critical values are needed -- but the characterizations are sufficiently explicit to yield finite-sample (Monte-Carlo) tests of the exogeneity hypothesis. The procedures so obtained are robust to weak identification, missing instruments or misspecified reduced forms, and can easily be adapted to allow for parametric non-Gaussian error distributions. We give a general invariance result (block triangular invariance) for exogeneity test statistics. This property yields a convenient exogeneity canonical form and a parsimonious reduction of the parameters on which power depends. In the extreme case where no structural parameter is identified, the distributions under the alternative hypothesis and the null hypothesis are identical, so the power function is flat, for all the exogeneity statistics. However, as soon as identification does not fail completely, this phenomenon typically disappears.


title: The multiresolution analysis of flow graphs

abstract: We introduce and prove basic results about several graph-theoretic notions relevant to the multiresolution analysis of flow graphs that represent the transfer of control in computer programs. We take a category-theoretical viewpoint to demonstrate that our definitions are natural and to motivate particular incarnations of related constructions.


title: Causality and Unitarity via the Tree-Loop Duality Relation

abstract: The tree-loop duality relation is used as a starting point to derive the constraints of causality and unitarity. Specifically, the Bogoliubov causality condition is ab initio derived at the individual graph level. It leads to a representation of a graph in terms of lower order cut graphs. Extracting the absorptive part gives then the general unitarity relation (Cutkosky rule). The derivation, being carried out directly in momentum space, holds for any local (polynomial) hermitian interaction vertices. This is in contrast to the technical difficulties arising from contact terms in the spacetime approach based on the largest time equation.


title: Natural and fishing mortalities affecting eastern sea garfish, Hyporhamphus australis, inferred from age-frequency data using hazard functions

abstract: Estimates of age-specific natural (M) and fishing (F) mortalities among economically important stocks are required to determine sustainable yields and, ultimately, facilitate effective resource management. Here we used hazard functions to estimate mortality rates for eastern sea garfish, Hyporhamphus australis, a pelagic species that forms the basis of an Australian commercial lampara-net fishery. Data describing annual (2004 to 2015) age frequencies (0-1 to 5-6 years), yield, effort (boat-days), and average weights at age were used to fit various stochastic models to estimate mortality rates by maximum likelihood. The model best supported by the data implied: (i) the escape of fish aged 0-1 years increased from approximately 90 to 97% as a result of a mandated increase in stretched mesh opening from 25 to 28 mm; (ii) full selectivity among older age groups; (iii) a constant M of 0.52 +- 0.06 per year; and (iv) a decline in F between 2004 and 2015. Recruitment and biomass were estimated to vary, but increased during the sampled period. The results reiterate the utility of hazard functions to estimate and partition mortality rates, and support traditional input controls designed to reduce both accounted and unaccounted F.


title: Interaction effects on galaxy pairs with GeminiGMOS-III: Stellar population synthesis

abstract: We present an observational study of the impacts of the interactions on the stellar population in a sample of galaxy pairs. Long-slit spectra in the wavelength range 3440-7300 {\AA} obtained with the Gemini Multi-Object Spectrograph (GMOS) at Gemini South for fifteen galaxies in nine close pairs were used. The spatial distributions of the stellar population contributions were obtained using the stellar population synthesis code STARLIGHT. Taking into account the different contributions to the emitted light, we found that most of the galaxies in our sample are dominated by the young/intermediate stellar populations. This result differs from the one derived for isolated galaxies where the old stellar population dominates the disc surface brightness. We interpreted such different behavior as being due to the effect of gas inflows along the disk of interacting galaxies on the star formation in a time scale of the order of about 2Gyr. We also found that, in general, the secondary galaxy of the pairs has a higher contribution of the young stellar population than the primary one. We compared the estimated values of the stellar and nebular extinctions derived from the synthesis method and the H{\alpha}/H\b{eta} emission-line ratio finding that the nebular extinctions are systematically higher than stellar ones by about a factor of 2. We did not find any correlation between nebular and stellar metallicities. We neither found a correlation between stellar metallicities and ages while a positive correlation between nebular metallicities and stellar ages was obtained, with the older regions being the most metal-rich.


title: Tailoring spin defects in diamond

abstract: Atomic-size spin defects in solids are unique quantum systems. Most applications require nanometer positioning accuracy, which is typically achieved by low energy ion implantation. So far, a drawback of this technique is the significant residual implantation-induced damage to the lattice, which strongly degrades the performance of spins in quantum applications. In this letter we show that the charge state of implantation-induced defects drastically influences the formation of lattice defects during thermal annealing. We demonstrate that charging of vacancies localized at e.g. individual nitrogen implantation sites suppresses the formation of vacancy complexes, resulting in a tenfold-improved spin coherence time of single nitrogen-vacancy (NV) centers in diamond. This has been achieved by confining implantation defects into the space charge layer of free carriers generated by a nanometer-thin boron-doped diamond structure. Besides, a twofold-improved yield of formation of NV centers is observed. By combining these results with numerical calculations, we arrive at a quantitative understanding of the formation and dynamics of the implanted spin defects. The presented results pave the way for improved engineering of diamond spin defect quantum devices and other solid-state quantum systems.


title: Kitaev Materials

abstract: In transition-metal compounds with partially filled $4d$ and $5d$ shells spin-orbit entanglement, electronic correlations, and crystal-field effects conspire to give rise to a variety of novel forms of topological quantum matter. This includes Kitaev materials -- a family of spin-orbit assisted Mott insulators, in which local, spin-orbit entangled $j=1/2$ moments form that are subject to strong bond-directional interactions. On a conceptual level, Kitaev materials attract much interest for their unconventional forms of magnetism, such as spin liquid physics in two- and three-dimensional lattice geometries or the formation of non-trivial spin textures. Experimentally, a number of Kitaev materials have been synthesized, which includes the honeycomb materials Na$_2$IrO$_3$, $\alpha$-Li$_2$IrO$_3$, and RuCl$_3$, the triangular materials Ba$_3$Ir$_x$Ti$_{3-x}$O$_9$, as well as the three-dimensional hyper-honeycomb and stripy-honeycomb materials $\beta$-Li$_2$IrO$_3$ and $\gamma$-Li$_2$IrO$_3$. These lecture notes provide a short review of the current status of the theoretical and experimental exploration of these Kitaev materials.


title: Inferring Properties of the ISM from Supernova Remnant Size Distributions

abstract: We model the size distribution of supernova remnants to infer the surrounding ISM density. Using simple, yet standard SNR evolution models, we find that the distribution of ambient densities is remarkably narrow; either the standard assumptions about SNR evolution are wrong, or observable SNRs are biased to a narrow range of ambient densities. We show that the size distributions are consistent with log-normal, which severely limits the number of model parameters in any SNR population synthesis model. Simple Monte Carlo simulations demonstrate that the size distribution is indistinguishable from log-normal when the SNR sample size is less than 600. This implies that these SNR distributions provide only information on the mean and variance, yielding additional information only when the sample size grows larger than $\sim{600}$ SNRs. To infer the parameters of the ambient density, we use Bayesian statistical inference under the assumption that SNR evolution is dominated by the Sedov phase. In particular, we use the SNR sizes and explosion energies to estimate the mean and variance of the ambient medium surrounding SNR progenitors. We find that the mean ISM particle density around our sample of SNRs is $\mu_{\log{n}} = -1.33$, in $\log_{10}$ of particles per cubic centimeter, with variance $\sigma^2_{\log{n}} = 0.49$. If interpreted at face value, this implies that most SNRs result from supernovae propagating in the warm, ionized medium. However, it is also likely that either SNR evolution is not dominated by the simple Sedov evolution or SNR samples are biased to the warm, ionized medium (WIM).


title: If you are not paying for it, you are the product: How much do advertisers pay to reach you?

abstract: Online advertising is progressively moving towards a programmatic model in which ads are matched to actual interests of individuals collected as they browse the web. Letting the huge debate around privacy aside, a very important question in this area, for which little is known, is: How much do advertisers pay to reach an individual? In this study, we develop a first of its kind methodology for computing exactly that -- the price paid for a web user by the ad ecosystem -- and we do that in real time. Our approach is based on tapping on the Real Time Bidding (RTB) protocol to collect cleartext and encrypted prices for winning bids paid by advertisers in order to place targeted ads. Our main technical contribution is a method for tallying winning bids even when they are encrypted. We achieve this by training a model using as ground truth prices obtained by running our own "probe" ad-campaigns. We design our methodology through a browser extension and a back-end server that provides it with fresh models for encrypted bids. We validate our methodology using a one year long trace of 1600 mobile users and demonstrate that it can estimate a user's advertising worth with more than 82% accuracy.


title: Stable explicit schemes for simulation of nonlinear moisture transfer in porous materials

abstract: Implicit schemes have been extensively used in building physics to compute the solution of moisture diffusion problems in porous materials for improving stability conditions. Nevertheless, these schemes require important sub-iterations when treating non-linear problems. To overcome this disadvantage, this paper explores the use of improved explicit schemes, such as Dufort-Frankel, Crank-Nicolson and hyperbolisation approaches. A first case study has been considered with the hypothesis of linear transfer. The Dufort-Frankel, Crank-Nicolson and hyperbolisation schemes were compared to the classical Euler explicit scheme and to a reference solution. Results have shown that the hyperbolisation scheme has a stability condition higher than the standard Courant-Friedrichs-Lewy (CFL) condition. The error of this schemes depends on the parameter \tau representing the hyperbolicity magnitude added into the equation. The Dufort-Frankel scheme has the advantages of being unconditionally stable and is preferable for non-linear transfer, which is the second case study. Results have shown the error is proportional to O(\Delta t). A modified Crank-Nicolson scheme has been proposed in order to avoid sub-iterations to treat the non-linearities at each time step. The main advantages of the Dufort-Frankel scheme are (i) to be twice faster than the Crank-Nicolson approach; (ii) to compute explicitly the solution at each time step; (iii) to be unconditionally stable and (iv) easier to parallelise on high-performance computer systems. Although the approach is unconditionally stable, the choice of the time discretisation $\Delta t$ remains an important issue to accurately represent the physical phenomena.