title: Astronomical random numbers for quantum foundations experiments

abstract: Photons from distant astronomical sources can be used as a classical source of randomness to improve fundamental tests of quantum nonlocality, wave-particle duality, and local realism through Bell's inequality and delayed-choice quantum eraser tests inspired by Wheeler's cosmic-scale Mach-Zehnder interferometer gedankenexperiment. Such sources of random numbers may also be useful for information-theoretic applications such as key distribution for quantum cryptography. Building on the design of an "astronomical random-number generator" developed for the recent "cosmic Bell" experiment [Handsteiner et al., Phys. Rev. Lett. 118, 060401 (2017)], in this paper we report on the design and characterization of a device that, with 20-nanosecond latency, outputs a bit based on whether the wavelength of an incoming photon is greater than or less than 700 nm. Using the one-meter telescope at the Jet Propulsion Laboratory (JPL) Table Mountain Observatory, we generated random bits from astronomical photons in both color channels from 50 stars of varying color and magnitude, and from 12 quasars with redshifts up to $z = 3.9$. With stars, we achieved bit rates of $\sim 1 \times 10^6$ Hz / m$^2$, limited by saturation for our single-photon detectors, and with quasars of magnitudes between 12.9 and 16, we achieved rates between $\sim 10^2$ and $2 \times 10^3$ Hz /m$^2$. For bright quasars, the resulting bitstreams exhibit sufficiently low amounts of statistical predictability as quantified by the mutual information. In addition, a sufficiently high fraction of bits generated are of true astronomical origin in order to address both the locality and freedom-of-choice loopholes when used to set the measurement settings in a test of the Bell-CHSH inequality.


title: A Note on Multiparty Communication Complexity and the Hales-Jewett Theorem

abstract: For integers $n$ and $k$, the density Hales-Jewett number $c_{n,k}$ is defined as the maximal size of a subset of $[k]^n$ that contains no combinatorial line. We show that for $k \ge 3$ the density Hales-Jewett number $c_{n,k}$ is equal to the maximal size of a cylinder intersection in the problem $Part_{n,k}$ of testing whether $k$ subsets of $[n]$ form a partition. It follows that the communication complexity, in the Number On the Forehead (NOF) model, of $Part_{n,k}$, is equal to the minimal size of a partition of $[k]^n$ into subsets that do not contain a combinatorial line. Thus, the bound in \cite{chattopadhyay2007languages} on $Part_{n,k}$ using the Hales-Jewett theorem is in fact tight, and the density Hales-Jewett number can be thought of as a quantity in communication complexity. This gives a new angle to this well studied quantity. As a simple application we prove a lower bound on $c_{n,k}$, similar to the lower bound in \cite{polymath2010moser} which is roughly $c_{n,k}/k^n \ge \exp(-O(\log n)^{1/\lceil \log_2 k\rceil})$. This lower bound follows from a protocol for $Part_{n,k}$. It is interesting to better understand the communication complexity of $Part_{n,k}$ as this will also lead to the better understanding of the Hales-Jewett number. The main purpose of this note is to motivate this study.


title: Photometry and spectroscopy of multiple populations along the AGB of NGC2808 and NGC6121 (M4)

abstract: We present a photometric and spectroscopic study of multiple populations along the asymptotic-giant branch (AGB) of the intermediate-metallicity globular clusters (GCs) NGC2808 and NGC6121 (M4). Chemical abundances of O, Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Fe, Co, Ni, Zn, Y, and Ce in AGB stars from high-resolution FLAMES+UVES@VLT spectra are reported for both clusters. Our spectroscopic results have been combined with multi-wavelength photometry from the HST UV survey of Galactic GCs and ground-based photometry, plus proper motions derived by combining stellar positions from ground-based images and Gaia DR1. Our analysis reveals that the AGBs of both clusters host multiple populations with different chemical composition. In M4 we have identified two main populations of stars with different Na/O content, lying on distinct AGBs in the mF438W vs. C_F275W,F336W,F438W and the V vs.C_U,B,I pseudo-CMDs. In the more massive and complex GC NGC2808 three groups of stars with different chemical abundances occupy different locations on the so-called "chromosome map" photometric diagram. The spectroscopic+photometric comparison of stellar populations along the AGB and the red giants of this GC suggests that the AGB hosts stellar populations with a range in helium abundances spanning from primordial up to high contents of Y~0.32. On the other hand, from our dataset, there is no evidence for stars with extreme helium abundance (Y~0.38) on the AGB, suggesting that the most He-rich stars of NGC2808 do not reach this phase.


title: Cosmological dynamics of magnetic Bianchi I in viable $f(R)$ models of gravity

abstract: Standard dynamical system analysis of Einstein-Maxwell equation in $f(R)$ theories is considered in this work. We investigate cosmological dynamics of a uniform magnetic field in the Orthogonal Spatially Homogeneous (OSH) Bianchi type I universe with viable $f(R)$ models of gravity. In this work, the $f(R) = R -\alpha R^n$ and $f(R) = \left( R^b - \Lambda\right)^c$ models are examined by using our dynamical system analysis. Our results show that both of two $f(R)$ models have a viable cosmological consequence identical to the analysis present in Ref.\cite{Amendola:2007nt} for the FLRW background. Contrary to Ref.\cite{Amendola:2007nt}, we discover in our models that there is an additional anisotropic and non-zero cosmological magnetic fields fixed point emerging before the present of the standard matter epoch. This means that the universe has initially isotropic stage with the intermediated epoch as the anisotropic background and it ends up with the isotropic late-time acceleration. The primordial magnetic fields play a crucial role of the shear evolutions obtained from these two models which have the same scaling of the cosmic time as $\sigma\sim t^{-\frac13}$, instead of $\sigma\sim t^{-1}$ for the absence of the primordial magnetic cases.


title: Modeling influenza-like illnesses through composite compartmental models

abstract: Epidemiological models for the spread of pathogens in a population are usually only able to describe a single pathogen. This makes their application unrealistic in cases where multiple pathogens with similar symptoms are spreading concurrently within the same population. Here we describe a method which makes possible the application of multiple single-strain models under minimal conditions. As such, our method provides a bridge between theoretical models of epidemiology and data-driven approaches for modeling of influenza and other similar viruses. Our model extends the Susceptible-Infected-Recovered model to higher dimensions, allowing the modeling of a population infected by multiple viruses. We further provide a method, based on an overcomplete dictionary of feasible realizations of SIR solutions, to blindly partition the time series representing the number of infected people in a population into individual components, each representing the effect of a single pathogen. We demonstrate the applicability of our proposed method on five years of seasonal influenza-like illness (ILI) rates, estimated from Twitter data. We demonstrate that our method describes, on average, 44\% of the variance in the ILI time series. The individual infectious components derived from our model are matched to known viral profiles in the populations, which we demonstrate matches that of independently collected epidemiological data. We further show that the basic reproductive numbers ($R0$) of the matched components are in range known for these pathogens. Our results suggest that the proposed method can be applied to other pathogens and geographies, providing a simple method for estimating the parameters of epidemics in a population.


title: A multilevel block building algorithm for fast modeling generalized separable systems

abstract: Data-driven modeling plays an increasingly important role in different areas of engineering. For most of existing methods, such as genetic programming (GP), the convergence speed might be too slow for large scale problems with a large number of variables. It has become the bottleneck of GP for practical applications. Fortunately, in many applications, the target models are separable in some sense. In this paper, we analyze different types of separability of some real-world engineering equations and establish a mathematical model of generalized separable system (GS system). In order to get the structure of the GS system, a multilevel block building (MBB) algorithm is proposed, in which the target model is decomposed into a number of blocks, further into minimal blocks and factors. Compare to the conventional GP, MBB can make large reductions to the search space. This makes MBB capable of modeling a complex system. The minimal blocks and factors are optimized and assembled with a global optimization search engine, low dimensional simplex evolution (LDSE). An extensive study between the proposed MBB and a state-of-the-art data-driven fitting tool, Eureqa, has been presented with several man-made problems, as well as some real-world problems. Test results indicate that the proposed method is more effective and efficient under all the investigated cases.


title: Proxy-SU(3) symmetry in heavy deformed nuclei

abstract: Background: Microscopic calculations of heavy nuclei face considerable difficulties due to the sizes of the matrices that need to be solved. Various approximation schemes have been invoked, for example by truncating the spaces, imposing seniority limits, or appealing to various symmetry schemes such as pseudo-SU(3). This paper proposes a new symmetry scheme also based on SU(3). This proxy-SU(3) can be applied to well-deformed nuclei, is simple to use, and can yield analytic predictions. Purpose: To present the new scheme and its microscopic motivation, and to test it using a Nilsson model calculation with the original shell model orbits and with the new proxy set. Method: We invoke an approximate, analytic, treatment of the Nilsson model, that allows the above vetting and yet is also transparent in understanding the approximations involved in the new proxy-SU(3). Results: It is found that the new scheme yields a Nilsson diagram for well-deformed nuclei that is very close to the original Nilsson diagram. The specific levels of approximation in the new scheme are also shown, for each major shell. Conclusions: The new proxy-SU(3) scheme is a good approximation to the full set of orbits in a major shell. Being able to replace a complex shell model calculation with a symmetry-based description now opens up the possibility to predict many properties of nuclei analytically and often in a parameter-free way. The new scheme works best for heavier nuclei, precisely where full microscopic calculations are most challenged. Some cases in which the new scheme can be used, often analytically, to make specific predictions, are shown in a subsequent paper.


title: Propensity score-based estimators with multiple error-prone covariates

abstract: Propensity score methods are an important tool to help reduce confounding in non-experimental studies. Most propensity score methods assume that covariates are measured without error. However, covariates are often measured with error, which leads to biased causal effect estimates if the true underlying covariates are the actual confounders. Although some studies have investigated the impact of a single mismeasured covariate on estimating a causal effect and proposed methods for handling the measurement error, almost no work exists investigating the case where multiple covariates are mismeasured. In this paper, we examine the consequences of multiple error-prone covariates when estimating causal effects using propensity score-based estimators via extensive simulation studies and real data analyses. We find that causal effect estimates are less biased when the propensity score model includes mismeasured covariates whose true underlying values are strongly correlated with each other. However, when the measurement \emph{errors} are correlated with each other, additional bias is introduced. In addition, it is beneficial to include correctly measured auxiliary variables that are correlated with confounders whose true underlying values are mismeasured in the propensity score model.


title: Generalized hydrodynamics of active polar suspensions

abstract: We utilize a generalized Irving-Kirkwood procedure to derive the hydrodynamic equations of an active matter suspension with internal structure and driven by internal torque. The internal structure and torque of the active Brownian particles give rise to a balance law for internal angular momentum density, making the hydrodynamic description a polar theory of continuum mechanics. We derive exact microscopic expressions for the stress tensor, couple stress tensor, internal energy density, and heat flux vector. Unlike passive matter, the symmetry of the stress tensor is broken explicitly due to active internal torque and the antisymmetric component drives the internal angular momentum density. These results provide a molecular basis to understand the transport characteristics and collectively provide a strategy to develop the theory of linear irreversible thermodynamics of active matter.


title: \textit{Ab Initio} Study of the Magnetic Behavior of Metal Hydrides: A Comparison with the Slater-Pauling Curve

abstract: We investigated the magnetic behavior of metal hydrides FeH$_{x}$, CoH$_{x}$ and NiH$_{x}$ for several concentrations of hydrogen ($x$) by using Density Functional Theory calculations. Several structural phases of the metallic host: bcc ($\alpha$), fcc ($\gamma$), hcp ($\varepsilon$), dhcp ($\varepsilon'$), tetragonal structure for FeH$_{x}$ and $\varepsilon$-$\gamma$ phases for CoH$_{x}$, were studied. We found that for CoH$_{x}$ and NiH$_{x}$ the magnetic moment ($m$) decreases regardless the concentration $x$. However, for FeH$_{x}$ systems, $m$ increases or decreases depending on the variation in $x$. In order to find a general trend for these changes of $m$ in magnetic metal hydrides, we compare our results with the Slater-Pauling curve for ferromagnetic metallic binary alloys. It is found that the $m$ of metal hydrides made of Fe, Co and Ni fits the shape of the Slater-Pauling curve as a function of $x$. Our results indicate that there are two main effects that determine the $m$ value due to hydrogenation: an increase of volume causes $m$ to increase, and the addition of an extra electron to the metal always causes it to decrease. We discuss these behaviors in detail.


title: The orbits of subdwarf-B + main-sequence binaries: III. The period - eccentricity distribution

abstract: The predicted orbital-period distribution of the subdwarf-B (sdB) population is bi-modal with a peak at short (< 10 days) and long (> 500 days) periods. Observationally, many short-period sdB systems are known, but only few wide sdB binaries have been studied in detail. Based on a long-term monitoring program the wide sdB sample has been increased, finding an unexpected correlation between the eccentricity and period. In this article we present the orbital solution and spectral analysis of four new systems, BD-7.5977, EC11031-1348, TYC2084-448-1 and TYC3871-835-1, and update the orbital solution of PG1104+243. Using the whole sample of wide sdBs, we aim at finding possible correlations between orbital and spectral properties, with as goal improving theoretical models of Roche-lobe overflow. High-resolution spectra were obtained to determine the radial velocities of both the sdB and MS components. Surface gravities and temperatures of both component were derived from photometric spectral-energy distributions. Spectral parameters of the cool companion were verified using the GSSP code. Furthermore the amount of accreted mass was estimated. Orbital parameters matching the earlier observed period-eccentricity relation were found for three systems, while TYC 2084-448-1 is found to have a lower eccentricity than expected from the period-eccentricity trend indicated by the other systems. Based on new observations, the orbit of PG 1104+243 has a small but significant eccentricity of 0.04 $\pm$ 0.02, matching other systems with similar periods. Furthermore, a correlation between accreted mass and orbital period was found, as well as a possible relation between the initial mass-ratio and the final period-eccentricity. The wide sdB-binary sample shows interesting possible correlations between orbital and spectral properties. However, a larger sample is necessary to statistically validate them.


title: Statistical Physics Models of Belief Dynamics: Theory and Empirical Tests

abstract: We build simple computational models of belief dynamics within the framework of discrete-spin statistical physics models, and explore how suitable they are for understanding and predicting real-world belief change on both the individual and group levels. We find that accurate modeling of real-world patterns requires attending to social interaction rules that people use, network structures in which they are embedded, distributions of initial beliefs and intrinsic preferences, and the relative importance of social information and intrinsic preferences. We demonstrate that these model parameters can be constrained by empirical measurement, and the resulting models can be used to investigate the mechanisms underlying belief dynamics in actual societies. We use data from two longitudinal studies of belief change, one on 80~individuals living in an MIT dorm during the 2008 presidential election season, and another on 94~participants recruited from Mechanical Turk during the 2016 presidential election primary season. We find that simple statistical physics-based models contain predictive value for real-world belief dynamics and enable empirical tests of different assumptions about the underlying network structure and the social interaction rules.


title: Revisiting the foundations of the quantum theory of atoms in molecules: Some open problems

abstract: In a series of papers in the last ten years, various aspects of the mathematical foundations of the quantum theory of atoms in molecules have been considered by this author and his coworkers in some details. Although these considerations answered part of the questions raised by some critics on the mathematical foundations of the quantum theory of atoms in molecules, however, new mathematical problems also emerged during these studies that were reviewed elsewhere [Int. J. Quantum Chem. 111, 4497 (2011)]. Beyond mathematical subtleties of the formalism that were the original motivation for initial exchanges and disputes, the questions raised by critics had a constructive effect and prompted the author to propose a novel extension of the theory, now called the multi-component quantum theory of atoms in molecules [Theor. Chem. Acc. 132, 1365 (2013)]. Taking this background into account, in this paper a new set of open problems is put forward that the author believes proper answers to these questions, may open new doors for future theoretical developments of the quantum theory of atoms in molecules. Accordingly, rather than emphasizing on the rigorous mathematical formulation, the practical motivations behind proposing these questions are discussed in detail and the relevant literature are reviewed while when possible, evidence and routes to answers are also provided. The author hopes that proposing these open questions as a compact package may motivate more mathematically oriented people to participate in future developments of the quantum theory of atoms in molecules and its multi-component version.


title: Meta-Learning for Resampling Recommendation Systems

abstract: One possible approach to tackle the class imbalance in classification tasks is to resample a training dataset, i.e., to drop some of its elements or to synthesize new ones. There exist several widely-used resampling methods. Recent research showed that the choice of resampling method significantly affects the quality of classification, which raises resampling selection problem. Exhaustive search for optimal resampling is time-consuming and hence it is of limited use. In this paper, we describe an alternative approach to the resampling selection. We follow the meta-learning concept to build resampling recommendation systems, i.e., algorithms recommending resampling for datasets on the basis of their properties.


title: How Retrocausality Helps

abstract: It has become increasingly apparent that a number of perplexing issues associated with the interpretation of quantum mechanics are more easily resolved once the notion of retrocausality is introduced. The aim here is to list and discuss various examples where a clear explanation has become available via this approach. In so doing, the intention is to highlight that this direction of research deserves more attention than it presently receives.


title: Sound Event Detection Using Spatial Features and Convolutional Recurrent Neural Network

abstract: This paper proposes to use low-level spatial features extracted from multichannel audio for sound event detection. We extend the convolutional recurrent neural network to handle more than one type of these multichannel features by learning from each of them separately in the initial stages. We show that instead of concatenating the features of each channel into a single feature vector the network learns sound events in multichannel audio better when they are presented as separate layers of a volume. Using the proposed spatial features over monaural features on the same network gives an absolute F-score improvement of 6.1% on the publicly available TUT-SED 2016 dataset and 2.7% on the TUT-SED 2009 dataset that is fifteen times larger.


title: Stacked Convolutional and Recurrent Neural Networks for Music Emotion Recognition

abstract: This paper studies the emotion recognition from musical tracks in the 2-dimensional valence-arousal (V-A) emotional space. We propose a method based on convolutional (CNN) and recurrent neural networks (RNN), having significantly fewer parameters compared with the state-of-the-art method for the same task. We utilize one CNN layer followed by two branches of RNNs trained separately for arousal and valence. The method was evaluated using the 'MediaEval2015 emotion in music' dataset. We achieved an RMSE of 0.202 for arousal and 0.268 for valence, which is the best result reported on this dataset.


title: Sound Event Detection in Multichannel Audio Using Spatial and Harmonic Features

abstract: In this paper, we propose the use of spatial and harmonic features in combination with long short term memory (LSTM) recurrent neural network (RNN) for automatic sound event detection (SED) task. Real life sound recordings typically have many overlapping sound events, making it hard to recognize with just mono channel audio. Human listeners have been successfully recognizing the mixture of overlapping sound events using pitch cues and exploiting the stereo (multichannel) audio signal available at their ears to spatially localize these events. Traditionally SED systems have only been using mono channel audio, motivated by the human listener we propose to extend them to use multichannel audio. The proposed SED system is compared against the state of the art mono channel method on the development subset of TUT sound events detection 2016 database. The usage of spatial and harmonic features are shown to improve the performance of SED.


title: Chiral Flat Bands: Existence, Engineering and Stability

abstract: We study flat bands in bipartite tight-binding networks with discrete translational invariance. Chiral flat bands with chiral symmetry eigenenergy E = 0 and host compact localized eigenstates for finite range hopping. For a bipartite network with a majority sublattice chiral flat bands emerge. We present a simple generating principle of chiral flat band networks and as a showcase add to the previously observed cases a number of new potentially realizable chiral flat bands in various lattice dimensions. Chiral symmetry respecting network perturbations - including disorder and synthetic magnetic fields - preserve both the flatband and the modified compact localized states. Chiral flatbands are spectrally protected by gaps, and pseudogaps in the presence of disorder due to Griffiths effects.


title: Generative-Discriminative Variational Model for Visual Recognition

abstract: The paradigm shift from shallow classifiers with hand-crafted features to end-to-end trainable deep learning models has shown significant improvements on supervised learning tasks. Despite the promising power of deep neural networks (DNN), how to alleviate overfitting during training has been a research topic of interest. In this paper, we present a Generative-Discriminative Variational Model (GDVM) for visual classification, in which we introduce a latent variable inferred from inputs for exhibiting generative abilities towards prediction. In other words, our GDVM casts the supervised learning task as a generative learning process, with data discrimination to be jointly exploited for improved classification. In our experiments, we consider the tasks of multi-class classification, multi-label classification, and zero-shot learning. We show that our GDVM performs favorably against the baselines or recent generative DNN models.


title: Adiabatic and nonadiabatic spin torques induced by spin-triplet supercurrent

abstract: We study spin transfer torques induced by a spin-triplet supercurrent in a magnet with the superconducting proximity effect. By a perturbative approach, we show that spin-triplet correlations realize new types of torques, which are analogous to the adiabatic and non-adiabatic ($\beta$) torques, without extrinsic spin-flip scattering. Remarkable advantages compared to conventional spin-transfer torques are highlighted in domain wall manipulation. Oscillatory motions of a domain wall do not occur for a small Gilbert damping, and the threshold current density to drive its motion becomes zero in the absence of extrinsic pinning potentials due to the nonadiabatic torque controlled by the triplet correlations.


title: Crystalline Soda Can Metamaterial exhibiting Graphene-like Dispersion at subwavelength scale

abstract: Graphene, a honeycomb lattice of carbon atoms ruled by tight-binding interaction, exhibits extraordinary electronic properties due to the presence of Dirac cones within its band structure. These intriguing singularities have naturally motivated the discovery of their classical analogues. In this work, we present a general and direct procedure to reproduce the peculiar physics of graphene within a very simple acoustic metamaterial: a double lattice of soda cans resonant at two different frequencies. The first triangular sub-lattice generates a bandgap at low frequency, which induces a tight-binding coupling between the resonant defects of the second Honeycomb one, hence allowing us to obtain a graphene-like band structure. We prove the relevance of this approach by showing that both numerical and experimental dispersion relations exhibit the requested Dirac cone. We also demonstrate the straightforward monitoring of the coupling strength within the crystal of resonant defects. This work shows that crystalline metamaterials are very promising candidates to investigate tantalizing solid-state physics phenomena with classical waves.


title: Heavy neutrinos from gluon fusion

abstract: Heavy neutrinos, a key prediction of many standard model extensions, remain some of the most searched-for objects at collider experiments. In this context, we revisit the premise that the gluon fusion production mechanism, $gg \to Z^*/h^* \to N\nu_\ell$, is phenomenologically irrelevant at the CERN LHC and report the impact of soft gluon corrections to the production cross section. We resum threshold logarithms up to next-to-next-to-next-to-leading logarithmic accuracy (N$^3$LL), thus capturing the dominant contributions to the inclusive cross section up to next-to-next-to-leading order (N$^2$LO). For $m_N > 150$ GeV and collider energies $\sqrt{s} = 7 - 100$ TeV, corrections to the Born rates span $+160$ to $+260\%$. At $\sqrt{s}$=14 TeV, the resummed channel is roughly equal in size to the widely-believed-to-be-dominant charged current Drell-Yan process and overtakes it outright at $\sqrt{s} \gtrsim 20-25$ TeV. Results are independent of the precise nature/mixing of $N$ and hold generically for other low-scale seesaws. Findings are also expected to hold for other exotic leptons and broken axial-vector currents, particularly as the $Z^*$ contribution identically reduces to that of a pseudoscalar.


title: On Algorithmic Universality in F-theory Compactifications

abstract: We study universality of geometric gauge sectors in the string landscape in the context of F-theory compactifications. A finite time construction algorithm is presented for $\frac43 \times 2.96 \times 10^{755}$ F-theory geometries that are connected by a network of topological transitions in a connected moduli space. High probability geometric assumptions uncover universal structures in the ensemble without explicitly constructing it. For example, non-Higgsable clusters of seven-branes with intricate gauge sectors occur with probability above $1-1.01\times 10^{-755}$, and the geometric gauge group rank is above $160$ with probability $.999995$. In the latter case there are at least $10$ $E_8$ factors, the structure of which fixes the gauge groups on certain nearby seven-branes. Visible sectors may arise from $E_6$ or $SU(3)$ seven-branes, which occur in certain random samples with probability $\simeq 1/200$.


title: Non-thermal photons and direct photodissociation of H2, HD and HeH+ in the chemistry of the primordial Universe

abstract: Non-thermal photons deriving from radiative transitions among the internal ladder of atoms and molecules are an important source of photons in addition to thermal and stellar sources in many astrophysical environments. In the present work the calculation of reaction rates for the direct photodissociation of some molecules relevant in early Universe chemistry is presented; in particular, the calculations include non-thermal photons deriving from the recombination of primordial hydrogen and helium atoms for the cases of H2, HD and HeH+. New effects on the fractional abundances of chemical species are investigated and the fits for the HeH+ photodissociation rates by thermal photons are provided.


title: Galactic googly: the rotation-metallicity bias in the inner stellar halo of the Milky Way

abstract: The first and second moments of stellar velocities encode important information about the formation history of the Galactic halo. However, due to the lack of tangential motion and inaccurate distances of the halo stars, the velocity moments in the Galactic halo have largely remained "known unknowns". Fortunately, our off-centric position within the Galaxy allows us to estimate these moments in the galacto-centric frame using the observed radial velocities of the stars alone. We use these velocities coupled with the Hierarchical Bayesian scheme, which allows easy marginalisation over the missing data (the proper-motion, and uncertainty-free distance and line-of-sight velocity), to measure the velocity dispersions, orbital anisotropy ($\beta$) and streaming motion ($v_{\rm rot}$) of the halo main-sequence turn-off (MSTO) and K-giant (KG) stars in the inner stellar halo (r $\lesssim 15$ kpc). We study the metallicity bias in kinematics of the halo stars and observe that the comparatively metal-rich ([Fe/H]$>-1.4$) and the metal-poor ([Fe/H]$\leq - 1.4$) MSTO samples show a clear systematic difference in $v_{\rm rot} \sim 20-40$ km s$^{-1}$, depending on how restrictive the spatial cuts to cull the disk contamination are. The bias is also detected in KG samples but with less certainty. Both MSTO and KG populations suggest that the inner stellar halo of the Galaxy is radially biased i.e. $\sigma_r>\sigma_\theta$ or $\sigma_\phi$ and $\beta \simeq 0.5$. The apparent metallicity contrariety in the rotation velocity among the halo sub-populations supports the co-existence of multiple populations in the galactic halo that may have formed through distinct formation scenarios, i.e. in-situ versus accretion.


title: GENESIS: New Self-Consistent Models of Exoplanetary Spectra

abstract: We are entering the era of high-precision and high-resolution spectroscopy of exoplanets. Such observations herald the need for robust self-consistent spectral models of exoplanetary atmospheres to investigate intricate atmospheric processes and to make observable predictions. Spectral models of plane-parallel exoplanetary atmospheres exist, mostly adapted from other astrophysical applications, with different levels of sophistication and accuracy. There is a growing need for a new generation of models custom-built for exoplanets and incorporating state-of-the-art numerical methods and opacities. The present work is a step in this direction. Here we introduce GENESIS, a plane-parallel, self-consistent, line-by-line exoplanetary atmospheric modelling code which includes (a) formal solution of radiative transfer using the Feautrier method, (b) radiative-convective equilibrium with temperature correction based on the Rybicki linearisation scheme, (c) latest absorption cross sections, and (d) internal flux and external irradiation, under the assumptions of hydrostatic equilibrium, local thermodynamic equilibrium and thermochemical equilibrium. We demonstrate the code here with cloud-free models of giant exoplanetary atmospheres over a range of equilibrium temperatures, metallicities, C/O ratios and spanning non-irradiated and irradiated planets, with and without thermal inversions. We provide the community with theoretical emergent spectra and pressure-temperature profiles over this range, along with those for several known hot Jupiters. The code can generate self-consistent spectra at high resolution and has the potential to be integrated into general circulation and non-equilibrium chemistry models as it is optimised for efficiency and convergence. GENESIS paves the way for high-fidelity remote sensing of exoplanetary atmospheres at high resolution with current and upcoming observations.


title: Constraining the Physical State of the Hot Gas Halos in NGC 4649 and NGC 5846

abstract: We present results of a joint \textit{Chandra}/\textit{XMM-Newton} analysis of the early-type galaxies NGC 4649 and NGC 5846 aimed at investigating differences between mass profiles derived from X-ray data and those from optical data, to probe the state of the hot ISM in these galaxies. If the hot ISM is at a given radius in hydrostatic equilibrium (HE) the X-ray data can be used to measure the total enclosed mass of the galaxy. Differences from optically-derived mass distributions therefore yield information about departures from HE in the hot halos. The X-ray mass profiles in different angular sectors of NGC 4649 are generally smooth with no significant azimuthal asymmetries within \(12\) kpc. Extrapolation of these profiles beyond this scale yields results consistent with the optical estimate. However, in the central region (\(r < 3\) kpc) the X-ray data underpredict the enclosed mass, when compared with the optical mass profiles. Consistent with previous results we estimate a non-thermal pressure component accounting for \(30\%\) of the gas pressure, likely linked to nuclear activity. In NGC 5846 the X-ray mass profiles show significant azimuthal asymmetries, especially in the NE direction. Comparison with optical mass profiles in this direction suggests significant departures from HE, consistent with bulk gas compression and decompression due to sloshing on \(\sim 15\) kpc scales; this effect disappears in the NW direction where the emission is smooth and extended. In this sector we find consistent X-ray and optical mass profiles, suggesting that the hot halo is not responding to strong non-gravitational forces.


title: Quantum spin chains with multiple dynamics

abstract: Many-body systems with multiple emergent time scales arise in various contexts, including classical critical systems, correlated quantum materials, and ultra-cold atoms. We investigate such non-trivial quantum dynamics in a new setting: a spin-1 bilinear-biquadratic chain. It has a solvable entangled groundstate, but a gapless excitation spectrum that is poorly understood. By using large-scale DMRG simulations, we find that the lowest excitations have a dynamical exponent $z$ that varies from 2 to 3.2 as we vary a coupling in the Hamiltonian. We find an additional gapless mode with a continuously varying exponent $2\leq z <2.7$, which establishes the presence of multiple dynamics. In order to explain these striking properties, we construct a continuum wavefunction for the groundstate, which correctly describes the correlations and entanglement properties. We also give a continuum parent Hamiltonian, but show that additional ingredients are needed to capture the excitations of the chain. By using an exact mapping to the non-equilibrium dynamics of a classical spin chain, we find that the large dynamical exponent is due to subdiffusive spin motion. Finally, we discuss the connections to other spin chains and to a family of quantum critical models in 2d.


title: Pairing versus phase coherence of doped holes in distinct quantum spin backgrounds

abstract: We examine the pairing structure of holes injected into two \emph{distinct} spin backgrounds: a short-range antiferromagnetic phase versus a symmetry protected topological phase. Based on density matrix renormalization group (DMRG) simulation, we find that although there is a strong binding between two holes in both phases, \emph{phase fluctuations} can significantly influence the pair-pair correlation depending on the spin-spin correlation in the background. Here the phase fluctuation is identified as an intrinsic string operator nonlocally controlled by the spins. We show that while the pairing amplitude is generally large, the coherent Cooper pairing can be substantially weakened by the phase fluctuation in the symmetry-protected topological phase, in contrast to the short-range antiferromagnetic phase. It provides an example of a non-BCS mechanism for pairing, in which the paring phase coherence is determined by the underlying spin state self-consistently, bearing an interesting resemblance to the pseudogap physics in the cuprate.


title: Charge density wave instabilities of type-II Weyl semimetals in a strong magnetic field

abstract: Shortly after the discovery of Weyl semimetals properties related to the topology of their bulk band structure have been observed, e.g. signatures of the chiral anomaly and Fermi arc surface states. These essentially single particle phenomena are well understood but whether interesting many-body effects due to interactions arise in Weyl systems remains much less explored. Here, we investigate the effect of interactions in a microscopic model of a type-II Weyl semimetal in a strong magnetic field. We identify a charge density wave (CDW) instability even for weak interactions stemming from the emergent nesting properties of the type-II Weyl Landau level dispersion. We map out the dependence of this CDW on magnetic field strength. Remarkably, as a function of decreasing temperature a cascade of CDW transitions emerges and we predict characteristic signatures for experiments.


title: Galaxy groups in the low-redshift Universe

abstract: We apply a halo-based group finder to four large redshift surveys, the 2MRS, 6dFGS, SDSS and 2dFGRS, to construct group catalogs in the low-redshift Universe. The group finder is based on that of Yang et al. but with an improved halo mass assignment so that it can be applied uniformly to various redshift surveys of galaxies. Halo masses are assigned to groups according to proxies based on the stellar mass/luminosity of member galaxies. The performances of the group finder in grouping galaxies according to common halos and in halo mass assignments are tested using realistic mock samples constructed from hydrodynamical simulations and empirical models of galaxy occupation in dark matter halos. Our group finder finds $\sim 94\%$ of the correct true member galaxies for $90-95\%$ of the groups in the mock samples; the halo masses assigned by the group finder are un-biased with respect to the true halo masses, and have a typical uncertainty of $\sim0.2\,{\rm dex}$. The properties of group catalogs constructed from the observational samples are described and compared with other similar catalogs in the literature.


title: The sequence of Compton dominance in blazars based on data from WISE and Fermi/LAT

abstract: The two-component broad-band spectral energy distributions of blazars were suggested to form a sequence in which (1) the peak frequency of the low-energy (synchrotron) component nu_syn is anticorrelated with the synchrotron luminosity L_syn, and (2) the luminosity ratio of the high-energy (inverse Compton) to synchrotron components q = L_IC / L_syn (Compton dominance) increases with L_syn from the BL Lac objects (BL Lacs) to the flat spectrum radio quasars (FSRQs). The Compton dominance parameter is an important probe of plasma magnetisation in the blazar zones within relativistic jets. We investigate a sample of blazars detected by WISE in the mid-infrared (MIR) band and by Fermi/LAT in the GeV gamma-ray band, with the focus on the distribution of luminosities and photon indices. Our findings are the following: (1) the MIR photon index Gamma_W12 is a useful probe of the blazar sequence, with the exception of low-luminosity BL Lacs that are most likely contaminated by their host galaxies (L_W1 ~ 10^44 erg/s and Gamma_W12 < 1); (2) Gamma_W12 is correlated with the gamma-ray photon index Gamma_1-100GeV, with the MIR luminosity L_W1, and with the Fermi/WISE Compton dominance q_FW = L_1GeV/L_W1; (3) a clean separation between FSRQs and BL Lacs can be seen in the parameter space of Gamma_W12 and q_FW; (4) the observed distribution of MIR luminosity L_W1 vs. Compton dominance q_FW for the entire sample of blazars can be modeled as a sequence of lepto-magnetic jet powers in the range (log_10 P_eB) in [42:45] with the preference for sub-equipartition magnetic fields P_B/P_e in [0.05:1], assuming fixed bulk Lorentz factor Gamma_j = 15, fixed jet opening angle Gamma_j Theta_j = 0.3, fixed radiative efficiency of jet electrons epsilon_em = 50%, and that external radiation luminosity scales like L_ext propto (P_eB)^(1.6) (parameter degeneracies are discussed).


title: Spin-Orbit Misalignment of Merging Black-Hole Binaries with Tertiary Companions

abstract: We study the effect of external companion on the orbital and spin evolution of merging black-hole (BH) binaries. An sufficiently close by and inclined companion can excite Lidov-Kozai (LK) eccentricity oscillations in the binary, thereby shortening its merger time. During such LK-enhanced orbital decay, the spin axis of the BH generally exhibits chaotic evolution, leading to a wide range ($0^\circ$-$180^\circ$) of final spin-orbit misalignment angle from an initially aligned configuration. For systems that do not experience eccentricity excitation, only modest ($\lesssim 20^\circ$) spin-orbit misalignment can be produced, and we derive an analytic expression for the final misalignment using the principle of adiabatic invariance. The spin-orbit misalignment directly impacts the gravitational waveform, and can be used to constrain the formation scenarios of BH binaries and dynamical influences of external companions.


title: A GMOS-N IFU study of the central H II region in the blue compact dwarf galaxy NGC 4449: Kinematics, Nebular Metallicity and Star-Formation

abstract: We use integral field spectroscopic (IFS) observations from the Gemini North Multi-Object Spectrograph (GMOS-N) to study the central H II region in a nearby blue compact dwarf (BCD) galaxy NGC 4449. The IFS data enable us to explore the variation of physical and chemical conditions of the star-forming region and the surrounding gas on spatial scales as small as 5.5 pc. Our kinematical analysis shows possible signatures of shock ionisation and shell structures in the surroundings of the star-forming region. The metallicity maps of the region, created using direct T$_e$ and indirect strong line methods (R$_{23}$, O3N2 and N2), do not show any chemical variation. From the integrated spectrum of the central H II region, we find a metallicity of 12 + log(O/H) = 7.88 $\pm$ 0.14 ($\sim$ 0.15$^{+0.06}_{-0.04}$ Z$_{\odot}$) using the direct method. Comparing the central H II region metallicity derived here with those of H II regions throughout this galaxy from previous studies, we find evidence of increasing metallicity with distance from the central nucleus. Such chemical inhomogeneities can be due to several mechanisms, including gas-loss via supernova blowout, galactic winds, or metal-poor gas accretion. However, we find that the localised area of decreased metallicity aligns spatially with the peak of star-forming activity in the galaxy, suggesting that gas-accretion may be at play here. Spatially-resolved IFS data for the entire galaxy is required to confirm the metallicity inhomogeneity found in this study, and determine its possible cause.


title: The intrinsic characteristics of galaxies on the SFR-stellar mass plane at 1.2<z<4: I. the correlation between stellar age, central density and position relative to the main sequence

abstract: We use the deep CANDELS observations in the GOODS North and South fields to revisit the correlations between stellar mass ($M_*$), star--formation rate (SFR) and morphology, and to introduce a fourth dimension, the mass-weighted stellar age, in galaxies at $1.2<z<4$. We do this by making new measures of $M_*$, $SFR$, and stellar age thanks to an improved SED fitting procedure that allows various star formation history for each galaxy. Like others, we find that the slope of the Main Sequence (MS) of star formation in the $(M_*;SFR)$ plane bends at high mass. We observe clear morphological differences among galaxies across the MS, which also correlate with stellar age. At all redshifts, galaxies that are quenching or quenched, and thus old, have high $\Sigma_1$ (the projected density within the central 1 kpc), while younger, star-forming galaxies span a much broader range of $\Sigma_1$, which includes the high values observed for quenched galaxies, but also extends to much lower values. As galaxies age and quench, the stellar age and the dispersion of $\Sigma_1$ for fixed values of $M_{*}$ shows two different regimes, one, at the low--mass end, where quenching might be driven by causes external to the galaxies; the other, at the high--mass end, where quenching is driven by internal causes, very likely the mass given the low scatter of $\Sigma_1$ (mass quenching). We suggest that the monotonic increase of central density as galaxies grow is one manifestation of a more general phenomenon of structural transformation that galaxies undergo as they evolve.


title: Emergent Dark Energy in Classical Channel Gravity with Matter

abstract: Motivated by the recent increased interest in energy non-conserving models in cosmology, we extend the analysis of the cosmological consequences of the Classical Channel Model of Gravity (CCG). This model is based on the classical-quantum interaction between a test particle and a metric (classical) and results in a theory with a modified Wheeler-deWitt equation that in turn leads to non conservation of energy. We show that CCG applied to a cosmological scenario with primordial matter leads to an emergent dark fluid that at late times behaves as a curvature term in the Friedmann equations, showing that the late time behaviour is always dominated by the vacuum solutions. We discuss possible observational constraints for this model and that, in its current formulation, CCG eludes any meaningful constraints from current observations.


title: Unification of gauge and Yukawa couplings

abstract: The unification of gauge and top Yukawa couplings is an attractive feature of gauge-Higgs unification models in extra-dimensions. This feature is usually considered difficult to obtain based on simple group theory analyses. We reconsider a minimal toy model including the renormalisation group running at one loop. Our results show that the gauge couplings unify asymptotically at high energies, and that this may result from the presence of an UV fixed point. The Yukawa coupling in our toy model is enhanced at low energies, showing that a genuine unification of gauge and Yukawa couplings may be achieved.


title: Holomorphic Classical Limit for Spin Effects in Gravitational and Electromagnetic Scattering

abstract: We provide universal expressions for the classical piece of the amplitude given by the graviton/photon exchange between massive particles of arbitrary spin, at both tree and one loop level. In the gravitational case this leads to higher order terms in the post-Newtonian expansion, which have been previously used in the binary inspiral problem. The expressions are obtained in terms of a contour integral that computes the Leading Singularity, which was recently shown to encode the relevant information up to one loop. The classical limit is performed along a holomorphic trajectory in the space of kinematics, such that the leading order is enough to extract arbitrarily high multipole corrections. These multipole interactions are given in terms of a recently proposed representation for massive particles of any spin by Arkani-Hamed et al. This explicitly shows universality of the multipole interactions in the effective potential with respect to the spin of the scattered particles. We perform the explicit match to standard EFT operators for $S=\frac{1}{2}$ and $S=1$. As a natural byproduct we obtain the classical pieces up to one loop for the bending of light.


title: Evolution of vorticity and enstrophy in the intracluster medium

abstract: Turbulence generated by large-scale motions during structure formation affects the evolution of the thermal and non-thermal components of the intracluster medium. As enstrophy is a measure of the magnitude of vorticity, we study the generation and evolution of turbulence by analysing the Lagrangian history of enstrophy. For this purpose we combine cosmological simulations carried out with the ENZO-code with our Lagrangian post-processing tool CRaTer. This way we are able to quantify the individual source terms of enstrophy in the course of the accretion of groups onto galaxy clusters. Here we focus on the redshift range from $z=1$ to $z=0$. Finally, we measure the rate of dissipation of turbulence and estimate the resulting amplification of intracluster magnetic fields. We find that compressive and baroclinic motions are the main sources of enstrophy, while stretching motions and dissipation affect most of the ensuing enstrophy evolution. The rate of turbulent dissipation is able to sustain the amplification of intracluster magnetic fields to observed levels.


title: Ionised Accretion Discs in Active Galactic Nuclei: The Effects of a Lamppost with a Variable Height

abstract: The X-ray emitting corona irradiates and ionises the surface of the inner accretion disc in Active Galactic Nuclei (AGNs). The ionisation parameter of the inner disc at a radius $r$ from the black hole, $\xi(r)$, can be used to infer information about the location of the corona. Here, a new formula is derived that predicts $\xi(r,h)$ for a disc irradiated by a X-ray source at a height $h$ above the black hole symmetry axis (i.e., a lamppost geometry). The equation is independent of the black hole mass and the X-ray spectrum, and accounts for the effects of gravitational light-bending on the ionisation state and a variable coronal dissipation factor. We predict a strong ionisation gradient across the surface of the inner disc that depends on the black hole spin and lamppost height. For a fixed $h$, the ionisation parameter is also expected to increase as $\lambda^3$, where $\lambda$ is the observed bolometric Eddington ratio of the AGN. Comparing this formula to the observed $\xi$-$\lambda$ relationship for Mrk 335 yields $h \propto \lambda^{0.5-0.6}$, supporting the view of a dynamic X-ray corona in AGNs.


title: The primordial environment of supermassive black holes (II): deep Y and J band images around the z=6.3 quasar SDSS J1030+0524

abstract: Many cosmological studies predict that early supermassive black holes (SMBHs) can only form in the most massive dark matter halos embedded within large scale structures marked by galaxy over-densities that may extend up to 10 physical Mpc. This scenario, however, has not been confirmed observationally, as the search for galaxy over-densities around high-z quasars has returned conflicting results. The field around the z=6.28 quasar SDSSJ1030+0524 (J1030) is unique for multi-band coverage and represents an excellent data legacy for studying the environment around a primordial SMBH. In this paper we present wide-area (25x25 arcmin), Y- and J-band imaging of the J1030 field obtained with the near infrared camera WIRCam at the Canada-France-Hawaii Telescope (CFHT). We built source catalogues in the Y- and J-band, and matched those with our photometric catalogue in the r, z, i bands presented in Morselli et al. (2014). We used these new infrared data together with H and K and Spitzer/IRAC data to refine our selection of Lyman Break Galaxies (LBGs), extending our selection criteria to galaxies in the range 25.2<zAB<25.7. We selected 21 robust high-z candidates in the J1030 field with photometric redshift around 6 and colors i-z>=1.3. We found a significant asymmetry in the distribution of the high-z galaxies in J1030, supporting the existence of a coherent large-scale structure around the quasar. We compared our results with those of Bowler et al. (2015), who adopted similar LBGs selection criteria, and estimated an over-density of galaxies in the field of delta = 2.4, which is significant at >4 sigma. The over-density value and its significance are higher than those found in Morselli et al. (2014), and we interpret this as evidence of an improved LBG selection.


title: Spatially Resolved Thermometry of Resistive Memory Devices

abstract: The operation of resistive and phase-change memory (RRAM and PCM) is controlled by highly localized self-heating effects, yet detailed studies of their temperature are rare due to challenges of nanoscale thermometry. Here we show that the combination of Raman thermometry and scanning thermal microscopy (SThM) can enable such measurements with high spatial resolution. We report temperature-dependent Raman spectra of HfO$_2$, TiO$_2$ and Ge$_2$Sb$_2$Te$_5$ (GST) films, and demonstrate direct measurements of temperature profiles in lateral PCM devices. Our measurements reveal that electrical and thermal interfaces dominate the operation of such devices, uncovering a thermal boundary resistance of 30 m$^2$K$^{-1}$GW$^{-1}$ at GST-SiO$_2$ interfaces and an effective thermopower 350 $\mu$V/K at GST-Pt interfaces. We also discuss possible pathways to apply Raman thermometry and SThM techniques to nanoscale and vertical resistive memory devices.


title: Multifield Stochastic Particle Production: Beyond a Maximum Entropy Ansatz

abstract: We explore non-adiabatic particle production for $N_{\rm f}$ coupled scalar fields in a time-dependent background with stochastically varying effective masses, cross-couplings and intervals between interactions. Under the assumption of weak scattering per interaction, we provide a framework for calculating the typical particle production rates after a large number of interactions. After setting up the framework, for analytic tractability, we consider interactions (effective masses and cross couplings) characterized by series of Dirac-delta functions in time with amplitudes and locations drawn from different distributions. Without assuming that the fields are statistically equivalent, we present closed form results (up to quadratures) for the asymptotic particle production rates for the $N_{\rm f}=1$ and $N_{\rm f}=2$ cases. We also present results for the general $N_{\rm f} >2$ case, but with more restrictive assumptions. We find agreement between our analytic results and direct numerical calculations of the total occupation number of the produced particles, with departures that can be explained in terms of violation of our assumptions. We elucidate the precise connection between the maximum entropy ansatz (MEA) used in Amin & Baumann (2015) and the underlying statistical distribution of the self and cross couplings. We provide and justify a simple to use (MEA-inspired) expression for the particle production rate, which agrees with our more detailed treatment when the parameters characterizing the effective mass and cross-couplings between fields are all comparable to each other. However, deviations are seen when some parameters differ significantly from others. We show that such deviations become negligible for a broad range of parameters when $N_{\rm f}\gg 1$.


title: Evaluation of the operatorial Q-system for non-compact super spin chains

abstract: We present an approach to evaluate the full operatorial Q-system of all $\mathfrak{u}(p,q|r+s)$-invariant spin chains with representations of Jordan-Schwinger type. In particular, this includes the super spin chain of planar $\mathcal{N}=4$ super Yang-Mills theory at one loop in the presence of a diagonal twist. Our method is based on the oscillator construction of Q-operators. The Q-operators are built as traces over Lax operators which are degenerate solutions of the Yang-Baxter equation. For non-compact representations these Lax operators may contain multiple infinite sums that conceal the form of the resulting functions. We determine these infinite sums and calculate the matrix elements of the lowest level Q-operators. Transforming the Lax operators corresponding to the Q-operators into a representation involving only finite sums allows us to take the supertrace and to obtain the explicit form of the Q-operators in terms of finite matrices for a given magnon sector. Imposing the functional relations, we then bootstrap the other Q-operators from those of the lowest level. We exemplify this approach for non-compact spin $-s$ spin chains and apply it to $\mathcal{N}=4$ SYM at the one-loop level using the BMN vacuum as an example.


title: Analytic predictions for nuclear shapes, the prolate dominance and the prolate-oblate shape transition in the proxy-SU(3) model

abstract: Using a new approximate analytic parameter-free proxy-SU(3) scheme, we make simple predictions of shape observables for deformed nuclei, namely gamma and beta deformation variables, the global feature of prolate dominance and the locus of the prolate-oblate shape transition. The predictions are compared with empirical results.


title: An exactly solvable quench protocol for integrable spin models

abstract: Quantum quenches in continuum field theory across critical points are known to display different scaling behaviours in different regimes of the quench rate. We extend these results to integrable lattice models such as the transverse field Ising model on a one-dimensional chain and the Kitaev model on a two-dimensional honeycomb lattice using a nonlinear quench protocol which allows for exact analytical solutions of the dynamics. Our quench protocol starts with a finite mass gap at early times and crosses a critical point or a critical region, and we study the behaviour of one point functions of the quenched operator at the critical point or in the critical region as a function of the quench rate. For quench rates slow compared to the initial mass gap, we find the expected Kibble-Zurek scaling. In contrast, for rates fast compared to the mass gap, but slow compared to the inverse lattice spacing, we find scaling behaviour similar to smooth fast continuum quenches. For quench rates of the same order of the lattice scale, the one point function saturates as a function of the rate, approaching the results of an abrupt quench. The presence of an extended critical surface in the Kitaev model leads to a variety of scaling exponents depending on the starting point and on the time where the operator is measured. We discuss the role of the amplitude of the quench in determining the extent of the slow (Kibble-Zurek) and fast quench regimes, and the onset of the saturation.


title: Critical behavior of the black hole / black string transition

abstract: We numerically construct static localized black holes in five and six spacetime dimensions which are solutions to Einstein's vacuum field equations with one compact periodic dimension. In particular, we investigate the critical regime in which the poles of the localized black hole are about to merge. A well adapted multi-domain pseudo-spectral scheme provides us with accurate results and enables us to investigate the phase diagram of those localized solutions within the critical regime, which goes far beyond previous results. We find that in this regime the phase diagram possesses a spiral structure adapting to the one recently found for non-uniform black strings. When approaching the common endpoint of both phases, the behavior of physical quantities is described by complex critical exponents giving rise to a discrete scaling symmetry. The numerically obtained values of the critical exponents agree remarkably well with those derived from the double-cone metric.


title: Efficient method for estimating the number of communities in a network

abstract: While there exist a wide range of effective methods for community detection in networks, most of them require one to know in advance how many communities one is looking for. Here we present a method for estimating the number of communities in a network using a combination of Bayesian inference with a novel prior and an efficient Monte Carlo sampling scheme. We test the method extensively on both real and computer-generated networks, showing that it performs accurately and consistently, even in cases where groups are widely varying in size or structure.


title: Adaptively enriched coarse space for the discontinuous Galerkin multiscale problems

abstract: In this paper, we propose a two-level overlapping additive Schwarz domain decomposition preconditioner for the symmetric interior penalty discontinuous Galerkin method for the second order elliptic boundary value problem with highly heterogeneous coefficients. A specific feature of this preconditioner is that it is based on adaptively enriching its coarse space with functions created by solving generalized eigenvalue problems on thin patches covering the subdomain interfaces. It is shown that the condition number of the underlined preconditioned system is independent of the contrast if an adequate number of functions are used to enrich the coarse space. Numerical results are provided to confirm this claim.


title: Improving Variational Auto-Encoders using convex combination linear Inverse Autoregressive Flow

abstract: In this paper, we propose a new volume-preserving flow and show that it performs similarly to the linear general normalizing flow. The idea is to enrich a linear Inverse Autoregressive Flow by introducing multiple lower-triangular matrices with ones on the diagonal and combining them using a convex combination. In the experimental studies on MNIST and Histopathology data we show that the proposed approach outperforms other volume-preserving flows and is competitive with current state-of-the-art linear normalizing flow.


title: Eigenvector Centrality Distribution for Characterization of Protein Allosteric Pathways

abstract: Determining the principal energy pathways for allosteric communication in biomolecules, that occur as a result of thermal motion, remains challenging due to the intrinsic complexity of the systems involved. Graph theory provides an approach for making sense of such complexity, where allosteric proteins can be represented as networks of amino acids. In this work, we establish the eigenvector centrality metric in terms of the mutual information, as a mean of elucidating the allosteric mechanism that regulates the enzymatic activity of proteins. Moreover, we propose a strategy to characterize the range of the physical interactions that underlie the allosteric process. In particular, the well known enzyme, imidazol glycerol phosphate synthase (IGPS), is utilized to test the proposed methodology. The eigenvector centrality measurement successfully describes the allosteric pathways of IGPS, and allows to pinpoint key amino acids in terms of their relevance in the momentum transfer process. The resulting insight can be utilized for refining the control of IGPS activity, widening the scope for its engineering. Furthermore, we propose a new centrality metric quantifying the relevance of the surroundings of each residue. In addition, the proposed technique is validated against experimental solution NMR measurements yielding fully consistent results. Overall, the methodologies proposed in the present work constitute a powerful and cost effective strategy to gain insight on the allosteric mechanism of proteins.


title: Delay Optimal Scheduling for Chunked Random Linear Network Coding Broadcast

abstract: We study the broadcast transmission of a single file to an arbitrary number of receivers using Random Linear Network Coding (RLNC) in a network with unreliable channels. Due to the increased computational complexity of the decoding process (especially for large files) we apply chunked RLNC (i.e. RLNC is applied within non-overlapping subsets of the file). In our work we show the optimality of the Least Received (LR) batch scheduling policy (which was introduced in our prior work) with regards to the expected file transfer completion time. Furthermore, we refine some of our earlier results, namely the expected file transfer completion time of the LR policy and the minimum achievable coding window size in the case of a user defined delay constraint. Finally, we experimentally evaluate a modification of the LR policy in a more realistic system setting with reduced feedback from the receivers.


title: Wigner's theorem on Grassmann spaces

abstract: Wigner's celebrated theorem, which is particularly important in the mathematical foundations of quantum mechanics, states that every bijective transformation on the set of all rank-one projections of a complex Hilbert space which preserves the transition probability is induced by a unitary or an antiunitary operator. This vital theorem has been generalised in various ways by several scientists. In 2001, Moln\'ar provided a natural generalisation, namely, he provided a characterisation of (not necessarily bijective) maps which act on the Grassmann space of all rank-$n$ projections and leave the system of Jordan principal angles invariant (see [20] and [17]). In this paper we give a very natural joint generalisation of Wigner's and Moln\'ar's theorems, namely, we prove a characterisation of all (not necessarily bijective) transformations on the Grassmann space which fix the quantity $\mathrm{tr} PQ$ (i.e.~the sum of the squares of cosines of principal angles) for every pair of rank-$n$ projections $P$ and $Q$.


title: Initial state qqg correlations as a background for the Chiral Magnetic Effect in collision of small systems

abstract: Motivated by understanding the background to Chiral Magnetic Effect in proton-nucleus collisions from first principles, we compute the three particle correlation in the projectile wave function. We extract the correlations between two quarks and one gluon in the framework of the Color Glass Condensate. This is related to the same-charge correlation of the conventional observable for the Chiral Magnetic Effect. We show that there are two different contributions to this correlation function. One contribution is rapidity-independent and as such can be identified with the pedestal; while the other displays rather strong rapidity dependence. The pedestal contribution and the rapidity-dependent contribution at large rapidity separation between the two quarks result in the negative same charge correlations, while at small rapidity separation the second contribution changes sign. We argue that the computed initial state correlations might be partially responsible for the experimentally observed signal in proton-nucleus collisions.


title: CoMaL Tracking: Tracking Points at the Object Boundaries

abstract: Traditional point tracking algorithms such as the KLT use local 2D information aggregation for feature detection and tracking, due to which their performance degrades at the object boundaries that separate multiple objects. Recently, CoMaL Features have been proposed that handle such a case. However, they proposed a simple tracking framework where the points are re-detected in each frame and matched. This is inefficient and may also lose many points that are not re-detected in the next frame. We propose a novel tracking algorithm to accurately and efficiently track CoMaL points. For this, the level line segment associated with the CoMaL points is matched to MSER segments in the next frame using shape-based matching and the matches are further filtered using texture-based matching. Experiments show improvements over a simple re-detect-and-match framework as well as KLT in terms of speed/accuracy on different real-world applications, especially at the object boundaries.


title: Low-shot learning with large-scale diffusion

abstract: This paper considers the problem of inferring image labels from images when only a few annotated examples are available at training time. This setup is often referred to as low-shot learning, where a standard approach is to re-train the last few layers of a convolutional neural network learned on separate classes for which training examples are abundant. We consider a semi-supervised setting based on a large collection of images to support label propagation. This is possible by leveraging the recent advances on large-scale similarity graph construction. We show that despite its conceptual simplicity, scaling label propagation up to hundred millions of images leads to state of the art accuracy in the low-shot learning regime.


title: K-polynomials of type A quiver orbit closures and lacing diagrams

abstract: This article contains an overview of the author's joint work with Allen Knutson and Jenna Rajchgot on $K$-polynomials of orbit closures for type $A$ quivers. It is written to an audience interested in interactions between representations of algebras, algebraic geometry, and commutative algebra. A few open problems resulting from the work are also explained.


title: Observational evidence that positive and negative AGN feedback depends on galaxy mass and jet power

abstract: Several studies support the existence of a link between the AGN and star formation activity. Radio jets have been argued to be an ideal mechanism for direct interaction between the AGN and the host galaxy. A drawback of previous surveys of AGN is that they are fundamentally limited by the degeneracy between redshift and luminosity in flux-density limited samples. To overcome this limitation, we present far-infrared Herschel observations of 74 radio-loud quasars (RLQs), 72 radio-quiet quasars (RQQs) and 27 radio galaxies (RGs), selected at 0.9<z<1.1 which span over two decades in optical luminosity. By decoupling luminosity from evolutionary effects, we investigate how the star formation rate (SFR) depends on AGN luminosity, radio-loudness and orientation. We find that: 1) the SFR shows a weak correlation with the bolometric luminosity for all AGN sub-samples, 2) the RLQs show a SFR excess of about a factor of 1.4 compared to the RQQs, matched in terms of black hole mass and bolometric luminosity, suggesting that either positive radio-jet feedback or radio AGN triggering are linked to star-formation triggering and 3) RGs have lower SFRs by a factor of 2.5 than the RLQ sub-sample with the same BH mass and bolometric luminosity. We suggest that there is some jet power threshold at which radio-jet feedback switches from enhancing star formation (by compressing gas) to suppressing it (by ejecting gas). This threshold depends on both galaxy mass and jet power.


title: Injective chromatic number of outerplanar graphs

abstract: An injective coloring of a graph is a vertex coloring where two vertices with common neighbor receive distinct colors. The minimum integer $k$ that $G$ has a $k-$injective coloring is called injective chromatic number of $G$ and denoted by $\chi_i(G)$. In this paper, the injective chromatic number of outerplanar graphs with maximum degree $\Delta$ and girth $g$ is studied. It is shown that for every outerplanar graph, $\chi_i(G)\leq \Delta+2$, and this bound is tight. Then, it is proved that for outerplanar graphs with $\Delta=3$, $\chi_i(G)\leq \Delta+1$ and the bound is tight for outerplanar graphs of girth three and $4$. Finally, it is proved that, the injective chromatic number of $2-$connected outerplanar graphs with $\Delta=3$, $g\geq 6$ and $\Delta\geq 4$, $g\geq 4$ is equal to $\Delta$.


title: The Semi-Hooperon: Gamma-ray and anti-proton excesses in the Galactic Center

abstract: A puzzling excess in gamma-rays at GeV energies has been observed in the center of our galaxy using Fermi-LAT data. Its origin is still unknown, but it is well fitted by Weakly Interacting Massive Particles (WIMPs) annihilations into quarks with a cross section around $10^{-26}{\rm cm^3 s^{-1}}$ with masses of $20-50$~GeV, scenario which is promptly revisited. An excess favoring similar WIMP properties has also been seen in anti-protons with AMS-02 data potentially coming from the Galactic Center as well. In this work, we explore the possibility of fitting these excesses in terms of semi-annihilating dark matter, dubbed as semi-Hooperon, with the process ${\rm WIMP\, WIMP \rightarrow WIMP\, X}$ being responsible for the gamma-ray excess, where X=h,Z. An interesting feature of semi-annihilations is the change in the relic density prediction compared to the standard case, and the possibility to alleviate stringent limits stemming from direct detection searches. Moreover, we discuss which models might give rise to a successful semi-Hooperon setup in the context of $\mathcal{Z}_3$,$\mathcal{Z}_4$ and extra "dark" gauge symmetries.


title: Learning to Extract Semantic Structure from Documents Using Multimodal Fully Convolutional Neural Network

abstract: We present an end-to-end, multimodal, fully convolutional network for extracting semantic structures from document images. We consider document semantic structure extraction as a pixel-wise segmentation task, and propose a unified model that classifies pixels based not only on their visual appearance, as in the traditional page segmentation task, but also on the content of underlying text. Moreover, we propose an efficient synthetic document generation process that we use to generate pretraining data for our network. Once the network is trained on a large set of synthetic documents, we fine-tune the network on unlabeled real documents using a semi-supervised approach. We systematically study the optimum network architecture and show that both our multimodal approach and the synthetic data pretraining significantly boost the performance.


abstract: Testing the simplifying assumption in high-dimensional vine copulas is a difficult task. Tests must be based on estimated observations and check constraints on high-dimensional distributions. So far, corresponding tests have been limited to single conditional copulas with a low-dimensional set of conditioning variables. We propose a novel testing procedure that is computationally feasible for high-dimensional data sets and that exhibits a power that decreases only slightly with the dimension. By discretizing the support of the conditioning variables and incorporating a penalty in the test statistic, we mitigate the curse of dimensionality by looking for the possibly strongest deviation from the simplifying assumption. The use of a decision tree renders the test computationally feasible for large dimensions. We derive the asymptotic distribution of the test and analyze its finite sample performance in an extensive simulation study. An application of the test to four real data sets is provided.


title: Flux growth in a horizontal configuration: an analogue to vapor transport growth

abstract: Flux growth of single crystals is normally performed in a vertical configuration with an upright refractory container holding the flux melt. At high temperatures, flux dissolves the charge forming a homogeneous solution before nucleation and growth of crystals take place under proper supersaturation generated by cooling or evaporating the flux. In this work, we report flux growth in a horizontal configuration with a temperature gradient along the horizontal axis: a liquid transport growth analogous to the vapor transport technique. In a typical liquid transport growth, the charge is kept at the hot end of the refractory container and the flux melt dissolves the charge and transfers it to the cold end. Once the concentration of charge is above the solubility limit at the cold end, the thermodynamically stable phase nucleates and grows. Compared to the vertical flux growth, the liquid transport growth can provide a large quantity of crystals in a single growth since the charge/flux ratio is not limited by the solubility limit at the growth temperature. This technique is complementary to the vertical flux growth and can be considered when a large amount of crystals are needed but the yield from the conventional vertical flux growth is limited. We applied this technique to the growth of IrSb$_3$, Mo$_3$Sb$_7$, MnBi from self flux, and the growth of FeSe, CrTe$_3$, NiPSe$_3$, FePSe$_3$, and InCuP$_2$S$_6$ from a halide flux.


title: Classification of finite $p$-groups by the size of their Schur multipliers

abstract: Let $d(G)$ be the minimum number of elements required to generated a group $G.$ For a group $G $ of order $p^n$ with derived subgroup of order $ p^k $ and $d(G) = d,$ we knew the order of the Schur multiplier of $G$ is bounded by $ p^{\frac{1}{2}(d-1)(n-k+2)+1}. $ In the current paper, we find the structure of all $p$-groups that attains the mentioned bound. Moreover, we show that all of them are capable.


title: Hirota bilinear equations for Painlev\'e transcendents

abstract: We present some observations on the tau-function for the fourth Painlev\'e equation. By considering a Hirota bilinear equation of order four for this tau-function, we describe the general form of the Taylor expansion around an arbitrary movable zero. The corresponding Taylor series for the tau-functions of the first and second Painlev\'e equations, as well as that for the Weierstrass sigma function, arise naturally as special cases, by setting certain parameters to zero.


title: Active Learning for Structured Prediction from Partially Labeled Data

abstract: We propose a general purpose active learning algorithm for structured prediction, gathering labeled data for training a model that outputs a set of related labels for an image or video. Active learning starts with a limited initial training set, then iterates querying a user for labels on unlabeled data and retraining the model. We propose a novel algorithm for selecting data for labeling, choosing examples to maximize expected information gain based on belief propagation inference. This is a general purpose method and can be applied to a variety of tasks or models. As a specific example we demonstrate this framework for learning to recognize human actions and group activities in video sequences. Experiments show that our proposed algorithm outperforms previous active learning methods and can achieve accuracy comparable to fully supervised methods while utilizing significantly less labeled data.


title: Some methods for constructing new operator monotone functions from old ones

abstract: We observe that if f is a continuous function on an interval I and x_0 \in I, then f is operator monotone if and only if the function (f(x) - f(x_0)/(x - x_0) is strongly operator convex. Then starting with an operator monotone function f_0, we construct a strongly operator convex function f_1, an (ordinary) operator convex function f_2, and then a new operator monotone function f_3. The process can be continued to obtain an infinite sequence which cycles between the three classes of functions. We also describe two other constructions, similar in spirit. We prove two lemmas which enable a treatment of those aspects of strong operator convexity needed for this paper which is more elementary than previous treatments. And we discuss the functions phi such that the composite phi \circ f is operator convex or strongly operator convex whenever f is strongly operator convex.


title: Energy-Efficient Hybrid Stochastic-Binary Neural Networks for Near-Sensor Computing

abstract: Recent advances in neural networks (NNs) exhibit unprecedented success at transforming large, unstructured data streams into compact higher-level semantic information for tasks such as handwriting recognition, image classification, and speech recognition. Ideally, systems would employ near-sensor computation to execute these tasks at sensor endpoints to maximize data reduction and minimize data movement. However, near- sensor computing presents its own set of challenges such as operating power constraints, energy budgets, and communication bandwidth capacities. In this paper, we propose a stochastic- binary hybrid design which splits the computation between the stochastic and binary domains for near-sensor NN applications. In addition, our design uses a new stochastic adder and multiplier that are significantly more accurate than existing adders and multipliers. We also show that retraining the binary portion of the NN computation can compensate for precision losses introduced by shorter stochastic bit-streams, allowing faster run times at minimal accuracy losses. Our evaluation shows that our hybrid stochastic-binary design can achieve 9.8x energy efficiency savings, and application-level accuracies within 0.05% compared to conventional all-binary designs.


title: Imaging Ferroelectric Domains via Charge Gradient Microscopy Enhanced by Principal Component Analysis

abstract: Local domain structures of ferroelectrics have been studied extensively using various modes of scanning probes at the nanoscale, including piezoresponse force microscopy (PFM) and Kelvin probe force microscopy (KPFM), though none of these techniques measure the polarization directly, and the fast formation kinetics of domains and screening charges cannot be captured by these quasi-static measurements. In this study, we used charge gradient microscopy (CGM) to image ferroelectric domains of lithium niobate based on current measured during fast scanning, and applied principal component analysis (PCA) to enhance the signal-to-noise ratio of noisy raw data. We found that the CGM signal increases linearly with the scan speed while decreases with the temperature under power-law, consistent with proposed imaging mechanisms of scraping and refilling of surface charges within domains, and polarization change across domain wall. We then, based on CGM mappings, estimated the spontaneous polarization and the density of surface charges with order of magnitude agreement with literature data. The study demonstrates that PCA is a powerful method in imaging analysis of scanning probe microscopy (SPM), with which quantitative analysis of noisy raw data becomes possible.


title: Khovanov spectra for tangles

abstract: We define stable homotopy refinements of Khovanov's arc algebras and tangle invariants.


title: Quantum Simulation of Competing Orders with Fermions in Quantum Optical Lattices

abstract: Ultracold Fermi atoms confined in optical lattices coupled to quantized modes of an optical cavity are an ideal scenario to engineer quantum simulators in the strongly interacting regime. The system has both short range and cavity induced long range interactions. We propose such a scheme to investigate the coexistence of superfluid pairing, density order and quantum domains having an- tiferromagnetic or density order in the Hubbard model in a high finesse optical cavity at T = 0. We demonstrate that those phases can be accessed by properly tuning the linear polarizer of an external pump beam via the cavity back-action effect, while modulating the system doping. This allows emulate the typical scenarios of analog strongly correlated electronic systems.

categories: cond-mat.quant-gas cond-mat.str-el cond-mat.supr-con physics.atom-ph quant-ph


title: Resampling Strategy in Sequential Monte Carlo for Constrained Sampling Problems

abstract: Sequential Monte Carlo (SMC) methods are a class of Monte Carlo methods that are used to obtain random samples of a high dimensional random variable in a sequential fashion. Many problems encountered in applications often involve different types of constraints. These constraints can make the problem much more challenging. In this paper, we formulate a general framework of using SMC for constrained sampling problems based on forward and backward pilot resampling strategies. We review some existing methods under the framework and develop several new algorithms. It is noted that all information observed or imposed on the underlying system can be viewed as constraints. Hence the approach outlined in this paper can be useful in many applications.


title: Connes' Embedding Problem and Winning Strategies for Quantum XOR Games

abstract: We consider quantum XOR games, defined in [11], from the perspective of unitary correlations defined in [7]. We show that Connes' embedding problem has a positive answer if and only if every quantum XOR game has entanglement bias equal to the commuting bias. In particular, the embedding problem is equivalent to determining whether every quantum XOR game $G$ with a winning strategy in the commuting model also has a winning strategy in the approximate finite-dimensional model.


title: On the postulation of lines and a fat line

abstract: In this note we show that the union of $r$ general lines and one fat line in ${\mathbb P}^3$ imposes independent conditions on forms of sufficiently high degree $d$, where the bound on $d$ is independent of the number of lines. This extends former results of Hartshorne and Hirschowitz on unions of general lines, and of Aladpoosh on unions of general lines and one double line.


title: Recognizing difference quotients of real functions

abstract: For a real function $f:[0,1]\to\mathbb{R}$, the difference quotient of $f$ is the function of two real variables $\operatorname{DQ}_f(a,b)=\dfrac{f(b)-f(a)}{b-a}$, which we view as defined on the triangle $\mathcal{T}=\{(a,b):0\leq a<b\leq1\}$. In this paper we investigate how to determine whether a given function of two variables $H(a,b)$ is the difference quotient of some real function $f(x)$. We develop three independent methods for recognizing such a function $H$ as a difference quotient, and corresponding methods for recovering the underlying function $f$ from $H$.


title: Direct Experimental Observation of the Gas Filamentation Effect using a Two-bunch X-ray FEL Beam

abstract: We report the experimental observation of the filamentation effect in gas devices designed for X-ray Free-electron Lasers. The measurements were carried out at the Linac Coherent Light Source on the X-ray Correlation Spectroscopy (XCS) instrument using a Two-bunch FEL beam at 6.5 keV with 122.5 ns separation passing through an Argon gas cell. The relative intensities of the two pulses of the Two-bunch beam were measured, after and before the gas cell, from the X-ray scattering off thin targets by using fast diodes with sufficient temporal resolution. It was found that the after-to-before ratio of the intensities of the second pulse was consistently higher than that of the first pulse, revealing lower effective attenuation of the gas cell due to the heating and subsequent gas density reduction in the beam path by the first pulse. This measurement is important in guiding the design and/or mitigating the adverse effect in gas devices for high repetition-rate FELs such as the LCLS-II and the European XFEL or other future high repetition-rate upgrade to existing FEL facilities


title: Sparse Wavelet Estimation in Quantile Regression with Multiple Functional Predictors

abstract: In this manuscript, we study quantile regression in partial functional linear model where response is scalar and predictors include both scalars and multiple functions. Wavelet basis are adopted to better approximate functional slopes while effectively detect local features. The sparse group lasso penalty is imposed to select important functional predictors while capture shared information among them. The estimation problem can be reformulated into a standard second-order cone program and then solved by an interior point method. We also give a novel algorithm by using alternating direction method of multipliers (ADMM) which was recently employed by many researchers in solving penalized quantile regression problems. The asymptotic properties such as the convergence rate and prediction error bound have been established. Simulations and a real data from ADHD-200 fMRI data are investigated to show the superiority of our proposed method.


title: Is a pole type singularity an alternative to inflation?

abstract: In this paper, we apply a method of reducing the dynamics of FRW cosmological models with the barotropic form of the equation of state to the dynamical system of the Newtonian type to detect the finite scale factor singularities and the finite-time singularities. In this approach all information concerning the dynamics of the system is contained in a diagram of the potential function $V(a)$ of the scale factor. Singularities of the finite scale factor manifest by poles of the potential function. In our approach the different types of singularities are represented by critical exponents in the power-law approximation of the potential. The classification can be given in terms of these exponents. We have found that the pole singularity can mimick an inflation epoch. We demonstrate that the cosmological singularities can be investigated in terms of the critical exponents of the potential function of the cosmological dynamical systems. We assume the general form of the model contains matter and some kind of dark energy which is parameterized by the potential. We distinguish singularities (by ansatz about the Lagrangian) of the pole type with the inflation and demonstrate that such a singularity can appear in the past.


title: Shadows of a Closed Curve

abstract: A shadow of a geometric object $A$ in a given direction $v$ is the orthogonal projection of $A$ on the hyperplane orthogonal to $v$. We show that any topological embedding of a circle into Euclidean $d$-space can have at most two shadows that are simple paths in linearly independent directions. The proof is topological and uses an analog of basic properties of degree of maps on a circle to relations on a circle. This extends a previous result which dealt with the case $d=3$.


title: Cover time for random walks on arbitrary complex networks

abstract: We present an analytical method for computing the mean cover time of a random walk process on arbitrary, complex networks. The cover time is defined as the time a random walker requires to visit every node in the network at least once. This quantity is particularly important for random search processes and target localization in network topologies. Based on the global mean first passage time of target nodes we derive an estimate for the cumulative distribution function of the cover time based on first passage time statistics. We show that our result can be applied to various model networks, including Erd\H{o}s-R\'enyi and Barab\'asi-Albert networks, as well as various real-world networks. Our results reveal an intimate link between first passage and cover time statistics in networks in which structurally induced temporal correlations decay quickly and offer a computationally efficient way for estimating cover times in network related applications.


title: Quasi-particles, thermodynamic consistency and the gap equation

abstract: The thermodynamic properties of superconducting electrons are usually studied by means of the quasi-particles distribution; but in this approach, the ground state energy and the dependence of the chemical potential on the electron density cannot be determined. In order to solve these problems, the thermodynamic potentials are derived by means of the Bogoliubov-Valatin formalism. The thermodynamic potentials can be obtained by computing the free energy of a gas of quasi-particles, whose energy spectrum is conditional on the gap function. However, the nontrivial dependence of the gap on the temperature jeopardises the validity of the standard thermodynamic relations. In this article it is shown how the thermodynamic consistency (i.e. the validity of the Maxwell relations) is recovered, and the correction terms to the quasi-particles potentials are computed. It is shown that the Bogoliubov-Valatin transformation avoids the problem of the thermodynamic consistency of the quasi-particle approach; in facts, the correct identification of the variables, which are associated with the quasi-particles, leads to a precise calculation of the quasi-particles vacuum energy and of the dependence of the chemical potential on the electron density. The stationarity condition for the grand potential coincides with the gap equation, which guarantees the thermodynamic consistency. The expressions of various thermodynamic potentials, as functions of the (T,V,N) variables, are produced in the low temperature limit; as a final check, a rederivation of the condensation energy is presented.

categories: cond-mat.supr-con cond-mat.stat-mech math-ph math.MP


title: Mind Your Credit: Assessing the Health of the Ripple Credit Network

abstract: The Ripple credit network has emerged as a payment backbone with key advantages for financial institutions and the remittance industry. Its path-based IOweYou (IOU) settlements across different (crypto)currencies conceptually distinguishes the Ripple blockchain from cryptocurrencies, and makes it highly suitable to an orthogonal yet vast set of applications in the remittance world for cross-border transactions and beyond. This work studies the structure and evolution of the Ripple network since its inception, and investigates its vulnerability to devilry attacks that affect the credit of linnet users' wallets. We find that about 13M USD are at risk in the current Ripple network due to inappropriate configuration of the rippling flag on credit links, facilitating undesired redistribution of credit across those links. Although the Ripple network has grown around a few highly connected hub (gateway) wallets that constitute the network's core and provide high liquidity to users, such a credit link distribution results in a user base of around 112,000 wallets that can be financially isolated by as few as 10 highly connected gateway wallets. Indeed, today about 4.9M USD cannot be withdrawn by their owners from the Ripple network due to PayRoutes, a gateway tagged as faulty by the Ripple community. Finally, we observe that stale exchange offers pose a real problem, and exchanges (market makers) have not always been vigilant about periodically updating their exchange offers according to current real-world exchange rates. For example, stale offers were used by 84 Ripple wallets to gain more than 4.5M USD from mid-July to mid-August 2017. Our findings should prompt the Ripple community to improve the health of the network by educating its users on increasing their connectivity, and by appropriately maintaining the credit limits, rippling flags, and exchange offers on their credit links.


title: New Factor Pairs for Factorizations of Lambert Series Generating Functions

abstract: We prove several new variants of the Lambert series factorization theorem established in the first article "Generating special arithmetic functions by Lambert series factorizations" by Merca and Schmidt (2017). Several characteristic examples of our new results are presented in the article to motivate the formulations of the generalized factorization theorems. Applications of these new factorization results include new identities involving the Euler partition function and the generalized sum-of-divisors functions, the M\"obius function, Euler's totient function, the Liouville lambda function, von Mangoldt's lambda function, and the Jordan totient function.


title: Galerkin approximations for the optimal control of nonlinear delay differential equations

abstract: Optimal control problems of nonlinear delay differential equations (DDEs) are considered for which we propose a general Galerkin approximation scheme built from Koornwinder polynomials. Error estimates for the resulting Galerkin-Koornwinder approximations to the optimal control and the value function, are derived for a broad class of cost functionals and nonlinear DDEs. The approach is illustrated on a delayed logistic equation set not far away from its Hopf bifurcation point in the parameter space. In this case, we show that low-dimensional controls for a standard quadratic cost functional can be efficiently computed from Galerkin-Koornwinder approximations to reduce at a nearly optimal cost the oscillation amplitude displayed by the DDE's solution. Optimal controls computed from the Pontryagin's maximum principle (PMP) and the Hamilton-Jacobi-Bellman equation (HJB) associated with the corresponding ODE systems, are shown to provide numerical solutions in good agreement. It is finally argued that the value function computed from the corresponding reduced HJB equation provides a good approximation of that obtained from the full HJB equation.