In this paper, we study the singularities of Feynman integrals using homological techniques. We a... more In this paper, we study the singularities of Feynman integrals using homological techniques. We analyse the Feynman integrals by compactifying the integration domain as well as the ambient space by embedding them in higher-dimensional space. In this compactified space the singularities occur due to the meeting of compactified propagators in the non-general position. The present analysis which had been previously used only for the singularities of second-type is further extended to study other kinds of singularities viz threshold, pseudo-threshold and anomalous threshold singularities. We will study various one-loop and two-loop examples and obtain their various singularities. We will also discuss how using the method we can also determine whether the singularities lie on the physical sheet or not.
Motivated by the foundational work of Tarasov, who pointed out that the algebraic relations of th... more Motivated by the foundational work of Tarasov, who pointed out that the algebraic relations of the type considered here can lead to functional reduction of Feynman integrals, we suitably modify the original method to be able to implement and automatize it and present a Mathematica package AlgRel.wl. The purpose of this package is to help derive the algebraic relations with arbitrary kinematic quantities, for the product of propagators. Under specific choices of the arbitrary parameters that appear in these relations, we can write the original integral with all massive propagators in general, as a sum of integrals which have fewer massive propagators. The resulting integrals are of reduced complexity for computational purposes. For the one-loop cases, with all different and non-zero masses, this would result in integrals with one massive propagator. We also devise a strategy so that the method can also be applied to higher-loop integrals. We demonstrate the procedure and the results obtained using the package for various oneloop and higher-loop examples. Due to the fact that the Feynman integrals are intimately related to the hypergeometric functions, a useful consequence of these algebraic relations is in deriving the sets of non-trivial reduction formulae. We present various such reduction formulae and further discuss how, more such formulae can be obtained than described here. The AlgRel.wl package and an example notebook Examples.nb can be found at GitHub.
We present new closed-form expressions for certain improper integrals of Mathematical Physics suc... more We present new closed-form expressions for certain improper integrals of Mathematical Physics such as Ising, Box, and Associated integrals. The techniques we employ here include (a) the Method of Brackets and its modifications and suitable extensions and (b) the evaluation of the resulting Mellin-Barnes representations via the recently discovered Conic Hull method. Analytic continuations of these series solutions are then produced using the automated method of Olsson. Thus, combining all the recent advances allows for closed-form solutions for the hitherto unknown B 3 (s) and related integrals in terms of multivariable hypergeometric functions. Along the way, we also discuss certain complications while using the Original Method of Brackets for these evaluations and how to rectify them. The interesting cases of C 5,k is also studied. It is not yet fully resolved for the reasons we discuss in this paper.
Considering the large q expansion of the Sachdev-Ye-Kitaev (SYK) model in the two-stage limit, we... more Considering the large q expansion of the Sachdev-Ye-Kitaev (SYK) model in the two-stage limit, we compute the Lanczos coefficients, Krylov complexity, and the higher Krylov cumulants in subleading order, along with the t/q effects. The Krylov complexity naturally describes the “size” of the distribution while the higher cumulants encode richer information. We further consider the double-scaled limit of SYKq at infinite temperature, where q ~ $$ \sqrt{N} $$ N . In such a limit, we find that the scrambling time shrinks to zero, and the Lanczos coefficients diverge. The growth of Krylov complexity appears to be “hyperfast”, which is previously conjectured to be associated with scrambling in de Sitter space.
We study the operator growth in open quantum systems with dephasing dissipation terms, extending ... more We study the operator growth in open quantum systems with dephasing dissipation terms, extending the Krylov complexity formalism of [1]. Our results are based on the study of the dissipative q-body Sachdev-Ye-Kitaev (SYKq) model, governed by the Markovian dynamics. We introduce a notion of “operator size concentration” which allows a diagrammatic and combinatorial proof of the asymptotic linear behavior of the two sets of Lanczos coefficients (an and bn) in the large q limit. Our results corroborate with the semi-analytics in finite q in the large N limit, and the numerical Arnoldi iteration in finite q and finite N limit. As a result, Krylov complexity exhibits exponential growth following a saturation at a time that grows logarithmically with the inverse dissipation strength. The growth of complexity is suppressed compared to the closed system results, yet it upper bounds the growth of the normalized out-of-time-ordered correlator (OTOC). We provide a plausible explanation of the ...
Considering the large q expansion of the Sachdev-Ye-Kitaev (SYK) model in the two-stage limit, we... more Considering the large q expansion of the Sachdev-Ye-Kitaev (SYK) model in the two-stage limit, we compute the Lanczos coefficients, Krylov complexity, and the higher Krylov cumulants in subleading order, along with the t/q effects. The Krylov complexity naturally describes the "size" of the distribution while the higher cumulants encode richer information. We further consider the double-scaled limit of SYK q at infinite temperature, where q ∼ √ N. In such a limit, we find that the scrambling time shrinks to zero, and the Lanczos coefficients diverge. The growth of Krylov complexity appears to be "hyperfast", which is previously conjectured to be associated with scrambling in de Sitter space.
We study the operator growth in open quantum systems with dephasing dissipation terms, extending ... more We study the operator growth in open quantum systems with dephasing dissipation terms, extending the Krylov complexity formalism of [1]. Our results are based on the study of the dissipative q-body Sachdev-Ye-Kitaev (SYK q) model, governed by the Markovian dynamics. We introduce a notion of "operator size concentration" which allows a diagrammatic and combinatorial proof of the asymptotic linear behavior of the two sets of Lanczos coefficients (a n and b n) in the large q limit. Our results corroborate with the semi-analytics in finite q in the large N limit, and the numerical Arnoldi iteration in finite q and finite N limit. As a result, Krylov complexity exhibits exponential growth following a saturation at a time that grows logarithmically with the inverse dissipation strength. The growth of complexity is suppressed compared to the closed system results, yet it upper bounds the growth of the normalized out-of-time-ordered correlator (OTOC). We provide a plausible explanation of the results from the dual gravitational side.
The analytic continuations (ACs) of the double variable Horn H 1 and H 5 functions have been deri... more The analytic continuations (ACs) of the double variable Horn H 1 and H 5 functions have been derived for the first time using the Mathematica package Olsson.wl. The corresponding region of convergences (ROCs) of the ACs are obtained using a companion package ROC2.wl. A Mathematica package HornH1H5.wl, containing all the derived ACs and the associated ROCs, along with a demonstration file of the same is publicly available in this URL : https://github.com/souvik5151/Horn H1 H5.git
The transformation theory of the Appell F2(a, b1, b2; c1, c2; x, y) double hypergeometric functio... more The transformation theory of the Appell F2(a, b1, b2; c1, c2; x, y) double hypergeometric function is used to obtain a set of series representations of F2 which provide an efficient way to evaluate F2 for real values of its arguments x and y and generic complex values of its parameters a, b1, b2, c1 and c2 (i.e. in the nonlogarithmic case). This study rests on a classical approach where the usual double series representation of F2 and other double hypergeometric series that appear in the intermediate steps of the calculations are written as infinite sums of one variable hypergeometric series, such as the Gauss 2F1 or the 3F2, various linear transformations of the latter being then applied to derive known and new formulas. Using the three well-known Euler transformations of F2 on these results allows us to obtain a total of 44 series which form the basis of the Mathematica package AppellF2, dedicated to the evaluation of F2. A brief description of the package and of the numerical analysis that we have performed to test it are also presented.
The Method of Brackets (MoB) is a technique used to compute definite integrals, that has its orig... more The Method of Brackets (MoB) is a technique used to compute definite integrals, that has its origin in the negative dimensional integration method. It was originally proposed for the evaluation of Feynman integrals for which, when applicable, it gives the results in terms of combinations of (multiple) series. We focus here on some of the limitations of MoB and address them by studying the Mellin-Barnes (MB) representation technique. There has been significant process recently in the study of the latter due to the development of a new computational approach based on conic hulls (see Phys. Rev. Lett. 127, 151601 (2021)). The comparison between the two methods helps to understand the limitations of the MoB, in particular when termwise divergent series appear. As a consequence, the MB technique is found to be superior over MoB for two major reasons: 1. the selection of the sets of series that form a series representation for a given integral follows, in the MB approach, from specific intersections of conic hulls, which, in contrast to MoB, does not need any convergence analysis of the involved series, and 2. MB can be used to evaluate resonant (i.e. logarithmic) cases where MoB fails due to the appearance of termwise divergent series. Furthermore, we show that the recently added Rule 5 of MoB naturally emerges as a consequence of the residue theorem in the context of MB.
Capacity of entanglement (CoE), an information-theoretic measure of entanglement, defined as the ... more Capacity of entanglement (CoE), an information-theoretic measure of entanglement, defined as the variance of modular Hamiltonian, is known to capture the deviation from the maximal entanglement. We derive an exact expression for the average eigenstate CoE in fermionic Gaussian states as a finite series, valid for arbitrary bi-partition of the total system. Further, we consider the complex SYK2 model, and we show that the average CoE in the thermodynamic limit upper bounds the average CoE for fermionic Gaussian states. In this limit, the variance of the average CoE becomes independent of the system size. Moreover, when the subsystem size is half of the total system, the leading volume-law coefficient approaches a value of π 2 /8 − 1. We identify this as a distinguishing feature between integrable and quantum-chaotic systems. We confirm our analytical results by numerical computations.
We present the Olsson.wl Mathematica package which aims to find linear transformations for some c... more We present the Olsson.wl Mathematica package which aims to find linear transformations for some classes of multivariable hypergeometric functions. It is based on a well-known method developed by P. O. M. Olsson in J. Math. Phys. 5, 420 (1964) in order to derive the analytic continuations of the Appell $F_1$ double hypergeometric series from the linear transformations of the Gauss $_2F_1$ hypergeometric function. We provide a brief description of Olsson's method and demonstrate the commands of the package, along with examples. We also provide a companion package, called ROC2.wl and dedicated to the derivation of the regions of convergence of double hypergeometric series. This package can be used independently of Olsson.wl.
In semi-classical systems, the exponential growth of the out-of-time-order correlator (OTOC) is b... more In semi-classical systems, the exponential growth of the out-of-time-order correlator (OTOC) is believed to be the hallmark of quantum chaos. However, on several occasions, it has been argued that, even in integrable systems, OTOC can grow exponentially due to the presence of unstable saddle points in the phase space. In this work, we probe such an integrable system exhibiting saddle-dominated scrambling through Krylov complexity and the associated Lanczos coefficients. In the realm of the universal operator growth hypothesis, we demonstrate that the Lanczos coefficients follow the linear growth, which ensures the exponential behavior of Krylov complexity at early times. The linear growth arises entirely due to the saddle, which dominates other phase-space points even away from itself. Our results reveal that the exponential growth of Krylov complexity can be observed in integrable systems with saddle-dominated scrambling and thus need not be associated with the presence of chaos.
Capacity of entanglement (CoE), an information-theoretic measure of entanglement, defined as the ... more Capacity of entanglement (CoE), an information-theoretic measure of entanglement, defined as the variance of modular Hamiltonian, is known to capture the deviation from the maximal entanglement. We derive an exact expression for the average eigenstate CoE in fermionic Gaussian states as a finite series, valid for arbitrary bi-partition of the total system. Further, we consider the complex SYK2 model, and we show that the average CoE in the thermodynamic limit upper bounds the average CoE for fermionic Gaussian states. In this limit, the variance of the average CoE becomes independent of the system size. Moreover, when the subsystem size is half of the total system, the leading volume-law coefficient approaches a value of π/8 − 1. We identify this as a distinguishing feature between integrable and quantum-chaotic systems. We confirm our analytical results by numerical computations.
We use the method of brackets to evaluate quadratic and quartic type integrals. We recall the ope... more We use the method of brackets to evaluate quadratic and quartic type integrals. We recall the operational rules of the method and give examples to illustrate its working. The method is then used to evaluate the quadratic type integrals which occur in entries 3.251.1,3,4 in the table of integrals by Gradshteyn and Ryzhik and obtain closed form expressions in terms of hypergeometric functions. The method is further used to evaluate the quartic integrals, entry 2.161.5 and 6 in the table. We also present generalization of both types of integrals with closed form expression in terms of hypergeometric functions.
In this paper, we study the singularities of Feynman integrals using homological techniques. We a... more In this paper, we study the singularities of Feynman integrals using homological techniques. We analyse the Feynman integrals by compactifying the integration domain as well as the ambient space by embedding them in higher-dimensional space. In this compactified space the singularities occur due to the meeting of compactified propagators in the non-general position. The present analysis which had been previously used only for the singularities of second-type is further extended to study other kinds of singularities viz threshold, pseudo-threshold and anomalous threshold singularities. We will study various one-loop and two-loop examples and obtain their various singularities. We will also discuss how using the method we can also determine whether the singularities lie on the physical sheet or not.
Motivated by the foundational work of Tarasov, who pointed out that the algebraic relations of th... more Motivated by the foundational work of Tarasov, who pointed out that the algebraic relations of the type considered here can lead to functional reduction of Feynman integrals, we suitably modify the original method to be able to implement and automatize it and present a Mathematica package AlgRel.wl. The purpose of this package is to help derive the algebraic relations with arbitrary kinematic quantities, for the product of propagators. Under specific choices of the arbitrary parameters that appear in these relations, we can write the original integral with all massive propagators in general, as a sum of integrals which have fewer massive propagators. The resulting integrals are of reduced complexity for computational purposes. For the one-loop cases, with all different and non-zero masses, this would result in integrals with one massive propagator. We also devise a strategy so that the method can also be applied to higher-loop integrals. We demonstrate the procedure and the results obtained using the package for various oneloop and higher-loop examples. Due to the fact that the Feynman integrals are intimately related to the hypergeometric functions, a useful consequence of these algebraic relations is in deriving the sets of non-trivial reduction formulae. We present various such reduction formulae and further discuss how, more such formulae can be obtained than described here. The AlgRel.wl package and an example notebook Examples.nb can be found at GitHub.
We present new closed-form expressions for certain improper integrals of Mathematical Physics suc... more We present new closed-form expressions for certain improper integrals of Mathematical Physics such as Ising, Box, and Associated integrals. The techniques we employ here include (a) the Method of Brackets and its modifications and suitable extensions and (b) the evaluation of the resulting Mellin-Barnes representations via the recently discovered Conic Hull method. Analytic continuations of these series solutions are then produced using the automated method of Olsson. Thus, combining all the recent advances allows for closed-form solutions for the hitherto unknown B 3 (s) and related integrals in terms of multivariable hypergeometric functions. Along the way, we also discuss certain complications while using the Original Method of Brackets for these evaluations and how to rectify them. The interesting cases of C 5,k is also studied. It is not yet fully resolved for the reasons we discuss in this paper.
Considering the large q expansion of the Sachdev-Ye-Kitaev (SYK) model in the two-stage limit, we... more Considering the large q expansion of the Sachdev-Ye-Kitaev (SYK) model in the two-stage limit, we compute the Lanczos coefficients, Krylov complexity, and the higher Krylov cumulants in subleading order, along with the t/q effects. The Krylov complexity naturally describes the “size” of the distribution while the higher cumulants encode richer information. We further consider the double-scaled limit of SYKq at infinite temperature, where q ~ $$ \sqrt{N} $$ N . In such a limit, we find that the scrambling time shrinks to zero, and the Lanczos coefficients diverge. The growth of Krylov complexity appears to be “hyperfast”, which is previously conjectured to be associated with scrambling in de Sitter space.
We study the operator growth in open quantum systems with dephasing dissipation terms, extending ... more We study the operator growth in open quantum systems with dephasing dissipation terms, extending the Krylov complexity formalism of [1]. Our results are based on the study of the dissipative q-body Sachdev-Ye-Kitaev (SYKq) model, governed by the Markovian dynamics. We introduce a notion of “operator size concentration” which allows a diagrammatic and combinatorial proof of the asymptotic linear behavior of the two sets of Lanczos coefficients (an and bn) in the large q limit. Our results corroborate with the semi-analytics in finite q in the large N limit, and the numerical Arnoldi iteration in finite q and finite N limit. As a result, Krylov complexity exhibits exponential growth following a saturation at a time that grows logarithmically with the inverse dissipation strength. The growth of complexity is suppressed compared to the closed system results, yet it upper bounds the growth of the normalized out-of-time-ordered correlator (OTOC). We provide a plausible explanation of the ...
Considering the large q expansion of the Sachdev-Ye-Kitaev (SYK) model in the two-stage limit, we... more Considering the large q expansion of the Sachdev-Ye-Kitaev (SYK) model in the two-stage limit, we compute the Lanczos coefficients, Krylov complexity, and the higher Krylov cumulants in subleading order, along with the t/q effects. The Krylov complexity naturally describes the "size" of the distribution while the higher cumulants encode richer information. We further consider the double-scaled limit of SYK q at infinite temperature, where q ∼ √ N. In such a limit, we find that the scrambling time shrinks to zero, and the Lanczos coefficients diverge. The growth of Krylov complexity appears to be "hyperfast", which is previously conjectured to be associated with scrambling in de Sitter space.
We study the operator growth in open quantum systems with dephasing dissipation terms, extending ... more We study the operator growth in open quantum systems with dephasing dissipation terms, extending the Krylov complexity formalism of [1]. Our results are based on the study of the dissipative q-body Sachdev-Ye-Kitaev (SYK q) model, governed by the Markovian dynamics. We introduce a notion of "operator size concentration" which allows a diagrammatic and combinatorial proof of the asymptotic linear behavior of the two sets of Lanczos coefficients (a n and b n) in the large q limit. Our results corroborate with the semi-analytics in finite q in the large N limit, and the numerical Arnoldi iteration in finite q and finite N limit. As a result, Krylov complexity exhibits exponential growth following a saturation at a time that grows logarithmically with the inverse dissipation strength. The growth of complexity is suppressed compared to the closed system results, yet it upper bounds the growth of the normalized out-of-time-ordered correlator (OTOC). We provide a plausible explanation of the results from the dual gravitational side.
The analytic continuations (ACs) of the double variable Horn H 1 and H 5 functions have been deri... more The analytic continuations (ACs) of the double variable Horn H 1 and H 5 functions have been derived for the first time using the Mathematica package Olsson.wl. The corresponding region of convergences (ROCs) of the ACs are obtained using a companion package ROC2.wl. A Mathematica package HornH1H5.wl, containing all the derived ACs and the associated ROCs, along with a demonstration file of the same is publicly available in this URL : https://github.com/souvik5151/Horn H1 H5.git
The transformation theory of the Appell F2(a, b1, b2; c1, c2; x, y) double hypergeometric functio... more The transformation theory of the Appell F2(a, b1, b2; c1, c2; x, y) double hypergeometric function is used to obtain a set of series representations of F2 which provide an efficient way to evaluate F2 for real values of its arguments x and y and generic complex values of its parameters a, b1, b2, c1 and c2 (i.e. in the nonlogarithmic case). This study rests on a classical approach where the usual double series representation of F2 and other double hypergeometric series that appear in the intermediate steps of the calculations are written as infinite sums of one variable hypergeometric series, such as the Gauss 2F1 or the 3F2, various linear transformations of the latter being then applied to derive known and new formulas. Using the three well-known Euler transformations of F2 on these results allows us to obtain a total of 44 series which form the basis of the Mathematica package AppellF2, dedicated to the evaluation of F2. A brief description of the package and of the numerical analysis that we have performed to test it are also presented.
The Method of Brackets (MoB) is a technique used to compute definite integrals, that has its orig... more The Method of Brackets (MoB) is a technique used to compute definite integrals, that has its origin in the negative dimensional integration method. It was originally proposed for the evaluation of Feynman integrals for which, when applicable, it gives the results in terms of combinations of (multiple) series. We focus here on some of the limitations of MoB and address them by studying the Mellin-Barnes (MB) representation technique. There has been significant process recently in the study of the latter due to the development of a new computational approach based on conic hulls (see Phys. Rev. Lett. 127, 151601 (2021)). The comparison between the two methods helps to understand the limitations of the MoB, in particular when termwise divergent series appear. As a consequence, the MB technique is found to be superior over MoB for two major reasons: 1. the selection of the sets of series that form a series representation for a given integral follows, in the MB approach, from specific intersections of conic hulls, which, in contrast to MoB, does not need any convergence analysis of the involved series, and 2. MB can be used to evaluate resonant (i.e. logarithmic) cases where MoB fails due to the appearance of termwise divergent series. Furthermore, we show that the recently added Rule 5 of MoB naturally emerges as a consequence of the residue theorem in the context of MB.
Capacity of entanglement (CoE), an information-theoretic measure of entanglement, defined as the ... more Capacity of entanglement (CoE), an information-theoretic measure of entanglement, defined as the variance of modular Hamiltonian, is known to capture the deviation from the maximal entanglement. We derive an exact expression for the average eigenstate CoE in fermionic Gaussian states as a finite series, valid for arbitrary bi-partition of the total system. Further, we consider the complex SYK2 model, and we show that the average CoE in the thermodynamic limit upper bounds the average CoE for fermionic Gaussian states. In this limit, the variance of the average CoE becomes independent of the system size. Moreover, when the subsystem size is half of the total system, the leading volume-law coefficient approaches a value of π 2 /8 − 1. We identify this as a distinguishing feature between integrable and quantum-chaotic systems. We confirm our analytical results by numerical computations.
We present the Olsson.wl Mathematica package which aims to find linear transformations for some c... more We present the Olsson.wl Mathematica package which aims to find linear transformations for some classes of multivariable hypergeometric functions. It is based on a well-known method developed by P. O. M. Olsson in J. Math. Phys. 5, 420 (1964) in order to derive the analytic continuations of the Appell $F_1$ double hypergeometric series from the linear transformations of the Gauss $_2F_1$ hypergeometric function. We provide a brief description of Olsson's method and demonstrate the commands of the package, along with examples. We also provide a companion package, called ROC2.wl and dedicated to the derivation of the regions of convergence of double hypergeometric series. This package can be used independently of Olsson.wl.
In semi-classical systems, the exponential growth of the out-of-time-order correlator (OTOC) is b... more In semi-classical systems, the exponential growth of the out-of-time-order correlator (OTOC) is believed to be the hallmark of quantum chaos. However, on several occasions, it has been argued that, even in integrable systems, OTOC can grow exponentially due to the presence of unstable saddle points in the phase space. In this work, we probe such an integrable system exhibiting saddle-dominated scrambling through Krylov complexity and the associated Lanczos coefficients. In the realm of the universal operator growth hypothesis, we demonstrate that the Lanczos coefficients follow the linear growth, which ensures the exponential behavior of Krylov complexity at early times. The linear growth arises entirely due to the saddle, which dominates other phase-space points even away from itself. Our results reveal that the exponential growth of Krylov complexity can be observed in integrable systems with saddle-dominated scrambling and thus need not be associated with the presence of chaos.
Capacity of entanglement (CoE), an information-theoretic measure of entanglement, defined as the ... more Capacity of entanglement (CoE), an information-theoretic measure of entanglement, defined as the variance of modular Hamiltonian, is known to capture the deviation from the maximal entanglement. We derive an exact expression for the average eigenstate CoE in fermionic Gaussian states as a finite series, valid for arbitrary bi-partition of the total system. Further, we consider the complex SYK2 model, and we show that the average CoE in the thermodynamic limit upper bounds the average CoE for fermionic Gaussian states. In this limit, the variance of the average CoE becomes independent of the system size. Moreover, when the subsystem size is half of the total system, the leading volume-law coefficient approaches a value of π/8 − 1. We identify this as a distinguishing feature between integrable and quantum-chaotic systems. We confirm our analytical results by numerical computations.
We use the method of brackets to evaluate quadratic and quartic type integrals. We recall the ope... more We use the method of brackets to evaluate quadratic and quartic type integrals. We recall the operational rules of the method and give examples to illustrate its working. The method is then used to evaluate the quadratic type integrals which occur in entries 3.251.1,3,4 in the table of integrals by Gradshteyn and Ryzhik and obtain closed form expressions in terms of hypergeometric functions. The method is further used to evaluate the quartic integrals, entry 2.161.5 and 6 in the table. We also present generalization of both types of integrals with closed form expression in terms of hypergeometric functions.
Uploads
Papers by tanay pathak