Kostas Orginos
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Kostas Orginos has a total of 283 co-authors. Here are the top co-authors by number of shared publications:
- William Detmold
Massachusetts Institute of Technology
26 shared publications - André Walker-Loud
Lawrence Berkeley National Laboratory
25 shared publications - Assumpta Parreno
Departament de Física Quàntica i Astrofísica and Institut de Ciències del Cosmos, Universitat de Barcelona, Martí i Franquès 1, E08028, Spain
18 shared publications - Savvas Zafeiropoulos
Aix Marseille Univ Université de Toulon CNRS CPT Marseille France
18 shared publications Martin J. Savage
Department of Physics, University of Washington, Seattle, Washington 98195-1560, USA
16 shared publications- Carleton DeTar
University of Utah
15 shared publications
Latest Publications
- Open Access
- 23 April 2026
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2026 (4) , 182
Gaussian processes for inferring parton distributions
A bstract: The extraction of parton distribution functions (PDFs) from experimental or lattice QCD data is an ill-posed inverse problem, where regularization strongly impacts both systematic uncertainties and the reliability of the results. We study a framework based on Gaussian Process Regression (GPR) to reconstruct PDFs from lattice QCD matrix elements. Within a Bayesian framework, Gaussian processes serve as flexible priors that encode uncertainties, correlations, and constraints without imposing rigid functional forms. We investigate a wide range of kernel choices, mean functions, and hyperparameter treatments. We quantify information gained from the data using the Kullback-Leibler divergence. Synthetic data tests demonstrate the consistency and robustness of the method. Our study establishes GPR as a systematic and non-parametric approach to PDF reconstruction, offering controlled uncertainty estimates and reduced model bias in lattice QCD analyses.
- 23 April 2026
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2026 (4) , 182
Open AccessGaussian processes for inferring parton distributions
A bstract: The extraction of parton distribution functions (PDFs) from experimental or lattice QCD data is an ill-posed inverse problem, where regularization strongly impacts both systematic uncertainties and the reliability of the results. We study a framework based on Gaussian Process Regression (GPR) to reconstruct PDFs from lattice QCD matrix elements. Within a Bayesian framework, Gaussian processes serve as flexible priors that encode uncertainties, correlations, and constraints without imposing rigid functional forms. We investigate a wide range of kernel choices, mean functions, and hyperparameter treatments. We quantify information gained from the data using the Kullback-Leibler divergence. Synthetic data tests demonstrate the consistency and robustness of the method. Our study establishes GPR as a systematic and non-parametric approach to PDF reconstruction, offering controlled uncertainty estimates and reduced model bias in lattice QCD analyses.
- Open Access
Accessing the Gluon Momentum Fraction of Nucleons through the Gradient Flow
We calculate the gluon momentum fraction of the nucleon using lattice quantum chromodynamics (QCD), with a nonperturbative renormalization technique based on the gradient flow. The gluon momentum fraction is determined on a single Wilson-clover ensemble using Nf = 2+1 flavors with pion mass 358 MeV and lattice spacing 0.094 fm. We employ the variational method to reduce excited-state contamination and apply the distillation framework to ensure a large operator basis. To reduce systematic uncertainties, we apply Bayesian model averaging to all fit procedures. We apply matching coefficients to the flow-time dependent lattice results to recover the gluon momentum fraction in the MS-scheme at 2 GeV. Our final result is _g(μ= 2 GeV) = 0.482(35), where we quote only statistical uncertainties.
Open AccessAccessing the Gluon Momentum Fraction of Nucleons through the Gradient Flow
We calculate the gluon momentum fraction of the nucleon using lattice quantum chromodynamics (QCD), with a nonperturbative renormalization technique based on the gradient flow. The gluon momentum fraction is determined on a single Wilson-clover ensemble using Nf = 2+1 flavors with pion mass 358 MeV and lattice spacing 0.094 fm. We employ the variational method to reduce excited-state contamination and apply the distillation framework to ensure a large operator basis. To reduce systematic uncertainties, we apply Bayesian model averaging to all fit procedures. We apply matching coefficients to the flow-time dependent lattice results to recover the gluon momentum fraction in the MS-scheme at 2 GeV. Our final result is _g(μ= 2 GeV) = 0.482(35), where we quote only statistical uncertainties.
- Open Access
The Spectrum and Scale Setting on 2+1-flavor NME Lattices
This paper describes the thirteen ensembles, named NME, generated with 2+1-flavor Wilson-clover fermions by the JLab/W\&M/LANL/MIT/Marseille collaborations, and presents an analysis of the meson and baryon spectrum, decay constants $f_π$ and $f_K$, flow scales $t_0$ and $w_0$, and time histories of the $Θ$ and Weinberg operators under gradient flow. Using these quantities, the physical point values of the two flow scales, ${t_0^{\rm Phy}}$ and ${w_0^{\rm Phy}}$, and the ratio $\mathop{f_K / f_π}^{\rm Phy}$ are determined. The masses of the octet and decuplet baryons are analyzed using both the next-to-leading order (NLO) and the next-next-to-leading order (NNLO) ansatz from heavy baryon chiral perturbation theory (HB$χ$PT). The NNLO fit to the octet baryons, $M_N$, $M_Σ$, $M_Λ$ and $M_Ξ$, is preferred while the corresponding fits to the decuplet Omega mass, $M_Ω$, are not distinguished. We also present a study of the autocorrelations in the data and show that there is no evidence, even at large flow time, of the freezing of the topological charge or the Weinberg three-gluon operator.Open AccessThe Spectrum and Scale Setting on 2+1-flavor NME Lattices
This paper describes the thirteen ensembles, named NME, generated with 2+1-flavor Wilson-clover fermions by the JLab/W\&M/LANL/MIT/Marseille collaborations, and presents an analysis of the meson and baryon spectrum, decay constants $f_π$ and $f_K$, flow scales $t_0$ and $w_0$, and time histories of the $Θ$ and Weinberg operators under gradient flow. Using these quantities, the physical point values of the two flow scales, ${t_0^{\rm Phy}}$ and ${w_0^{\rm Phy}}$, and the ratio $\mathop{f_K / f_π}^{\rm Phy}$ are determined. The masses of the octet and decuplet baryons are analyzed using both the next-to-leading order (NLO) and the next-next-to-leading order (NNLO) ansatz from heavy baryon chiral perturbation theory (HB$χ$PT). The NNLO fit to the octet baryons, $M_N$, $M_Σ$, $M_Λ$ and $M_Ξ$, is preferred while the corresponding fits to the decuplet Omega mass, $M_Ω$, are not distinguished. We also present a study of the autocorrelations in the data and show that there is no evidence, even at large flow time, of the freezing of the topological charge or the Weinberg three-gluon operator. - Open Access
- 3 November 2025
- Published by American Physical Society (APS) in Physical Review Letters
- Vol. 135 (19) , 191901
Window Observables for Benchmarking Parton Distribution Functions
Global analysis of collider and fixed-target experimental data and calculations from lattice quantum chromodynamics (QCD) are used to gain complementary information on the structure of hadrons. We propose novel “window observables” that allow for higher precision cross-validation between the different approaches, a critical step for studies that wish to combine the datasets. Global analyses are limited by the kinematic regions accessible to experiment, particularly in a range of Bjorken- x , and lattice QCD calculations also have limitations requiring extrapolations to obtain the parton distributions. We provide two different window observables that can be defined within a region of x where extrapolations and interpolations in global analyses remain reliable and where lattice QCD results retain sensitivity and precision.
- 3 November 2025
- Published by American Physical Society (APS) in Physical Review Letters
- Vol. 135 (19) , 191901
Open AccessWindow Observables for Benchmarking Parton Distribution Functions
Global analysis of collider and fixed-target experimental data and calculations from lattice quantum chromodynamics (QCD) are used to gain complementary information on the structure of hadrons. We propose novel “window observables” that allow for higher precision cross-validation between the different approaches, a critical step for studies that wish to combine the datasets. Global analyses are limited by the kinematic regions accessible to experiment, particularly in a range of Bjorken- x , and lattice QCD calculations also have limitations requiring extrapolations to obtain the parton distributions. We provide two different window observables that can be defined within a region of x where extrapolations and interpolations in global analyses remain reliable and where lattice QCD results retain sensitivity and precision.
- Open Access
- 28 February 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 111 (3)
Simple nonparametric reconstruction of parton distributions from limited Fourier information
Some calculations of parton distributions from first principles only give access to a limited range of Fourier modes of the function to reconstruct. We present a physically motivated procedure to regularize the inverse integral problem using a Gaussian process as a Bayesian prior. We propose to fix the hyperparameters of the prior in a meaningful physical fashion, offering a simple implementation, great numerical efficiency, and allowing us to understand and keep control easily of the uncertainty of the reconstruction.
- 28 February 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 111 (3)
Open AccessSimple nonparametric reconstruction of parton distributions from limited Fourier information
Some calculations of parton distributions from first principles only give access to a limited range of Fourier modes of the function to reconstruct. We present a physically motivated procedure to regularize the inverse integral problem using a Gaussian process as a Bayesian prior. We propose to fix the hyperparameters of the prior in a meaningful physical fashion, offering a simple implementation, great numerical efficiency, and allowing us to understand and keep control easily of the uncertainty of the reconstruction.
- Open Access
- 21 August 2024
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2024 (8) , 1-60
Towards unpolarized GPDs from pseudo-distributions
A bstract: We present an exploration of the unpolarized isovector proton generalized parton distributions (GPDs) H u − d ( x , ξ , t ) and E u − d ( x , ξ , t ) in the pseudo-distribution formalism using distillation. Taking advantage of the large kinematic coverage made possible by this approach, we present results on the moments of GPDs up to the order x 3 — including their skewness dependence — at a pion mass m π = 358 MeV and a lattice spacing a = 0 . 094 fm.
- 21 August 2024
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2024 (8) , 1-60
Open AccessTowards unpolarized GPDs from pseudo-distributions
A bstract: We present an exploration of the unpolarized isovector proton generalized parton distributions (GPDs) H u − d ( x , ξ , t ) and E u − d ( x , ξ , t ) in the pseudo-distribution formalism using distillation. Taking advantage of the large kinematic coverage made possible by this approach, we present results on the moments of GPDs up to the order x 3 — including their skewness dependence — at a pion mass m π = 358 MeV and a lattice spacing a = 0 . 094 fm.
- Open Access
- 2 May 2024
- Published by Sissa Medialab Srl
- Vol. 453, 300
Extracting the Pion Distribution Amplitude from Lattice QCD through Pseudo-Distributions
The Light-Cone Distribution Amplitude (LCDA) encodes the non-perturbative information of the leading Fock component of the hadron wave function, therefore required for processes including exclusive hadron production. As the Pseudo-Nambu-Goldstone boson of QCD, the nonperturbative structure of the pion is of particular interest. Progress on the Lattice QCD calculation of the pion LCDA on ${\cal O}(a)$-improved Wilson fermion ensembles at several lattice spacings is presented. Excited-state systematics are taken into account within a Bayesian Model Averaging framework. A Renormalization-Group-Invariant (RGI) ratio of matrix elements is formed for further extraction of the pion LCDA.
- 2 May 2024
- Published by Sissa Medialab Srl
- Vol. 453, 300
Open AccessExtracting the Pion Distribution Amplitude from Lattice QCD through Pseudo-Distributions
The Light-Cone Distribution Amplitude (LCDA) encodes the non-perturbative information of the leading Fock component of the hadron wave function, therefore required for processes including exclusive hadron production. As the Pseudo-Nambu-Goldstone boson of QCD, the nonperturbative structure of the pion is of particular interest. Progress on the Lattice QCD calculation of the pion LCDA on ${\cal O}(a)$-improved Wilson fermion ensembles at several lattice spacings is presented. Excited-state systematics are taken into account within a Bayesian Model Averaging framework. A Renormalization-Group-Invariant (RGI) ratio of matrix elements is formed for further extraction of the pion LCDA.
- Open Access
- 11 April 2024
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2024 (4) , 1-57
Evolution of parton distribution functions in the short-distance factorization scheme
Lattice QCD offers the possibility of computing parton distributions from first principles, although not in the usual $$ \overline{MS} $$ factorization scheme. Calculations are therefore matched to $$ \overline{MS} $$ using a perturbative procedure which is the source of significant uncertainty within the currently accessible kinematics. We present the possibility of computing the z 2 evolution of non-singlet pseudo-parton distribution functions within the short factorization scheme in a numerically improvable way. The goal is to have tools to evolve a calculation to a scale where perturbative uncertainties are less pronounced. We compare a numerical extraction of the evolution operator from lattice data to the computation of z 2 dependence in perturbation theory. Finally, we discuss how this numerical work may be extended to address the two-scale problem that arises when the Ioffe time range must be made large to extend the reach of the calculation of the pseudo-PDF to smaller values of the momentum fraction.- 11 April 2024
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2024 (4) , 1-57
Open AccessEvolution of parton distribution functions in the short-distance factorization scheme
Lattice QCD offers the possibility of computing parton distributions from first principles, although not in the usual $$ \overline{MS} $$ factorization scheme. Calculations are therefore matched to $$ \overline{MS} $$ using a perturbative procedure which is the source of significant uncertainty within the currently accessible kinematics. We present the possibility of computing the z 2 evolution of non-singlet pseudo-parton distribution functions within the short factorization scheme in a numerically improvable way. The goal is to have tools to evolve a calculation to a scale where perturbative uncertainties are less pronounced. We compare a numerical extraction of the evolution operator from lattice data to the computation of z 2 dependence in perturbation theory. Finally, we discuss how this numerical work may be extended to address the two-scale problem that arises when the Ioffe time range must be made large to extend the reach of the calculation of the pseudo-PDF to smaller values of the momentum fraction. - Open Access
Trivializing Flow in 2D O(3) sigma model
The two dimensional O(3) sigma model, just as quantum chromodynamics, is an asymptotically free theory with a mass gap. Therefore, it is an interesting and simple toy model to investigate algorithms for Markov Chain Monte Carlo simulations of quantum chromodynamics. In this talk, we discuss the construction of a trivializing map, a field transformation from a given theory to a trivial one, through a suitably chosen gradient flow. An analytic solution for the generating functional of this trivializing flow has been obtained by a perturbative expansion in the flow time. Utilizing this solution allows for new approaches to be considered when proposing updates for a Markov Chain Monte Carlo algorithm.
Open AccessTrivializing Flow in 2D O(3) sigma model
The two dimensional O(3) sigma model, just as quantum chromodynamics, is an asymptotically free theory with a mass gap. Therefore, it is an interesting and simple toy model to investigate algorithms for Markov Chain Monte Carlo simulations of quantum chromodynamics. In this talk, we discuss the construction of a trivializing map, a field transformation from a given theory to a trivial one, through a suitably chosen gradient flow. An analytic solution for the generating functional of this trivializing flow has been obtained by a perturbative expansion in the flow time. Utilizing this solution allows for new approaches to be considered when proposing updates for a Markov Chain Monte Carlo algorithm.
- Open Access
Extracting the Pion Distribution Amplitude from Lattice QCD through Pseudo-Distributions
The Light-Cone Distribution Amplitude (LCDA) encodes the non-perturbative information of the leading Fock component of the hadron wave function, therefore required for processes including exclusive hadron production. As the pseudo-Nambu-Goldstone boson of QCD, the nonperturbative structure of the pion is of particular interest. Progress on the Lattice QCD calculation of the pion LCDA on ${\cal O}(a)$-improved Wilson fermion ensembles at several lattice spacings is presented. Excited-state systematics are taken into account within a Bayesian Model Averaging framework. A Renormalization-Group-Invariant (RGI) ratio of matrix elements is formed for further extraction of the pion LCDA.
Open AccessExtracting the Pion Distribution Amplitude from Lattice QCD through Pseudo-Distributions
The Light-Cone Distribution Amplitude (LCDA) encodes the non-perturbative information of the leading Fock component of the hadron wave function, therefore required for processes including exclusive hadron production. As the pseudo-Nambu-Goldstone boson of QCD, the nonperturbative structure of the pion is of particular interest. Progress on the Lattice QCD calculation of the pion LCDA on ${\cal O}(a)$-improved Wilson fermion ensembles at several lattice spacings is presented. Excited-state systematics are taken into account within a Bayesian Model Averaging framework. A Renormalization-Group-Invariant (RGI) ratio of matrix elements is formed for further extraction of the pion LCDA.