孪生素数（英语：twin prime），也称为孪生质数、双生质数，是指一对素数，它们之间相差2。例如3和5，5和7，11和13，10016957和10016959等等都是孪生素数。

关于孪生素数有著名的孪生素数猜想，即是否存在无穷多对孪生素数。这是数论中未解决的一个重要问题。哈代-李特尔伍德猜想是孪生素数猜想的一个增强形式，猜测孪生素数的分布与素数定理中描述的素数分布规律相类似。

与之相关的，两者相差为1的素数对只有 (2, 3)；两者相差为3的素数对只有 (2, 5)。

素数在自然数中的分布是不规则的。欧几里得在他的著作《几何原本》中首次证明了素数有无穷多个。十九世纪后，素数定理的证明给出了素数在自然数中大致的分布情况。根据素数定理，在前${\displaystyle N}$ 个自然数里孪生素数的个数。那么

其中的常数${\displaystyle C_{twin}}$是所谓的孪生素数常数：

其中的${\displaystyle p}$表示素数。

哈代和李特尔伍德的猜测实际上是存在已久的孪生素数猜想的加强版。孪生素数猜想是指“孪生素数有无穷多个”。这个猜想至今仍未被证明。然而，哈代和李特尔伍德的猜测并不是需要建立在孪生素数猜想成立的前提上。很多时候，对于无法证明的命题，数学家会尝试证明比它更强或更为广泛的命题，从而解决原来的命题。例如数学家安德鲁·怀尔斯就是证明了比费马最后猜想更广泛的命题，从而完成了费马最后猜想的证明。

2013年5月14日，《自然》杂志报道，数学家张益唐证明存在无穷多个素数对相差（上界）都小于7000万。论文已被《数学年刊》（Annals of Mathematics）接受。截至2014年10月9日 (2014-10-09), 素数对之差被缩小为${\displaystyle \le <!--\; \le \; -->246}$。

孪生素数猜想也可以用另一种形式表达：

1920年代，通过使用著名的筛理论（英语：Sieve theory），基于埃拉托斯特尼筛法的理论），挪威的维果·布朗（英语：Viggo Brun）证明了2能表示成两个最多有9个素数因子的数的差。这个结论已经有些近似于孪生素数猜想了。可以看到，只要将这个证明中的“最多有9个素数因子的数”改进到“最多有1个素数因子的数”，就可以证明孪生素数猜想了。利用同样的方法，布朗证明了所有偶数都能表达成两个最多有9个素数因子的数的和，也就是所谓的“9+9”。这个思路被不少数学家沿用，1966年陈景润利用筛法证明了“1+2”。基于陈景润的工作，也可以证明2有无限多种方法表示成一个素数和一个最多有两个素数因子的数的差。

布朗的另一个结论，是发现所有孪生素数的倒数之和收敛，即收敛到布朗常数${\displaystyle B_2}$：

${\displaystyle B_2}$的值大约在1.9与2之间。与之相对的，所有素数的倒数之和是发散的。由于孪生素数的倒数之和收敛，所以无法依此证明孪生素数有无限个。

布朗还发现了孪生素数数量的一个上限。他证明了：

也就是说，当${\displaystyle x}$足够大的时候，小于${\displaystyle x}$的孪生素数的数量比起小于${\displaystyle x}$的素数的数量是可以忽略不计的。1987年的一个结果改进了这个上限：

其中${\displaystyle C^{\prime }}$是一个常数。1998年上限中的7.1被改进为6.833。

孪生素数还必须满足一些必要的条件，比如：

统计分析所有小于${\displaystyle 4.35\times <!--\; \times \; -->{10}^{15}}$的孪生素数，可以得到小于${\displaystyle x}$的素数对的个数是${\displaystyle \frac{x\cdot <!--\; \cdot \; -->f(x)}{(\log <!--\; \; -->x{)}^2}}$。当${\displaystyle x}$较小时，${\displaystyle f(x)}$大约为 1.7， 当${\displaystyle x}$较大时大约为 1.3。这个值和${\displaystyle 2C_{twin}=1.3203236\dots <!--\; \dots \; -->}$相近。

孪生素数的概念可以扩展到多元组，即由多个间隔为2的素数构成的序列。由于三个相邻奇数总有一个能被3整除，不可能是素数，因此 (3, 5, 7) 是唯一的孪生素数三元组。而且由于更多元素构成的孪生素数多元组必定包含三元组的结构，因此多于三个元素的孪生素数多元组不存在。

以下的多项式时由维也纳大学数学系教授克里斯多夫·巴萨（Christoph Baxa）提出的，基于丢番图不定方程理论。

其中有二十六个不定量${\displaystyle a,b,c,\cdots <!--\; \cdots \; -->,y,z}$。当这二十六个变量取遍所有的自然数的时候，这个多项式的取值中正数的部分就会取遍所有孪生素数对${\displaystyle (p,p+2)}$中的${\displaystyle p}$。

意大利作家保罗·裘唐诺的小说《素数的孤独》即是以孪生素数现象，比喻故事中相爱的男女主角牵挛乖隔的处境。