Source author record

Leonid Petrov

Leonid Petrov appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

42works
16topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

42 published item(s)

preprint2022arXiv

Irreversible Markov Dynamics and Hydrodynamics for KPZ States in the Stochastic Six Vertex Model

We introduce a family of Markov growth processes on discrete height functions defined on the 2-dimensional square lattice. Each height function corresponds to a configuration of the six vertex model on the infinite square lattice. We focus on the stochastic six vertex model corresponding to a particular two-parameter family of weights within the ferroelectric regime. It is believed (and partially proven, see Aggarwal, arXiv:2004.13272) that the stochastic six vertex model displays nontrivial pure (i.e., translation invariant and ergodic) Gibbs states of two types, KPZ and liquid. These phases have very different long-range correlation structure. The Markov processes we construct preserve the KPZ pure states in the full plane. We also show that the same processes put on the torus preserve arbitrary Gibbs measures for generic six vertex weights (not necessarily in the ferroelectric regime). Our dynamics arise naturally from the Yang-Baxter equation for the six vertex model via its bijectivisation, a technique first used in Bufetov-Petrov (arXiv:1712.04584). The dynamics we construct are irreversible; in particular, the height function has a nonzero average drift. In each KPZ pure state, we explicitly compute the average drift (also known as the current) as a function of the slope. We use this to analyze the hydrodynamics of a non-stationary version of our process acting on quarter plane stochastic six vertex configurations. The fixed-time limit shapes in the quarter plane model were obtained in Borodin-Corwin-Gorin (arXiv:1407.6729).

preprint2021arXiv

Mapping TASEP back in time

We obtain a new relation between the distributions $μ_t$ at different times $t\ge 0$ of the continuous-time TASEP (Totally Asymmetric Simple Exclusion Process) started from the step initial configuration. Namely, we present a continuous-time Markov process with local interactions and particle-dependent rates which maps the TASEP distributions $μ_t$ backwards in time. Under the backwards process, particles jump to the left, and the dynamics can be viewed as a version of the discrete-space Hammersley process. Combined with the forward TASEP evolution, this leads to a stationary Markov dynamics preserving $μ_t$ which in turn brings new identities for expectations with respect to $μ_t$. The construction of the backwards dynamics is based on Markov maps interchanging parameters of Schur processes, and is motivated by bijectivizations of the Yang-Baxter equation. We also present a number of corollaries, extensions, and open questions arising from our constructions.

preprint2021arXiv

Parameter Permutation Symmetry in Particle Systems and Random Polymers

Many integrable stochastic particle systems in one space dimension (such as TASEP - totally asymmetric simple exclusion process - and its various deformations, with a notable exception of ASEP) remain integrable when we equip each particle $x_i$ with its own jump rate parameter $ν_i$. It is a consequence of integrability that the distribution of each particle $x_n(t)$ in a system started from the step initial configuration depends on the parameters $ν_j$, $j\le n$, in a symmetric way. A transposition $ν_n \leftrightarrow ν_{n+1}$ of the parameters thus affects only the distribution of $x_n(t)$. For $q$-Hahn TASEP and its degenerations ($q$-TASEP and directed beta polymer) we realize the transposition $ν_n \leftrightarrow ν_{n+1}$ as an explicit Markov swap operator acting on the single particle $x_n(t)$. For beta polymer, the swap operator can be interpreted as a simple modification of the lattice on which the polymer is considered. Our main tools are Markov duality and contour integral formulas for joint moments. In particular, our constructions lead to a continuous time Markov process $\mathsf{Q}^{(\mathsf{t})}$ preserving the time $\mathsf{t}$ distribution of the $q$-TASEP (with step initial configuration, where $\mathsf{t}\in \mathbb{R}_{>0}$ is fixed). The dual system is a certain transient modification of the stochastic $q$-Boson system. We identify asymptotic survival probabilities of this transient process with $q$-moments of the $q$-TASEP, and use this to show the convergence of the process $\mathsf{Q}^{(\mathsf{t})}$ with arbitrary initial data to its stationary distribution. Setting $q=0$, we recover the results about the usual TASEP established recently in [arXiv:1907.09155] by a different approach based on Gibbs ensembles of interlacing particles in two dimensions.

preprint2020arXiv

Influence of the Galactic gravitational field on the positional accuracy of extragalactic sources. II Observational appearances and detectability

We consider a possibility of detecting the jitter effect of apparent celestial positions of distant sources due to local fluctuations of the Galaxy gravitational field. It is proposed to observe two samples of extragalactic sources (target and control) in different sky directions using the high-precision radio interferometry. It is shown that on a scale of ~2 years, it is possible to detect a systematic increase in the standard deviation of measured arc lengths of pairs of target sources compared to the control ones at the $3σ$-level if the accuracy of differential astrometric observations is around 10 $μ$as. For the current state-of-the-art accuracy of 30 $μ$as achieved at the KVN or VERA interferometers, which have shorter baselines in comparison with VLBI, the target and control samples will differ only at the 2$σ$-level on the scale of 10 years. To achieve the $3σ$-level on this time interval, it is necessary to improve the accuracy up to ~20 $μ$as. Other possible effects that can also affect the arc length measurements between two sources are discussed, and an observational strategy to minimize them is suggested.

preprint2020arXiv

PushTASEP in inhomogeneous space

We consider the PushTASEP (pushing totally asymmetric simple exclusion process, also sometimes called long-range TASEP) with the step initial configuration evolving in an inhomogeneous space. That is, the rate of each particle's jump depends on the location of this particle. We match the distribution of the height function of this PushTASEP with Schur processes. Using this matching and determinantal structure of Schur processes, we obtain limit shape and fluctuation results which are typical for stochastic particle systems in the Kardar-Parisi-Zhang universality class. PushTASEP is a close relative of the usual TASEP. In inhomogeneous space the former is integrable, while the integrability of the latter is not known.

preprint2020arXiv

Spin q-Whittaker polynomials and deformed quantum Toda

Spin $q$-Whittaker symmetric polynomials labeled by partitions $λ$ were recently introduced by Borodin and Wheeler (arXiv:1701.06292) in the context of integrable $\mathfrak{sl}_2$ vertex models. They are a one-parameter deformation of the $t=0$ Macdonald polynomials. We present a new, more convenient modification of spin $q$-Whittaker polynomials and find two Macdonald type $q$-difference operators acting diagonally in these polynomials with eigenvalues, respectively, $q^{-λ_1}$ and $q^{λ_N}$ (where $λ$ is the polynomial's label). We study probability measures on interlacing arrays based on spin $q$-Whittaker polynomials, and match their observables with known stochastic particle systems such as the $q$-Hahn TASEP. In a scaling limit as $q\nearrow 1$, spin $q$-Whittaker polynomials turn into a new one-parameter deformation of the $\mathfrak{gl}_n$ Whittaker functions. The rescaled Pieri type rule gives rise to a one-parameter deformation of the quantum Toda Hamiltonian. The deformed Hamiltonian acts diagonally on our new spin Whittaker functions. On the stochastic side, as $q\nearrow 1$ we discover a multilevel extension of the beta polymer model of Barraquand and Corwin (arXiv:1503.04117), and relate it to spin Whittaker functions.

preprint2020arXiv

VLBI20-30: a scientific roadmap for the next decade -- The future of the European VLBI Network

This white paper describes the science case for Very Long Baseline Interferometry (VLBI) and provides suggestions towards upgrade paths for the European VLBI Network (EVN). The EVN is a distributed long-baseline radio interferometric array, that operates at the very forefront of astronomical research. Recent results, together with the new science possibilities outlined in this vision document, demonstrate the EVN's potential to generate new and exciting results that will transform our view of the cosmos. Together with e-MERLIN, the EVN provides a range of baseline lengths that permit unique studies of faint radio sources to be made over a wide range of spatial scales. The science cases are reviewed in six chapters that cover the following broad areas: cosmology, galaxy formation and evolution, innermost regions of active galactic nuclei, explosive phenomena and transients, stars and stellar masers in the Milky Way, celestial reference frames and space applications. The document concludes with identifying the synergies with other radio, as well as multi-band/multi-messenger instruments, and provide the recommendations for future improvements. The appendices briefly describe other radio VLBI arrays, the technological framework for EVN developments, and a selection of spectral lines of astrophysical interest below 100 GHz. The document includes a glossary for non-specialists, and a list of acronyms at the end.

preprint2019arXiv

Generalizations of TASEP in discrete and continuous inhomogeneous space

We investigate a rich new class of exactly solvable particle systems generalizing the Totally Asymmetric Simple Exclusion Process (TASEP). Our particle systems can be thought of as new exactly solvable examples of tandem queues, directed first- or last-passage percolation models, or Robinson-Schensted-Knuth type systems with random input. One of the novel features of the particle systems is the presence of spatial inhomogeneity which can lead to the formation of traffic jams. For systems with special step-like initial data, we find explicit limit shapes, describe hydrodynamic evolution, and obtain asymptotic fluctuation results which put the systems into the Kardar-Parisi-Zhang universality class. At a critical scaling around a traffic jam in the continuous space TASEP, we observe deformations of the Tracy-Widom distribution and the extended Airy kernel, revealing the finer structure of this novel type of phase transitions. A homogeneous version of a discrete space system we consider is a one-parameter deformation of the geometric last-passage percolation, and we obtain extensions of the limit shape parabola and the corresponding asymptotic fluctuation results. The exact solvability and asymptotic behavior results are powered by a new nontrivial connection to Schur measures and processes.

preprint2016arXiv

Higher spin six vertex model and symmetric rational functions

We consider a fully inhomogeneous stochastic higher spin six vertex model in a quadrant. For this model we derive concise integral representations for multi-point q-moments of the height function and for the q-correlation functions. At least in the case of the step initial condition, our formulas degenerate in appropriate limits to many known formulas of such type for integrable probabilistic systems in the (1+1)d KPZ universality class, including the stochastic six vertex model, ASEP, various q-TASEPs, and associated zero range processes. Our arguments are largely based on properties of a family of symmetric rational functions which can be defined as partition functions of the inhomogeneous higher spin six vertex model for suitable domains. In the homogeneous case, such functions were previously studied in http://arxiv.org/abs/1410.0976; they also generalize classical Hall-Littlewood and Schur polynomials. A key role is played by Cauchy-like summation identities for these functions, which are obtained as a direct corollary of the Yang-Baxter equation for the higher spin six vertex model.

preprint2016arXiv

Lectures on Integrable probability: Stochastic vertex models and symmetric functions

We consider a homogeneous stochastic higher spin six vertex model in a quadrant. For this model we derive concise integral representations for multi-point q-moments of the height function and for the q-correlation functions. At least in the case of the step initial condition, our formulas degenerate in appropriate limits to many known formulas of such type for integrable probabilistic systems in the (1+1)d KPZ universality class, including the stochastic six vertex model, ASEP, various q-TASEPs, and associated zero range processes. Our arguments are largely based on properties of a family of symmetric rational functions (introduced in arXiv:1410.0976) that can be defined as partition functions of the higher spin six vertex model for suitable domains; they generalize classical Hall-Littlewood and Schur polynomials. A key role is played by Cauchy-like summation identities for these functions, which are obtained as a direct corollary of the Yang-Baxter equation for the higher spin six vertex model. These are lecture notes for a course given by A.B. at the Ecole de Physique des Houches in July of 2015. All the results and proofs presented here generalize to the setting of the fully inhomogeneous higher spin six vertex model, see arXiv:1601.05770 for a detailed exposition of the inhomogeneous case.

preprint2016arXiv

q-randomized Robinson-Schensted-Knuth correspondences and random polymers

We introduce and study q-randomized Robinson-Schensted-Knuth (RSK) correspondences which interpolate between the classical (q=0) and geometric (q->1) RSK correspondences (the latter ones are sometimes also called tropical). For 0<q<1 our correspondences are randomized, i.e., the result of an insertion is a certain probability distribution on semistandard Young tableaux. Because of this randomness, we use the language of discrete time Markov dynamics on two-dimensional interlacing particle arrays (these arrays are in a natural bijection with semistandard tableaux). Our dynamics act nicely on a certain class of probability measures on arrays, namely, on q-Whittaker processes (which are t=0 versions of Macdonald processes). We present four Markov dynamics which for q=0 reduce to the classical row or column RSK correspondences applied to a random input matrix with independent geometric or Bernoulli entries. Our new two-dimensional discrete time dynamics generalize and extend several known constructions: (1) The discrete time q-TASEPs arise as one-dimensional marginals of our "column" dynamics. In a similar way, our "row" dynamics lead to discrete time q-PushTASEPs - new integrable particle systems in the Kardar-Parisi-Zhang universality class. We employ these new one-dimensional discrete time systems to establish a Fredholm determinantal formula for the two-sided continuous time q-PushASEP conjectured by Corwin-Petrov (2013). (2) In a certain Poisson-type limit (from discrete to continuous time), our two-dimensional dynamics reduce to the q-randomized column and row Robinson-Schensted correspondences introduced by O'Connell-Pei (2012) and Borodin-Petrov (2013), respectively. (3) In a scaling limit as q->1, two of our four dynamics on interlacing arrays turn into the geometric RSK correspondences associated with log-Gamma or strict-weak directed random polymers.

preprint2016arXiv

Stochastic higher spin six vertex model and q-TASEPs

We present two new connections between the inhomogeneous stochastic higher spin six vertex model in a quadrant and integrable stochastic systems from the Macdonald processes hierarchy. First, we show how Macdonald $q$-difference operators with $t=0$ (an algebraic tool crucial for studying the corresponding Macdonald processes) can be utilized to get $q$-moments of the height function $\mathfrak{h}$ in the higher spin six vertex model first computed in arXiv:1601.05770 using Bethe ansatz. This result in particular implies that for the vertex model with the step Bernoulli boundary condition, the value of $\mathfrak{h}$ at an arbitrary point $(N+1,T)\in\mathbb{Z}_{\ge2}\times\mathbb{Z}_{\ge1}$ has the same distribution as the last component $λ_N$ of a random partition under a specific $t=0$ Macdonald measure. On the other hand, it is known that $\mathbf{x}_N:=λ_N-N$ can be identified with the location of the $N$th particle in a certain discrete time $q$-TASEP started from the step initial configuration. The second construction we present is a coupling of this $q$-TASEP and the higher spin six vertex model (with the step Bernoulli boundary condition) along time-like paths providing an independent probabilistic explanation of the equality of $\mathfrak{h}(N+1,T)$ and $\mathbf{x}_N+N$ in distribution. Combined with the identification of averages of observables between the stochastic higher spin six vertex model and Schur measures (which are $t=q$ Macdonald measures) obtained recently in arXiv:1608.01553, this produces GUE Tracy--Widom asymptotics for a discrete time $q$-TASEP with the step initial configuration and special jump parameters.

preprint2015arXiv

1FGL J1417.7-4407: A likely gamma-ray bright binary with a massive neutron star and a giant secondary

We present multiwavelength observations of the persistent Fermi-LAT unidentified gamma-ray source 1FGL J1417.7-4407, showing it is likely to be associated with a newly discovered X-ray binary containing a massive neutron star (nearly 2 M_sun) and a ~ 0.35 M_sun giant secondary with a 5.4 day period. SOAR optical spectroscopy at a range of orbital phases reveals variable double-peaked H-alpha emission, consistent with the presence of an accretion disk. The lack of radio emission and evidence for a disk suggests the gamma-ray emission is unlikely to originate in a pulsar magnetosphere, but could instead be associated with a pulsar wind, relativistic jet, or could be due to synchrotron self-Compton at the disk--magnetosphere boundary. Assuming a wind or jet, the high ratio of gamma-ray to X-ray luminosity (~ 20) suggests efficient production of gamma-rays, perhaps due to the giant companion. The system appears to be a low-mass X-ray binary that has not yet completed the pulsar recycling process. This system is a good candidate to monitor for a future transition between accretion-powered and rotational-powered states, but in the context of a giant secondary.

preprint2015arXiv

Asymptotics of uniformly random lozenge tilings of polygons. Gaussian free field

We study large-scale height fluctuations of random stepped surfaces corresponding to uniformly random lozenge tilings of polygons on the triangular lattice. For a class of polygons (which allows arbitrarily large number of sides), we show that these fluctuations are asymptotically governed by a Gaussian free (massless) field. This complements the similar result obtained in Kenyon [Comm. Math. Phys. 281 (2008) 675-709] about tilings of regions without frozen facets of the limit shape. In our asymptotic analysis we use the explicit double contour integral formula for the determinantal correlation kernel of the model obtained previously in Petrov [Asymptotics of random lozenge tilings via Gelfand-Tsetlin schemes (2012) Preprint].

preprint2015arXiv

First scientific VLBI observations using New Zealand 30 metre radio telescope WARK30M

We report the results of a successful 24 hour 6.7 GHz VLBI experiment using the 30 meter radio telescope WARK30M near Warkworth, New Zealand, recently converted from a radio telecommunications antenna, and two radio telescopes located in Australia: Hobart 26-m and Ceduna 30-m. The geocentric position of WARK30M is determined with a 100 mm uncertainty for the vertical component and 10 mm for the horizontal components. We report correlated flux densities at 6.7 GHz of 175 radio sources associated with Fermi gamma-ray sources. A parsec scale emission from the radio source 1031-837 is detected, and its association with the gamma-ray object 2FGL J1032.9-8401 is established with a high likelihood ratio. We conclude that the new Pacific area radio telescope WARK30M is ready to operate for scientific projects.

preprint2015arXiv

Law of Large Numbers for Infinite Random Matrices over a Finite Field

Asymptotic representation theory of general linear groups GL(n,q) over a finite field leads to studying probability measures ρon the group U of all infinite uni-uppertriangular matrices over F_q, with the condition that ρis invariant under conjugations by arbitrary infinite matrices. Such probability measures form an infinite-dimensional simplex, and the description of its extreme points (in other words, ergodic measures ρ) was conjectured by Kerov in connection with nonnegative specializations of Hall-Littlewood symmetric functions. Vershik and Kerov also conjectured the following Law of Large Numbers. Consider an n by n diagonal submatrix of the infinite random matrix drawn from an ergodic measure coming from the Kerov's conjectural classification. The sizes of Jordan blocks of the submatrix can be interpreted as a (random) partition of n, or, equivalently, as a (random) Young diagram λ(n) with n boxes. Then, as n goes to infinity, the rows and columns of λ(n) have almost sure limiting frequencies corresponding to parameters of this ergodic measure. Our main result is the proof of this Law of Large Numbers. We achieve it by analyzing a new randomized Robinson-Schensted-Knuth (RSK) insertion algorithm which samples random Young diagrams λ(n) coming from ergodic measures. The probability weights of these Young diagrams are expressed in terms of Hall-Littlewood symmetric functions. Our insertion algorithm is a modified and extended version of a recent construction by Borodin and the second author (arXiv:1305.5501). On the other hand, our randomized RSK insertion generalizes a version of the RSK insertion introduced by Vershik and Kerov (1986) in connection with asymptotic representation theory of symmetric groups (which is governed by nonnegative specializations of Schur symmetric functions).

preprint2015arXiv

Modeling of path delay in the neutral atmosphere: a paradigm shift

Computation of propagation effects in the neutral atmosphere, namely path delay, extinction, and bending angle is a trivial task provided the 4D state of the atmosphere is known. Unfortunately, the mixing ratio of water vapor is highly variable and it cannot be deduced from surface measurements. That fact led to a paradigm that considers path delay and extinction in the atmosphere as a~priori unknown quantities that have to be evaluated from the radio astronomy data themselves. Development of our ability to model the atmosphere and to digest humongous outputs of these models that took place over the course of the 21st century changed the game. Using the publicly available output of operational numerical weather model GEOS run by NASA, we are in a position to compute path delay through the neutral atmosphere for any station and for any epoch from 1979 through now with accuracy of 45 ps * cosec elevation. We are in a position to compute extinction with accuracy better than 10 pro cents. We are in a position to do it routinely, in a similar way how we update apparent star positions for precession and nutation. Moreover, we are in a position to do it now. As a demonstration of current capabilities, I have computed time series of path delays for all radiotelecopes that I was aware of (220 sites) since 1979 with a step 3-6 hours. Results of the validation tests are presented. A new paradigm of data analysis assumes that we know the atmosphere propagation effects a priori with the accuracy higher that one could deduce them from radio astronomy observations.

preprint2015arXiv

New Associations of Gamma-Ray Sources from the Fermi Second Source Catalog

We present the results of an all-sky radio survey between 5 and 9 GHz of the fields surrounding all unassociated gamma-ray objects listed in the Fermi Large Area Telescope Second Source Catalog (2FGL). The goal of these observations is to find all new gamma-ray AGN associations with radio sources >10 mJy at 8 GHz. We observed with the Very Large Array and the Australia Telescope Compact Array the areas around unassociated sources, providing localizations of weak radio point sources found in 2FGL fields at arcmin scales. Then we followed-up a subset of those with the Very Long Baseline and the Long Baseline Arrays to confirm detections of radio emission on parsec-scales. We quantified association probabilities based on known statistics of source counts and assuming a uniform distribution of background sources. In total we found 865 radio sources at arcsec scales as candidates for association and detected 95 of 170 selected for follow-up observations at milliarcsecond resolution. Based on this we obtained firm associations for 76 previously unknown gamma-ray AGN. Comparison of these new AGN associations with the predictions from using the WISE color-color diagram shows that half of the associations are missed. We found that 129 out of 588 observed gamma-ray sources at arcmin scales not a single radio continuum source was detected above our sensitivity limit within the 3-sigma gamma-ray localization. These "empty" fields were found to be particularly concentrated at low Galactic latitudes. The nature of these Galactic gamma-ray emitters is not yet determined.

preprint2015arXiv

The International Mass Loading Service

The International Mass Loading Service computes four loadings: a) atmospheric pressure loading; b) land water storage loading; c) oceanic tidal loading; and d) non-tidal oceanic loading. The service provides to users the mass loading time series in three forms: 1) pre-computed time series for a list of 849 space geodesy stations; 2) pre-computed time series on the global 1deg x 1deg grid; and 3) on-demand Internet service for a list of stations and a time range specified by the user. The loading displacements are provided for the time period from 1979.01.01 through present, updated on an hourly basis, and have latencies 8-20 hours.

preprint2015arXiv

The q-PushASEP: A New Integrable Model for Traffic in 1+1 Dimension

We introduce a new interacting (stochastic) particle system q-PushASEP which interpolates between the q-TASEP introduced by Borodin and Corwin (see arXiv:1111.4408, and also arXiv:1207.5035; arXiv:1305.2972; arXiv:1212.6716) and the q-PushTASEP introduced recently by Borodin and Petrov (arXiv:1305.5501). In the q-PushASEP, particles can jump to the left or to the right, and there is a certain partially asymmetric pushing mechanism present. This particle system has a nice interpretation as a model of traffic on a one-lane highway in which cars are able to accelerate or slow down. Using the quantum many body system approach, we explicitly compute the expectations of a large family of observables for this system in terms of nested contour integrals. We also discuss relevant Fredholm determinantal formulas for the distribution of the location of each particle, and connections of the model with a certain two-sided version of Macdonald processes and with the semi-discrete stochastic heat equation.

preprint2013arXiv

Asymptotics of Random Lozenge Tilings via Gelfand-Tsetlin Schemes

A Gelfand-Tsetlin scheme of depth N is a triangular array with m integers at level m, m=1,...,N, subject to certain interlacing constraints. We study the ensemble of uniformly random Gelfand-Tsetlin schemes with arbitrary fixed N-th row. We obtain an explicit double contour integral expression for the determinantal correlation kernel of this ensemble (and also of its q-deformation). This provides new tools for asymptotic analysis of uniformly random lozenge tilings of polygons on the triangular lattice; or, equivalently, of random stepped surfaces. We work with a class of polygons which allows arbitrarily large number of sides. We show that the local limit behavior of random tilings (as all dimensions of the polygon grow) is directed by ergodic translation invariant Gibbs measures. The slopes of these measures coincide with the ones of tangent planes to the corresponding limit shapes described by Kenyon and Okounkov in arXiv:math-ph/0507007. We also prove that at the edge of the limit shape, the asymptotic behavior of random tilings is given by the Airy process. In particular, our results cover the most investigated case of random boxed plane partitions (when the polygon is a hexagon).

preprint2013arXiv

Australia Telescope Compact Array observations of Fermi unassociated sources

We report results of the first phase of observations with the Australia Telescope Compact Array (ATCA) at 5 and 9 GHz of the fields around 411 gamma-ray sources with declinations < +10 deg detected by Fermi but marked as unassociated in the 2FGL catalogue. We have detected 424 sources with flux densities in a range of 2 mJy to 6 Jy that lie within the 99 per cent localisation uncertainty of 283 gamma-ray sources. Of these, 146 objects were detected in both the 5 and 9 GHz bands. We found 84 sources in our sample with a spectral index flatter than -0.5. The majority of detected sources are weaker than 100 mJy and for this reason were not found in previous surveys. Approximately 1/3 of our sample, 128 objects, have the probability of being associated by more than 10 times than the probability of being a background source found in the vicinity of a gamma-ray object by chance. We present the catalogue of positions of these sources, estimates of their flux densities and spectral indices where available.

preprint2013arXiv

Nearest neighbor Markov dynamics on Macdonald processes

Macdonald processes are certain probability measures on two-dimensional arrays of interlacing particles introduced by Borodin and Corwin (arXiv:1111.4408 [math.PR]). They are defined in terms of nonnegative specializations of the Macdonald symmetric functions and depend on two parameters (q,t), where 0<= q, t < 1. Our main result is a classification of continuous time, nearest neighbor Markov dynamics on the space of interlacing arrays that act nicely on Macdonald processes. The classification unites known examples of such dynamics and also yields many new ones. When t = 0, one dynamics leads to a new integrable interacting particle system on the one-dimensional lattice, which is a q-deformation of the PushTASEP (= long-range TASEP). When q = t, the Macdonald processes become the Schur processes of Okounkov and Reshetikhin (arXiv:math/0107056 [math.CO]). In this degeneration, we discover new Robinson--Schensted-type correspondences between words and pairs of Young tableaux that govern some of our dynamics.

preprint2013arXiv

The catalogue of positions of optically bright extragalactic radio sources OBRS-2

It is anticipated that future space-born missions, such as Gaia, will be able to determine in optical domain positions of more than 100,000 bright quasars with sub-mas accuracies that are comparable to very long baseline interferometry (VLBI) accuracies. Comparisons of coordinate systems from space-born missions and from VLBI will be very important, first for investigation of possible systematic errors, second for investigation of possible shift between centroids of radio and optical emissions in active galaxy nuclea. In order to make such a comparison more robust, a program of densification of the grid of radio sources detectable with both VLBI and Gaia was launched in 2006. In the second observing campaign a set of 290 objects from the list of 398 compact extragalactic radio sources with declinations greater -10 deg was observed with the VLBA+EVN in 2010-2011 with the primary goal of producing their images with milliarcsecond resolution. These sources are brighter than 18 magnitude at V band. In this paper coordinates of observed sources have been derived with milliarcsecond accuracies from analysis of these VLBI observations following the method of absolute astrometry and their images were produced. The catalogue of positions of 295 target sources and estimates of their correlated flux densities at 2.2 and 8.4 GHz is presented. The accuracies of source coordinates are in a range of 2 to 200 mas, with the median 3.2 mas.

preprint2012arXiv

Early science with Korean VLBI network: the QCAL-1 43GHz calibrator survey

This paper presents the catalog of correlated flux densities in three ranges of baseline projection lengths of 637 sources from a 43 GHz (Q-band) survey observed with the Korean VLBI Network. Of them, 623 sources have not been observed before at Q-band with VLBI. The goal of this work in the early science phase of the new VLBI array is twofold: to evaluate the performance of the new instrument that operates in a frequency range of 22-129 GHz and to build a list of objects that can be used as targets and as calibrators. We have observed the list of 799 target sources with declinations down to -40 degrees. Among them, 724 were observed before with VLBI at 22 GHz and had correlated flux densities greater than 200 mJy. The overall detection rate is 78%. The detection limit, defined as the minimum flux density for a source to be detected with 90% probability in a single observation, was in a range of 115-180 mJy depending on declination. However, some sources as weak as 70 mJy have been detected. Of 623 detected sources, 33 objects are detected for the first time in VLBI mode. We determined their coordinates with the median formal uncertainty 20 mas. The results of this work set the basis for future efforts to build the complete flux-limited sample of extragalactic sources at frequencies 22 GHz and higher at 3/4 of the celestial sphere.

preprint2012arXiv

The Boundary of the Gelfand-Tsetlin Graph: New Proof of Borodin-Olshanski's Formula, and its q-analogue

In the recent paper [arXiv:1109.1412], Borodin and Olshanski have presented a novel proof of the celebrated Edrei-Voiculescu theorem which describes the boundary of the Gelfand-Tsetlin graph as a region in an infinite-dimensional coordinate space. This graph encodes branching of irreducible characters of finite-dimensional unitary groups. Points of the boundary of the Gelfand-Tsetlin graph can be identified with finite indecomposable (= extreme) characters of the infinite-dimensional unitary group. An equivalent description identifies the boundary with the set of doubly infinite totally nonnegative sequences. A principal ingredient of Borodin-Olshanski's proof is a new explicit determinantal formula for the number of semi-standard Young tableaux of a given skew shape (or of Gelfand-Tsetlin schemes of trapezoidal shape). We present a simpler and more direct derivation of that formula using the Cauchy-Binet summation involving the inverse Vandermonde matrix. We also obtain a q-generalization of that formula, namely, a new explicit determinantal formula for arbitrary q-specializations of skew Schur polynomials. Its particular case is related to the q-Gelfand-Tsetlin graph and q-Toeplitz matrices introduced and studied by Gorin [arXiv:1011.1769].

preprint2011arXiv

First geodetic observations using new VLBI stations ASKAP-29 and WARK12M

We report the results of a successful 7 hour 1.4 GHz VLBI experiment using two new stations, ASKAP-29 located in Western Australia and WARK12M located on the North Island of New Zealand. This was the first geodetic VLBI observing session with the participation of these new stations. We have determined the positions of ASKAP-29 and WARK12M. Random errors on position estimates are 150-200 mm for the vertical component and 40-50 mm for the horizontal component. Systematic errors caused by the unmodeled ionosphere path delay may reach 1.3 m for the vertical component.

preprint2011arXiv

On Measures on Partitions Arising in Harmonic Analysis for Linear and Projective Characters of the Infinite Symmetric Group

The z-measures on partitions originated from the problem of harmonic analysis of linear representations of the infinite symmetric group in the works of Kerov, Olshanski and Vershik (1993, 2004). A similar family corresponding to projective representations was introduced by Borodin (1997). The latter measures live on strict partitions (i.e., partitions with distinct parts), and the z-measures are supported by all partitions. In this note we describe some combinatorial relations between these two families of measures using the well-known doubling of shifted Young diagrams.

preprint2011arXiv

Pfaffian Stochastic Dynamics of Strict Partitions

We study a family of continuous time Markov jump processes on strict partitions (partitions with distinct parts) preserving the distributions introduced by Borodin (1997) in connection with projective representations of the infinite symmetric group. The one-dimensional distributions of the processes (i.e., the Borodin's measures) have determinantal structure. We express the dynamical correlation functions of the processes in terms of certain Pfaffians and give explicit formulas for both the static and dynamical correlation kernels using the Gauss hypergeometric function. Moreover, we are able to express our correlation kernels (both static and dynamical) through those of the z-measures on partitions obtained previously by Borodin and Olshanski in a series of papers. The results about the fixed time case were announced in the author's note arXiv:1002.2714. A part of the present paper contains proofs of those results.

preprint2011arXiv

Precise absolute astrometry from the VLBA imaging and polarimetry survey at 5 GHz

We present in this paper accurate positions of 857 sources derived from the astrometric analysis of 16 eleven-hour experiments from the Very Long Baseline Array imaging and polarimetry survey at 5 GHz (VIPS). Among observed sources, positions of 430 objects were not determined before at a milliarcsecond level of accuracy. For 95% of the sources the uncertainty of their positions range from 0.3 to 0.9 mas, with the median value of 0.5 mas. This estimate of accuracy is substantiated by the comparison of positions of 386 sources that were previously observed in astrometric programs simultaneously at 2.3/8.6 GHz. Surprisingly, the ionosphere contribution to group delay was adequately modeled with the use of the total electron contents maps derived from GPS observations and only marginally affected estimates of source coordinates.

preprint2011arXiv

sl(2) Operators and Markov Processes on Branching Graphs

We present a unified approach to various examples of Markov dynamics on partitions studied by Borodin, Olshanski, Fulman, and the author. Our technique generalizes the Kerov's operators first appeared in [Okounkov, arXiv:math/0002135], and also stems from the study of duality of graded graphs in [Fomin, 1994]. Our main object is a countable branching graph carrying an sl(2,C)-module of a special kind. Using this structure, we introduce distinguished probability measures on the floors of the graph, and define two related types of Markov dynamics associated with these measures. We study spectral properties of the dynamics, and our main result is the explicit description of eigenfunctions of the Markov generator of one of the processes. For the Young graph our approach reconstructs the z-measures on partitions and the associated dynamics studied by Borodin and Olshanski [arXiv:math-ph/0409075, arXiv:0706.1034]. The generator of the dynamics of [arXiv:math-ph/0409075] is diagonal in the basis of the Meixner symmetric functions introduced recently by Olshanski [arXiv:1009.2037, arXiv:1103.5848]. We give new proofs to some of the results of these two papers. Other graphs to which our technique is applicable include the Pascal triangle, the Kingman graph (with the two-parameter Poisson-Dirichlet measures), the Schur graph and the general Young graph with Jack edge multiplicities.

preprint2011arXiv

The catalogue of positions of optically bright extragalactic radio sources OBRS-1

It is expected that the European Space Agency mission Gaia will make possible to determine coordinates in the optical domain of more than 500000 quasars. In 2006, a radio astrometry project was launched with the overall goal to make comparison of coordinate systems derived from future space-born astrometry instruments with the coordinate system constructed from analysis of the global very long baseline interferometry (VLBI) more robust. Investigation of their rotation, zonal errors, and the non-alignment of the radio and optical positions caused by both radio and optical structures are important for validation of both techniques. In order to support these studies, the densification of the list of compact extragalactic objects that are bright in both radio and optical ranges is desirable. A set of 105 objects from the list of 398 compact extragalactic radio sources with declination > -10 deg was observed with the VLBA+EVN with the primary goal of producing their images with milliarcsecond resolution. These sources are brighter than 18 magnitude at V band, and they were previously detected at the European VLBI network. In this paper coordinates of observed sources have been derived with milliarcsecond accuracies from analysis of these VLBI observations following the method of absolute astrometry. The catalogue of positions of 105 target sources is presented. The accuracies of sources coordinates are in the range of 0.3 to 7 mas, with the median 1.1 mas.

preprint2011arXiv

The EVN Galactic Plane Survey - EGaPS

I present a catalogue of positions and correlated flux densities of 109 compact extragalactic radio sources in the Galactic plane determined from analysis of a 48 hour VLBI experiment at 22 GHz with the European VLBI Network. The median position uncertainty is 9 mas. The correlated flux densities of detected sources are in the range of 20 to 300 mJy. In addition to target sources, nine water masers have been detected, two of them new. I derived position of masers with accuracies 30 to 200 mas and determined velocities of maser components and their correlated flux densities. The catalogue and supporting material is available at http://astrogeo.org/egaps.

preprint2011arXiv

The KCAL VERA 22 GHz calibrator survey

We observed at 22 GHz with the VLBI array VERA a sample of 1536 sources with correlated flux densities brighter than 200 mJy at 8 GHz. One half of target sources has been detected. The detection limit was around 200 mJy. We derived the correlated flux densities of 877 detected sources in three ranges of projected baseline lengths. The objective of these observations was to determine the suitability of given sources as phase calibrators for dual-beam and phase-referencing observations at high frequencies. Preliminary results indicate that the number of compact extragalactic sources at 22 GHz brighter than a given correlated flux density level is twice less than at 8 GHz.

preprint2011arXiv

The LBA Calibrator Survey of southern compact extragalactic radio sources - LCS1

We present a catalogue of positions and correlated flux densities of 410 flat-spectrum, compact extragalactic radio sources, previously detected in the AT20G survey. The catalogue spans the declination range -90deg, -40deg and was constructed from four 24 hour VLBI observing sessions with the Australian Long Baseline Array made at 8.3 GHz. The detection rate in these experiments is 97%. The median uncertainty of source positions is 2.6 mas, the median correlated flux density at baseline projections lengths longer than 1000 km is 0.14 Jy. The goal of this work is 1) to provide a pool of sources with positions known at the milliarcsecond level of accuracy that are needed for phase referencing observations, for geodetic VLBI, and for space navigation; 2) to extend the complete flux-limited sample of compact extragalactic sources to the southern hemisphere; and 3) to investigate parsec-scale properties of high-frequency selected sources from the AT20G survey. As a result of the campaign, the number of compact radio sources with declinations < -40deg detectable with VLBI with measured correlated flux densities and positions known with the milliarcsec level of accuracies increased by a factor of 3.5. The catalogue and supporting material is available at http://astrogeo.org/lcs1 .

preprint2010arXiv

Random Strict Partitions and Determinantal Point Processes

In this note we present new examples of determinantal point processes with infinitely many particles. The particles live on the half-lattice {1,2,...} or on the open half-line (0,+\infty). The main result is the computation of the correlation kernels. They have integrable form and are expressed through the Euler gamma function (the lattice case) and the classical Whittaker functions (the continuous case). Our processes are obtained via a limit transition from a model of random strict partitions introduced by Borodin (1997) in connection with the problem of harmonic analysis for projective characters of the infinite symmetric group.

preprint2009arXiv

Precise geodesy with the Very Long Baseline Array

We report on a program of geodetic measurements between 1994 and 2007 which used the Very Long Baseline Array and up to 10 globally distributed antennas. One of the goals of this program was to monitor positions of the array at a 1 millimeter level of accuracy and to tie the VLBA into the International Terrestrial Reference Frame. We describe the analysis of these data and report several interesting geophysical results including measured station displacements due to crustal motion, earthquakes, and antenna tilt. In terms of both formal errors and observed scatter, these sessions are among the very best geodetic VLBI experiments.

preprint2009arXiv

Random Walks on Strict Partitions

We consider a certain sequence of random walks. The state space of the n-th random walk is the set of all strict partitions of n (that is, partitions without equal parts). We prove that, as n goes to infinity, these random walks converge to a continuous-time Markov process. The state space of this process is the infinite-dimensional simplex consisting of all nonincreasing infinite sequences of nonnegative numbers with sum less than or equal to one. The main result about the limit process is the expression of its the pre-generator as a formal second order differential operator in a polynomial algebra. Of separate interest is the generalization of Kerov interlacing coordinates to the case of shifted Young diagrams.

preprint2008arXiv

A Two-Parameter Family of Infinite-Dimensional Diffusions in the Kingman Simplex

The aim of the paper is to introduce a two-parameter family of infinite-dimensional diffusion processes X(alpha,theta) related to Pitman's two-parameter Poisson-Dirichlet distributions PD(alpha,theta). The diffusions X(alpha,theta) are obtained in a scaling limit transition from certain finite Markov chains on partitions of natural numbers. The state space of X(alpha,theta) is an infinite-dimensional simplex called the Kingman simplex. In the special case when parameter alpha vanishes, our finite Markov chains are similar to Moran-type model in population genetics, and our diffusion processes reduce to the infinitely-many-neutral-alleles diffusion model studied by Ethier and Kurtz (1981). Our main results extend those of Ethier and Kurtz to the two-parameter case and are as follows: The Poisson-Dirichlet distribution PD(alpha,theta) is a unique stationary distribution for the corresponding process X(alpha,theta); the process is ergodic and reversible; the spectrum of its generator is explicitly described. The general two-parameter case seems to fall outside the setting of models of population genetics, and our approach differs in some aspects from that of Ethier and Kurtz. We also consider the case of degenerate series of parameters alpha and theta and conclude that the diffusions in finite-dimensional simplexes studied by Ethier and Kurtz (1981) arise as a special case of our two-parameter family of diffusions.

preprint2007arXiv

On observability of the free core nutation

Neither astronomical technique, including VLBI, can measure nutation directly. Estimates of parameters of the nutation model are produced by solving the LSQ problem of adjusting millions parameters using estimates of group delay. The choice of the mathematical model for nutation used in the estimation process of analysis of group delays affects our ability to interpret the results. Ignoring these subtleties and using parameters of the nutation model either in the form of time series, or in the form of empirical expansion as "VLBI measurement of nutation", opens a room for misinterpretation and mistakes. Detailed analysis of the problem reveals that the separation of forced nutations, atmospheric nutations, ocean nutations, and the retrograde free core nutation requires invoking some hypotheses, and beyond a specific level becomes uncertain. This sets a limit of our ability to make an inference about the free core nutation.