Showing 1–50 of 52 results for author: Román, P

Generative AI in Cybersecurity

Authors: Shivani Metta, Isaac Chang, Jack Parker, Michael P. Roman, Arturo F. Ehuan

Abstract: The dawn of Generative Artificial Intelligence (GAI), characterized by advanced models such as Generative Pre-trained Transformers (GPT) and other Large Language Models (LLMs), has been pivotal in reshaping the field of data analysis, pattern recognition, and decision-making processes. This surge in GAI technology has ushered in not only innovative opportunities for data processing and automation… ▽ More The dawn of Generative Artificial Intelligence (GAI), characterized by advanced models such as Generative Pre-trained Transformers (GPT) and other Large Language Models (LLMs), has been pivotal in reshaping the field of data analysis, pattern recognition, and decision-making processes. This surge in GAI technology has ushered in not only innovative opportunities for data processing and automation but has also introduced significant cybersecurity challenges. As GAI rapidly progresses, it outstrips the current pace of cybersecurity protocols and regulatory frameworks, leading to a paradox wherein the same innovations meant to safeguard digital infrastructures also enhance the arsenal available to cyber criminals. These adversaries, adept at swiftly integrating and exploiting emerging technologies, may utilize GAI to develop malware that is both more covert and adaptable, thus complicating traditional cybersecurity efforts. The acceleration of GAI presents an ambiguous frontier for cybersecurity experts, offering potent tools for threat detection and response, while concurrently providing cyber attackers with the means to engineer more intricate and potent malware. Through the joint efforts of Duke Pratt School of Engineering, Coalfire, and Safebreach, this research undertakes a meticulous analysis of how malicious agents are exploiting GAI to augment their attack strategies, emphasizing a critical issue for the integrity of future cybersecurity initiatives. The study highlights the critical need for organizations to proactively identify and develop more complex defensive strategies to counter the sophisticated employment of GAI in malware creation. △ Less

Submitted 2 May, 2024; originally announced May 2024.

TMM$-$Sim: A Versatile Tool for Optical Simulation of Thin$-$Film Solar Cells

Authors: Leandro Benatto, Omar Mesquita, Kaike R. M. Pachecoand Lucimara S. Roman, Marlus Koehler, Rodrigo B. Capaz, Graziâni Candiotto

Abstract: The Transfer Matrix Method (TMM) has become a prominent tool for the optical simulation of thin$-$film solar cells, particularly among researchers specializing in organic semiconductors and perovskite materials. As the commercial viability of these solar cells continues to advance, driven by rapid developments in materials and production processes, the importance of optical simulation has grown si… ▽ More The Transfer Matrix Method (TMM) has become a prominent tool for the optical simulation of thin$-$film solar cells, particularly among researchers specializing in organic semiconductors and perovskite materials. As the commercial viability of these solar cells continues to advance, driven by rapid developments in materials and production processes, the importance of optical simulation has grown significantly. By leveraging optical simulation, researchers can gain profound insights into photovoltaic phenomena, empowering the implementation of device optimization strategies to achieve enhanced performance. However, existing TMM$-$based packages exhibit limitations, such as requiring programming expertise, licensing fees, or lack of support for bilayer device simulation. In response to these gaps and challenges, we present the TMM Simulator (TMM$-$Sim), an intuitive and user$-$friendly tool to calculate essential photovoltaic parameters, including the optical electric field profile, exciton generation profile, fraction of light absorbed per layer, photocurrent, external quantum efficiency, internal quantum efficiency, and parasitic losses. An additional advantage of TMM$-$Sim lies in its capacity to generate outcomes suitable as input parameters for electro$-$optical device simulations. In this work, we offer a comprehensive guide, outlining a step$-$by$-$step process to use TMM$-$Sim, and provide a thorough analysis of the results. TMM$-$Sim is freely available, accessible through our web server (nanocalc.org), or downloadable from the TMM$-$Sim repository (for \textit{Unix}, \textit{Windows}, and \textit{macOS}) on \textit{GitHub}. With its user$-$friendly interface and powerful capabilities, TMM$-$Sim aims to facilitate and accelerate research in thin$-$film solar cells, fostering advancements in renewable energy technologies. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Comments: 9 pages, 8 figures

Report number: COMPHY-D-23-00846

Quantum simulation of entanglement dynamics in a quantum processor

Authors: C. Inzulza, S. Saavedra-Pino, F. Albarrán-Arriagada, P. Roman, J. C. Retamal

Abstract: We implement a five-qubit protocol in IBM quantum processors to get entanglement dynamics in a two qubit system in the presence of an environment. Specifically, two qubits represent the main system, another two qubits the environment, and an additional qubit is used as an auxiliary qubit to perform the quantum entanglement estimation. We focus on measuring, in this superconducting quantum processo… ▽ More We implement a five-qubit protocol in IBM quantum processors to get entanglement dynamics in a two qubit system in the presence of an environment. Specifically, two qubits represent the main system, another two qubits the environment, and an additional qubit is used as an auxiliary qubit to perform the quantum entanglement estimation. We focus on measuring, in this superconducting quantum processor, the sudden death and sudden birth of entanglement. We obtain the quantum entanglement evolution of the main system qubits and the environment qubits as the average of $N=10$ independent experiments in the same quantum device, observing that the noisy nature of current quantum processors produce a shift on times signaling sudden death o sudden birth of entanglement. This work takes relevance showing the usefulness of current noisy quantum devices to test fundamental concepts in quantum information. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: 6 pages, and 12 figures

Matrix exceptional Laguerre polynomials

Authors: Erik Koelink, Lucía Morey, Pablo Román

Abstract: We give an analog of exceptional polynomials in the matrix valued setting by considering suitable factorizations of a given second order differential operator and performing Darboux transformations. Orthogonality and density of the exceptional sequence is discussed in detail. We give an example of matrix valued exceptional Laguerre polynomials of arbitrary size. We give an analog of exceptional polynomials in the matrix valued setting by considering suitable factorizations of a given second order differential operator and performing Darboux transformations. Orthogonality and density of the exceptional sequence is discussed in detail. We give an example of matrix valued exceptional Laguerre polynomials of arbitrary size. △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: 28 pages, 1 figure

MSC Class: 33C45; 33C47

Chop Chop: Byzantine Atomic Broadcast to the Network Limit

Authors: Martina Camaioni, Rachid Guerraoui, Matteo Monti, Pierre-Louis Roman, Manuel Vidigueira, Gauthier Voron

Abstract: At the heart of state machine replication, the celebrated technique enabling decentralized and secure universal computation, lies Atomic Broadcast, a fundamental communication primitive that orders, authenticates, and deduplicates messages. This paper presents Chop Chop, a Byzantine Atomic Broadcast system that amortizes the cost of ordering, authenticating and deduplicating messages, achieving "l… ▽ More At the heart of state machine replication, the celebrated technique enabling decentralized and secure universal computation, lies Atomic Broadcast, a fundamental communication primitive that orders, authenticates, and deduplicates messages. This paper presents Chop Chop, a Byzantine Atomic Broadcast system that amortizes the cost of ordering, authenticating and deduplicating messages, achieving "line rate" (i.e., closely matching the complexity of a protocol that does not ensure any ordering, authentication or Byzantine resilience) even when processing messages as small as 8 bytes. Chop Chop attains this performance by means of a new form of batching we call distillation. A distilled batch is a set of messages that are fast to authenticate and deduplicate, as well as order. Batches are distilled using a novel interactive mechanism involving brokers, an untrusted layer of facilitating processes between clients and servers. In a geo-distributed deployment of 64 medium-sized servers, with clients situated cross-cloud, Chop Chop processes 43,600,000 messages per second with an average latency of 3.6 seconds. Under the same conditions, state-of-the-art alternatives offer two orders of magnitude less throughput for the same latency. We showcase three simple Chop Chop applications: a Payment system, an Auction house and a "Pixel war" game, respectively achieving 32, 2.3 and 35 million operations per second. △ Less

Submitted 14 April, 2023; originally announced April 2023.

Lie algebras of differential operators for Matrix valued Laguerre type polynomials

Authors: Andrea L. Gallo, Pablo Román

Abstract: We study algebras of differential and difference operators acting on matrix valued orthogonal polynomials (MVOPs) with respect to a weight matrix of the form $W^{(νにゅー)}_φふぁい(x) = x^νにゅーe^{-φふぁい(x)} W^{(νにゅー)}_{pol}(x)$, where $νにゅー>0$, $W^{(νにゅー)}_{pol}(x)$ is certain matrix valued polynomial and $φふぁい$ an entire function. We introduce a pair differential operators $\mathcal{D}$, $\mathcal{D}^{\dagger}$ which are mutual… ▽ More We study algebras of differential and difference operators acting on matrix valued orthogonal polynomials (MVOPs) with respect to a weight matrix of the form $W^{(νにゅー)}_φふぁい(x) = x^νにゅーe^{-φふぁい(x)} W^{(νにゅー)}_{pol}(x)$, where $νにゅー>0$, $W^{(νにゅー)}_{pol}(x)$ is certain matrix valued polynomial and $φふぁい$ an entire function. We introduce a pair differential operators $\mathcal{D}$, $\mathcal{D}^{\dagger}$ which are mutually adjoint with respect to the matrix inner product induced by $W^{(νにゅー)}_φふぁい(x)$. We prove that the Lie algebra generated by $\mathcal{D}$ and $\mathcal{D}^{\dagger}$ is finite dimensional if and only if $φふぁい$ is a polynomial, giving a partial answer to a problem by M. Ismail. In the case $φふぁい$ polynomial, we describe the structure of this Lie algebra. The case $φふぁい(x)=x$, is discussed in detail. We derive difference and differential relations for the MVOPs. We give explicit expressions for the entries of the MVOPs in terms of classical Laguerre and Dual Hahn polynomials. △ Less

Submitted 12 March, 2023; originally announced March 2023.

Non Abelian Toda-type equations and matrix valued orthogonal polynomials

Authors: Alfredo Deaño, Lucía Morey, Pablo Román

Abstract: In this paper, we study parameter deformations of matrix valued orthogonal polynomials (MVOPs). These deformations are built on the use of certain matrix valued operators which are symmetric with respect to the matrix valued inner product defined by the orthogonality weight. We show that the recurrence coefficients associated with these operators satisfy generalizations of the non-Abelian lattice… ▽ More In this paper, we study parameter deformations of matrix valued orthogonal polynomials (MVOPs). These deformations are built on the use of certain matrix valued operators which are symmetric with respect to the matrix valued inner product defined by the orthogonality weight. We show that the recurrence coefficients associated with these operators satisfy generalizations of the non-Abelian lattice equations. We provide a Lax pair formulation for these equations, and an example of deformed Hermite-type matrix valued polynomials is discussed in detail. △ Less

Submitted 17 March, 2024; v1 submitted 28 February, 2023; originally announced February 2023.

Comments: 20 pages, small typos corrected

MSC Class: 37K10; 33C47; 33C45

SurferMonkey: A Decentralized Anonymous Blockchain Intercommunication System via Zero Knowledge Proofs

Authors: Miguel Díaz Montiel, Rachid Guerraoui, Pierre-Louis Roman

Abstract: Blockchain intercommunication systems enable the exchanges of messages between blockchains. This interoperability promotes innovation, unlocks liquidity and access to assets. However, blockchains are isolated systems that originally were not designed for interoperability. This makes cross-chain communication, or bridges for short, insecure by nature. More precisely, cross-chain systems face securi… ▽ More Blockchain intercommunication systems enable the exchanges of messages between blockchains. This interoperability promotes innovation, unlocks liquidity and access to assets. However, blockchains are isolated systems that originally were not designed for interoperability. This makes cross-chain communication, or bridges for short, insecure by nature. More precisely, cross-chain systems face security challenges in terms of selfish rational players such as maximal extractable value (MEV) and censorship. We propose to solve these challenges using zero knowledge proofs (ZKPs) for cross-chain communication. Securing cross-chain communication is remarkably more complex than securing single-chain events as such a system must preserve user security against both on- and off-chain analysis. To achieve this goal, we propose the following pair of contributions: the DACT protocol and the SurferMonkey infrastructure that supports the DACT protocol. The DACT protocol is a global solution for the anonymity and security challenges of agnostic blockchain intercommunication. DACT breaks on- and off-chain analysis thanks to the use of ZKPs. SurferMonkey is a decentralized infrastructure that implements DACT in practice. Since SurferMonkey works at the blockchain application layer, any decentralized application (dApp) can use SurferMonkey to send any type of message to a dApp on another blockchain. With SurferMonkey, users can neither be censored nor be exposed to MEV. By applying decentralized proactive security, we obtain resilience against selfish rational players, and raise the security bar against cyberattacks. We have implemented a proof of concept (PoC) of SurferMonkey by reverse engineering Tornado Cash and by applying IDEN3 ZKP circuits. SurferMonkey enables new usecases, ranging from anonymous voting and gaming, to a new phase of anonymous decentralized finance (aDeFi). △ Less

Submitted 24 October, 2022; originally announced October 2022.

Design and implementation of a device based on an off-axis parabolic mirror to perform luminescence experiments in a scanning tunneling microscope

Authors: Ricardo Javier Peña Román, Yves Auad, Lucas Grasso, Lazaro A Padilha, Fernando Alvarez, Ingrid David Barcelos, Mathieu Kociak, Luiz Fernando Zagonel

Abstract: We present the design, implementation, and illustrative results of a light collection/injection strategy based on an off-axis parabolic mirror collector for a low-temperature Scanning Tunneling Microscope (STM). This device allows us to perform STM induced Light Emission (STM-LE) and Cathodoluminescence (STM-CL) experiments and in situ Photoluminescence (PL) and Raman spectroscopy as complementary… ▽ More We present the design, implementation, and illustrative results of a light collection/injection strategy based on an off-axis parabolic mirror collector for a low-temperature Scanning Tunneling Microscope (STM). This device allows us to perform STM induced Light Emission (STM-LE) and Cathodoluminescence (STM-CL) experiments and in situ Photoluminescence (PL) and Raman spectroscopy as complementary techniques. Considering the Étendue conservation and using an off-axis parabolic mirror, it is possible to design a light collection and injection system that displays 72% of collection efficiency (considering the hemisphere above the sample surface) while maintaining high spectral resolution and minimizing signal loss. The performance of the STM is tested by atomically resolved images and scanning tunneling spectroscopy results on standard sample surfaces. The capabilities of our system are demonstrated by performing STM-LE on metallic surfaces and two-dimensional semiconducting samples, observing both plasmonic and excitonic emissions. In addition, we carried out in situ PL measurements on semiconducting monolayers and quantum dots and in situ Raman on graphite and hexagonal boron nitride (h-BN) samples. Additionally, STM-CL and PL were obtained on monolayer h-BN gathering luminescence spectra that are typically associated with intragap states related to carbon defects. The results show that the flexible and efficient light injection and collection device based on an off-axis parabolic mirror is a powerful tool to study several types of nanostructures with multiple spectroscopic techniques in correlation with their morphology at the atomic scale and electronic structure. △ Less

Submitted 10 October, 2022; originally announced October 2022.

Comments: 19 pages, 14 Figures

Journal ref: Review of Scientific Instruments 93, 043704 (2022);

Tunneling-current-induced local excitonic luminescence in p-doped WSe$_2$ monolayers

Authors: Ricardo Javier Peña Román, Yves Auad, Lucas Grasso, Fernando Alvarez, Ingrid David Barcelos, Luiz Fernando Zagonel

Abstract: We have studied the excitonic properties of exfoliated tungsten diselenide (WSe$_2$) monolayers transferred to gold substrates using the tunneling current in a Scanning Tunneling Microscope (STM) operated in air to excite the light emission locally. In obtained spectra, emission energies are independent of the applied bias voltage and resemble photoluminescence (PL) results, indicating that, in bo… ▽ More We have studied the excitonic properties of exfoliated tungsten diselenide (WSe$_2$) monolayers transferred to gold substrates using the tunneling current in a Scanning Tunneling Microscope (STM) operated in air to excite the light emission locally. In obtained spectra, emission energies are independent of the applied bias voltage and resemble photoluminescence (PL) results, indicating that, in both cases, the light emission is due to neutral and charged exciton recombination. Interestingly, the electron injection rate, that is, the tunneling current, can be used to control the ratio of charged to neutral exciton emission. The obtained quantum yield in the transition metal dichalcogenide (TMD) is $~5x10^{-7}$ photons per electron. The proposed excitation mechanism is the direct injection of carriers into the conduction band. The monolayer WSe$_2$ presents bright and dark defects spotted by STM images performed under UHV. STS confirms the sample as p-doped, possibly as a net result of the observed defects. The presence of an interfacial water layer decouples the monolayer from the gold support and allows excitonic emission from the WSe$_2$ monolayer. The creation of a water layer is an inherent feature of the sample transferring process due to the ubiquitous air moisture. Consequently, vacuum thermal annealing, which removes the water layer, quenches excitonic luminescence from the TMD. The tunneling current can locally displace water molecules leading to excitonic emission quenching and to plasmonic emission due to the gold substrate. The present findings extend the use and the understanding of STM induced light emission (STM-LE) on semiconducting TMDs to probe exciton emission and dynamics with high spatial resolution. △ Less

Submitted 4 October, 2022; originally announced October 2022.

Comments: Includes main text and supplementary material, 22 pages

Journal ref: Nanoscale, 2020,12, 13460-13470

Asymptotics of matrix valued orthogonal polynomials on $[-1,1]$

Authors: Alfredo Deaño, Arno B. J. Kuijlaars, Pablo Román

Abstract: We analyze the large degree asymptotic behavior of matrix valued orthogonal polynomials (MVOPs), with a weight that consists of a Jacobi scalar factor and a matrix part. Using the Riemann-Hilbert formulation for MVOPs and the Deift-Zhou method of steepest descent, we obtain asymptotic expansions for the MVOPs as the degree tends to infinity, in different regions of the complex plane (outside the i… ▽ More We analyze the large degree asymptotic behavior of matrix valued orthogonal polynomials (MVOPs), with a weight that consists of a Jacobi scalar factor and a matrix part. Using the Riemann-Hilbert formulation for MVOPs and the Deift-Zhou method of steepest descent, we obtain asymptotic expansions for the MVOPs as the degree tends to infinity, in different regions of the complex plane (outside the interval of orthogonality, on the interval away from the endpoints and in neighborhoods of the endpoints), as well as for the matrix coefficients in the three-term recurrence relation for these MVOPs. The asymptotic analysis follows the work of Kuijlaars, McLaughlin, Van Assche and Vanlessen on scalar Jacobi-type orthogonal polynomials, but it also requires several different factorizations of the matrix part of the weight, in terms of eigenvalues/eigenvectors and using a matrix Szegő function. We illustrate the results with two main examples, MVOPs of Jacobi and Gegenbauer type, coming from group theory. △ Less

Submitted 10 April, 2023; v1 submitted 3 October, 2022; originally announced October 2022.

Comments: 48 pages, 8 figures. Some typos corrected in the last version

Duality and difference operators for matrix valued discrete polynomials on the nonnegative integers

Authors: Bruno Eijsvoogel, Lucía Morey, Pablo Román

Abstract: In this paper we introduce a notion of duality for matrix valued orthogonal polynomials with respect to a measure supported on the nonnegative integers. We show that the dual families are closely related to certain difference operators acting on the matrix orthogonal polynomials. These operators belong to the so called Fourier algebras, which play a key role in the construction of the families.… ▽ More In this paper we introduce a notion of duality for matrix valued orthogonal polynomials with respect to a measure supported on the nonnegative integers. We show that the dual families are closely related to certain difference operators acting on the matrix orthogonal polynomials. These operators belong to the so called Fourier algebras, which play a key role in the construction of the families. In order to illustrate duality, we describe a family of Charlier type matrix orthogonal polynomials with explicit shift operators which allow us to find explicit formulas for three term recurrences, difference operators and square norms. These are the essential ingredients for the construction of different dual families. △ Less

Submitted 25 October, 2021; originally announced October 2021.

Comments: 54 pages, 1 figure

MSC Class: 42C05; 33C45; 33C47; 33E30

Band gap measurements of monolayer h-BN and insights into carbon-related point defects

Authors: Ricardo Javier Peña Román, Fábio J R Costa Costa, Alberto Zobelli, Christine Elias, Pierre Valvin, Guillaume Cassabois, Bernard Gil, Alex Summerfield, Tin S Cheng, Christopher J Mellor, Peter H Beton, Sergei V Novikov, Luiz F Zagonel

Abstract: Being a flexible wide band gap semiconductor, hexagonal boron nitride (h-BN) has great potential for technological applications like efficient deep ultraviolet light sources, building block for two-dimensional heterostructures and room temperature single photon emitters in the ultraviolet and visible spectral range. To enable such applications, it is mandatory to reach a better understanding of th… ▽ More Being a flexible wide band gap semiconductor, hexagonal boron nitride (h-BN) has great potential for technological applications like efficient deep ultraviolet light sources, building block for two-dimensional heterostructures and room temperature single photon emitters in the ultraviolet and visible spectral range. To enable such applications, it is mandatory to reach a better understanding of the electronic and optical properties of h-BN and the impact of various structural defects. Despite the large efforts in the last years, aspects such as the electronic band gap value, the exciton binding energy and the effect of point defects remained elusive, particularly when considering a single monolayer. Here, we directly measured the density of states of a single monolayer of h-BN epitaxially grown on highly oriented pyrolytic graphite, by performing low temperature scanning tunneling microscopy (STM) and spectroscopy (STS). The observed h-BN electronic band gap on defect-free regions is $(6.8\pm0.2)$ eV. Using optical spectroscopy to obtain the h-BN optical band gap, the exciton binding energy is determined as being of $(0.7\pm0.2)$ eV. In addition, the locally excited cathodoluminescence and photoluminescence show complex spectra that are typically associated to intragap states related to carbon defects. Moreover, in some regions of the monolayer h-BN we identify, using STM, point defects which have intragap electronic levels around 2.0 eV below the Fermi level. △ Less

Submitted 16 July, 2021; originally announced July 2021.

Comments: 50 Pages, 8 Figures, 100+ references

Journal ref: 2D Material 8 044001 (2021)

Application of the Adams-Bashfort-Mowlton Method to the Numerical Study of Linear Fractional Oscillators Models

Authors: Parovik Roman

Abstract: The paper presents a numerical analysis of the class of mathematical models of linear fractional oscillators, which is the Cauchy problem for a differential equation with derivatives of fractional orders in the sense of Gerasimov-Caputo. A method based on an explicit nonlocal finite-difference scheme (ENFDS) and the Adams-Bashfort-Moulton (ABM) method is considered a numerical analysis tool. An an… ▽ More The paper presents a numerical analysis of the class of mathematical models of linear fractional oscillators, which is the Cauchy problem for a differential equation with derivatives of fractional orders in the sense of Gerasimov-Caputo. A method based on an explicit nonlocal finite-difference scheme (ENFDS) and the Adams-Bashfort-Moulton (ABM) method is considered a numerical analysis tool. An analysis of the errors of the methods is carried out, and it is shown that the ABM method is more accurate and converges faster to an exact solution than the ENFDS method. △ Less

Submitted 15 May, 2021; originally announced May 2021.

MSC Class: 65L05 (Primary); 34A08 (Secondary) ACM Class: G.1.7

Let Me Know What to Ask: Interrogative-Word-Aware Question Generation

Authors: Junmo Kang, Haritz Puerto San Roman, Sung-Hyon Myaeng

Abstract: Question Generation (QG) is a Natural Language Processing (NLP) task that aids advances in Question Answering (QA) and conversational assistants. Existing models focus on generating a question based on a text and possibly the answer to the generated question. They need to determine the type of interrogative word to be generated while having to pay attention to the grammar and vocabulary of the que… ▽ More Question Generation (QG) is a Natural Language Processing (NLP) task that aids advances in Question Answering (QA) and conversational assistants. Existing models focus on generating a question based on a text and possibly the answer to the generated question. They need to determine the type of interrogative word to be generated while having to pay attention to the grammar and vocabulary of the question. In this work, we propose Interrogative-Word-Aware Question Generation (IWAQG), a pipelined system composed of two modules: an interrogative word classifier and a QG model. The first module predicts the interrogative word that is provided to the second module to create the question. Owing to an increased recall of deciding the interrogative words to be used for the generated questions, the proposed model achieves new state-of-the-art results on the task of QG in SQuAD, improving from 46.58 to 47.69 in BLEU-1, 17.55 to 18.53 in BLEU-4, 21.24 to 22.33 in METEOR, and from 44.53 to 46.94 in ROUGE-L. △ Less

Submitted 30 October, 2019; originally announced October 2019.

Comments: Accepted at 2nd Workshop on Machine Reading for Question Answering (MRQA), EMNLP 2019

Ladder relations for a class of matrix valued orthogonal polynomials

Authors: Alfredo Deaño, Bruno Eijsvoogel, Pablo Román

Abstract: Using the theory introduced by Casper and Yakimov, we investigate the structure of algebras of differential and difference operators acting on matrix valued orthogonal polynomials (MVOPs) on $\mathbb{R}$, and we derive algebraic and differential relations for these MVOPs. A particular case of importance is that of MVOPs with respect to a matrix weight of the form… ▽ More Using the theory introduced by Casper and Yakimov, we investigate the structure of algebras of differential and difference operators acting on matrix valued orthogonal polynomials (MVOPs) on $\mathbb{R}$, and we derive algebraic and differential relations for these MVOPs. A particular case of importance is that of MVOPs with respect to a matrix weight of the form $W(x)=e^{-v(x)}e^{xA} e^{xA^\ast}$ on the real line, where $v$ is a scalar polynomial of even degree with positive leading coefficient and $A$ is a constant matrix. △ Less

Submitted 16 November, 2020; v1 submitted 17 July, 2019; originally announced July 2019.

Matrix valued Hermite polynomials, Burchnall formulas and non-abelian Toda lattice

Authors: Mourad E. H. Ismail, Erik Koelink, Pablo Román

Abstract: A general family of matrix valued Hermite type orthogonal polynomials is introduced and studied in detail by deriving Pearson equations for the weight and matrix valued differential equations for these matrix polynomials. This is used to derive Rodrigues formulas, explicit formulas for the squared norm and to give an explicit expression of the matrix entries as well to derive a connection formula… ▽ More A general family of matrix valued Hermite type orthogonal polynomials is introduced and studied in detail by deriving Pearson equations for the weight and matrix valued differential equations for these matrix polynomials. This is used to derive Rodrigues formulas, explicit formulas for the squared norm and to give an explicit expression of the matrix entries as well to derive a connection formula for the matrix polynomials of Hermite type. We derive matrix valued analogues of Burchnall formulas in operational form as well explicit expansions for the matrix valued Hermite type orthogonal polynomials as well as for previously introduced matrix valued Gegenbauer type orthogonal polynomials. The Burchnall approach gives two descriptions of the matrix valued orthogonal polynomials for the Toda modification of the matrix weight for the Hermite setting. In particular, we obtain a non-trivial solution to the non-abelian Toda lattice equations. △ Less

Submitted 17 November, 2018; originally announced November 2018.

Comments: 28 pages

Journal ref: Advances in Applied Mathematics 110 (2019), 235-269

Matrix valued Laguerre polynomials

Authors: Erik Koelink, Pablo Román

Abstract: Matrix valued Laguerre polynomials are introduced via a matrix weight function involving several degrees of freedom using the matrix nature. Under suitable conditions on the parameters the matrix weight function satisfies matrix Pearson equations, which allow to introduce shift operators for these polynomials. The shift operators lead to explicit expressions for the structures of these matrix valu… ▽ More Matrix valued Laguerre polynomials are introduced via a matrix weight function involving several degrees of freedom using the matrix nature. Under suitable conditions on the parameters the matrix weight function satisfies matrix Pearson equations, which allow to introduce shift operators for these polynomials. The shift operators lead to explicit expressions for the structures of these matrix valued Laguerre polynomials, such as a Rodrigues formula, the coefficients in the three-term recurrence, differential operators, and expansion formulas. △ Less

Submitted 15 November, 2018; originally announced November 2018.

Comments: 20 pages, to appear in Positivity and Noncommutative Analysis Festschrift in Honour of Ben de Pagter (eds. G. Buskes, M. de Jeu, P. Dodds, A. Schep, F. Sukochev, J. van Neerven and A. Wickstead)

Journal ref: pp. 295-320 in "Positivity and Noncommutative Analysis Festschrift in Honour of Ben de Pagter'' (eds. G. Buskes, M. de Jeu, P. Dodds, A. Schep, F. Sukochev, J. van Neerven and A. Wickstead), Birkhäuser, 2019

An inner warp in the DoAr 44 T Tauri transition disk

Authors: Simon Casassus, Henning Avenhaus, Sebastian Perez, Victor Navarro, Miguel Carcamo, Sebastian Marino, Lucas Cieza, Sascha P. Quanz, Felipe Alarcon, Alice Zurlo, Axel Osses, Fernando R. Rannou, Pablo E. Romaan, Marcelo Barraza

Abstract: Optical/IR images of transition disks (TDs) have revealed deep intensity decrements in the rings of HAeBes HD142527 and HD100453, that can be interpreted as shadowing from sharply tilted inner disks, such that the outer disks are directly exposed to stellar light. Here we report similar dips in SPHERE+IRDIS differential polarized imaging (DPI) of TTauri DoAr44. With a fairly axially symmetric ring… ▽ More Optical/IR images of transition disks (TDs) have revealed deep intensity decrements in the rings of HAeBes HD142527 and HD100453, that can be interpreted as shadowing from sharply tilted inner disks, such that the outer disks are directly exposed to stellar light. Here we report similar dips in SPHERE+IRDIS differential polarized imaging (DPI) of TTauri DoAr44. With a fairly axially symmetric ring in the submm radio continuum, DoAr44 is likely also a warped system. We constrain the warp geometry by comparing radiative transfer predictions with the DPI data in H band (Q_φふぁい(H)) and with a re-processing of archival 336GHz ALMA observations. The observed DPI shadows have coincident radio counterparts, but the intensity drops are much deeper in Q_φふぁい(H) (~88%), compared to the shallow drops at 336GHz (~24%). Radiative transfer predictions with an inner disk tilt of ~30+-5deg approximately account for the observations. ALMA long-baseline observations should allow the observation of the warped gas kinematics inside the cavity of DoAr44. △ Less

Submitted 6 April, 2018; originally announced April 2018.

Comments: 11 pages, 10 figures, 3 tables, accepted for publication in MNRAS

Dietcoin: shortcutting the Bitcoin verification process for your smartphone

Authors: Davide Frey, Marc X. Makkes, Pierre-Louis Roman, François Taïani, Spyros Voulgaris

Abstract: Blockchains have a storage scalability issue. Their size is not bounded and they grow indefinitely as time passes. As of August 2017, the Bitcoin blockchain is about 120 GiB big while it was only 75 GiB in August 2016. To benefit from Bitcoin full security model, a bootstrapping node has to download and verify the entirety of the 120 GiB. This poses a challenge for low-resource devices such as sma… ▽ More Blockchains have a storage scalability issue. Their size is not bounded and they grow indefinitely as time passes. As of August 2017, the Bitcoin blockchain is about 120 GiB big while it was only 75 GiB in August 2016. To benefit from Bitcoin full security model, a bootstrapping node has to download and verify the entirety of the 120 GiB. This poses a challenge for low-resource devices such as smartphones. Thankfully, an alternative exists for such devices which consists of downloading and verifying just the header of each block. This partial block verification enables devices to reduce their bandwidth requirements from 120 GiB to 35 MiB. However, this drastic decrease comes with a safety cost implied by a partial block verification. In this work, we enable low-resource devices to fully verify subchains of blocks without having to pay the onerous price of a full chain download and verification; a few additional MiB of bandwidth suffice. To do so, we propose the design of diet nodes that can securely query full nodes for shards of the UTXO set, which is needed to perform full block verification and can otherwise only be built by sequentially parsing the chain. △ Less

Submitted 28 March, 2018; originally announced March 2018.

Comments: HAL link: https://hal.inria.fr/hal-01743995v1

Report number: RR-9162

Generalized Burchnall-Type Identities for Orthogonal Polynomials and Expansions

Authors: Mourad E. H. Ismail, Erik Koelink, Pablo Román

Abstract: Burchnall's method to invert the Feldheim-Watson linearization formula for the Hermite polynomials is extended to all polynomial families in the Askey-scheme and its $q$-analogue. The resulting expansion formulas are made explicit for several families corresponding to measures with infinite support, including the Wilson and Askey-Wilson polynomials. An integrated version gives the possibility to g… ▽ More Burchnall's method to invert the Feldheim-Watson linearization formula for the Hermite polynomials is extended to all polynomial families in the Askey-scheme and its $q$-analogue. The resulting expansion formulas are made explicit for several families corresponding to measures with infinite support, including the Wilson and Askey-Wilson polynomials. An integrated version gives the possibility to give alternate expression for orthogonal polynomials with respect to a modified weight. This gives expansions for polynomials, such as Hermite, Laguerre, Meixner, Charlier, Meixner-Pollaczek and big $q$-Jacobi polynomials and big $q$-Laguerre polynomials. We show that one can find expansions for the orthogonal polynomials corresponding to the Toda-modification of the weight for the classical polynomials that correspond to known explicit solutions for the Toda lattice, i.e., for Hermite, Laguerre, Charlier, Meixner, Meixner-Pollaczek and Krawtchouk polynomials. △ Less

Submitted 17 July, 2018; v1 submitted 26 February, 2018; originally announced February 2018.

Journal ref: SIGMA 14 (2018), 072, 24 pages

Matrix elements of irreducible representations of $\mathrm{SU}(n+1)\times\mathrm{SU}(n+1)$ and multivariable matrix-valued orthogonal polynomials

Authors: Erik Koelink, Maarten van Pruijssen, Pablo Román

Abstract: In Part 1 we study the spherical functions on compact symmetric pairs of arbitrary rank under a suitable multiplicity freeness assumption and additional conditions on the branching rules. The spherical functions are taking values in the spaces of linear operators of a finite dimensional representation of the subgroup, so the spherical functions are matrix-valued. Under these assumptions these func… ▽ More In Part 1 we study the spherical functions on compact symmetric pairs of arbitrary rank under a suitable multiplicity freeness assumption and additional conditions on the branching rules. The spherical functions are taking values in the spaces of linear operators of a finite dimensional representation of the subgroup, so the spherical functions are matrix-valued. Under these assumptions these functions can be described in terms of matrix-valued orthogonal polynomials in several variables, where the number of variables is the rank of the compact symmetric pair. Moreover, these polynomials are uniquely determined as simultaneous eigenfunctions of a commutative algebra of differential operators. In Part 2 we verify that the group case $\mathrm{SU}(n+1)$ meets all the conditions that we impose in Part 1. For any $k\in\mathbb{N}_{0}$ we obtain families of orthogonal polynomials in $n$ variables with values in the $N\times N$-matrices, where $N=\binom{n+k}{k}$. The case $k=0$ leads to the classical Heckman-Opdam polynomials of type $A_{n}$ with geometric parameter. For $k=1$ we obtain the most complete results. In this case we give an explicit expression of the matrix weight, which we show to be irreducible whenever $n\ge2$. We also give explicit expressions of the spherical functions that determine the matrix weight for $k=1$. These expressions are used to calculate the spherical functions that determine the matrix weight for general $k$ up to invertible upper-triangular matrices. This generalizes and gives a new proof of a formula originally obtained by Koornwinder for the case $n=1$. The commuting family of differential operators that have the matrix-valued polynomials as simultaneous eigenfunctions contains an element of order one. We give explicit formulas for differential operators of order one and two for $(n,k)$ equal to $(2,1)$ and $(3,1)$. △ Less

Submitted 6 June, 2017; originally announced June 2017.

Comments: 45 pages, 2 figures

Multi-GPU maximum entropy image synthesis for radio astronomy

Authors: M. Cárcamo, P. Román, S. Casassus, V. Moral, F. R. Rannou

Abstract: The maximum entropy method (MEM) is a well known deconvolution technique in radio-interferometry. This method solves a non-linear optimization problem with an entropy regularization term. Other heuristics such as CLEAN are faster but highly user dependent. Nevertheless, MEM has the following advantages: it is unsupervised, it has a statistical basis, it has a better resolution and better image qua… ▽ More The maximum entropy method (MEM) is a well known deconvolution technique in radio-interferometry. This method solves a non-linear optimization problem with an entropy regularization term. Other heuristics such as CLEAN are faster but highly user dependent. Nevertheless, MEM has the following advantages: it is unsupervised, it has a statistical basis, it has a better resolution and better image quality under certain conditions. This work presents a high performance GPU version of non-gridding MEM, which is tested using real and simulated data. We propose a single-GPU and a multi-GPU implementation for single and multi-spectral data, respectively. We also make use of the Peer-to-Peer and Unified Virtual Addressing features of newer GPUs which allows to exploit transparently and efficiently multiple GPUs. Several ALMA data sets are used to demonstrate the effectiveness in imaging and to evaluate GPU performance. The results show that a speedup from 1000 to 5000 times faster than a sequential version can be achieved, depending on data and image size. This allows to reconstruct the HD142527 CO(6-5) short baseline data set in 2.1 minutes, instead of 2.5 days that takes a sequential version on CPU. △ Less

Submitted 6 November, 2017; v1 submitted 8 March, 2017; originally announced March 2017.

Comments: 11 pages, 13 figures

Journal ref: M. Cárcamo, P.E. Román, S. Casassus, V. Moral, F.R. Rannou, Multi-GPU maximum entropy image synthesis for radio astronomy, Astronomy and Computing, Volume 22, 2018, Pages 16-27, ISSN 2213-1337

Spherical functions approach to sums of random Hermitian matrices

Authors: Arno B. J. Kuijlaars, Pablo Román

Abstract: We present an approach to sums of random Hermitian matrices via the theory of spherical functions for the Gelfand pair $(\mathrm{U}(n) \ltimes \mathrm{Herm}(n), \mathrm{U}(n))$. It is inspired by a similar approach of Kieburg and Kösters for products of random matrices. The spherical functions have determinantal expressions because of the Harish-Chandra/Itzykson-Zuber integral formula. It leads to… ▽ More We present an approach to sums of random Hermitian matrices via the theory of spherical functions for the Gelfand pair $(\mathrm{U}(n) \ltimes \mathrm{Herm}(n), \mathrm{U}(n))$. It is inspired by a similar approach of Kieburg and Kösters for products of random matrices. The spherical functions have determinantal expressions because of the Harish-Chandra/Itzykson-Zuber integral formula. It leads to remarkably simple expressions for the spherical transform and its inverse. The spherical transform is applied to sums of unitarily invariant random matrices from polynomial ensembles and the subclass of polynomial ensembles of derivative type (in the additive sense), which turns out to be closed under addition. We finally present additional detailed calculations for the sum with a random matrix from a Laguerre Unitary Ensemble. △ Less

Submitted 27 November, 2016; originally announced November 2016.

Comments: 25 pages

Journal ref: International Mathematics Research Notices 2019 (4) (2019), 1005-1029

Deformation of matrix-valued orthogonal polynomials related to Gelfand pairs

Authors: Maarten van Pruijssen, Pablo Román

Abstract: In this paper we present a method to obtain deformations of families of matrix-valued orthogonal polynomials that are associated to the representation theory of compact Gelfand pairs. These polynomials have the Sturm-Liouville property in the sense that they are simultaneous eigenfunctions of a symmetric second order differential operator and we deform this operator accordingly so that the deforme… ▽ More In this paper we present a method to obtain deformations of families of matrix-valued orthogonal polynomials that are associated to the representation theory of compact Gelfand pairs. These polynomials have the Sturm-Liouville property in the sense that they are simultaneous eigenfunctions of a symmetric second order differential operator and we deform this operator accordingly so that the deformed families also have the Sturm-Liouville property. Our strategy is to deform the system of spherical functions that is related to the matrix-valued orthogonal polynomials and then check that the polynomial structure is respected by the deformation. Crucial in these considerations is the full spherical function $Ψぷさい_{0}$, which relates the spherical functions to the polynomials. We prove an explicit formula for $Ψぷさい_{0}$ in terms of Krawtchouk polynomials for the Gelfand pair $(\mathrm{SU}(2)\times\mathrm{SU}(2),\mathrm{diag}(\mathrm{SU}(2)))$. For the matrix-valued orthogonal polynomials associated to this pair, a deformation was already available by different methods and we show that our method gives same results using explicit knowledge of $Ψぷさい_{0}$. Furthermore we apply our method to some of the examples of size $2\times2$ for more general Gelfand pairs. We prove that the families related to the groups $\mathrm{SU}(n)$ are deformations of one another. On the other hand, the families associated to the symplectic groups $\mathrm{Sp}(n)$ give rise to a new family with an extra free parameter. △ Less

Submitted 4 October, 2016; originally announced October 2016.

Comments: 19 pages

Some bivariate stochastic models arising from group representation theory

Authors: Manuel D. de la Iglesia, Pablo Román

Abstract: The aim of this paper is to study some continuous-time bivariate Markov processes arising from group representation theory. The first component (level) can be either discrete (quasi-birth-and-death processes) or continuous (switching diffusion processes), while the second component (phase) will always be discrete and finite. The infinitesimal operators of these processes will be now matrix-valued… ▽ More The aim of this paper is to study some continuous-time bivariate Markov processes arising from group representation theory. The first component (level) can be either discrete (quasi-birth-and-death processes) or continuous (switching diffusion processes), while the second component (phase) will always be discrete and finite. The infinitesimal operators of these processes will be now matrix-valued (either a block tridiagonal matrix or a matrix-valued second-order differential operator). The matrix-valued spherical functions associated to the compact symmetric pair $(\mathrm{SU}(2)\times \mathrm{SU}(2), \mathrm{diag} \, \mathrm{SU}(2))$ will be eigenfunctions of these infinitesimal operators, so we can perform spectral analysis and study directly some probabilistic aspects of these processes. Among the models we study there will be rational extensions of the one-server queue and Wright-Fisher models involving only mutation effects. △ Less

Submitted 4 October, 2016; originally announced October 2016.

Comments: 22 pages, 7 figures

Saliency Driven Object recognition in egocentric videos with deep CNN

Authors: Philippe Pérez de San Roman, Jenny Benois-Pineau, Jean-Philippe Domenger, Florent Paclet, Daniel Cataert, Aymar de Rugy

Abstract: The problem of object recognition in natural scenes has been recently successfully addressed with Deep Convolutional Neuronal Networks giving a significant break-through in recognition scores. The computational efficiency of Deep CNNs as a function of their depth, allows for their use in real-time applications. One of the key issues here is to reduce the number of windows selected from images to b… ▽ More The problem of object recognition in natural scenes has been recently successfully addressed with Deep Convolutional Neuronal Networks giving a significant break-through in recognition scores. The computational efficiency of Deep CNNs as a function of their depth, allows for their use in real-time applications. One of the key issues here is to reduce the number of windows selected from images to be submitted to a Deep CNN. This is usually solved by preliminary segmentation and selection of specific windows, having outstanding "objectiveness" or other value of indicators of possible location of objects. In this paper we propose a Deep CNN approach and the general framework for recognition of objects in a real-time scenario and in an egocentric perspective. Here the window of interest is built on the basis of visual attention map computed over gaze fixations measured by a glass-worn eye-tracker. The application of this set-up is an interactive user-friendly environment for upper-limb amputees. Vision has to help the subject to control his worn neuro-prosthesis in case of a small amount of remaining muscles when the EMG control becomes unefficient. The recognition results on a specifically recorded corpus of 151 videos with simple geometrical objects show the mAP of 64,6\% and the computational time at the generalization lower than a time of a visual fixation on the object-of-interest. △ Less

Submitted 23 June, 2016; originally announced June 2016.

Comments: 20 pages, 8 figures, 3 tables, Submitted to the Journal of Computer Vision and Image Understanding

Branching rules for finite-dimensional $\mathcal{U}_q(\mathfrak{su}(3))$-representations with respect to a right coideal subalgebra

Authors: Noud Aldenhoven, Erik Koelink, Pablo Román

Abstract: We consider the quantum symmetric pair $(\mathcal{U}_q(\mathfrak{su}(3)), \mathcal{B})$ where $\mathcal{B}$ is a right coideal subalgebra. We prove that all finite-dimensional irreducible representations of $\mathcal{B}$ are weight representations and are characterised by their highest weight and dimension. We show that the restriction of a finite-dimensional irreducible representation of… ▽ More We consider the quantum symmetric pair $(\mathcal{U}_q(\mathfrak{su}(3)), \mathcal{B})$ where $\mathcal{B}$ is a right coideal subalgebra. We prove that all finite-dimensional irreducible representations of $\mathcal{B}$ are weight representations and are characterised by their highest weight and dimension. We show that the restriction of a finite-dimensional irreducible representation of $\mathcal{U}_q(\mathfrak{su}(3))$ to $\mathcal{B}$ decomposes multiplicity free into irreducible representations of $\mathcal{B}$. Furthermore we give explicit expressions for the highest weight vectors in this decomposition in terms of dual $q$-Krawtchouk polynomials. △ Less

Submitted 25 January, 2016; originally announced January 2016.

Comments: 21 pages, 3 figures

Orthogonal vs. Non-Orthogonal Reducibility of Matrix-Valued Measures

Authors: Erik Koelink, Pablo Román

Abstract: A matrix-valued measure $Θしーた$ reduces to measures of smaller size if there exists a constant invertible matrix $M$ such that $MΘしーたM^*$ is block diagonal. Equivalently, the real vector space ${\mathscr A}$ of all matrices $T$ such that $TΘしーた(X)=Θしーた(X) T^*$ for any Borel set $X$ is non-trivial. If the subspace $A_h$ of self-adjoints elements in the commutant algebra $A$ of $Θしーた$ is non-trivial, then $Θしーた$ is re… ▽ More A matrix-valued measure $Θしーた$ reduces to measures of smaller size if there exists a constant invertible matrix $M$ such that $MΘしーたM^*$ is block diagonal. Equivalently, the real vector space ${\mathscr A}$ of all matrices $T$ such that $TΘしーた(X)=Θしーた(X) T^*$ for any Borel set $X$ is non-trivial. If the subspace $A_h$ of self-adjoints elements in the commutant algebra $A$ of $Θしーた$ is non-trivial, then $Θしーた$ is reducible via a unitary matrix. In this paper we prove that ${\mathscr A}$ is $*$-invariant if and only if $A_h={\mathscr A}$, i.e., every reduction of $Θしーた$ can be performed via a unitary matrix. The motivation for this paper comes from families of matrix-valued polynomials related to the group ${\rm SU}(2)\times {\rm SU}(2)$ and its quantum analogue. In both cases the commutant algebra $A=A_h\oplus iA_h$ is of dimension two and the matrix-valued measures reduce unitarily into a $2\times 2$ block diagonal matrix. Here we show that there is no further non-unitary reduction. △ Less

Submitted 23 January, 2016; v1 submitted 21 September, 2015; originally announced September 2015.

Journal ref: SIGMA 12 (2016), 008, 9 pages

Matrix-valued orthogonal polynomials related to the quantum analogue of $(SU(2) \times SU(2), \text{diag})$

Authors: Noud Aldenhoven, Erik Koelink, Pablo Román

Abstract: Matrix-valued spherical functions related to the quantum symmetric pair for the quantum analogue of $(SU(2) \times SU(2), \text{diag})$ are introduced and studied in detail. The quantum symmetric pair is given in terms of a quantised universal enveloping algebra with a coideal subalgebra. The matrix-valued spherical functions give rise to matrix-valued orthogonal polynomials, which are matrix-valu… ▽ More Matrix-valued spherical functions related to the quantum symmetric pair for the quantum analogue of $(SU(2) \times SU(2), \text{diag})$ are introduced and studied in detail. The quantum symmetric pair is given in terms of a quantised universal enveloping algebra with a coideal subalgebra. The matrix-valued spherical functions give rise to matrix-valued orthogonal polynomials, which are matrix-valued analogues of a subfamily of Askey-Wilson polynomials. For these matrix-valued orthogonal polynomials a number of properties are derived using this quantum group interpretation: the orthogonality relations from the Schur orthogonality relations, the three-term recurrence relation and the structure of the weight matrix in terms of Chebyshev polynomials from tensor product decompositions, the matrix-valued Askey-Wilson type $q$-difference operators from the action of the Casimir elements. A more analytic study of the weight gives an explicit LDU-decomposition in terms of continuous $q$-ultraspherical polynomials. The LDU-decomposition gives the possibility to find explicit expressions of the matrix entries of the matrix-valued orthogonal polynomials in terms of continuous $q$-ultraspherical polynomials and $q$-Racah polynomials. △ Less

Submitted 18 February, 2021; v1 submitted 13 July, 2015; originally announced July 2015.

Comments: 42 pages, accompanying sage-file available through author's website. V2: correction of Remark 6.5 with proof, references updated

MSC Class: 33D80; 33D45

Journal ref: Ramanujan Journal (2017) 43, 243-311

A compact concentration of large grains in the HD142527 protoplanetary dust trap

Authors: Simon Casassus, Chris Wright, Sebastian Marino, Sarah T. Maddison, Al Wootten, Pablo Roman, Sebastian Perez, Paola Pinilla, Mark Wyatt, Victor Moral, Francois Menard, Valentin Christiaens, Lucas Cieza, Gerrit van der Plas

Abstract: A pathway to the formation of planetesimals, and eventually giant planets, may occur in concentrations of dust grains trapped in pressure maxima. Dramatic crescent-shaped dust concentrations have been seen in recent radio images at sub-mm wavelengths. These disk asymmetries could represent the initial phases of planet formation in the dust trap scenario, provided that grain sizes are spatially seg… ▽ More A pathway to the formation of planetesimals, and eventually giant planets, may occur in concentrations of dust grains trapped in pressure maxima. Dramatic crescent-shaped dust concentrations have been seen in recent radio images at sub-mm wavelengths. These disk asymmetries could represent the initial phases of planet formation in the dust trap scenario, provided that grain sizes are spatially segregated. A testable prediction of azimuthal dust trapping is that progressively larger grains should be more sharply confined and furthermore the trapped grains should follow a distribution that is markedly different from the gas. However, gas tracers such as CO and the infrared emission from small grains are both very optically thick where the submm continuum originates, so observations have been unable to test the trapping predictions or to identify compact concentrations of larger grains required for planet formation by core-accretion. Here we report multifrequency observations of HD142527, from 34GHz to 700GHz, that reveal a compact concentration of ~cm-sized grains, with a few Earth masses, embedded in a large-scale crescent of ~mm-sized particles. The emission peaks at wavelengths shorter than ~1mm are optically thick and trace the temperature structure resulting from shadows cast by the inner regions. Given this temperature structure, we infer that the largest dust grains are concentrated in the 34 GHz clump. We conclude that dust trapping is efficient for approximately cm-sized grains and leads to enhanced concentrations, while the smaller grains largely reflect the gas distribution. △ Less

Submitted 28 May, 2015; originally announced May 2015.

Comments: Submitted to ApJ

Journal ref: 2015, ApJ, 812, 126

Accretion kinematics through the warped transition disk in HD142527 from resolved CO(6-5) observations

Authors: Simon Casassus, Sebastian Marino, Sebastian Perez, Pablo Roman, Alex Dunhill, Philip Armitage, Jorge Cuadra, Al Wootten, Gerrit van der Plas, Lucas Cieza, Victor Moral, Valentin Christiaens, Matias Montesinos

Abstract: The finding of residual gas in the large central cavity of the HD142527 disk motivates questions on the origin of its non-Keplerian kinematics, and possible connections with planet formation. We aim to understand the physical structure that underlies the intra-cavity gaseous flows, guided by new molecular-line data in CO(6-5) with unprecedented angular resolutions. Given the warped structure infer… ▽ More The finding of residual gas in the large central cavity of the HD142527 disk motivates questions on the origin of its non-Keplerian kinematics, and possible connections with planet formation. We aim to understand the physical structure that underlies the intra-cavity gaseous flows, guided by new molecular-line data in CO(6-5) with unprecedented angular resolutions. Given the warped structure inferred from the identification of scattered-light shadows cast on the outer disk, the kinematics are consistent, to first order, with axisymmetric accretion onto the inner disk occurring at all azimuth. A steady-state accretion profile, fixed at the stellar accretion rate, explains the depth of the cavity as traced in CO isotopologues. The abrupt warp and evidence for near free-fall radial flows in HD 142527 resemble theoretical models for disk tearing, which could be driven by the reported low mass companion, whose orbit may be contained in the plane of the inner disk. The companion's high inclination with respect to the massive outer disk could drive Kozai oscillations over long time-scales; high-eccentricity periods may perhaps account for the large cavity. While shadowing by the tilted disk could imprint an azimuthal modulation in the molecular-line maps, further observations are required to ascertain the significance of azimuthal structure in the density field inside the cavity of HD142527. △ Less

Submitted 28 May, 2015; originally announced May 2015.

Comments: Submitted to ApJ

Planet formation signposts: observability of circumplanetary disks via gas kinematics

Authors: Sebastian Perez, Alex Dunhill, Simon Casassus, Pablo Roman, Judit Szulágyi, Christian Flores, Sebastian Marino, Matias Montesinos

Abstract: The identification of on-going planet formation requires the finest angular resolutions and deepest sensitivities in observations inspired by state-of-the-art numerical simulations. Hydrodynamic simulations of planet-disk interactions predict the formation of circumplanetary disks (CPDs) around accreting planetary cores. These CPDs have eluded unequivocal detection -their identification requires p… ▽ More The identification of on-going planet formation requires the finest angular resolutions and deepest sensitivities in observations inspired by state-of-the-art numerical simulations. Hydrodynamic simulations of planet-disk interactions predict the formation of circumplanetary disks (CPDs) around accreting planetary cores. These CPDs have eluded unequivocal detection -their identification requires predictions in CPD tracers. In this work, we aim to assess the observability of embedded CPDs with ALMA as features imprinted in the gas kinematics. We use 3D Smooth Particle Hydrodynamic (SPH) simulations of CPDs around 1 and 5 M_Jup planets at large stellocentric radii, in locally isothermal and adiabatic disks. The simulations are then connected with 3D radiative transfer for predictions in CO isotopologues. Observability is assessed by corrupting with realistic long baseline phase noise extracted from the recent HL Tau ALMA data. We find that the presence of a CPD produces distinct signposts: 1) compact emission separated in velocity from the overall circumstellar disk's Keplerian pattern, 2) a strong impact on the velocity pattern when the Doppler shifted line emission sweeps across the CPD location, and 3) a local increase in the velocity dispersion. We test our predictions with a simulation tailored for HD 100546 -which has a reported protoplanet candidate. We find that the CPDs are detectable in all 3 signposts with ALMA Cycle 3 capabilities for both 1 and 5 M_Jup protoplanets, when embedded in an isothermal disk. △ Less

Submitted 3 November, 2015; v1 submitted 26 May, 2015; originally announced May 2015.

Comments: 6 pages, 4 figures. Published in ApJL

Journal ref: The Astrophysical Journal Letters, Volume 811, Issue 1, article id. L5, 6 pp. (2015)

Compact dust concentration in the MWC 758 protoplanetary disk

Authors: S. Marino, S. Casassus, S. Perez, W. Lyra, P. E. Roman, H. Avenhaus, C. M. Wright, S. T. Maddison

Abstract: The formation of planetesimals requires that primordial dust grains grow from micron- to km-sized bodies. Dust traps caused by gas pressure maxima have been proposed as regions where grains can concentrate and grow fast enough to form planetesimals, before radially migrating onto the star. We report new VLA Ka & Ku observations of the protoplanetary disk around the Herbig Ae/Be star MWC 758. The K… ▽ More The formation of planetesimals requires that primordial dust grains grow from micron- to km-sized bodies. Dust traps caused by gas pressure maxima have been proposed as regions where grains can concentrate and grow fast enough to form planetesimals, before radially migrating onto the star. We report new VLA Ka & Ku observations of the protoplanetary disk around the Herbig Ae/Be star MWC 758. The Ka image shows a compact emission region in the outer disk indicating a strong concentration of big dust grains. Tracing smaller grains, archival ALMA data in band 7 continuum shows extended disk emission with an intensity maximum to the north-west of the central star, which matches the VLA clump position. The compactness of the Ka emission is expected in the context of dust trapping, as big grains are trapped more easily than smaller grains in gas pressure maxima. We develop a non-axisymmetric parametric model inspired by a steady state vortex solution with parameters adequately selected to reproduce the observations, including the spectral energy distribution. Finally, we compare the radio continuum with SPHERE scattered light data. The ALMA continuum spatially coincides with a spiral-like feature seen in scattered light, while the VLA clump is offset from the scattered light maximum. Moreover, the ALMA map shows a decrement that matches a region devoid of scattered polarised emission. Continuum observations at a different wavelength are necessary to conclude if the VLA-ALMA difference is an opacity or a real dust segregation. △ Less

Submitted 2 October, 2015; v1 submitted 25 May, 2015; originally announced May 2015.

Comments: Accepted for publication in ApJ

Asymptotics of orthogonal polynomials generated by a Geronimus perturbation of the Laguerre measure

Authors: Alfredo Deaño, Edmundo J. Huertas, Pablo Román

Abstract: This paper deals with monic orthogonal polynomials generated by a Geronimus canonical spectral transformation of the Laguerre classical measure: \[ \frac{1}{x-c}x^{αあるふぁ}e^{-x}dx+Nδでるた(x-c), \] for $x\in[0,\infty)$, $αあるふぁ>-1$, a free parameter $N\in \mathbb{R}_{+}$ and a shift $c<0$. We analyze the asymptotic behavior (both strong and relative to classical Laguerre polynomials) of these orthogonal polynomia… ▽ More This paper deals with monic orthogonal polynomials generated by a Geronimus canonical spectral transformation of the Laguerre classical measure: \[ \frac{1}{x-c}x^{αあるふぁ}e^{-x}dx+Nδでるた(x-c), \] for $x\in[0,\infty)$, $αあるふぁ>-1$, a free parameter $N\in \mathbb{R}_{+}$ and a shift $c<0$. We analyze the asymptotic behavior (both strong and relative to classical Laguerre polynomials) of these orthogonal polynomials as $n$ tends to infinity. △ Less

Submitted 22 April, 2015; originally announced April 2015.

Comments: 17 pages, no figures

MSC Class: 33C45; 41A60; 33C15

Gas inside the 97 au cavity around the transition disk Sz\,91

Authors: H. Canovas, M. R. Schreiber, C. Cáceres, F. Ménard, C. Pinte, G. S. Mathews, L. Cieza, S. Casassus, A. Hales, J. P. Williams, P. Román, A. Hardy

Abstract: We present ALMA (Cycle 0) band-6 and band-3 observations of the transition disk Sz\,91. The disk inclination and position angle are determined to be $i=49.5\degr\pm3.5\degr$ and $\mathrm{PA}=18.2\degr\pm3.5\degr$ and the dusty and gaseous disk are detected up to $\sim220$ au and $\sim400$ au from the star, respectively. Most importantly, our continuum observations indicate that the cavity size in… ▽ More We present ALMA (Cycle 0) band-6 and band-3 observations of the transition disk Sz\,91. The disk inclination and position angle are determined to be $i=49.5\degr\pm3.5\degr$ and $\mathrm{PA}=18.2\degr\pm3.5\degr$ and the dusty and gaseous disk are detected up to $\sim220$ au and $\sim400$ au from the star, respectively. Most importantly, our continuum observations indicate that the cavity size in the mm-sized dust distribution must be $\sim97$ au in radius, the largest cavity observed around a T Tauri star. Our data clearly confirms the presence of \co(2-1) well inside the dust cavity. Based on these observational constrains we developed a disk model that simultaneously accounts for the \co and continuum observations (i.e., gaseous and dusty disk). According to our model, most of the millimeter emission comes from a ring located between 97 and 140 au. We also find that the dust cavity is divided into an innermost region largely depleted of dust particles ranging from the dust sublimation radius up to 85 au, and a second, moderately dust-depleted region, extending from 85 to 97 au. The extremely large size of the dust cavity, the presence of gas and small dust particles within the cavity and the accretion rate of Sz\,91 are consistent with the formation of multiple (giant) planets. △ Less

Submitted 16 March, 2015; originally announced March 2015.

Comments: Accepted for publication in ApJ

CO gas inside the protoplanetary disk cavity in HD 142527: disk structure from ALMA

Authors: Sebastian Perez, S. Casassus, F. Ménard, P. Roman, G. van der Plas, L. Cieza, C. Pinte, V. Christiaens, A. S. Hales

Abstract: Inner cavities and annular gaps in circumstellar disks are possible signposts of giant planet formation. The young star HD 142527 hosts a massive protoplanetary disk with a large cavity that extends up to 140 au from the central star, as seen in continuum images at infrared and millimeter wavelengths. Estimates of the survival of gas inside disk cavities are needed to discriminate between clearing… ▽ More Inner cavities and annular gaps in circumstellar disks are possible signposts of giant planet formation. The young star HD 142527 hosts a massive protoplanetary disk with a large cavity that extends up to 140 au from the central star, as seen in continuum images at infrared and millimeter wavelengths. Estimates of the survival of gas inside disk cavities are needed to discriminate between clearing scenarios. We present a spatially and spectrally resolved carbon monoxide isotopologue observations of the gas-rich disk HD 142527, in the J=2-1 line of 12CO, 13CO and C18O, obtained with the Atacama Large Millimeter Array (ALMA). We detect emission coming from inside the dust-depleted cavity in all three isotopologues. Based on our analysis of the gas in the dust cavity, the 12CO emission is optically thick, while 13CO and C18O emission are both optically thin. The total mass of residual gas inside the cavity is about 1.5-2 Jupiter masses. We model the gas with an axisymmetric disk model. Our best fit model shows that the cavity radius is much smaller in CO than it is in millimeter continuum and scattered light observations, with a gas cavity that does not extend beyond 105 au (at 3-sigma). The gap wall at its outer edge is diffuse and smooth in the gas distribution, while in dust continuum it is manifestly sharper. The inclination angle, as estimated from the high velocity channel maps, is 28+/-0.5 degrees, higher than in previous estimates, assuming a fix central star mass of 2.2 Solar masses. △ Less

Submitted 29 October, 2014; originally announced October 2014.

Comments: 13 pages, 10 figures. Accepted for publication in ApJ Main Journal

Asymptotic behavior and zero distribution of polynomials orthogonal with respect to Bessel functions

Authors: Alfredo Deaño, Arno B. J. Kuijlaars, Pablo Román

Abstract: We consider polynomials $P_n$ orthogonal with respect to the weight $J_νにゅー$ on $[0,\infty)$, where $J_νにゅー$ is the Bessel function of order $νにゅー$. Asheim and Huybrechs considered these polynomials in connection with complex Gaussian quadrature for oscillatory integrals. They observed that the zeros are complex and accumulate as $n \to \infty$ near the vertical line $\textrm{Re}\, z = \frac{νにゅーπぱい}{2}$. We pr… ▽ More We consider polynomials $P_n$ orthogonal with respect to the weight $J_νにゅー$ on $[0,\infty)$, where $J_νにゅー$ is the Bessel function of order $νにゅー$. Asheim and Huybrechs considered these polynomials in connection with complex Gaussian quadrature for oscillatory integrals. They observed that the zeros are complex and accumulate as $n \to \infty$ near the vertical line $\textrm{Re}\, z = \frac{νにゅーπぱい}{2}$. We prove this fact for the case $0 \leq νにゅー\leq 1/2$ from strong asymptotic formulas that we derive for the polynomials $P_n$ in the complex plane. Our main tool is the Riemann-Hilbert problem for orthogonal polynomials, suitably modified to cover the present situation, and the Deift-Zhou steepest descent method. A major part of the work is devoted to the construction of a local parametrix at the origin, for which we give an existence proof that only works for $νにゅー\leq 1/2$. △ Less

Submitted 4 June, 2014; originally announced June 2014.

Comments: 42 pages, 5 figures

MSC Class: 33C47; 34M50; 30E15; 33C10

Journal ref: Constructive Approximation 43 (2016), 153-196

Matrix-valued Gegenbauer polynomials

Authors: Erik Koelink, Ana M. de los Rios, Pablo Roman

Abstract: We introduce matrix-valued weight functions of arbitrary size, which are analogues of the weight function for the Gegenbauer or ultraspherical polynomials for the parameter $νにゅー>0$. The LDU-decomposition of the weight is explicitly given in terms of Gegenbauer polynomials. We establish a matrix-valued Pearson equation for these matrix weights leading to explicit shift operators relating the weights… ▽ More We introduce matrix-valued weight functions of arbitrary size, which are analogues of the weight function for the Gegenbauer or ultraspherical polynomials for the parameter $νにゅー>0$. The LDU-decomposition of the weight is explicitly given in terms of Gegenbauer polynomials. We establish a matrix-valued Pearson equation for these matrix weights leading to explicit shift operators relating the weights for parameter $νにゅー$ and $νにゅー+1$. The matrix coefficients of the Pearson equation are obtained using a special matrix-valued differential operator in a commutative algebra of symmetric differential operators. The corresponding orthogonal polynomials are the matrix-valued Gegenbauer polynomials which are eigenfunctions for the symmetric matrix-valued differential operators. Using the shift operators we find the squared norm and we establish a simple Rodrigues formula. The three-term recurrence relation is obtained explicitly using the shift operators as well. We give an explicit non-trivial expression for the matrix entries of the matrix-valued Gegenbauer polynomials in terms of scalar-valued Gegenbauer and Racah polynomials using the LDU-decomposition and differential operators. The case $νにゅー=1$ reduces to the case of matrix-valued Chebyshev polynomials previously obtained using group theoretic considerations. △ Less

Submitted 14 April, 2016; v1 submitted 12 March, 2014; originally announced March 2014.

Comments: 25 pages, accompanying Maple worksheet available at http://www.math.ru.nl/~koelink/publist-ro.html. Proofs have been simplified by exploiting shift operators extensively. The proofs of the symmetry of some matrix-valued differential operators have been simplified using the LDU-decomposition of the weight. The proof of the matrix-valued Pearson equation has been simplified using the shift operator

Matrix Valued Classical Pairs Related to Compact Gelfand Pairs of Rank One

Authors: Maarten van Pruijssen, Pablo Román

Abstract: We present a method to obtain infinitely many examples of pairs $(W,D)$ consisting of a matrix weight $W$ in one variable and a symmetric second-order differential operator $D$. The method is based on a uniform construction of matrix valued polynomials starting from compact Gelfand pairs $(G,K)$ of rank one and a suitable irreducible $K$-representation. The heart of the construction is the existen… ▽ More We present a method to obtain infinitely many examples of pairs $(W,D)$ consisting of a matrix weight $W$ in one variable and a symmetric second-order differential operator $D$. The method is based on a uniform construction of matrix valued polynomials starting from compact Gelfand pairs $(G,K)$ of rank one and a suitable irreducible $K$-representation. The heart of the construction is the existence of a suitable base change $Ψぷさい_{0}$. We analyze the base change and derive several properties. The most important one is that $Ψぷさい_{0}$ satisfies a first-order differential equation which enables us to compute the radial part of the Casimir operator of the group $G$ as soon as we have an explicit expression for $Ψぷさい_{0}$. The weight $W$ is also determined by $Ψぷさい_{0}$. We provide an algorithm to calculate $Ψぷさい_{0}$ explicitly. For the pair $(\mathrm{USp}(2n),\mathrm{USp}(2n-2)\times\mathrm{USp}(2))$ we have implemented the algorithm in GAP so that individual pairs $(W,D)$ can be calculated explicitly. Finally we classify the Gelfand pairs $(G,K)$ and the $K$-representations that yield pairs $(W,D)$ of size $2\times2$ and we provide explicit expressions for most of these cases. △ Less

Submitted 20 December, 2014; v1 submitted 23 December, 2013; originally announced December 2013.

MSC Class: 22E46; 33C47

Journal ref: SIGMA 10 (2014), 113, 28 pages

Observations of gas flows inside a protoplanetary gap

Authors: Simon Casassus, Gerrit van der Plas, Sebastian Perez M., William R. F. Dent, Ed Fomalont, Janis Hagelberg, Antonio Hales, Andrés Jordán, Dimitri Mawet, Francois Ménard, Al Wootten, David Wilner, A. Meredith Hughes, Matthias R. Schreiber, Julien H. Girard, Barbara Ercolano, Hector Canovas, Pablo E. Román, Vachail Salinas

Abstract: Gaseous giant planet formation is thought to occur in the first few million years following stellar birth. Models predict that giant planet formation carves a deep gap in the dust component (shallower in the gas). Infrared observations of the disk around the young star HD142527, at ~140pc, found an inner disk ~10AUえーゆー in radius, surrounded by a particularly large gap, with a disrupted outer disk beyo… ▽ More Gaseous giant planet formation is thought to occur in the first few million years following stellar birth. Models predict that giant planet formation carves a deep gap in the dust component (shallower in the gas). Infrared observations of the disk around the young star HD142527, at ~140pc, found an inner disk ~10AUえーゆー in radius, surrounded by a particularly large gap, with a disrupted outer disk beyond 140AUえーゆー, indicative of a perturbing planetary-mass body at ~90 AUえーゆー. From radio observations, the bulk mass is molecular and lies in the outer disk, whose continuum emission has a horseshoe morphology. The vigorous stellar accretion rate would deplete the inner disk in less than a year, so in order to sustain the observed accretion, matter must flow from the outer-disk into the cavity and cross the gap. In dynamical models, the putative protoplanets channel outer-disk material into gap-crossing bridges that feed stellar accretion through the inner disk. Here we report observations with the Atacama Large Millimetre Array (ALMA) that reveal diffuse CO gas inside the gap, with denser HCO+ gas along gap-crossing filaments, and that confirm the horseshoe morphology of the outer disk. The estimated flow rate of the gas is in the range 7E-9 to 2E-7 Msun/yr, which is sufficient to maintain accretion onto the star at the present rate. △ Less

Submitted 26 May, 2013; originally announced May 2013.

Journal ref: 2013, Nature, 493, 191

Matrix-valued orthogonal polynomials related to (SU(2)\times SU(2),diag), II

Authors: Erik Koelink, Maarten van Pruijssen, Pablo Roman

Abstract: In a previous paper we have introduced matrix-valued analogues of the Chebyshev polynomials by studying matrix-valued spherical functions on SU(2)\times SU(2). In particular the matrix-size of the polynomials is arbitrarily large. The matrix-valued orthogonal polynomials and the corresponding weight function are studied. In particular, we calculate the LDU-decomposition of the weight where the mat… ▽ More In a previous paper we have introduced matrix-valued analogues of the Chebyshev polynomials by studying matrix-valued spherical functions on SU(2)\times SU(2). In particular the matrix-size of the polynomials is arbitrarily large. The matrix-valued orthogonal polynomials and the corresponding weight function are studied. In particular, we calculate the LDU-decomposition of the weight where the matrix entries of L are given in terms of Gegenbauer polynomials. The monic matrix-valued orthogonal polynomials P_n are expressed in terms of Tirao's matrix-valued hypergeometric function using the matrix-valued differential operator of first and second order to which the P_n's are eigenfunctions. From this result we obtain an explicit formula for coefficients in the three-term recurrence relation satisfied by the polynomials P_n. These differential operators are also crucial in expressing the matrix entries of P_nL as a product of a Racah and a Gegenbauer polynomial. We also present a group theoretic derivation of the matrix-valued differential operators by considering the Casimir operators corresponding to SU(2)\times SU(2). △ Less

Submitted 25 October, 2012; v1 submitted 29 February, 2012; originally announced March 2012.

Comments: 35 pages, sequel to arXiv:1012.2719, incorporating referee's comments (including change in title)

Journal ref: Publ. RIMS Kyoto 49, no. 2, (2013), 271-312

Non-intersecting squared Bessel paths with one positive starting and ending point

Authors: Steven Delvaux, Arno B. J. Kuijlaars, Pablo Román, Lun Zhang

Abstract: We consider a model of $n$ non-intersecting squared Bessel processes with one starting point $a>0$ at time t=0 and one ending point $b>0$ at time $t=T$. After proper scaling, the paths fill out a region in the $tx$-plane. Depending on the value of the product $ab$ the region may come to the hard edge at 0, or not. We formulate a vector equilibrium problem for this model, which is defined for three… ▽ More We consider a model of $n$ non-intersecting squared Bessel processes with one starting point $a>0$ at time t=0 and one ending point $b>0$ at time $t=T$. After proper scaling, the paths fill out a region in the $tx$-plane. Depending on the value of the product $ab$ the region may come to the hard edge at 0, or not. We formulate a vector equilibrium problem for this model, which is defined for three measures, with upper constraints on the first and third measures and an external field on the second measure. It is shown that the limiting mean distribution of the paths at time $t$ is given by the second component of the vector that minimizes this vector equilibrium problem. The proof is based on a steepest descent analysis for a $4 \times 4$ matrix valued Riemann-Hilbert problem which characterizes the correlation kernel of the paths at time $t$. We also discuss the precise locations of the phase transitions. △ Less

Submitted 13 May, 2011; v1 submitted 12 May, 2011; originally announced May 2011.

Comments: 51 pages, 9 figures

Matrix Valued Orthogonal Polynomials related to (SU(2) x SU(2),diag)

Authors: Erik Koelink, Maarten van Pruijssen, Pablo Roman

Abstract: The matrix-valued spherical functions for the pair (K x K, K), K=SU(2), are studied. By restriction to the subgroup A the matrix-valued spherical functions are diagonal. For suitable set of representations we take these diagonals into a matrix-valued function, which are the full spherical functions. Their orthogonality is a consequence of the Schur orthogonality relations. From the full spherical… ▽ More The matrix-valued spherical functions for the pair (K x K, K), K=SU(2), are studied. By restriction to the subgroup A the matrix-valued spherical functions are diagonal. For suitable set of representations we take these diagonals into a matrix-valued function, which are the full spherical functions. Their orthogonality is a consequence of the Schur orthogonality relations. From the full spherical functions we obtain matrix-valued orthogonal polynomials of arbitrary size, and they satisfy a three-term recurrence relation which follows by considering tensor product decompositions. An explicit expression for the weight and the complete block-diagonalization of the matrix-valued orthogonal polynomials is obtained. From the explicit expression we obtain right-hand sided differential operators of first and second order for which the matrix-valued orthogonal polynomials are eigenfunctions. We study the low-dimensional cases explicitly, and for these cases additional results, such as the Rodrigues formula and being eigenfunctions to first order differential-difference and second order differential operators, are obtained. △ Less

Submitted 18 January, 2011; v1 submitted 13 December, 2010; originally announced December 2010.

Comments: 33 pages, 4 figures

Journal ref: Int Math Res Notices (2012) 2012 (24): 5673-5730

Asymptotic zero distribution of multiple orthogonal polynomials associated with Macdonald functions

Authors: Lun Zhang, Pablo Román

Abstract: We study the asymptotic zero distribution of type II multiple orthogonal polynomials associated with two Macdonald functions (modified Bessel functions of the second kind). Based on the four-term recurrence relation, it is shown that, after proper scaling, the sequence of normalized zero counting measures converges weakly to the first component of a vector of two measures which satisfies a vector… ▽ More We study the asymptotic zero distribution of type II multiple orthogonal polynomials associated with two Macdonald functions (modified Bessel functions of the second kind). Based on the four-term recurrence relation, it is shown that, after proper scaling, the sequence of normalized zero counting measures converges weakly to the first component of a vector of two measures which satisfies a vector equilibrium problem with two external fields. We also give the explicit formula for the equilibrium vector in terms of solutions of an algebraic equation. △ Less

Submitted 2 February, 2011; v1 submitted 24 March, 2010; originally announced March 2010.

Comments: 22 pages, 4 figures

Journal ref: Journal of Approximation Theory 163 (2011), 143--162

Recurrence relations and vector equilibrium problems arising from a model of non-intersecting squared Bessel paths

Authors: A. B. J. Kuijlaars, P. Román

Abstract: In this paper we consider the model of $n$ non-intersecting squared Bessel processes with parameter $αあるふぁ$, in the confluent case where all particles start, at time $t=0$, at the same positive value $x=a$, remain positive, and end, at time $T=t$, at the position $x=0$. The positions of the paths have a limiting mean density as $n\to\infty$ which is characterized by a vector equilibrium problem. We… ▽ More In this paper we consider the model of $n$ non-intersecting squared Bessel processes with parameter $αあるふぁ$, in the confluent case where all particles start, at time $t=0$, at the same positive value $x=a$, remain positive, and end, at time $T=t$, at the position $x=0$. The positions of the paths have a limiting mean density as $n\to\infty$ which is characterized by a vector equilibrium problem. We show how to obtain this equilibrium problem from different considerations involving the recurrence relations for multiple orthogonal polynomials associated with the modified Bessel functions. We also extend the situation by rescaling the parameter $αあるふぁ$, letting it increase proportionally to $n$ as $n$ increases. In this case we also analyze the recurrence relation and obtain a vector equilibrium problem for it. △ Less

Submitted 19 November, 2009; originally announced November 2009.

Comments: 28 pages, 10 figures

Synchronization of Coupled Anizochronous Auto-Oscillating Systems

Authors: A. P. Kuznetsov, Ju. P. Roman

Abstract: The particular properties of synchronization are discussed for coupled auto-oscillating systems, which are characterized by non-quadratic law of potential dependence on the coordinate. In particular, structure of the parameter plane (frequency mismatch - coupling value) is considered for coupled van der Pol - Duffing oscillators. The arrangement of synchronization tongues and the particular prop… ▽ More The particular properties of synchronization are discussed for coupled auto-oscillating systems, which are characterized by non-quadratic law of potential dependence on the coordinate. In particular, structure of the parameter plane (frequency mismatch - coupling value) is considered for coupled van der Pol - Duffing oscillators. The arrangement of synchronization tongues and the particular properties of their internal structure in the parameter space are revealed. The features of attractors in the phase space are discussed. △ Less

Submitted 16 October, 2008; originally announced October 2008.

Comments: 6 pages, 5 figures

Journal ref: NPCS 2009, Vol. 12, No. 1, pp.54-60

Synchronization of coupled anizochronous oscillators

Authors: A. P. Kuznetsov, J. P. Roman

Abstract: The particular properties of synchronization are discussed for the system of coupled van der Pol - Duffing oscillators. The arrangement of synchronization tongues and the particular properties of their internal structure in the parameter space are revealed. The features of attractors in the phase space and in the Poincare section are considered. The particular properties of synchronization are discussed for the system of coupled van der Pol - Duffing oscillators. The arrangement of synchronization tongues and the particular properties of their internal structure in the parameter space are revealed. The features of attractors in the phase space and in the Poincare section are considered. △ Less

Submitted 9 June, 2008; originally announced June 2008.

Comments: 13 pages, 7 figures

Properties of synchronization in the systems of non-identical coupled van der Pol and van der Pol - Duffing oscillators. Broadband synchronization

Authors: Alexander P. Kuznetsov, Julia P. Roman

Abstract: The particular properties of dynamics are discussed for the dissipatively coupled van der Pol oscillators, non-identical in values of parameters controlling the Hopf bifurcation. Possibility of a special synchronization regime in an infinitively long band between oscillation death and quasiperiodic areas is shown for such system. Features of the bifurcation picture are discussed for different va… ▽ More The particular properties of dynamics are discussed for the dissipatively coupled van der Pol oscillators, non-identical in values of parameters controlling the Hopf bifurcation. Possibility of a special synchronization regime in an infinitively long band between oscillation death and quasiperiodic areas is shown for such system. Features of the bifurcation picture are discussed for different values of the control parameters and for the case of additional Duffing-type nonlinearity. Analysis of the abridged equations is presented. △ Less

Submitted 8 May, 2008; originally announced May 2008.

Comments: 19 pages, 9 figures

The generalized matrix valued hypergeometric equation

Authors: P. Roman, S. Simondi

Abstract: The matrix valued analog of the Euler's hypergeometric differential equation was introduced by Tirao in \cite{T2}. This equation arises in the study of matrix valued spherical functions and in the theory of matrix valued orthogonal polynomials. The goal of this paper is to extend naturally the number of parameters of Tirao's equation in order to get a generalized matrix valued hypergeometric equ… ▽ More The matrix valued analog of the Euler's hypergeometric differential equation was introduced by Tirao in \cite{T2}. This equation arises in the study of matrix valued spherical functions and in the theory of matrix valued orthogonal polynomials. The goal of this paper is to extend naturally the number of parameters of Tirao's equation in order to get a generalized matrix valued hypergeometric equation. We take advantage of the tools and strategies developed in \cite{T2} to identify the corresponding matrix hypergeometric functions ${}_nF_m$. We prove that, if n=m+1, this functions are analytic for |z|<1 and we give a necesary condition for the convergence on the unit circle |z|=1. △ Less

Submitted 19 October, 2007; v1 submitted 2 August, 2007; originally announced August 2007.

MSC Class: 33C20; 33C80