{"id":1013,"date":"2026-01-15T03:09:15","date_gmt":"2026-01-15T03:09:15","guid":{"rendered":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/?p=1013"},"modified":"2026-01-15T03:10:43","modified_gmt":"2026-01-15T03:10:43","slug":"tail-sensitive-convergence-guarantees-for-unadjusted-hamiltonian-monte-carlo","status":"publish","type":"post","link":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/tail-sensitive-convergence-guarantees-for-unadjusted-hamiltonian-monte-carlo\/","title":{"rendered":"Tail-Sensitive Convergence Guarantees for Unadjusted Hamiltonian Monte Carlo"},"content":{"rendered":"<p class=\"p1\">I\u2019m pleased to share a new paper with <strong>Siddharth Mitra<\/strong> and <strong>Andre Wibisono<\/strong> (both at Yale), titled:<\/p>\n<p class=\"p1\"><i><a href=\"https:\/\/arxiv.org\/abs\/2601.09019\">Tail-Sensitive KL and R\u00e9nyi Convergence of Unadjusted Hamiltonian Monte Carlo via One-Shot Couplings.<\/a> <\/i><\/p>\n<p class=\"p1\">This work develops a framework for tail-sensitive convergence analysis of unadjusted Hamiltonian Monte Carlo (HMC) in Kullback\u2013Leibler (KL) and R\u00e9nyi divergence. These divergences measure <span class=\"s1\"><b>relative density mismatch<\/b><\/span> and therefore govern Metropolis acceptance behavior: small errors in the tails\u2014where the target density is tiny\u2014can dominate acceptance probabilities, even when total variation distance is small.<\/p>\n<hr \/>\n<h3><b>Why tail-sensitive divergences matter<\/b><\/h3>\n<p class=\"p1\">Most theoretical analyses of Markov chain mixing time focus on metrics such as total variation distance or Wasserstein distance. While these are useful in many settings, they do not directly control relative density mismatch in the tails of the target distribution. For Metropolis-adjusted Markov chains, including Metropolis-adjusted HMC, acceptance probabilities depend on <i>ratios<\/i> of densities rather than additive differences in probability mass.<\/p>\n<p class=\"p1\">As a result, good performance requires controlling tail behavior. Divergences such as KL and R\u00e9nyi, which measure multiplicative discrepancies between densities, are therefore the natural metrics in this context.<\/p>\n<hr \/>\n<h3><b>What we prove<\/b><\/h3>\n<p class=\"p1\">Our main results provide a framework for lifting Wasserstein convergence guarantees (often easier to establish) into KL and R\u00e9nyi convergence bounds for unadjusted HMC. This is achieved by introducing and analyzing one-shot couplings that establish a regularization property of the unadjusted HMC transition kernel.<\/p>\n<p class=\"p1\">Concretely, we show that:<\/p>\n<ul>\n<li>\n<p class=\"p1\">A single step of unadjusted HMC has a smoothing effect in strong divergences, transforming even rough initial distributions into ones that are closer in KL or R\u00e9nyi divergence to the invariant distribution.<\/p>\n<\/li>\n<li>\n<p class=\"p1\">From this regularization, we derive mixing-time bounds and asymptotic bias estimates for unadjusted HMC in both KL and R\u00e9nyi divergence.<\/p>\n<\/li>\n<li>\n<p class=\"p1\">These bounds quantify relative density mismatch and yield principled warm-start guarantees for subsequent Metropolis-adjusted chains.<\/p>\n<\/li>\n<\/ul>\n<p class=\"p1\">A key ingredient in the analysis is the use of one-shot couplings. When combined with existing Wasserstein contractivity results (for example, under strong log-concavity assumptions), this yields end-to-end convergence estimates in divergences that directly govern algorithmic performance.<\/p>\n<hr \/>\n<h3><b>Practical implications<\/b><\/h3>\n<p class=\"p4\"><b>Warm starts for adjusted samplers.<\/b><b><\/b><\/p>\n<p class=\"p1\">A primary application of our results is the generation of rigorous warm starts for Metropolis-adjusted Markov chains. In practice, it is often difficult to initialize directly from a warm-start initial distribution. Unadjusted HMC can instead be run to reach a distribution that is close in KL or R\u00e9nyi divergence, after which a Metropolis-adjusted chain can be executed efficiently with good acceptance rates.<\/p>\n<p class=\"p4\"><b>Discretization bias.<\/b><\/p>\n<p class=\"p1\">By working explicitly in KL and R\u00e9nyi divergence, our results clarify how discretization bias (the discrepancy between the invariant distribution of an unadjusted sampler and the true target) affects strong notions of convergence. This provides a more complete picture of how unadjusted integrators behave in practice, beyond what is visible in weaker norms.<\/p>\n<hr \/>\n<h3><b>Broader relevance<\/b><\/h3>\n<p class=\"p1\">Tail-sensitive divergences also play an important role in:<\/p>\n<ul>\n<li>\n<p class=\"p1\">information contraction along Markov chains,<\/p>\n<\/li>\n<li>\n<p class=\"p1\">high-dimensional Bayesian inference,<\/p>\n<\/li>\n<li>\n<p class=\"p1\">R\u00e9nyi-divergence-based analyses in differential privacy.<\/p>\n<\/li>\n<\/ul>\n<hr \/>\n<h3><b>On the analysis<\/b><\/h3>\n<p class=\"p1\">Achieving these results required a delicate probabilistic analysis. It is fair to say that this project became a technical tour de force. Siddharth Mitra took the lead on the key probabilistic ideas, developing new arguments that substantially extend earlier total-variation convergence results for unadjusted HMC. The shift from additive to tail-sensitive, multiplicative control was challenging, but ultimately one of the most rewarding aspects of the project.<\/p>\n","protected":false},"excerpt":{"rendered":"<p>I\u2019m pleased to share a new paper with Siddharth Mitra and Andre Wibisono (both at Yale), titled: Tail-Sensitive KL and R\u00e9nyi Convergence of Unadjusted Hamiltonian Monte Carlo via One-Shot Couplings. &hellip; <a href=\"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/tail-sensitive-convergence-guarantees-for-unadjusted-hamiltonian-monte-carlo\/\" class=\"\">Read More<\/a><\/p>\n","protected":false},"author":2614,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"_acf_changed":false,"footnotes":""},"categories":[8],"tags":[],"class_list":["post-1013","post","type-post","status-publish","format-standard","hentry","category-preprints"],"acf":[],"yoast_head":"<!-- This site is optimized with the Yoast SEO plugin v23.5 - https:\/\/yoast.com\/wordpress\/plugins\/seo\/ -->\n<title>Tail-Sensitive Convergence Guarantees for Unadjusted Hamiltonian Monte Carlo - Nawaf Bou-Rabee<\/title>\n<meta name=\"robots\" content=\"index, follow, max-snippet:-1, max-image-preview:large, max-video-preview:-1\" \/>\n<link rel=\"canonical\" href=\"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/tail-sensitive-convergence-guarantees-for-unadjusted-hamiltonian-monte-carlo\/\" \/>\n<meta property=\"og:locale\" content=\"en_US\" \/>\n<meta property=\"og:type\" content=\"article\" \/>\n<meta property=\"og:title\" content=\"Tail-Sensitive Convergence Guarantees for Unadjusted Hamiltonian Monte Carlo - Nawaf Bou-Rabee\" \/>\n<meta property=\"og:description\" content=\"I\u2019m pleased to share a new paper with Siddharth Mitra and Andre Wibisono (both at Yale), titled: Tail-Sensitive KL and R\u00e9nyi Convergence of Unadjusted Hamiltonian Monte Carlo via One-Shot Couplings. &hellip; Read More\" \/>\n<meta property=\"og:url\" content=\"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/tail-sensitive-convergence-guarantees-for-unadjusted-hamiltonian-monte-carlo\/\" \/>\n<meta property=\"og:site_name\" content=\"Nawaf Bou-Rabee\" \/>\n<meta property=\"article:published_time\" content=\"2026-01-15T03:09:15+00:00\" \/>\n<meta property=\"article:modified_time\" content=\"2026-01-15T03:10:43+00:00\" \/>\n<meta name=\"author\" content=\"Nawaf Bou-Rabee\" \/>\n<meta name=\"twitter:card\" content=\"summary_large_image\" \/>\n<meta name=\"twitter:label1\" content=\"Written by\" \/>\n\t<meta name=\"twitter:data1\" content=\"Nawaf Bou-Rabee\" \/>\n\t<meta name=\"twitter:label2\" content=\"Est. reading time\" \/>\n\t<meta name=\"twitter:data2\" content=\"3 minutes\" \/>\n<script type=\"application\/ld+json\" class=\"yoast-schema-graph\">{\"@context\":\"https:\/\/schema.org\",\"@graph\":[{\"@type\":\"WebPage\",\"@id\":\"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/tail-sensitive-convergence-guarantees-for-unadjusted-hamiltonian-monte-carlo\/\",\"url\":\"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/tail-sensitive-convergence-guarantees-for-unadjusted-hamiltonian-monte-carlo\/\",\"name\":\"Tail-Sensitive Convergence Guarantees for Unadjusted Hamiltonian Monte Carlo - Nawaf Bou-Rabee\",\"isPartOf\":{\"@id\":\"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/#website\"},\"datePublished\":\"2026-01-15T03:09:15+00:00\",\"dateModified\":\"2026-01-15T03:10:43+00:00\",\"author\":{\"@id\":\"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/#\/schema\/person\/f9a7363ddc14a33eef77b520728dbde0\"},\"breadcrumb\":{\"@id\":\"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/tail-sensitive-convergence-guarantees-for-unadjusted-hamiltonian-monte-carlo\/#breadcrumb\"},\"inLanguage\":\"en-US\",\"potentialAction\":[{\"@type\":\"ReadAction\",\"target\":[\"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/tail-sensitive-convergence-guarantees-for-unadjusted-hamiltonian-monte-carlo\/\"]}]},{\"@type\":\"BreadcrumbList\",\"@id\":\"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/tail-sensitive-convergence-guarantees-for-unadjusted-hamiltonian-monte-carlo\/#breadcrumb\",\"itemListElement\":[{\"@type\":\"ListItem\",\"position\":1,\"name\":\"Home\",\"item\":\"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/\"},{\"@type\":\"ListItem\",\"position\":2,\"name\":\"Tail-Sensitive Convergence Guarantees for Unadjusted Hamiltonian Monte Carlo\"}]},{\"@type\":\"WebSite\",\"@id\":\"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/#website\",\"url\":\"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/\",\"name\":\"Nawaf Bou-Rabee\",\"description\":\"\",\"potentialAction\":[{\"@type\":\"SearchAction\",\"target\":{\"@type\":\"EntryPoint\",\"urlTemplate\":\"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/?s={search_term_string}\"},\"query-input\":{\"@type\":\"PropertyValueSpecification\",\"valueRequired\":true,\"valueName\":\"search_term_string\"}}],\"inLanguage\":\"en-US\"},{\"@type\":\"Person\",\"@id\":\"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/#\/schema\/person\/f9a7363ddc14a33eef77b520728dbde0\",\"name\":\"Nawaf Bou-Rabee\",\"image\":{\"@type\":\"ImageObject\",\"inLanguage\":\"en-US\",\"@id\":\"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/#\/schema\/person\/image\/\",\"url\":\"https:\/\/secure.gravatar.com\/avatar\/f3f33e8b5e6e2430997aea3dfec1a454?s=96&d=mm&r=g\",\"contentUrl\":\"https:\/\/secure.gravatar.com\/avatar\/f3f33e8b5e6e2430997aea3dfec1a454?s=96&d=mm&r=g\",\"caption\":\"Nawaf Bou-Rabee\"},\"description\":\"Math Professor at Rutgers\",\"sameAs\":[\"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/\"],\"url\":\"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/author\/nb361\/\"}]}<\/script>\n<!-- \/ Yoast SEO plugin. -->","yoast_head_json":{"title":"Tail-Sensitive Convergence Guarantees for Unadjusted Hamiltonian Monte Carlo - Nawaf Bou-Rabee","robots":{"index":"index","follow":"follow","max-snippet":"max-snippet:-1","max-image-preview":"max-image-preview:large","max-video-preview":"max-video-preview:-1"},"canonical":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/tail-sensitive-convergence-guarantees-for-unadjusted-hamiltonian-monte-carlo\/","og_locale":"en_US","og_type":"article","og_title":"Tail-Sensitive Convergence Guarantees for Unadjusted Hamiltonian Monte Carlo - Nawaf Bou-Rabee","og_description":"I\u2019m pleased to share a new paper with Siddharth Mitra and Andre Wibisono (both at Yale), titled: Tail-Sensitive KL and R\u00e9nyi Convergence of Unadjusted Hamiltonian Monte Carlo via One-Shot Couplings. &hellip; Read More","og_url":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/tail-sensitive-convergence-guarantees-for-unadjusted-hamiltonian-monte-carlo\/","og_site_name":"Nawaf Bou-Rabee","article_published_time":"2026-01-15T03:09:15+00:00","article_modified_time":"2026-01-15T03:10:43+00:00","author":"Nawaf Bou-Rabee","twitter_card":"summary_large_image","twitter_misc":{"Written by":"Nawaf Bou-Rabee","Est. reading time":"3 minutes"},"schema":{"@context":"https:\/\/schema.org","@graph":[{"@type":"WebPage","@id":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/tail-sensitive-convergence-guarantees-for-unadjusted-hamiltonian-monte-carlo\/","url":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/tail-sensitive-convergence-guarantees-for-unadjusted-hamiltonian-monte-carlo\/","name":"Tail-Sensitive Convergence Guarantees for Unadjusted Hamiltonian Monte Carlo - Nawaf Bou-Rabee","isPartOf":{"@id":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/#website"},"datePublished":"2026-01-15T03:09:15+00:00","dateModified":"2026-01-15T03:10:43+00:00","author":{"@id":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/#\/schema\/person\/f9a7363ddc14a33eef77b520728dbde0"},"breadcrumb":{"@id":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/tail-sensitive-convergence-guarantees-for-unadjusted-hamiltonian-monte-carlo\/#breadcrumb"},"inLanguage":"en-US","potentialAction":[{"@type":"ReadAction","target":["https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/tail-sensitive-convergence-guarantees-for-unadjusted-hamiltonian-monte-carlo\/"]}]},{"@type":"BreadcrumbList","@id":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/tail-sensitive-convergence-guarantees-for-unadjusted-hamiltonian-monte-carlo\/#breadcrumb","itemListElement":[{"@type":"ListItem","position":1,"name":"Home","item":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/"},{"@type":"ListItem","position":2,"name":"Tail-Sensitive Convergence Guarantees for Unadjusted Hamiltonian Monte Carlo"}]},{"@type":"WebSite","@id":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/#website","url":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/","name":"Nawaf Bou-Rabee","description":"","potentialAction":[{"@type":"SearchAction","target":{"@type":"EntryPoint","urlTemplate":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/?s={search_term_string}"},"query-input":{"@type":"PropertyValueSpecification","valueRequired":true,"valueName":"search_term_string"}}],"inLanguage":"en-US"},{"@type":"Person","@id":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/#\/schema\/person\/f9a7363ddc14a33eef77b520728dbde0","name":"Nawaf Bou-Rabee","image":{"@type":"ImageObject","inLanguage":"en-US","@id":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/#\/schema\/person\/image\/","url":"https:\/\/secure.gravatar.com\/avatar\/f3f33e8b5e6e2430997aea3dfec1a454?s=96&d=mm&r=g","contentUrl":"https:\/\/secure.gravatar.com\/avatar\/f3f33e8b5e6e2430997aea3dfec1a454?s=96&d=mm&r=g","caption":"Nawaf Bou-Rabee"},"description":"Math Professor at Rutgers","sameAs":["https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/"],"url":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/author\/nb361\/"}]}},"_links":{"self":[{"href":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/wp-json\/wp\/v2\/posts\/1013"}],"collection":[{"href":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/wp-json\/wp\/v2\/users\/2614"}],"replies":[{"embeddable":true,"href":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/wp-json\/wp\/v2\/comments?post=1013"}],"version-history":[{"count":1,"href":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/wp-json\/wp\/v2\/posts\/1013\/revisions"}],"predecessor-version":[{"id":1014,"href":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/wp-json\/wp\/v2\/posts\/1013\/revisions\/1014"}],"wp:attachment":[{"href":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/wp-json\/wp\/v2\/media?parent=1013"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/wp-json\/wp\/v2\/categories?post=1013"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/sites.rutgers.edu\/nawaf-bou-rabee\/wp-json\/wp\/v2\/tags?post=1013"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}