title: Some new formulas for Appell series over finite fields

abstract: In 1987 Greene introduced the notion of the finite field analogue of hypergeometric series. In this paper we give a finite field analogue of Appell series and obtain some transformation and reduction formulas. We also establish the generating functions for Appell series over finite fields.


title: Directional Total Generalized Variation Regularization

abstract: In inverse problems, prior information and a priori-based regularization techniques play important roles. In this paper, we focus on image restoration problems, especially on restoring images whose texture mainly follow one direction. In order to incorporate the directional information, we propose a new directional total generalized variation (DTGV) functional, which is based on total generalized variation (TGV) by Bredies \textit{et al}. [SIAM J. Imaging Sci., 3 (2010)]. After studying the mathematical properties of DTGV, we utilize it as regularizer and propose the L$^2$-DTGV variational model for solving image restoration problems. Due to the requirement of the directional information in DTGV, we give a direction estimation algorithm, and then apply a primal-dual algorithm to solve the minimization problem. Experimental results show the effectiveness of the proposed method for restoring the directional images. In comparison with isotropic regularizers like total variation and TGV, the improvement of texture preservation and noise removal is significant.


title: Unsupervised Image-to-Image Translation with Generative Adversarial Networks

abstract: It's useful to automatically transform an image from its original form to some synthetic form (style, partial contents, etc.), while keeping the original structure or semantics. We define this requirement as the "image-to-image translation" problem, and propose a general approach to achieve it, based on deep convolutional and conditional generative adversarial networks (GANs), which has gained a phenomenal success to learn mapping images from noise input since 2014. In this work, we develop a two step (unsupervised) learning method to translate images between different domains by using unlabeled images without specifying any correspondence between them, so that to avoid the cost of acquiring labeled data. Compared with prior works, we demonstrated the capacity of generality in our model, by which variance of translations can be conduct by a single type of model. Such capability is desirable in applications like bidirectional translation


title: A VLT-ULTRACAM study of the fast optical quasi-periodic oscillations in the polar V834 Centauri

abstract: Quasi-periodic oscillations (QPOs) of a few seconds have been detected in some Polars, the synchronised subclass of cataclysmic systems containing a strongly magnetised white dwarf (WD) which accretes matter from a red dwarf companion. The QPOs are thought to be related to instabilities of a shock formed in the accretion column, close to the WD photosphere above the impact region. We present optical observations of the polar V834 Centauri performed with the fast ULTRACAM camera mounted on the ESO-VLT simultaneously in three filters (u', He II 4686A, r') to study these oscillations and characterise their properties along the orbit. Fast Fourier transforms and wavelet analysis have been performed and the mean frequency, rms amplitude, and coherence of the QPOs are derived; a detailed inspection of individual pulses has also been performed. The observations confirm the probable ubiquity of the QPOs for this source at all epochs when the source is in a high state, with observed mean amplitude of 2.1 percent (r'), 1.5 percent (He II), and 0.6 percent (u'). Trains of oscillations are clearly observed in the r' light curve and can be mimicked by a superposition of damped sinusoids with various parameters. The QPO energy distribution is comparable to that of the cyclotron flux, consistent for the r' and He II filters but requiring a significant dilution in the u' filter. New 1D hydrodynamical simulations of shock instabilities, adapted to the physical parameters of V834 Cen, can account for the optical QPO amplitude and X-ray upper limit assuming a cross section of the accretion column in the range (4-5) E14 cm2. However, the predicted frequency is larger than the observed one by an order of magnitude. This shortcoming indicates that the QPO generation is more complex than that produced in a homogeneous column and calls for a more realistic 3D treatment of the accretion flow in future modelling.


title: Two-Higgs-doublet model of type-II confronted with the LHC run-I and run-II data

abstract: We examine the parameter space of two-Higgs-doublet model of type-II after imposing the relevant theoretical and experimental constraints from the precision electroweak data, $B$-meson decays, the LHC run-I and run-II data. We find that the searches for Higgs via the $\tau\bar{\tau},~WW,~ZZ,\gamma\gamma,~hh,~hZ,~HZ,~AZ$ channels can give strong constraints on the CP-odd Higgs $A$ and heavy CP-even Higgs $H$, and the parameter space excluded by each channel is respectively carved out in detail. The surviving samples are discussed in two different of regions: (i) In the SM-like coupling region of the 125 GeV Higgs, $m_A$ is allowed to be as low as 350 GeV, and $\tan\beta$ is imposed a strong upper limit. $m_H$ is allowed to be as low as 200 GeV for the proper $\tan\beta$, $\sin(\beta-\alpha)$ and $m_A$, but is required to be larger than 300 GeV for $m_A=700$ GeV. (ii) In the wrong sign Yukawa coupling of the 125 GeV, $\tan\beta$ and $\sin(\beta-\alpha)$ can be imposed upper limits by $b\bar{b}\to A/H\to \tau\bar{\tau}$ channel and lower limits by the $A\to hZ$ channel. $m_A$ and $m_H$ are respectively allowed to be as low as 60 GeV and 200 GeV, but 320 GeV $< m_A<$ 500 GeV is excluded for $m_H=$ 700 GeV.


title: Investigation of optimal control problems governed by a time-dependent Kohn-Sham model

abstract: Many application models in quantum physics and chemistry require to control multi-electron systems to achieve a desired target configuration. This challenging task appears possible in the framework of time-dependent density functional theory (TDDFT) that allows to describe these systems while avoiding the high dimensionality resulting from the multi-particle Schr\"{o}dinger equation. For this purpose, the theory and numerical solution of optimal control problems governed by a Kohn-Sham TDDFT model are investigated, considering different objectives and a bilinear control mechanism. Existence of optimal control solutions and their characterization as solutions to Kohn-Sham TDDFT optimality systems are discussed. To validate this control framework, a time-splitting discretization of the optimality systems and a nonlinear conjugate gradient scheme are implemented. Results of numerical experiments demonstrate the computational capability of the proposed control approach.


title: A stochastic root finding approach: The Homotopy Analysis Method applied to Dyson-Schwinger Equations

abstract: We present the construction and stochastic summation of rooted-tree diagrams, based on the expansion of a root finding algorithm applied to the Dyson-Schwinger equations (DSEs). The mathematical formulation shows superior convergence properties compared to the bold diagrammatic Monte Carlo approach and the developed algorithm allows one to tackle generic high-dimensional integral equations, to avoid the curse of dealing explicitly with high-dimensional objects and to access non-perturbative regimes. The sign problem remains the limiting factor, but it is not found to be worse than in other approaches. We illustrate the method for $\phi^4$ theory but note that it applies in principle to any model.


title: Recursive Marginal Quantization of Higher-Order Schemes

abstract: Quantization techniques have been applied in many challenging finance applications, including pricing claims with path dependence and early exercise features, stochastic optimal control, filtering problems and efficient calibration of large derivative books. Recursive Marginal Quantization of the Euler scheme has recently been proposed as an efficient numerical method for evaluating functionals of solutions of stochastic differential equations. This method involves recursively quantizing the conditional marginals of the discrete-time Euler approximation of the underlying process. By generalizing this approach, we show that it is possible to perform recursive marginal quantization for two higher-order schemes: the Milstein scheme and a simplified weak order 2.0 scheme. As part of this generalization a simple matrix formulation is presented, allowing efficient implementation. We further extend the applicability of recursive marginal quantization by showing how absorption and reflection at the zero boundary may be incorporated, when this is necessary. To illustrate the improved accuracy of the higher order schemes, various computations are performed using geometric Brownian motion and its generalization, the constant elasticity of variance model. For both processes, we show numerical evidence of improved weak order convergence and we compare the marginal distributions implied by the three schemes to the known analytical distributions. By pricing European, Bermudan and Barrier options, further evidence of improved accuracy of the higher order schemes is demonstrated.


title: A Devil's Advocate against Termination of Direct Recursion

abstract: A devil's advocate is one who argues against a claim, not as a committed opponent but in order to determine the validity of the claim. We are interested in a devil's advocate that argues against termination of a program. He does so by producing a maleficent program that can cause the non-termination of the original program. By inspecting and running the malicious program, one may gain insight into the potential reasons for non-termination and produce counterexamples for termination. We introduce our method using the concurrent programming language Constraint Handling Rules (CHR). Like in other declarative languages, non-termination occurs through unbounded recursion. Given a self-recursive rule, we automatically generate one or more devil's rules from it. The construction of the devil's rules is straightforward and involves no guessing. The devil's rules can be simple. For example, they are non-recursive for rules with single recursion. We show that the devil's rules are maximally vicious in the following sense: For any program that contains the self-recursive rule and for any infinite computation through that rule in that program, there is a corresponding infinite computation with the recursive rule and the devil's rules alone. In that case, the malicious rules serve as a finite witness for non-termination. On the other hand, if the devil's rules do not exhibit an infinite computation, the recursive rule is unconditionally terminating. We also identify cases where the static analysis of the devil's rule decides termination or non-termination of the recursive rule.


title: Reconstructing the ideal results of a perturbed analog quantum simulator

abstract: Well-controlled quantum systems can potentially be used as quantum simulators. However, a quantum simulator is inevitably perturbed by coupling to additional degrees of freedom. This constitutes a major roadblock to useful quantum simulations. So far there are only limited means to understand the effect of perturbation on the results of quantum simulation. Here, we present a method which, in certain circumstances, allows for the reconstruction of the ideal result from measurements on a perturbed quantum simulator. We consider extracting the value of the correlator $\langle\hat{O}^i(t) \hat{O}^j(0)\rangle$ from the simulated system, where $\hat{O}^i$ are the operators which couple the system to its environment. The ideal correlator can be straightforwardly reconstructed by using statistical knowledge of the environment, if any $n$-time correlator of operators $\hat{O}^i$ of the ideal system can be written as products of two-time correlators. We give an approach to verify the validity of this assumption experimentally by additional measurements on the perturbed quantum simulator. The proposed method can allow for reliable quantum simulations with systems subjected to environmental noise without adding an overhead to the quantum system.


title: Differential forms for fractal subspaces and finite energy coordinates

abstract: This paper introduces a notion of differential forms on closed, potentially fractal, subsets of Euclidean space by defining pointwise cotangent spaces using the restriction of $C^1$ functions to this set. Aspects of cohomology are developed: it is shown that the differential forms are a Banach algebra and it is possible to integrate these forms along rectifiable paths. These definitions are connected to the theory of differential forms on Dirichlet spaces by considering fractals with finite energy coordinates. In this situation, the $C^1$ differential forms project onto the space of Dirichlet differential forms. Further, it is shown that if the intrinsic metric of a Dirichlet form is a length space, then the image of any rectifiable path through a finite energy coordinate sequence is also rectifiable. The example of the harmonic Sierpinski gasket is worked out in detail.


title: Forming disc galaxies in major mergers II. The central mass concentration problem and a comparison of GADGET3 with GIZMO

abstract: Context: In a series of papers, we study the major merger of two disk galaxies in order to establish whether or not such a merger can produce a disc galaxy. Aims: Our aim here is to describe in detail the technical aspects of our numerical experiments. Methods: We discuss the initial conditions of our major merger, which consist of two protogalaxies on a collision orbit. We show that such merger simulations can produce a non-realistic central mass concentration, and we propose simple, parametric, AGN-like feedback as a solution to this problem. Our AGN-like feedback algorithm is very simple: at each time-step we take all particles whose local volume density is above a given threshold value and increase their temperature to a preset value. We also compare the GADGET3 and GIZMO codes, by applying both of them to the same initial conditions. Results: We show that the evolution of isolated protogalaxies resembles the evolution of disk galaxies, thus arguing that our protogalaxies are well suited for our merger simulations. We demonstrate that the problem with the unphysical central mass concentration in our merger simulations is further aggravated when we increase the resolution. We show that our AGN-like feedback removes this non-physical central mass concentration, and thus allows the formation of realistic bars. Note that our AGN-like feedback mainly affects the central region of a model, without significantly modifying the rest of the galaxy. We demonstrate that, in the context of our kind of simulation, GADGET3 gives results which are very similar to those obtained with the PSPH (density independent SPH) flavor of GIZMO. Moreover, in the examples we tried, the differences between the results of the two flavors of GIZMO, namely PSPH, and MFM (mesh-less algorithm) are similar to and, in some comparisons, larger than the differences between the results of GADGET3 and PSPH.


title: On The Rank Of Congruent Elliptic Curves

abstract: In this paper, $p$ and $q$ are two different odd primes. First, We construct the congruent elliptic curves corresponding to $p$, $2p$, $pq$, and $2pq,$ then, in the cases of congruent numbers, we determine the rank of the corresponding congruent elliptic curves.


title: On the sums of many biquadrates in two different ways

abstract: The beautiful quartic Diophantine equation $A^4+hB^4=C^4+hD^4$, where $h$ is a fixed arbitrary positive integer, has been studied by some mathematicians for many years. Although Choudhry, Gerardin and Piezas presented solutions of this equation for many values of $h$, the solutions were not known for arbitrary positive integer values of $h$. In a separate paper (see the arxiv), the authors completely solved the equation for arbitrary values of $h$, and worked out many examples for different values of $h$, in particular for the values which has not already been given a solution. Our method, give rise to infinitely many solutions and also infinitely many parametric solutions for the equation for arbitrary rational values of $h$. In the present paper, we use the above solutions as well as a simple idea to show that how some numbers can be written as the sums of two, three, four, five, or more biquadrates in two different ways. In particular we give examples for the sums of $2$, $3$, $\cdots$, and $10$, biquadrates expressed in two different ways.


title: Pion Distribution Amplitude and Quasi-Distributions

abstract: We extend our analysis of quasi-distributions onto the pion distribution amplitude. Using the formalism of parton virtuality distribution amplitudes (VDAs), we establish a connection between the pion transverse momentum dependent distribution amplitude (TMDA) $\Psi (x, k_\perp^2)$ and the pion quasi-distribution amplitude (QDA) $Q_\pi (y,p_3)$. We build models for the QDAs from the VDA-based models for soft TMDAs, and analyze the $p_3$ dependence of the resulting QDAs. As there are many models claimed to describe the primordial shape of the pion DA, we present the $p_3$-evolution patterns for models producing some popular proposals: Chernyak-Zhitnitsky, flat and asymptotic DAs. Our results may be used as a guide for future studies of the pion distribution amplitude on the lattice using the quasi-distribution approach.


title: Scattering above energy norm of a focusing size-dependent log energy-supercritical Schrodinger equation with radial data below ground state

abstract: Given $n \in \{ 3,4,5 \}$ and $k > 1$ (resp. $\frac{4}{3} > k > 1$) if $n \in \{ 3,4 \}$ (resp. $n=5$), we prove scattering of the radial $\tilde{H}^{k}:= \dot{H}^{k}(\mathbb{R}^{n}) \cap \dot{H}^{1}(\mathbb{R}^{n})$ solutions of a focusing size-dependent log energy-supercritical Schrodinger equation for critical energies below that of the ground states, and for critical potential energies below that of the ground states.


title: Spectral Projector-Based Graph Fourier Transforms

abstract: The paper presents the graph Fourier transform (GFT) of a signal in terms of its spectral decomposition over the Jordan subspaces of the graph adjacency matrix $A$. This representation is unique and coordinate free, and it leads to unambiguous definition of the spectral components ("harmonics") of a graph signal. This is particularly meaningful when $A$ has repeated eigenvalues, and it is very useful when $A$ is defective or not diagonalizable (as it may be the case with directed graphs). Many real world large sparse graphs have defective adjacency matrices. We present properties of the GFT and show it to satisfy a generalized Parseval inequality and to admit a total variation ordering of the spectral components. We express the GFT in terms of spectral projectors and present an illustrative example for a real world large urban traffic dataset.


title: Strategies for quantum computing molecular energies using the unitary coupled cluster ansatz

abstract: The variational quantum eigensolver (VQE) algorithm combines the ability of quantum computers to efficiently compute expectation values with a classical optimization routine in order to approximate ground state energies of quantum systems. In this paper, we study the application of VQE to the simulation of molecular energies using the unitary coupled cluster (UCC) ansatz. We introduce new strategies to reduce the circuit depth for the implementation of UCC and improve the optimization of the wavefunction based on efficient classical approximations of the cluster amplitudes. Additionally, we propose an analytical method to compute the energy gradient that reduces the sampling cost for gradient estimation by several orders of magnitude compared to numerical gradients. We illustrate our methodology with numerical simulations for a system of four hydrogen atoms that exhibit strong correlation and show that the circuit depth of VQE using a UCC ansatz can be reduced without introducing significant loss of accuracy in the final wavefunctions and energies.


title: Von Neumann Regular Cellular Automata

abstract: For any group $G$ and any set $A$, a cellular automaton (CA) is a transformation of the configuration space $A^G$ defined via a finite memory set and a local function. Let $\text{CA}(G;A)$ be the monoid of all CA over $A^G$. In this paper, we investigate a generalisation of the inverse of a CA from the semigroup-theoretic perspective. An element $\tau \in \text{CA}(G;A)$ is von Neumann regular (or simply regular) if there exists $\sigma \in \text{CA}(G;A)$ such that $\tau \circ \sigma \circ \tau = \tau$ and $\sigma \circ \tau \circ \sigma = \sigma$, where $\circ$ is the composition of functions. Such an element $\sigma$ is called a generalised inverse of $\tau$. The monoid $\text{CA}(G;A)$ itself is regular if all its elements are regular. We establish that $\text{CA}(G;A)$ is regular if and only if $\vert G \vert = 1$ or $\vert A \vert = 1$, and we characterise all regular elements in $\text{CA}(G;A)$ when $G$ and $A$ are both finite. Furthermore, we study regular linear CA when $A= V$ is a vector space over a field $\mathbb{F}$; in particular, we show that every regular linear CA is invertible when $G$ is torsion-free elementary amenable (e.g. when $G=\mathbb{Z}^d, \ d \in \mathbb{N}$) and $V=\mathbb{F}$, and that every linear CA is regular when $V$ is finite-dimensional and $G$ is locally finite with $\text{Char}(\mathbb{F}) \nmid o(g)$ for all $g \in G$.


title: Clique colourings of geometric graphs

abstract: A clique colouring of a graph is a colouring of the vertices such that no maximal clique is monochromatic (ignoring isolated vertices). The least number of colours in such a colouring is the clique chromatic number. Given $n$ points $x_1, \ldots,x_n$ in the plane, and a threshold $r>0$, the corresponding geometric graph has vertex set $\{v_1,\ldots,v_n\}$, and distinct $v_i$ and $v_j$ are adjacent when the Euclidean distance between $x_i$ and $x_j$ is at most $r$. We investigate the clique chromatic number of such graphs. We first show that the clique chromatic number is at most 9 for any geometric graph in the plane, and briefly consider geometric graphs in higher dimensions. Then we study the asymptotic behaviour of the clique chromatic number for the random geometric graph $RG$ in the plane, where $n$ random points are independently and uniformly distributed in a suitable square. We see that as $r$ increases from 0, with high probability the clique chromatic number is 1 for very small $r$, then 2 for small $r$, then at least 3 for larger $r$, and finally drops back to 2.


title: Limited individual attention and online virality of low-quality information

abstract: Social media are massive marketplaces where ideas and news compete for our attention. Previous studies have shown that quality is not a necessary condition for online virality and that knowledge about peer choices can distort the relationship between quality and popularity. However, these results do not explain the viral spread of low-quality information, such as the digital misinformation that threatens our democracy. We investigate quality discrimination in a stylized model of online social network, where individual agents prefer quality information, but have behavioral limitations in managing a heavy flow of information. We measure the relationship between the quality of an idea and its likelihood to become prevalent at the system level. We find that both information overload and limited attention contribute to a degradation in the market's discriminative power. A good tradeoff between discriminative power and diversity of information is possible according to the model. However, calibration with empirical data characterizing information load and finite attention in real social media reveals a weak correlation between quality and popularity of information. In these realistic conditions, the model predicts that high-quality information has little advantage over low-quality information.


title: Classification of the Second Minimal Odd Periodic Orbits in the Sharkovskii Ordering

abstract: This paper presents full classification of second minimal odd periodic orbits of a continuous endomorphisms on the real line. A $(2k+1)$-periodic orbit ($k\geq 3$) is called second minimal for the map $f$, if $2k-1$ is a minimal period of $f$ in the Sharkovskii ordering. We prove that there are $4k-3$ types of second minimal $(2k+1)$-orbits, each characterized with unique cyclic permutation and directed graph of transitions with accuracy up to inverses.


title: Coupling quantum emitters in 2D materials with tapered fibers

abstract: Realization of integrated photonic circuits on a single chip requires controlled manipulation and integration of solid-state quantum emitters with nanophotonic components. Previous works focused on emitters embedded in a three-dimensional crystals -- such as nanodiamonds or quantum dots. In contrast, in this work we demonstrate coupling of a single emitter in a two-dimensional (2D) material, namely hexagonal boron nitride (hBN), with a tapered optical fiber and find a collection efficiency of the system is found to be 10~\%. Furthermore, due to the single dipole character of the emitter, we were able to analyse the angular emission pattern of the coupled system via back focal plane imaging. The good coupling efficiency to the tapered fiber even allows excitation and detection in a fully fiber coupled way yielding a true integrated system. Our results provide evidence of the feasibility to efficiently integrate quantum emitters in 2D materials with photonic structures.


title: Sensitivity of the NEXT-100 detector to neutrinoless double beta decay

abstract: A high pressure xenon gas time projection chamber with electroluminescent amplification (EL HPGXe TPC) searching for the neutrinoless double beta ($0\nu\beta\beta$) decay offers: excellent energy resolution ($0.5-0.7\%$ FWHM at the $Q_{\beta\beta}$), by amplifying the ionization signal with electroluminescent light, and tracking capabilities, as demonstrated by the NEXT collaboration using two kg-scale prototypes. The NEXT collaboration is building an EL HPGXe TPC capable of holding 100 kg (NEXT-100) of xenon isotopically enriched in ${{}^{136}\rm Xe}$. The installation and commissioning of the NEXT-100 detector at the Laboratorio Subterr\'aneo de Canfranc (LSC) is planned for 2018. The current estimated background level for the NEXT-100 detector is of $4\times10^{-4}$ counts/keV-kg-yr or less in the energy region of interest. Assuming an energy resolution of 0.75$\%$ FWHM at the $Q_{\beta\beta}$ and a $0\nu\beta\beta$ signal efficiency of about 28$\%$, this gives an expected sensitivity (at 90$\%$ CL) to the $0\nu\beta\beta$ decay half life of $T^{0\nu}_{1/2}>6.0\times10^{25}$ yr for an exposure of 275 kg yr. A first phase of the NEXT experiment, called NEW, is currently being commissioned at the LSC. The NEW detector is a scale 1:2 in size (1:10 in mass) of the NEXT-100 detector using the same materials and photosensors and will be used to perform a characterization of the $0\nu\beta\beta$ backgrounds and a measurement of the standard double beta decay with neutrinos (${2\nu\beta\beta}$). An 8 sigma significance for the ${2\nu\beta\beta}$ signal in the NEW detector has been estimated for a 100-day run.


title: Evidences against cuspy dark matter halos in large galaxies

abstract: We develop and apply new techniques in order to uncover galaxy rotation curves (RC) systematics. Considering that an ideal dark matter (DM) profile should yield RCs that have no bias towards any particular radius, we find that the Burkert DM profile satisfies the test, while the Navarro-Frenk-While (NFW) profile has a tendency of better fitting the region between one and two disc scale lengths than the inner disc scale length region. Our sample indicates that this behaviour happens to more than 75% of the galaxies fitted with an NFW halo. Also, this tendency does not weaken by considering "large" galaxies, for instance those with $M_*\gtrsim 10^{10} M_\odot$. Besides the tests on the homogeneity of the fits, we also use a sample of 62 galaxies of diverse types to perform tests on the quality of the overall fit of each galaxy, and to search for correlations with stellar mass, gas mass and the disc scale length. In particular, we find that only 13 galaxies are better fitted by the NFW halo; and that even for the galaxies with $M_* \gtrsim 10^{10} M_\odot$ the Burkert profile either fits as good as, or better than, the NFW profile. This result is relevant since different baryonic effects important for the smaller galaxies, like supernova feedback and dynamical friction from baryonic clumps, indicate that at such large stellar masses the NFW profile should be preferred over the Burkert profile. Hence, our results either suggest a new baryonic effect or a change of the dark matter physics.


title: Optomechanically-induced chiral transport of phonons in one dimension

abstract: Non-reciprocal devices, with one-way transport properties, form a key component for isolating and controlling light in photonic systems. Optomechanical systems have emerged as a potential platform for optical non-reciprocity, due to ability of a pump laser to break time and parity symmetry in the system. Here we consider how the non-reciprocal behavior of light can also impact the transport of sound in optomechanical devices. We focus on the case of a quasi one dimensional optical ring resonator with many mechanical modes coupled to light via the acousto-optic effect. The addition of disorder leads to finite diffusion for phonon transport in the material, largely due to elastic backscattering between clockwise and counter-clockwise phonons. We show that a laser pump field, along with the assumption of high quality-factor, sideband-resolved optical resonances, suppresses the effects of disorder and leads to the emergence of chiral diffusion, with direction-dependent diffusion emerging in a bandwidth similar to the phase-matching bandwidth for Brillouin scattering. A simple diagrammatic theory connects the observation of reduced mechanical linewidths directly to the associated phonon diffusion properties, and helps explain recent experimental results.


title: Superconductivity and Dirac Fermions in 112-phase Pnictides

abstract: This article reviews the status of current research on the 112-phase of pnictides. The 112-phase has gained augmented attention due to the recent discovery of high-temperature superconductivity in $\cl$ with a maximum critical temperature $\tc\sim$ 47\,K upon Sb substitution. The structural, magnetic, and electronic properties of $\cl$ bear some similarities with other superconducting pnictide phases, however, the different valence states of the pnictogen and the presence of a metallic spacer layer are unique features of the 112-system. Low-temperature superconductivity which coexists with antiferromagnetic order was observed in transition metal (Ni, Pd) deficient 112-compounds like $\cn$, $\lpb$, $\lps$, $\lns$. Besides superconductivity, the presence of naturally occurring anisotropic Dirac Fermionic states were observed in the layered 112-compounds $\smb$, $\cmb$, $\lab$ which are of significant interest for future nanoelectronics as an alternative to graphene. In these compounds, the linear energy dispersion resulted in a high magnetoresistance that stayed unsaturated even at the highest applied magnetic fields. Here, we describe various 112-type materials systems combining experimental results and theoretical predictions to stimulate further research on this less well-known member of the pnictide family.


title: Carrier mobility and scattering lifetime in electric double-layer gated few-layer graphene

abstract: We fabricate electric double-layer field-effect transistor (EDL-FET) devices on mechanically exfoliated few-layer graphene. We exploit the large capacitance of a polymeric electrolyte to study the transport properties of three, four and five-layer samples under a large induced surface charge density both above and below the glass transition temperature of the polymer. We find that the carrier mobility shows a strong asymmetry between the hole and electron doping regime. We then employ ab-initio density functional theory (DFT) calculations to determine the average scattering lifetime from the experimental data. We explain its peculiar dependence on the carrier density in terms of the specific properties of the electrolyte we used in our experiments.


title: Detailed analysis of Rouse mode and dynamic scattering function of highly entangled polymer melts in equilibrium

abstract: We present large-scale molecular dynamics simulations for a coarse-grained model of polymer melts in equilibrium. From detailed Rouse mode analysis we show that the time-dependent relaxation of the autocorrelation function (ACF) of modes $p$ can be well described by the effective stretched exponential function due to the crossover from Rouse to reptation regime. The ACF is independent of chain sizes $N$ for $N/p<N_e$ ($N_e$ is the entanglement length), and there exists a minimum of the stretching exponent as $N/p \rightarrow N_e$. As $N/p$ increases, we verify the crossover scaling behavior of the effective relaxation time $\tau_{{\rm eff},p}$ from the Rouse regime to the reptation regime. We have also provided evidence that the incoherent dynamic scattering function follows the same crossover scaling behavior of the mean square displacement of monomers at the corresponding characteristic time scales. The decay of the coherent dynamic scattering function is slowed down and a plateau develops as chain sizes increase at the intermediate time and wave length scales. The tube diameter extracted from the coherent dynamic scattering function is equivalent to the previous estimate from the mean square displacement of monomers.


title: Detailed investigation of the duration of inflation in loop quantum cosmology for a Bianchi-I universe with different inflaton potentials and initial conditions

abstract: There is a wide consensus on the correct dynamics of the background in loop quantum cosmology. In this article we make a systematic investigation of the duration of inflation by varying what we think to be the most important "unknowns" of the model: the way to set initial conditions, the amount of shear at the bounce and the shape of the inflaton potential.


title: What are the visual features underlying human versus machine vision?

abstract: Although Deep Convolutional Networks (DCNs) are approaching the accuracy of human observers at object recognition, it is unknown whether they leverage similar visual representations to achieve this performance. To address this, we introduce Clicktionary, a web-based game for identifying visual features used by human observers during object recognition. Importance maps derived from the game are consistent across participants and uncorrelated with image saliency measures. These results suggest that Clicktionary identifies image regions that are meaningful and diagnostic for object recognition but different than those driving eye movements. Surprisingly, Clicktionary importance maps are only weakly correlated with relevance maps derived from DCNs trained for object recognition. Our study demonstrates that the narrowing gap between the object recognition accuracy of human observers and DCNs obscures distinct visual strategies used by each to achieve this performance.


title: Formation of Massive Rocky Exomoons by Giant Impact

abstract: The formation of satellites is thought to be a natural by-product of planet formation in our Solar System, and thus, moons of extrasolar planets (exomoons) may be abundant in extrasolar planetary systems, as well. Exomoons have yet to be discovered. However, moons larger than 0.1 Earth masses can be detected and characterized using current transit techniques. Here, we show that collisions between rocky planets with masses between a quarter to ten Earth masses can create impact-generated debris disks that could accrete into moons. Collisions between like-sized objects, at oblique impact angles, and velocities near escape speed create disks massive enough to form satellites that are dynamically stable against planetary tides. Impacts of this type onto a superearth between 2 to 7 Earth masses can launch into orbit enough mass to create a satellite large enough to be detected in Kepler transit data. Impact velocity is a crucial controlling factor on disk mass, which has been overlooked in all prior studies of moon formation via planetary collision.


title: Induced interactions in the BCS-BEC crossover of two-dimensional Fermi gases with Rashba spin-orbit coupling

abstract: We investigate the Gorkov--Melik-Barkhudarov (GM) correction to superfluid transition temperature in two-dimensional Fermi gases with Rashba spin-orbit coupling (SOC) across the SOC-driven BCS-BEC crossover. In the calculation of the induced interaction, we find that the spin-component mixing due to SOC can induce both of the conventional screening and additional antiscreening contributions that interplay significantly in the strong SOC regime. While the GM correction generally lowers the estimate of transition temperature, it turns out that at a fixed weak interaction, the correction effect exhibits a crossover behavior where the ratio between the estimates without and with the correction first decreases with SOC and then becomes insensitive to SOC when it goes into the strong SOC regime. We demonstrate the applicability of the GM correction by comparing the zero-temperature condensate fraction with the recent quantum Monte Carlo results.


title: Iteration of Quadratic Polynomials Over Finite Fields

abstract: For a finite field of odd cardinality $q$, we show that the sequence of iterates of $aX^2+c$, starting at $0$, always recurs after $O(q/\log\log q)$ steps. For $X^2+1$ the same is true for any starting value. We suggest that the traditional "Birthday Paradox" model is inappropriate for iterates of $X^3+c$, when $q$ is 2 mod 3.


title: Multiset Combinatorial Batch Codes

abstract: Batch codes, first introduced by Ishai, Kushilevitz, Ostrovsky, and Sahai, mimic a distributed storage of a set of $n$ data items on $m$ servers, in such a way that any batch of $k$ data items can be retrieved by reading at most some $t$ symbols from each server. Combinatorial batch codes, are replication-based batch codes in which each server stores a subset of the data items. In this paper, we propose a generalization of combinatorial batch codes, called multiset combinatorial batch codes (MCBC), in which $n$ data items are stored in $m$ servers, such that any multiset request of $k$ items, where any item is requested at most $r$ times, can be retrieved by reading at most $t$ items from each server. The setup of this new family of codes is motivated by recent work on codes which enable high availability and parallel reads in distributed storage systems. The main problem under this paradigm is to minimize the number of items stored in the servers, given the values of $n,m,k,r,t$, which is denoted by $N(n,k,m,t;r)$. We first give a necessary and sufficient condition for the existence of MCBCs. Then, we present several bounds on $N(n,k,m,t;r)$ and constructions of MCBCs. In particular, we determine the value of $N(n,k,m,1;r)$ for any $n\geq \left\lfloor\frac{k-1}{r}\right\rfloor{m\choose k-1}-(m-k+1)A(m,4,k-2)$, where $A(m,4,k-2)$ is the maximum size of a binary constant weight code of length $m$, distance four and weight $k-2$. We also determine the exact value of $N(n,k,m,1;r)$ when $r\in\{k,k-1\}$ or $k=m$.


title: Orthogonal polynomials of several variables

abstract: Preliminary version of Chapter 2 in the book "Encyclopedia of Special functions: The Askey-Bateman Project, Vol. 2: Multivariate special functions", T. H. Koornwinder and J. V. Stokman (eds.), Cambridge University Press, 2021.


title: Distributed Estimation of Dynamic Fields over Multi-agent Networks

abstract: This work presents distributed algorithms for estimation of time-varying random fields over multi-agent/sensor networks. A network of sensors makes sparse and noisy local measurements of the dynamic field. Each sensor aims to obtain unbiased distributed estimates of the entire field with bounded mean-squared error (MSE) based on its own local observations and its neighbors' estimates. This work develops three novel distributed estimators: Pseudo-Innovations Kalman Filter (PIKF), Distributed Information Kalman Filter (DIKF) and Consensus+Innovations Kalman Filter (CIKF). We design the gain matrices such that the estimators achieve unbiased estimates with bounded MSE under minimal assumptions on the local observation and network communication models. This work establishes trade-offs between these three distributed estimators and demonstrates how they outperform existing solutions. We validate our results through extensive numerical evaluations.


title: On the Feasibility of Malware Authorship Attribution

abstract: There are many occasions in which the security community is interested to discover the authorship of malware binaries, either for digital forensics analysis of malware corpora or for thwarting live threats of malware invasion. Such a discovery of authorship might be possible due to stylistic features inherent to software codes written by human programmers. Existing studies of authorship attribution of general purpose software mainly focus on source code, which is typically based on the style of programs and environment. However, those features critically depend on the availability of the program source code, which is usually not the case when dealing with malware binaries. Such program binaries often do not retain many semantic or stylistic features due to the compilation process. Therefore, authorship attribution in the domain of malware binaries based on features and styles that will survive the compilation process is challenging. This paper provides the state of the art in this literature. Further, we analyze the features involved in those techniques. By using a case study, we identify features that can survive the compilation process. Finally, we analyze existing works on binary authorship attribution and study their applicability to real malware binaries.


title: A proof of Kurdyka's conjecture on arc-analytic functions

abstract: We prove a conjecture of Kurdyka stating that every arc-symmetric semialgebraic set is precisely the zero locus of an arc-analytic semialgebraic function. This implies, in particular, that arc-symmetric semialgebraic sets are in one-to-one correspondence with radical ideals of the ring of arc-analytic semialgebraic functions.


title: Evidence of a substellar companion around a very young T Tauri star

abstract: We present results from a Near Infrared multi-epoch spectroscopic campaign to detect a young low-mass companion to a T Tauri star. AS 205A is a late-type dwarf (~K5) of ~1 Msun that belongs to a triple system. Independent photometric surveys discovered that AS 205A has two distinct periods (P1=6.78 and P2=24.78 days) detected in the light curve that persist over several years. Period P1 seems to be linked to the axial-rotation of the star and is caused by the presence of cool surface spots. Period P2 is correlated with the modulation in AS 205A brightness (V) and red color (V-R), consistent with a gravitating object within the accretion disk. We here derive precise Near Infrared radial velocities to investigate the origin of period P2 which is predicted to correspond to a cool source in a Keplerian orbit with a semi-major axis of ~0.17 AU positioned close to the inner disk radius of 0.14 AU. The radial velocity variations of AS 205A were found to have a period of P ~ 24.84 days and a semi-amplitude of 1.529 kms-1. This result closely resembles the P2 period in past photometric observations (P ~ 24.78 days). The analysis of the cross-correlation function bisector has shown no correlation with the radial velocity modulations, strongly suggesting that the period is not controlled by stellar rotation. Additional activity indicators should however be explored in future surveys. Taking this into account we found that the presence of a substellar companion is the explanation that best fits the results. We derived an orbital solution for AS 205A and found evidence of a m2sini 19.25 MJup object in an orbit with moderate eccentricity of e~0.34. If confirmed with future observations, preferably using a multiwavelength survey approach, this companion could provide interesting constraints on brown dwarf and planetary formation models.


title: H-infinity Filtering for Cloud-Aided Semi-active Suspension with Delayed Information

abstract: This chapter presents an H-infinity filtering framework for cloud-aided semiactive suspension system with time-varying delays. In this system, road profile information is downloaded from a cloud database to facilitate onboard estimation of suspension states. Time-varying data transmission delays are considered and assumed to be bounded. A quarter-car linear suspension model is used and an H-infinity filter is designed with both onboard sensor measurements and delayed road profile information from the cloud. The filter design procedure is designed based on linear matrix inequalities (LMIs). Numerical simulation results are reported that illustrates the fusion of cloud-based and on-board information that can be achieved in Vehicleto- Cloud-to-Vehicle (V2C2V) implementation.


title: Lattice symmetries, spectral topology and opto-electronic properties of graphene-like materials

abstract: The topology of the band structure, which is determined by the lattice symmetries, has a strong influence on the transport properties. Here we consider an anisotropic honeycomb lattice and study the effect of a continuously deformed band structure on the optical conductivity and on diffusion due to quantum fluctuations. In contrast to the behavior at an isotropic node we find super- and subdiffusion for the anisotropic node. The spectral saddle points create van Hove singularities in the optical conductivity, which could be used to characterize the spectral properties experimentally.


title: On the minimum value of sum-Balaban index

abstract: We consider extremal values of sum-Balaban index among graphs on $n$ vertices. We determine that the upper bound for the minimum value of the sum-Balaban index is at most $4.47934$ when $n$ goes to infinity. For small values of $n$ we determine the extremal graphs and we observe that they are similar to dumbbell graphs, in most cases having one extra edge added to the corresponding extreme for the usual Balaban index. We show that in the class of balanced dumbbell graphs, those with clique sizes $\sqrt[4]{\sqrt 2\log\big(1+\sqrt 2\big)}\sqrt n+o(\sqrt n)$ have asymptotically the smallest value of sum-Balaban index. We pose several conjectures and problems regarding this topic.


title: Searching for compact radio sources associated to UCHII regions

abstract: Ultra-Compact (UC)HII regions represent a very early stage of massive star formation whose structure and evolution are not yet fully understood. Interferometric observations in recent years show that some UCHII regions have associated compact sources of uncertain nature. Based on this, we carried out VLA 1.3 cm observations in the A configuration of selected UCHII regions in order to report additional cases of compact sources embedded in UCHII regions. From the observations, we find 13 compact sources associated to 9 UCHII regions. Although we cannot establish an unambiguous nature for the newly detected sources, we assess some of their observational properties. According to the results, we can distinguish between two types of compact sources. One type corresponds to sources that probably are deeply embedded in the dense ionized gas of the UCHII region. These sources are being photo-evaporated by the exciting star of the region and will last for 10$^4-10^5$ yr. They may play a crucial role in the evolution of the UCHII region as the photo-evaporated material could replenish the expanding plasma and might provide a solution to the so-called lifetime problem for these regions. The second type of compact sources is not associated with the densest ionized gas of the region. A few of these sources appear resolved and may be photo-evaporating objects such as those of the first type but with significantly lower mass depletion rates. The rest of sources of this second type appear unresolved and their properties are varied. We speculate on the similarity between the sources of the second type and those of the Orion population of radio sources.


title: See the Glass Half Full: Reasoning about Liquid Containers, their Volume and Content

abstract: Humans have rich understanding of liquid containers and their contents; for example, we can effortlessly pour water from a pitcher to a cup. Doing so requires estimating the volume of the cup, approximating the amount of water in the pitcher, and predicting the behavior of water when we tilt the pitcher. Very little attention in computer vision has been made to liquids and their containers. In this paper, we study liquid containers and their contents, and propose methods to estimate the volume of containers, approximate the amount of liquid in them, and perform comparative volume estimations all from a single RGB image. Furthermore, we show the results of the proposed model for predicting the behavior of liquids inside containers when one tilts the containers. We also introduce a new dataset of Containers Of liQuid contEnt (COQE) that contains more than 5,000 images of 10,000 liquid containers in context labelled with volume, amount of content, bounding box annotation, and corresponding similar 3D CAD models.


title: CCD time-series photometry of variable stars in globular clusters and the metallicity dependence of the horizontal branch luminosity

abstract: We describe and summarize the findings from our CCD time-series photometry of globular clusters (GCs) program and the use of difference image analysis (DIA) in the extraction of very precise light curves even in the crowded central regions down to $V\sim 19$ mag. We have discovered approximately 250 variable stars in a sample of 23 selected GCs and have, therefore, updated the census of variables in each system and in some cases we have actually completed it down to a certain magnitude. The absolute magnitude and [Fe/H] for each individual RR Lyrae is obtained via the Fourier decomposition of the light curve. An average of these parameters leads to the distance and metallicity of the host GCs. We have also calibrated the P-L relation for SX Phe stars which enables an independent calculation of the cluster distance. We present the mean [Fe/H], $M_V$ and distance for a group of selected GCs based exclusively on the RR Lyrae light curve Fourier decomposition technique and set on a rather unprecedented homogeneous scale. We also discuss the luminosity dependence of the horizontal branch (HB) via the $M_V$-[Fe/H] relation. We find that this relation should be considered separately for the RRab and RRc stars.


title: Towards End-to-End Speech Recognition with Deep Convolutional Neural Networks

abstract: Convolutional Neural Networks (CNNs) are effective models for reducing spectral variations and modeling spectral correlations in acoustic features for automatic speech recognition (ASR). Hybrid speech recognition systems incorporating CNNs with Hidden Markov Models/Gaussian Mixture Models (HMMs/GMMs) have achieved the state-of-the-art in various benchmarks. Meanwhile, Connectionist Temporal Classification (CTC) with Recurrent Neural Networks (RNNs), which is proposed for labeling unsegmented sequences, makes it feasible to train an end-to-end speech recognition system instead of hybrid settings. However, RNNs are computationally expensive and sometimes difficult to train. In this paper, inspired by the advantages of both CNNs and the CTC approach, we propose an end-to-end speech framework for sequence labeling, by combining hierarchical CNNs with CTC directly without recurrent connections. By evaluating the approach on the TIMIT phoneme recognition task, we show that the proposed model is not only computationally efficient, but also competitive with the existing baseline systems. Moreover, we argue that CNNs have the capability to model temporal correlations with appropriate context information.


title: One-dimensional von K\'arm\'an models for elastic ribbons

abstract: By means of a variational approach we rigorously deduce three one-dimensional models for elastic ribbons from the theory of von K\'arm\'an plates, passing to the limit as the width of the plate goes to zero. The one-dimensional model found starting from the "linearized" von K\'arm\'an energy corresponds to that of a linearly elastic beam that can twist but can deform in just one plane; while the model found from the von K\'arm\'an energy is a non-linear model that comprises stretching, bendings, and twisting. The "constrained" von K\'arm\'an energy, instead, leads to a new Sadowsky type of model.


title: Semianalytical calculation of the zonal-flow oscillation frequency in stellarators

abstract: Due to their capability to reduce turbulent transport in magnetized plasmas, understanding the dynamics of zonal flows is an important problem in the fusion programme. Since the pioneering work by Rosenbluth and Hinton in axisymmetric tokamaks, it is known that studying the linear and collisionless relaxation of zonal flow perturbations gives valuable information and physical insight. Recently, the problem has been investigated in stellarators and it has been found that in these devices the relaxation process exhibits a characteristic feature: a damped oscillation. The frequency of this oscillation might be a relevant parameter in the regulation of turbulent transport, and therefore its efficient and accurate calculation is important. Although an analytical expression can be derived for the frequency, its numerical evaluation is not simple and has not been exploited systematically so far. Here, a numerical method for its evaluation is considered, and the results are compared with those obtained by calculating the frequency from gyrokinetic simulations. This "semianalytical" approach for the determination of the zonal-flow frequency reveals accurate and faster than the one based on gyrokinetic simulations.


title: Measurements of Solar Differential Rotation and Meridional Circulation from Tracking of Photospheric Magnetic Features

abstract: Long-lived rotational and meridional flows are important ingredients of the solar cycle. Magnetic field images have typically been used to measure these flows on the solar surface by cross-correlating thin longitudinal strips or square patches across sufficiently long time gaps. Here, I use one month of SDO/HMI line-of-sight magnetic field observations, combined with the SWAMIS magnetic feature tracking algorithm to measure the motion of individual features in these magnetograms. By controlling for perturbations due to short-lived flows and due to false motions from feature interactions, I effectively isolate the long-lived flows traced by the magnetic features. This allows me to produce high-fidelity differential rotation measurements with well-characterized variances and covariances of the fit parameters.I find a sidereal rotational profile of $(14.296\pm0.006)+(-1.847\pm0.056)\sin^{2}b+(-2.615\pm0.093)\sin^{4}b$, with units of $\textrm{ deg d}^{-1}$, and a large covariance $\sigma_{BC}^{2}=-4.87\times10^{-3}(\textrm{ deg d}^{-1})^{2}$. I also produce medium-fidelity measurements of the much weaker meridional flow that is broadly consistent with previous results. This measurement shows a peak flow of $16.7\pm0.6\text{ m s}^{-1}$ at latitude $b=45^\circ$ but is insufficiently characterized at higher latitudes to ascertain whether the chosen functional form $2\cos b\sin b$ is appropriate. This work shows that measuring the motions of individual features in photospheric magnetograms can produce high precision results in relatively short time spans, and suggests that high resolution non-longitudinally averaged photospheric velocity residual measurements could be produced to compare with coronal results, and to provide other diagnostics of the solar dynamo.


title: The binegativity of two qubits

abstract: We introduce an entanglement-related quantity that we call the binegativity. Based on numerical evidence, we conjecture that the binegativity is an entanglement measure for two-qubit states. The binegativity is compared to the concurrence and negativity, the only known entanglement measures that can be computed analytically for all two-qubit states. The ordering of entangled two-qubit mixed states induced by the binegativity is distinct from the ordering determined by previously known measures, yielding possible new insights into the structure of entangled two-qubit mixed states.


title: Dynamic fracture of a discrete media under moving load

abstract: Most of the research concerting crack propagation in discrete media is concerned with specific types of external loading: displacements on the boundaries, or constant energy fluxes or feeding waves originating from infinity. In this paper the action of a moving load is analysed on the simplest lattice model: a thin strip, where the fault propagating in its middle portion as the result of the moving force acting on the destroyed part of the structure. We study both analytically and numerically how the load amplitude and its velocity influence the possible solution, and specifically the way the fracture process reaches its steady-state regime. We present the relation between the possible steady-state crack speed and the loading parameters, as well as the energy release rate. In particular, we show that there exists a class of loading regime corresponding to each point on the energy-speed diagram (and thus determine the same limiting steady-state regime). The phenomenon of the forbidden regimes is discussed in detail, from both the points of view of force and energy. For a sufficiently anisotropic structure, we find a stable steady-state propagation corresponding to the slow crack. Numerical simulations reveal various ways by which the process approaches - or fails to approach - the steady-state regime. The results extend our understanding of fracture processes in discrete structures, and reveal some new questions that should be addressed.


title: Low Mass X-Ray Binaries in the Inner Galaxy: Implications for Millisecond Pulsars and the GeV Excess

abstract: If millisecond pulsars (MSPs) are responsible for the excess gamma-ray emission observed from the region surrounding the Galactic Center, the same region should also contain a large population of low-mass X-ray binaries (LMXBs). In this study, we compile and utilize a sizable catalog of LMXBs observed in the the Milky Way's globular cluster system and in the Inner Galaxy, as well as the gamma-ray emission observed from globular clusters, to estimate the flux of gamma rays predicted from MSPs in the Inner Galaxy. From this comparison, we conclude that only up to $\sim$4-23% of the observed gamma-ray excess is likely to originate from MSPs. This result is consistent with, and more robust than, previous estimates which utilized smaller samples of both globular clusters and LMXBs. If MSPs had been responsible for the entirety of the observed excess, INTEGRAL should have detected $\sim$$10^3$ LMXBs from within a $10^{\circ}$ radius around the Galactic Center, whereas only 42 LMXBs (and 46 additional LMXB candidates) have been observed.

categories: astro-ph.HE astro-ph.CO astro-ph.GA hep-ph


title: On the accurate computation of the true contact-area in mechanical contact of random rough surfaces

abstract: We introduce a corrective function to compensate errors in contact area computations coming from mesh discretization. The correction is based on geometrical arguments and requires only one additional quantity to be computed: the length of contact/non-contact interfaces. The new technique enables us to evaluate accurately the true contact area using a coarse mesh for which the shortest wavelength in the surface spectrum reaches the grid size. The validity of the approach is demonstrated for surfaces with different fractal dimensions and different spectral content using a properly designed mesh convergence test. In addition, we use a topology preserving smoothing technique to adjust the morphology of contact clusters obtained with a coarse grid.


title: Suzaku spectroscopy of the neutron star transient 4U 1608-52 during its outburst decay

abstract: We test the proposed 3-component spectral model for neutron star low mass X-ray binaries using broad-band X-ray data. We have analysed 4 X-ray spectra (0.8-30 keV) obtained with Suzaku during the 2010 outburst of 4U 1608-52, which have allowed us to perform a comprehensive spectral study covering all the classical spectral states. We use a thermally Comptonized continuum component to account for the hard emission, as well as two thermal components to constrain the accretion disc and neutron star surface contributions. We find that the proposed combination of multicolor disc, single-temperature black body and Comptonization components successfully reproduces the data from soft to hard states. In the soft state, our study supports the neutron star surface (or boundary layer) as the dominant source for the Comptonization seed photons yielding the observed weak hard emission, while in the hard state both solutions, either the disc or the neutron star surface, are equally favoured. The obtained spectral parameters as well as the spectral/timing correlations are comparable to those observed in accreting black holes, which support the idea that black hole and neutron star low mass X-ray binaries undergo a similar state evolution during their accretion episodes.


title: Digital Advertising Traffic Operation: Flow Management Analysis

abstract: In a Web Advertising Traffic Operation the Trafficking Routing Problem (TRP) consists in scheduling the management of Web Advertising (Adv) campaign between Trafficking campaigns in the most efficient way to oversee and manage relationship with partners and internal teams, managing expectations through integration and post-launch in order to ensure success for every stakeholders involved. For our own interest we did that independent research projects also through specific innovative tasks validate towards average working time declared on "specification required" by the main worldwide industry leading Advertising Agency. We present a Mixed Integer Linear Programming (MILP) formulation for end-to-end management of campaign workflow along a predetermined path and generalize it to include alternative path to oversee and manage detail-oriented relationship with partners and internal teams to achieve the goals above mentioned. To meet clients' KPIs, we consider an objective function that includes the punctuality indicators (the average waiting time and completion times) but also the main punctuality indicators (the average delay and the on time performance). Then we investigate their analytical relationships in the advertising domain through experiments based on real data from a Traffic Office. We show that the classic punctuality indicators are in contradiction with the task of reducing waiting times. We propose new indicators used for a synthesize analysis on projects or process changes for the wider team that are more sustainable, but also more relevant for stakeholders. We also show that the flow of a campaign (adv-ways) is the main bottleneck of a Traffic Office and that alternate paths cannot improve the performance indicators.


title: Constraints on Bygone Nucleosynthesis of Accreting Neutron Stars

abstract: Nuclear burning near the surface of an accreting neutron star produces ashes that, when compressed deeper by further accretion, alter the star's thermal and compositional structure. Bygone nucleosynthesis can be constrained by the impact of compressed ashes on the thermal relaxation of quiescent neutron star transients. In particular, Urca cooling nuclei pairs in nuclear burning ashes, which cool the neutron star crust via neutrino emission from electron-capture/beta-decay cycles, provide signatures of prior nuclear burning over the ~century timescales it takes to accrete to the electron-capture depth of the strongest cooling pairs. Using crust cooling models of the accreting neutron star transient MAXI J0556-332, we show that this source likely lacked Type I X-ray bursts and superbursts >120 years ago. Reduced nuclear physics uncertainties in rp-process reaction rates and electron-capture ft-values for low-lying transitions will improve nucleosynthesis constraints using this technique.


title: The 'twin paradox' in relativistic rigid motion

abstract: Relativistic rigid motion suggests a new version for the so-called `twin paradox', comparing the ages of two astronauts on a very long spaceship. Although there is always an instantaneous inertial frame in which the whole spaceship, being rigid, is simultaneously at rest, the twins' ages, measured as the proper-times along their individual world lines, are different when they are located at remote parts of the spaceship. The age, or proper-time, difference depends on the distance at rest between the astronauts and the rapidity difference between start to end. The relation of the age difference with the relative Doppler shift of light signals transmitted between the astronauts, and implications for the possibility to assign common age (proper-time) to complex, spatially extended, relativistic systems, are also discussed. The condition for simultaneous arrival of light signals emitted simultaneously from the opposite ends of a rigidly accelerating spaceship is resolved.


title: The graphs of join-semilattices and the shape of congruence lattices of particle lattices

abstract: We attach to each $\langle 0, \vee \rangle$-semilattice a graph $\boldsymbol{G}_{\boldsymbol{S}}$ whose vertices are join-irreducible elements of $\boldsymbol{S}$ and whose edges correspond to the reflexive dependency relation. We study properties of the graph $\boldsymbol{G}_{\boldsymbol{S}}$ both when $\boldsymbol{S}$ is a join-semilattice and when it is a lattice. We call a $\langle 0, \vee \rangle$-semilattice $\boldsymbol{S}$ particle provided that the set of its join-irreducible elements join-generates $\boldsymbol{S}$ and it satisfies DCC. We prove that the congruence lattice of a particle lattice is anti-isomorphic to the lattice of hereditary subsets of the corresponding graph that are closed in a certain zero-dimensional topology. Thus we extend the result known for principally chain finite lattices.


title: Charge Coupled Device detectors with high quantum efficiency at UV wavelengths

abstract: We report on multilayer high efficiency antireflection coating (ARC) design and development for use at UV wavelengths on CCDs and other Si-based detectors. We have previously demonstrated a set of single-layer coatings, which achieve >50% quantum efficiency (QE) in four bands from 130 to 300 nm. We now present multilayer coating designs that significantly outperform our previous work between 195 and 215 nm. Using up to 11 layers, we present several model designs to reach QE above 80%. We also demonstrate the successful performance of 5 and 11 layer ARCs on silicon and fused silica substrates. Finally, we present a five-layer coat- ing deposited onto a thinned, delta-doped CCD and demonstrate external QE greater than 60% between 202 and 208 nm, with a peak of 67.6% at 206 nm.


title: The COSMOS2015 galaxy stellar mass function: 13 billion years of stellar mass assembly in 10 snapshots

abstract: We measure the stellar mass function (SMF) of galaxies in the COSMOS field up to $z\sim6$. We select them in the near-IR bands of the COSMOS2015 catalogue, which includes ultra-deep photometry from UltraVISTA-DR2, SPLASH, and Subaru/Hyper-SuprimeCam. At $z>2.5$ we use new precise photometric redshifts with error $\sigma_z=0.03(1+z)$ and an outlier fraction of $12\%$, estimated by means of the unique spectroscopic sample of COSMOS. The increased exposure time in the DR2, along with our panchromatic detection strategy, allow us to improve the stellar mass completeness at high $z$ with respect to previous UltraVISTA catalogues. We also identify passive galaxies through a robust colour-colour selection, extending their SMF estimate up to $z=4$. Our work provides a comprehensive view of galaxy stellar mass assembly between $z=0.1$ and 6, for the first time using consistent estimates across the entire redshift range. We fit these measurements with a Schechter function, correcting for Eddington bias. We compare the SMF fit with the halo mass function predicted from $\Lambda$CDM simulations. We find that at $z>3$ both functions decline with a similar slope in the high-mass end. This feature could be explained assuming that the mechanisms that quench star formation in massive haloes become less effective at high redshift; however further work needs to be done to confirm this scenario. Concerning the SMF low-mass end, it shows a progressive steepening as moving towards higher redshifts, with $\alpha$ decreasing from $-1.47_{-0.02}^{+0.02}$ at $z\simeq0.1$ to $-2.11_{-0.13}^{+0.30}$ at $z\simeq5$. This slope depends on the characterisation of the observational uncertainties, which is crucial to properly remove the Eddington bias. We show that there is currently no consensus on the method to quantify such errors: different error models result in different best-fit Schechter parameters. [Abridged]


title: Cosmology and time dependent parameters induced by misaligned light scalar

abstract: We consider a scenario where time dependence on physical parameters is introduced by the misalignment of an ultra-light scalar field. The initial VEV of such field at the early time remains a constant until Hubble becomes comparable to its mass. Interesting cosmological consequences are considered. Light sterile neutrinos hinted by terrestrial neutrino experiments are studied as a benchmark model. We show the BBN constraints can be easily avoided in this scenario, even if reheating temperature is high. The scalar can be naturally light in spite of its couplings to other fields. Parameters of sterile neutrino may remain changing with time nowadays. This can further relax the tension from the recent IceCube constraints.


title: Colour gradients in cluster ellipticals at z~1.4: the hidden content of the galaxy central regions

abstract: We present F775W-F850LP (rest-frame UV-U) and F850LP-F160W (rest-frame U-R) colour gradients for a sample of 17 elliptical galaxies morphologically selected in the cluster XMMU J2235.3-2557 at $z$=1.39. We detected significant negative (redder inwards) U-R colour gradients in $\sim$70 per cent of the galaxies and flat gradients for the remaining ones. On the other hand, the UV-U gradients are significant positive (bluer inwards) for $\sim$80 per cent of the galaxies and flat for the remaining ones. Using stellar population synthesis models, we found that the behaviour of the two colour gradients cannot be simultaneously explained by a radial variation of age, metallicity and/or dust. The observed U-R gradients are consistent with a metallicity gradient (mean value $\nabla_{Z} =-0.4$) in agreement with the one observed in the local elliptical galaxies. The positive UV-U gradients cannot be explained with age or metallicity variations and imply an excess of UV emission towards the galaxies' central regions.This excess calls into question mechanisms able to efficiently produce UV emission. The data require either steady weak star formation (< 1 M$_\odot$ yr$^{-1}$) or an He-rich population in the cores of these galaxies in order to simultaneously reproduce both the colour gradients. On the contrary, the presence of a QSO cannot account for the observed UV excess on its own. We discuss these hypotheses on the basis of current observations and available models.


title: Constraints on Supersymmetric Dark Matter for Heavy Scalar Superpartners

abstract: We study the constraints on neutralino dark matter in minimal low energy supersymmetry models and the case of heavy lepton and quark scalar superpartners. For values of the Higgsino and gaugino mass parameters of the order of the weak scale, direct detection experiments are already putting strong bounds on models in which the dominant interactions between the dark matter candidates and nuclei are governed by Higgs boson exchange processes, particularly for positive values of the Higgsino mass parameter mu. For negative values of mu, there can be destructive interference between the amplitudes associated with the exchange of the standard CP-even Higgs boson and the exchange of the non-standard one. This leads to specific regions of parameter space which are consistent with the current experimental constraints and a thermal origin of the observed relic density. In this article we study the current experimental constraints on these scenarios, as well as the future experimental probes, using a combination of direct and indirect dark matter detection and heavy Higgs and electroweakino searches at hadron colliders.


title: Effective Temperatures of Cataclysmic Variable White Dwarfs as a Probe of their Evolution

abstract: We present HST spectroscopy for 45 cataclysmic variables (CVs), observed with HST/COS and HST/STIS. For 36 CVs, the white dwarf is recognisable through its broad Ly$\alpha$ absorption profile and we measure the white dwarf effective temperatures ($T_{\mathrm{eff}}$) by fitting the HST data assuming $\log\,g=8.35$, which corresponds to the average mass for CV white dwarfs ($\simeq\,0.8\,\mathrm{M}_\odot$). Our results nearly double the number of CV white dwarfs with an accurate temperature measurement. We find that CVs above the period gap have, on average, higher temperatures ($\langle T_{\mathrm{eff}} \rangle \simeq 23\,000\,$K) and exhibit much more scatter compared to those below the gap ($\langle T_{\mathrm{eff}} \rangle \simeq 15\,000\,$K). While this behaviour broadly agrees with theoretical predictions, some discrepancies are present: (i) all our new measurements above the gap are characterised by lower temperatures ($T_{\mathrm{eff}} \simeq 16\,000 - 26\,000\,$K) than predicted by the present day CV population models ($T_{\mathrm{eff}} \simeq 38\,000 - 43\,000\,$K); (ii) our results below the gap are not clustered in the predicted narrow track and exhibit in particular a relatively large spread near the period minimum, which may point to some shortcomings in the CV evolutionary models. Finally, in the standard model of CV evolution, reaching the minimum period, CVs are expected to evolve back towards longer periods with mean accretion rates $\dot{M}\lesssim 2 \times 10^{-11}\,\mathrm{M}_\odot\,\mathrm{yr}^{-1}$, corresponding to $T_\mathrm{eff}\lesssim 11\,500\,$K. We do not unambiguously identify any such system in our survey, suggesting that this major component of the predicted CV population still remains elusive to observations.


title: Fast Weak Lensing Simulations with Halo Model

abstract: Full ray-tracing maps of gravitational lensing, constructed from N-Body simulations, represent a fundamental tool to interpret present and future weak lensing data. However the limitation of computational resources and storage capabilities severely restrict the number of realizations that can be performed in order to accurately sample both the cosmic shear models and covariance matrices. In this paper we present a halo model formalism for weak gravitational lensing that alleviates these issues by producing weak-lensing mocks at a reduced computational cost. Our model takes as input the halo population within a desired light-cone and the linear power spectrum of the underlined cosmological model. We examine the contribution given by the presence of substructures within haloes to the cosmic shear power spectrum and quantify it to the percent level. Our method allows us to reconstruct high-resolution convergence maps, for any desired source redshifts, of light-cones that realistically trace the matter density distribution in the universe, account for masked area and sample selections. We compare our analysis on the same large scale structures constructed using ray-tracing techniques and find very good agreements both in the linear and non-linear regimes up to few percent levels. The accuracy and speed of our method demonstrate the potential of our halo model for weak lensing statistics and the possibility to generate a large sample of convergence maps for different cosmological models as needed for the analysis of large galaxy redshift surveys.


title: Uncovering the Spatiotemporal Patterns of Collective Social Activity

abstract: Social media users and microbloggers post about a wide variety of (off-line) collective social activities as they participate in them, ranging from concerts and sporting events to political rallies and civil protests. In this context, people who take part in the same collective social activity often post closely related content from nearby locations at similar times, resulting in distinctive spatiotemporal patterns. Can we automatically detect these patterns and thus provide insights into the associated activities? In this paper, we propose a modeling framework for clustering streaming spatiotemporal data, the Spatial Dirichlet Hawkes Process (SDHP), which allows us to automatically uncover a wide variety of spatiotemporal patterns of collective social activity from geolocated online traces. Moreover, we develop an efficient, online inference algorithm based on Sequential Monte Carlo that scales to millions of geolocated posts. Experiments on synthetic data and real data gathered from Twitter show that our framework can recover a wide variety of meaningful social activity patterns in terms of both content and spatiotemporal dynamics, that it yields interesting insights about these patterns, and that it can be used to estimate the location from where a tweet was posted.


title: The accreted stellar halo as a window on halo assembly in L* galaxies

abstract: Theory and observations agree that the accreted stellar halos (ASHs) of Milky Way-like galaxies display significant scatter. I take advantage of this stochasticity to invert the link between halo assembly history (HAH) and ASH, using mock ASHs corresponding to 750 $\Lambda$CDM HAHs, sharing a final virial mass of $M_{h}(z=0)=10^{12.25}M_\odot$. Hosts with poor/rich ASHs assemble following orthogonal growth-patterns. Hosts with rich ASHs experience accretion events (AEs) with high virial mass ratios (HVMRs, $M_s/M_h\gtrsim 0.1$) at $0.5\lesssim z_{infall}\lesssim1.5$, in a phase of fast growth. This maximizes the accreted stellar mass under the condition these satellites are disrupted by $z=0$. At similar times, hosts with poor ASHs grow slowly through minor mergers, with only very recent HVMR AEs: this results in a globally more abundant satellite population and in distinctive surviving massive satellites (stellar mass $\log M_{s,*}/M_\odot\gtrsim 9$). Several properties of the Milky Way are in agreement with the predictions of this framework for hosts with poor, concentrated ASHs, including: i) the recent infall of Sagittarius and Magellanic Clouds, ii) the likely higher-than-average concentration of its dark halo, iii) the signatures of fast chemical enrichment of a sizable fraction of its halo stellar populations.


title: Detection of Rotational Spectral Variation on the M-type asteroid (16) Psyche

abstract: The asteroid (16) Psyche is of scientific interest because it contains ~ 1% of the total mass of the asteroid belt and is thought to be the remnant metallic core of a protoplanet. Radar observations have indicated the significant presence of metal on the surface with a small percentage of silicates. Prior ground-based observations showed rotational variations in the near-infrared (NIR) spectra and radar albedo of this asteroid. However, no comprehensive study that combines multi-wavelength data has been conducted so far. Here we present rotationally resolved NIR spectra (0.7-2.5 microns) of (16) Psyche obtained with the NASA Infrared Telescope Facility. These data have been combined with shape models of the asteroid for each rotation phase. Spectral band parameters extracted from the NIR spectra show that the pyroxene band center varies from ~ 0.92 to 0.94 microns. Band center values were used to calculate the pyroxene chemistry of the asteroid, whose average value was found to be Fs30En65Wo5. Variations in the band depth were also observed, with values ranging from 1.0 to 1.5%. Using a new laboratory spectral calibration we estimated an average orthopyroxene content of 6+/-1%. The mass-deficit region of Psyche, which exhibits the highest radar albedo, also shows the highest value for spectral slope and the minimum band depth. The spectral characteristics of Psyche suggest that its parent body did not have the typical structure expected for a differentiated body or that the sequence of events that led to its current state was more complex than previously thought.


title: ZOMG II: Does the halo assembly history influence central galaxies and gas accretion?

abstract: The growth-rate and the internal dynamics of galaxy-sized dark-matter haloes depend on their location within the cosmic web. Haloes that sit at the nodes grow in mass till the present time and are dominated by radial orbits. Conversely, haloes embedded in prominent filaments do not change much in size and are dominated by tangential orbits. Using zoom hydrodynamical simulations including star formation and feedback, we study how gas accretes onto these different classes of objects that, for simplicity, we dub 'accreting' and 'stalled' haloes. We find that all haloes get a fresh supply of newly accreted gas in their inner regions, although this slowly decreases with time, in particular for the stalled haloes. The inflow of new gas is always higher than (but comparable with) that of recycled material. Overall, the cold-gas fraction increases (decreases) with time for the accreting (stalled) haloes. In all cases, a stellar disc and a bulge form at the centre of the simulated haloes. The total stellar mass is in excellent agreement with expectations based on the abundance-matching technique. Many properties of the central galaxies do not seem to correlate with the large-scale environment in which the haloes reside. However, there are two notable exceptions that characterise stalled haloes with respect to their accreting counterparts: i) the galaxy disc contains much older stellar populations; ii) its vertical scale-height is larger by a factor of two or more. This thickening is likely due to the heating of the long-lived discs by mergers and close flybys.


title: Non-interacting central site model: localization and logarithmic entanglement growth

abstract: We investigate the stationary and dynamical behavior of an Anderson localized chain coupled to a single central bound state. The coupling to the central site partially dilutes the Anderson localized peak towards the nearly resonant sites. In particular, the number of resonantly coupled sites remains finite in the thermodynamic limit. This is further supported by a multifractal analysis of eigenstates that shows the frozen spectrum of fractal dimension, which is characteristic for localized phases in models with power-law hopping. Although the well-known Fano-resonance problem is seemingly similar to our system, it fails to describe it because of the absence of level repulsion within the energy spectrum. For weak coupling strengths to the central site, we identify a regime with a logarithmic in time transport of particles and information.


title: A counter-image to the gravitational arc in Abell 1201: Evidence for IMF variations, or a $10^{10}$M$_{\odot}$ black hole?

abstract: Abell 1201 is a massive galaxy cluster at z=0.169 with a brightest cluster galaxy (BCG) that acts as a gravitational lens to a background source at z=0.451. The lensing configuration is unusual, with a single bright arc formed at small radius $\sim$2 arcsec), where stars and dark matter are both expected to contribute substantially to the total lensing mass. Here, we present deep spectroscopic observations of the Abell 1201 BCG with MUSE, which reveal emission lines from a faint counter-image, opposite to the main arc, at a radius of 0.6 arcsec. We explore models in which the lensing mass is described by a combination of stellar mass and a standard dark-matter halo. The counter-image is not predicted in such models, unless the dark-matter component is negligible, which would imply an extremely heavy stellar initial mass function (IMF) in this galaxy. We consider two modifications to the model which can produce the observed configuration without resorting to extreme IMFs. Imposing a radial gradient in the stellar mass-to-light ratio, $\Upsilon$, can generate a counter-image close to the observed position if $\Upsilon$ increases by $\gtrsim$60 per cent within the inner $\sim$1 arcsec (e.g. variation from a Milky-Way-like to a Salpeter-like IMF). Alternatively, the counter-image can be produced by introducing a central super-massive black hole. The required mass is $M_{\rm BH}$ = (1.3$\pm$0.6)$\times$10$^{10}$ M$_\odot$, which is comparable to the largest black holes known to date, several of which are also hosted by BCGs. We comment on future observations which promise to distinguish between these alternatives.


title: Reconciling the diffuse Galactic gamma-ray and the cosmic ray spectra

abstract: Most of the diffuse Galactic GeV gamma-ray emission is produced via collisions of cosmic ray (CR) protons with ISM protons. As such the observed spectra of the gamma-rays and the CRs should be strongly linked. Recent observations of Fermi-LAT exhibit a hardening of the gamma-ray spectrum at around a hundred GeV, between the Sagittarius and Carina tangents, and a further hardening at a few degrees above and below the Galactic plane. However, standard CR propagation models that assume a time independent source distribution and a location independent diffusion cannot give rise to a spatially dependent CR (and hence gamma-ray) spectral slopes. Here we consider a dynamic spiral arm model in which the distribution of CR sources is concentrated in the (dynamic) spiral arms, and we study the effects of this model on the $\pi^0$-decay produced gamma-ray spectra. Within this model, near the Galactic arms the observed gamma-ray spectral slope is not trivially related to the CR injection spectrum and energy dependence of the diffusion coefficient. We find unique signatures that agree with the Fermi-LAT observations. This model also provides a physical explanation for the difference between the local CR spectral slope and the CR slope inferred from the average gamma-ray spectrum.


title: The Planetary Accretion Shock: I. Framework for Radiation-hydrodynamical Simulations and First Results

abstract: The key aspect determining the post-formation luminosity of gas giants has long been considered to be the energetics of the accretion shock at the planetary surface. We use 1D radiation-hydrodynamical simulations to study the radiative loss efficiency and to obtain post-shock temperatures and pressures and thus entropies. The efficiency is defined as the fraction of the total incoming energy flux which escapes the system (roughly the Hill sphere), taking into account the energy recycling which occurs ahead of the shock in a radiative precursor. We focus here on a constant equation of state to isolate the shock physics but use constant and tabulated opacities. While robust quantitative results will require a self-consistent treatment including hydrogen dissocation and ionization, the results show the correct qualitative behavior and can be understood semi-analytically. The shock is found to be isothermal and supercritical for a range of conditions relevant to core accretion (CA), with Mach numbers greater than ca. 3. Across the shock, the entropy decreases significantly, by a few entropy units (k_B/baryon). While nearly 100 percent of the incoming kinetic energy is converted to radiation locally, the efficiencies are found to be as low as roughly 40 percent, implying that a meaningful fraction of the total accretion energy is brought into the planet. For realistic parameter combinations in the CA scenario, a non-zero fraction of the luminosity always escapes the system. This luminosity could explain, at least in part, recent observations in the LkCa 15 and HD 100546 systems.


title: Physical Properties of Molecular Clouds at 2 parsec Resolution in the Low-Metallicity Dwarf Galaxy NGC 6822 and the Milky Way

abstract: We present the ALMA survey of CO(2-1) emission from the 1/5 solar metallicity, Local Group dwarf galaxy NGC 6822. We achieve high (0.9 arcsec ~ 2 pc) spatial resolution while covering large area: four 250 pc x 250 pc regions that encompass ~2/3 of NGC 6822's star formation. In these regions, we resolve ~150 compact CO clumps that have small radii (~2-3 pc), narrow line width (~1 km/s), and low filling factor across the galaxy. This is consistent with other recent studies of low metallicity galaxies, but here shown with a 15 times larger sample. At parsec scales, CO emission correlates with 8 micron emission better than with 24 micron emission and anti-correlates with Halpha, so that PAH emission may be an effective tracer of molecular gas at low metallicity. The properties of the CO clumps resemble those of similar-size structures in Galactic clouds except of slightly lower surface brightness and CO-to-H2 ratio ~1-2 times the Galactic value. The clumps exist inside larger atomic-molecular complexes with masses typical for giant molecular cloud. Using dust to trace H2 for the entire complex, we find CO-to-H2 to be ~20-25 times the Galactic value, but with strong dependence on spatial scale and variations between complexes that may track their evolutionary state. The H2-to-HI ratio is low globally and only mildly above unity within the complexes. The SFR-to-H2 ratio is ~3-5 times higher in the complexes than in massive disk galaxies, but after accounting for the bias from targeting star-forming regions, we conclude that the global molecular gas depletion time may be as long as in massive disk galaxies.


title: First Light: exploring the Spectra of High-Redshift Galaxies in the Renaissance Simulations

abstract: We present synthetic observations for the first generations of galaxies in the Universe and make predictions for future deep field observations for redshifts greater than 6. Due to the strong impact of nebular emission lines and the relatively compact scale of HII regions, high resolution cosmological simulations and a robust suite of analysis tools are required to properly simulate spectra. We created a software pipeline consisting of FSPS, Hyperion, Cloudy and our own tools to generate synthetic IR observations from a fully three-dimensional arrangement of gas, dust, and stars. Our prescription allows us to include emission lines for a complete chemical network and tackle the effect of dust extinction and scattering in the various lines of sight. We provide spectra, 2-D binned photon imagery for both HST and JWST IR filters, luminosity relationships, and emission line strengths for a large sample of high redshift galaxies in the Renaissance Simulations. Our resulting synthetic spectra show high variability between galactic halos with a strong dependence on stellar mass, metallicity, gas mass fraction, and formation history. Halos with the lowest stellar mass have the greatest variability in [OIII]/H$\beta$, [OIII] and CIII] while halos with higher masses are seen to show consistency in their spectra and [OIII] equivalent widths (EW) between 1\AA\ and 10\AA. Viewing angle accounted for three-fold difference in flux due to the presence of ionized gas channels in a halo. Furthermore, JWST color plots show a discernible relationship between redshift, color, and mean stellar age.


title: The effect of ram pressure on the molecular gas of galaxies: three case studies in the Virgo cluster

abstract: We present 12CO (2-1) data of three Virgo spirals - NGC 4330, NGC 4402 and NGC 4522 obtained using the Submillimeter Array. These three galaxies show clear evidence of ram pressure stripping due to the cluster medium as found in previous HI imaging studies. Using high-resolution CO data, we investigate how the properties of the inner molecular gas disc change while a galaxy is undergoing HI stripping in the cluster. At given sensitivity limits, we do not find any clear signs of molecular gas stripping. However, both its morphology and kinematics appear to be quite disturbed as those of HI. Morphological peculiarities present in the molecular and atomic gas are closely related with each other, suggesting that molecular gas can be also affected by strong ICM pressure even if it is not stripped. CO is found to be modestly enhanced along the upstream sides in these galaxies, which may change the local star formation activity in the disc. Indeed, the distribution of H$\alpha$ emission, a tracer of recent star formation, well coincides with that of the molecular gas, revealing enhancements near the local CO peak or along the CO compression. FUV and H$\alpha$ share some properties in common, but FUV is always more extended than CO/H$\alpha$ in the three galaxies, implying that the star-forming disc is rapidly shrinking as the molecular gas properties have changed. We discuss how ICM pressure affects dense molecular gas and hence star formation properties while diffuse atomic gas is being removed from a galaxy.


title: Vestigial nematicity from spin and/or charge order in the cuprates

abstract: Nematic order has manifested itself in a variety of materials in the cuprate family. We propose an effective field theory of a layered system with incommensurate, intertwined spin- and charge-density wave (SDW and CDW) orders, each of which consists of two components related by $C_4$ rotations. Using a variational method (which is exact in a large $N$ limit), we study the development of nematicity from partially melting those density waves by either increasing temperature or adding quenched disorder. As temperature decreases we first find a transition to a nematic phase, but depending on the range of parameters (e.g. doping concentration) the strongest fluctuations associated with this phase reflect either proximate SDW or CDW order. We also discuss the changes in parameters that can account for the differences in the SDW-CDW interplay between the (214) family and the other hole-doped cuprates.


title: Equation of State Effects on Gravitational Waves from Rotating Core Collapse

abstract: Gravitational waves (GWs) generated by axisymmetric rotating collapse, bounce, and early postbounce phases of a galactic core-collapse supernova will be detectable by current-generation gravitational wave observatories. Since these GWs are emitted from the quadrupole-deformed nuclear-density core, they may encode information on the uncertain nuclear equation of state (EOS). We examine the effects of the nuclear EOS on GWs from rotating core collapse and carry out 1824 axisymmetric general-relativistic hydrodynamic simulations that cover a parameter space of 98 different rotation profiles and 18 different EOS. We show that the bounce GW signal is largely independent of the EOS and sensitive primarily to the ratio of rotational to gravitational energy, and at high rotation rates, to the degree of differential rotation. The GW frequency of postbounce core oscillations shows stronger EOS dependence that can be parameterized by the core's EOS-dependent dynamical frequency $\sqrt{G\bar{\rho}_c}$. We find that the ratio of the peak frequency to the dynamical frequency follows a universal trend that is obeyed by all EOS and rotation profiles and that indicates that the nature of the core oscillations changes when the rotation rate exceeds the dynamical frequency. We find that differences in the treatments of low-density nonuniform nuclear matter, of the transition from nonuniform to uniform nuclear matter, and in the description of nuclear matter up to around twice saturation density can mildly affect the GW signal. We find that approximations and uncertainties in electron capture rates can lead to variations in the GW signal that are of comparable magnitude to those due to different nuclear EOS. This emphasizes the need for reliable nuclear electron capture rates and for self-consistent multi-dimensional neutrino radiation-hydrodynamic simulations of rotating core collapse.


title: The local spiral arm in the LAMOST-Gaia common stars?

abstract: Using the LAMOST-Gaia common stars, we demonstrate that the in-plane velocity field for the nearby young stars are significantly different from that for the old ones. For the young stars, the probably perturbed velocities similar to the old population are mostly removed from the velocity maps in the $X$--$Y$ plane. The residual velocity field shows that the young stars consistently move along $Y$ with faster $v_\phi$ at the trailing side of the local arm, while at the leading side, they move slower in azimuth direction. At both sides, the young stars averagely move inward with $v_R$ of $-5\sim-3$ km s$^{-1}$. The divergence of the velocity in $Y$ direction implies that the young stars are associated with a density wave nearby the local arm. We therefore suggest that the young stars may reflect the formation of the local spiral arm by correlating themselves with a density wave. The range of the age for the young stars is around 2 Gyr, which is sensible since the transient spiral arm can sustain that long. We also point out that alternative explanations of the peculiar velocity field for the young population cannot be ruled out only from this observed data.