"It feels different when blindfolded": Developing social empathy through inclusive designs in STEM

Tang Wee Teo; Aik Ling Tan; Ban Heng Choy; Tang Wee Teo; Aik Ling Tan; Ban Heng Choy

doi:10.3934/steme.2024002

STEM Education

2024, Volume 4, Issue 1: 20-34. doi: 10.3934/steme.2024002

Previous Article Next Article

Research article

"It feels different when blindfolded": Developing social empathy through inclusive designs in STEM

1.
Natural Sciences & Science Education, National Institute of Education, Nanyang Technological University, 1 Nanyang Walk, Singapore; tangwee.teo@nie.edu.sg, aikling.tan@nie.edu.sg
2.
Mathematics and Mathematics Education, National Institute of Education, Nanyang Technological University, 1 Nanyang Walk, Singapore; banheng.choy@nie.edu.sg

Academic Editor: Wei Li

Received: 23 October 2023 Revised: 02 January 2024 Published: 23 January 2024

Narratives about STEM (science, technology, engineering, and mathematics) education are strongly connected with conversations about developing learners' humanistic knowledge and their ability to listen with understanding and empathy. This is challenging because learners need to find resonance through first-hand contextual experiences with the issues at hand. In this paper, we describe and discuss an activity that was enacted to actively engage 74 teachers from Thailand in constructing a prototype cutting device for the blind to prepare food themselves. This activity underscores important considerations for inclusive design and offers affordances for teachers to develop their learners' inclusive mindsets. Findings were generated from voice recordings of reflections and written reflections collected after the activity. We highlight the importance of creating opportunities for learning to listen and resonate with others' experiences and argue that such STEM activities can offer a platform for learners to develop humanistic qualities such as social empathy.

Keywords:

Citation: Tang Wee Teo, Aik Ling Tan, Ban Heng Choy. 'It feels different when blindfolded': Developing social empathy through inclusive designs in STEM[J]. STEM Education, 2024, 4(1): 20-34. doi: 10.3934/steme.2024002

Related Papers:

[1]	Kwang-Wu Chen, Fu-Yao Yang . Infinite series involving harmonic numbers and reciprocal of binomial coefficients. AIMS Mathematics, 2024, 9(7): 16885-16900. doi: 10.3934/math.2024820
[2]	Zhi-Hong Sun . Supercongruences involving Apéry-like numbers and binomial coefficients. AIMS Mathematics, 2022, 7(2): 2729-2781. doi: 10.3934/math.2022153
[3]	Shelby Kilmer, Songfeng Zheng . A generalized alternating harmonic series. AIMS Mathematics, 2021, 6(12): 13480-13487. doi: 10.3934/math.2021781
[4]	Wenbin Zhang, Yong Zhang, Jiachen Wu . A Ramanujan-type supercongruence related to harmonic numbers. AIMS Mathematics, 2025, 10(5): 11444-11464. doi: 10.3934/math.2025521
[5]	Xuqiong Luo . A D-N alternating algorithm for exterior 3-D problem with ellipsoidal artificial boundary. AIMS Mathematics, 2022, 7(1): 455-466. doi: 10.3934/math.2022029
[6]	Jiaye Lin . Evaluations of some Euler-type series via powers of the arcsin function. AIMS Mathematics, 2025, 10(4): 8116-8130. doi: 10.3934/math.2025372
[7]	Ling Zhu . Asymptotic expansion of a finite sum involving harmonic numbers. AIMS Mathematics, 2021, 6(3): 2756-2763. doi: 10.3934/math.2021168
[8]	Andrea Ratto . Higher order energy functionals and the Chen-Maeta conjecture. AIMS Mathematics, 2020, 5(2): 1089-1104. doi: 10.3934/math.2020076
[9]	Aleksa Srdanov . Invariants in partition classes. AIMS Mathematics, 2020, 5(6): 6233-6243. doi: 10.3934/math.2020401
[10]	Kholood M. Alsager, Sheza M. El-Deeb, Ala Amourah, Jongsuk Ro . Some results for the family of holomorphic functions associated with the Babalola operator and combination binomial series. AIMS Mathematics, 2024, 9(10): 29370-29385. doi: 10.3934/math.20241423

Abstract

1. Introduction and outline

Let $\mathbb{N}$ be the set of natural numbers with $\mathbb{N}_0 = \mathbb{N}\cup\{0\}$ . For $n \in \mathbb{N}_0$ , the following harmonic-like numbers are given by partial sums:

● The classical harmonic number and its quadratic counterpart

$\mathbf{H}_{{n}} = \sum\limits_{k = 0}^{n-1}\frac1{k+1} \quad\text{and}\quad \mathbf{H}_{{n}}^{ \langle{{{2}}}\rangle} = \sum\limits_{k = 0}^{n-1}\frac1{(k+1)^{2}}.$

● Skew harmonic number and its quadratic counterpart

$\mathbf{O}_{{n}} = \sum\limits_{k = 0}^{n-1}\frac1{2k+1} \quad\text{and}\quad \mathbf{O}_{{n}}^{ \langle{{{2}}}\rangle} = \sum\limits_{k = 0}^{n-1}\frac1{(2k+1)^{2}}.$

● Alternating harmonic number and its quadratic counterpart

$\bar{ \mathbf{H}}_{{n}} = \sum\limits_{k = 0}^{n-1}\frac{(-1)^k}{k+1} \quad\text{and}\quad \bar{ \mathbf{H}}_{{n}}^{ \langle{{{2}}}\rangle} = \sum\limits_{k = 0}^{n-1}\frac{(-1)^k}{(k+1)^{2}}.$

Among them, there are the following simple but useful relations:

$\begin{align} & \mathbf{H}_{{2n}} = \mathbf{O}_{{n}}+\frac{ \mathbf{H}_{{n}}}{2}, &\qquad& \mathbf{H}_{{2n}}^{ \langle{{{2}}}\rangle} = \mathbf{O}_{{n}}^{ \langle{{{2}}}\rangle}+\frac{ \mathbf{H}_{{n}}^{ \langle{{{2}}}\rangle}}{4};\\ &\bar{ \mathbf{H}}_{{2n}} = \mathbf{O}_{{n}}-\frac{ \mathbf{H}_{{n}}}{2}, &&\bar{ \mathbf{H}}_{{2n}}^{ \langle{{{2}}}\rangle} = \mathbf{O}_{{n}}^{ \langle{{{2}}}\rangle}-\frac{ \mathbf{H}_{{n}}^{ \langle{{{2}}}\rangle}}{4}. \end{align}$

Harmonic numbers and their variants have been studied since the distant past and are involved in diverse fields (cf. ^{[1,5,6,7,11,17,18]}) such as analysis of algorithms in computer science, various topics of number theory, and combinatorial analysis. In this paper, we aim at presenting several algebraic identities about infinite series involving harmonic-like numbers, the binomial coefficient $\binom{3n}{n}$ in numerators, and a free variable " $y$ ".

For $n\in\mathbb{{N}}_0$ and an indeterminate $x$ , the shifted factorials (cf. [3, §1.1])) are usually defined by

$(x)_0 = 1 \quad\text{and}\quad (x)_n = x(x+1)\cdots(x+n-1) \quad\text{for}\quad n\in\mathbb{{N}}.$

Denote by $[x^m]\phi(x)$ the coefficient of $x^m$ in the formal power series $\phi(x)$ . Then we can express harmonic-like numbers (see ^[11,18]) by extracting the initial coefficients of $x$ from the quotients of shifted factorials as below:

$\begin{align} &[x]\frac{(1+x)_n}{n!} = \mathbf{H}_{{n}}, &\qquad&[x^2]\frac{(1+x)_n}{n!} = \frac{ \mathbf{H}_{{n}}^{{2}}- \mathbf{H}_{{n}}^{ \langle{{{2}}}\rangle}}{2}, \\ &[x]\frac{n!}{(1-x)_n} = \mathbf{H}_{{n}}, &&[x^2]\frac{n!}{(1-x)_n} = \frac{ \mathbf{H}_{{n}}^{{2}}+ \mathbf{H}_{{n}}^{ \langle{{{2}}}\rangle}}{2}, \\ &[y]\frac{(\frac12+y)_n}{(\frac12)_n} = 2 \mathbf{O}_{{n}}, &&[y^2]\frac{(\frac12+y)_n}{(\frac12)_n} = 2( \mathbf{O}_{{n}}^{{2}}- \mathbf{O}_{{n}}^{ \langle{{{2}}}\rangle}), \\ &[y]\frac{(\frac12)_n}{(\frac12-y)_n} = 2 \mathbf{O}_{{n}}, &&[y^2]\frac{(\frac12)_n}{(\frac12-y)_n} = 2( \mathbf{O}_{{n}}^{{2}}+ \mathbf{O}_{{n}}^{ \langle{{{2}}}\rangle}). \end{align}$

There exist many infinite series identities involving central binomial coefficients (see ^{[12,14,19,25]}). However, the closed form evaluations for the known series with $\binom{3n}{n}$ are rare. In a recent paper, Sun ^[23] reported several remarkable infinite series identities, detected by numerical experimentation. Among them, the following three (see [23, Eqs (2.12), (2.13), and (2.15)]) seem rather challenging:

$\begin{align} & \sum\limits_{n = 1}^{\infty}\Big(\frac{2}{27}\Big)^n \binom{3n}{n} \mathbf{H}_{{2n}} = \frac{1+\sqrt3}2\Big\{\sqrt3\ln2-3\ln(3-\sqrt3)\Big\}, \end{align}$

(1.1)

$\begin{align} & \sum\limits_{n = 1}^{\infty}\Big(\frac2{27}\Big)^n \binom{3n}{n} \mathbf{H}_{{3n}} = \frac{1+\sqrt3}4\ln\Big(\frac{4(7+4\sqrt3)}3\Big)-\sqrt3\ln2, \end{align}$

(1.2)

$\begin{align} & \sum\limits_{n = 1}^{\infty}\Big(\frac{3+\sqrt5}{54}\Big)^n \binom{3n}{n}\Big\{ \mathbf{H}_{{3n}}- \mathbf{H}_{{2n}}\Big\} = \frac{1+\sqrt5}{2}\ln\Big(\frac{9-3\sqrt5}{2}\Big). \end{align}$

(1.3)

When making efforts to confirm these three formulae, two substantial novelties emerged. First, not only the harmonic numbers $\mathbf{H}_{{3n}}$ and $\mathbf{H}_{{2n}}$ are separable in the third series, but also independent series involving $\mathbf{H}_{{n}}$ and $\mathbf{O}_{{n}}$ can be evaluated explicitly in closed form (cf. Corollaries 6 and 10). Another surprising novelty (see Theorems 1 and 5) lies in the fact that these series are particular cases of more general algebraic identities containing a free variable " $y$ ", corresponding to the convergent rates of these numerical series. To our knowledge, most of the infinite series identities presented in this paper have not appeared previously in the literature.

The rest of the paper is organized as follows. As preliminaries, we shall evaluate, in closed form, three series containing binomial coefficient $\binom{3n}{n}$ , harmonic number $\mathbf{H}_{{n}}$ and a free variable " $y$ " by integrating Lambert's series. Then in Sections 3–5, three classes of infinite series with harmonic-like numbers will be examined by applying the "coefficient extraction" method (see ^[8,13,24]) to cubic transformations for the $_3{F_2}$ -series, that result in several algebraic identities about these series and confirm or refine, as special numerical examples, three conjectured identities made recently by Sun ^[23]. Finally, the paper will end with Section 6, where concluding remarks and further problems are provided.

2. Lambert series weighted by harmonic numbers

In classical analysis and enumerative combinatorics, Lambert's binomial series are well-known (see Riordan [, §4.5] and ^[9,10,16]), which can be reproduced as follows: Let $T$ and $\tau$ be the two variables related by the equation $\tau = T/(1+T)^{ \beta}$ . Then

$\begin{equation} { (1+T)^{ \alpha} = \sum\limits_{k = 0}^{\infty} \frac{ \alpha}{ \alpha+k \beta} \binom{ \alpha+k \beta}{k}\tau^k \quad\text{and}\quad \frac{(1+T)^{ \alpha+1}}{1+T-T \beta} = \sum\limits_{k = 0}^{\infty} \binom{ \alpha+k \beta}{k}\tau^k.} \end{equation}$

(2.1)

As preliminaries, we are going to evaluate, in closed form, three series containing ternary coefficient $\binom{3n}{n}$ and a free variable " $y$ " by integrating across the above two equations, that will be crucial for demonstrating infinite series identities involving harmonic numbers in the subsequent sections.

2.1. Series with $\mathbf{H}_{{n}}$

For " $\alpha = 0$ and $\beta = 3$ ", performing the replacement $\tau\to{xy}/{(1+y)^3}$ in Lambert's second series, we can state the resulting equality as

$\begin{equation} { \phi(T): = \frac{1+T}{1-2T} = \sum\limits_{n = 0}^{\infty} \binom{3n}{n}\frac{(xy)^n}{(1+y)^{3n}}: \qquad {x = \frac{T(1+y)^3}{y(1+T)^3}}.} \end{equation}$

(2.2)

Since the rightmost function is monotone around $T = 0$ , it determines an implicit inverse function $T = T(x)$ in the neighborhood of $x = 0$ with particular values " $T(0) = 0$ and $T(1) = y$ ".

By expressing the harmonic number $\mathbf{H}_{{n}}$ in terms of the integral

$\mathbf{H}_{{n}} = \int_0^1\frac{1-x^n}{1-x}dx,$

we can insert $\mathbf{H}_{{n}}$ in the series (2.2) as follows:

$\begin{align*} \sum\limits_{n = 0}^{\infty} \binom{3n}{n}\frac{ \mathbf{H}_{{n}}\; y^n}{(1+y)^{3n}} & = \int_0^1\frac{\phi(T(1))-\phi(T(x))}{1-x}dx\\ & = \frac{(1+y)^3}{y}\int_0^y\frac{1-2T}{(1+T)^4} \frac{\frac{1+y}{1-2y}-\frac{1+T}{1-2T}}{1-\frac{T(1+y)^3}{y(1+T)^3}}dT\\ & = (1+y)^3\int_0^y\frac{1-\frac{(1-2T)(1+y)}{(1+T)(1-2y)}}{T(1+y)^3-y(1+T)^3}dT, \end{align*}$

where the change of variable $x = \dfrac{T(1+y)^3}{y(1+T)^3}$ has been made so that the integral becomes explicit.

For the sake of brevity, the notation " $\rho$ " will be fixed, throughout the paper, for the algebraic number

$\begin{equation} {\rho: = \frac32(\sqrt2+1)^{1/3}-\frac32(\sqrt2-1)^{1/3}-1\approx-0.105892543.} \end{equation}$

(2.3)

The convergence region of the preceding series is determined by

$-1 < \frac{27y}{4(1+y)^3} < 1, \quad \rho < y < \frac12.$

This can be illustrated graphically as below:

Figure 1. Range of convergence.

DownLoad: Full-Size Img PowerPoint

By decomposing the integrand into partial fractions

$\begin{align*} \frac{1-\frac{(1-2T)(1+y)}{(1+T)(1-2y)}}{T(1+y)^3-y(1+T)^3} & = \frac{3}{(1+T) (1+y)^2 (1-2 y)} -\frac3{2(1+y)^2(1-2y)\sqrt{4+y}}\\ &\times\Bigg\{\frac{\sqrt{4+y}+\sqrt{y}} {T+\frac{\sqrt{y}(3+y)-(1+y) \sqrt{4+y}}{2\sqrt{y}}} +\frac{\sqrt{4+y}-\sqrt{y}} {T+\frac{\sqrt{y}(3+y)+(1+y)\sqrt{4+y}}{2\sqrt{y}}}\Bigg\}, \end{align*}$

we can explicitly compute the integral

$\begin{align*} \int_0^y\frac{1-\frac{(1-2T)(1+y)}{(1+T)(1-2y)}}{T(1+y)^3-y(1+T)^3}dT & = \frac{3\ln(1+y)}{(1+y)^2(1-2y)}\\ &\; +\frac{3(\sqrt{y}-\sqrt{4+y})}{2(1+y)^2(1-2y)\sqrt{4+y}} \ln\Big(\frac{(1+y)(2-2y-y^2+y\sqrt{y^2+4y})}{2}\Big)\\ &\; -\frac{3(\sqrt{y}+\sqrt{4+y})}{2(1+y)^2(1-2y)\sqrt{4+y}} \ln\Big(\frac{(1+y)(2-2y-y^2-y\sqrt{y^2+4y})}{2}\Big). \end{align*}$

By cancelling the three terms containing $\ln(1+y)$ and then making substitution, we find the closed formula for the series as in the theorem below:

Theorem 1 ( $\rho < y < 1/2$ ).

$\begin{align*} \sum\limits_{n = 1}^{\infty} \binom{3n}{n}\frac{ \mathbf{H}_{{n}}\; y^n}{(1+y)^{3n}} & = (1+y)^3\int_0^y\frac{1-\frac{(1-2T)(1+y)}{(1+T)(1-2y)}}{T(1+y)^3-y(1+T)^3}dT\\ & = \frac{3(1+y)\ln(1-2y)}{2(2y-1)} +\frac{3(1+y)\sqrt{y}}{2(1-2y)\sqrt{4+y}} \ln\Big(\frac{2-y(2+y-\sqrt{y^2+4y})}{2-y(2+y+\sqrt{y^2+4y})}\Big). \end{align*}$

2.2. Series with $\mathbf{H}_{{n}}/(1+2n)$

For " $\alpha = 1$ and $\beta = 3$ ", performing the replacement $\tau\to{xy}/{(1+y)^3}$ , we can rewrite Lambert's first series as

$\begin{equation} { \psi(T): = 1+T = \sum\limits_{n = 0}^{\infty} \binom{3n}{n}\frac{(xy)^n}{(1+2n)(1+y)^{3n}}: \qquad {x = \frac{T(1+y)^3}{y(1+T)^3}}.} \end{equation}$

(2.4)

Analogously, by putting the harmonic number $\mathbf{H}_{{n}}$ inside the series (2.4), we have

$\begin{align*} \sum\limits_{n = 0}^{\infty} \binom{3n}{n}\frac{ \mathbf{H}_{{n}}\; y^n}{(1+2n)(1+y)^{3n}} & = \int_0^1\frac{\psi(T(1))-\psi(T(x))}{1-x}dx \quad \boxed{x = \frac{T(1+y)^3}{y(1+T)^3}}\\ & = \frac{(1+y)^3}{y}\int_0^y \frac{(y-T)(1-2T)\; dT}{(1+T)^4\big\{1-\frac{T(1+y)^3}{y(1+T)^3}\big\}}\\ & = \int_0^y\frac{(1+y)^3(T-y)(1-2T)\; dT}{(1+T)\big\{T(1+y)^3-y(1+T)^3\big\}}\\ & = \int_0^y\frac{(1+y)^3(1-2T)\; dT}{(1+T)(1-T^2y-Ty^2-3Ty)}. \end{align*}$

By decomposing the integrand into partial fractions

$\begin{align*} \frac{(1+y)^3(1-2T)\; dT}{(1+T)(1-T^2y-Ty^2-3Ty)} = &3\frac{1+y}{1+T} -\frac{(1+y)(3\sqrt{y}+\sqrt{4+y})}{2\sqrt{y} \big(T+\frac{(3+y)\sqrt{y}+(1+y)\sqrt{4+y}}{2\sqrt{y}}\big)}\\ -&\frac{(1+y)(3\sqrt{y}-\sqrt{4+y})}{2\sqrt{y} \big(T+\frac{(3+y)\sqrt{y}-(1+y)\sqrt{4+y}}{2\sqrt{y}}\big)}, \end{align*}$

we can explicitly compute the integral

$\begin{align*} \int_0^y\frac{(1+y)^3(1-2T)\; dT}{(1+T)(1-T^2y-Ty^2-3Ty)} & = 3(1+y)\ln(1+y)\\ &-\frac{(1+y)(3\sqrt{y}+\sqrt{4+y})}{2\sqrt{y}} \ln\bigg(\frac{(1+y)(3\sqrt{y}+\sqrt{4+y})} {(3+y)\sqrt{y}+(1+y)\sqrt{4+y}}\bigg)\\ &-\frac{(1+y)(3\sqrt{y}-\sqrt{4+y})}{2\sqrt{y}} \ln\bigg(\frac{(1+y)(3\sqrt{y}-\sqrt{4+y})} {(3+y)\sqrt{y}-(1+y)\sqrt{4+y}}\bigg), \end{align*}$

which is simplified into the closed formula as in the theorem below:

Theorem 2 ( $\rho < y < 1/2$ ).

$\begin{align*} \sum\limits_{n = 1}^{\infty} \binom{3n}{n}\frac{ \mathbf{H}_{{n}}\; y^n}{(1+2n)(1+y)^{3n}} = &\int_0^y\frac{(1+y)^3(1-2T)\; dT}{(1+T)(1-T^2y-Ty^2-3Ty)}\\ = &\frac{(1+y)\sqrt{4+y}}{2\sqrt{y}} \ln\Big(\frac{2-y(2+y+\sqrt{y^2+4y})}{2-y(2+y-\sqrt{y^2+4y})}\Big) -\frac32(1+y)\ln(1-2y). \end{align*}$

2.3. Series with $\mathbf{H}_{{n}}/n$

Alternatively, rewriting Lambert's first series

$\frac{(1+T)^{ \alpha}-1}{ \alpha} = \sum\limits_{n = 1}^{\infty} \frac{1}{ \alpha+n \beta} \binom{ \alpha+n \beta}{n}\tau^n: \; \quad{\tau = T/(1+T)^{ \beta}}$

and letting " $\alpha\to0$ and $\beta = 3$ ", we deduce the identity

$\begin{equation} { { \sigma(T)}: = 3\ln(1+T) = \sum\limits_{n = 1}^{\infty} \binom{3n}{n}\frac{(xy)^n}{n(1+y)^{3n}}: \quad {x = \frac{T(1+y)^3}{y(1+T)^3}}.} \end{equation}$

(2.5)

By inserting $\mathbf{H}_{{n}}$ in the above series, we obtain an integral expression

$\begin{align*} \sum\limits_{n = 1}^{\infty} \binom{3n}{n}\frac{ \mathbf{H}_{{n}}(xy)^n}{n(1+y)^{3n}} & = \int_0^1\frac{ \sigma\big(T(1)\big)- \sigma\big(T(x)\big)}{1-x}dx \qquad\boxed{x\to\frac{T(1+y)^3}{y(1+T)^3}}\\ & = 3\int_0^y\ln\Big(\frac{1+T}{1+y}\Big) \frac{(1-2 T) (1+y)^3}{(1+T) (T-y)(1-T^2 y-T y^2-3 T y)}dT\\ & = 3\int_0^y\ln\Big(\frac{1+T}{1+y}\Big)\bigg\{\frac{1}{T-y}-\frac{3}{1+T} +\frac{y^2+3 y+2 T y}{T^2 y+T y^2+3 T y-1}\bigg\}dT. \end{align*}$

In view of the partial fraction decomposition

$\begin{align*} \frac{y^2+3 y+2 T y}{T^2 y+T y^2+3 T y-1} = \frac1{T+\dfrac{\sqrt{y}(3+y)+(1+y)\sqrt{4+y}}{2\sqrt{y}}\big)} +\frac1{T+\dfrac{\sqrt{y}(3+y)-(1+y)\sqrt{4+y}}{2\sqrt{y}}\big)}, \end{align*}$

it suffices to evaluate, for a parameter $\lambda$ subject to $T+ \lambda\not = 0$ for $T\in(0, y)$ , the parametric integral

$\begin{align*} J( \lambda, y) = \int_0^y\ln\Big(\frac{1+T}{1+y}\Big)\frac{dT}{T+ \lambda} & = \mathrm{Li}_2\Big(\frac{1+y}{1- \lambda}\Big) - \mathrm{Li}_2\Big(\frac{1}{1- \lambda}\Big) +\ln\Big(\frac{ \lambda}{ \lambda-1}\Big)\ln(1+y), \end{align*}$

with particular cases

$J(1, y) = -\frac12\ln^2(1+y) \quad\text{and}\quad J(-y, y) = \frac{\pi^2}6 +\ln(1+y)\ln\Big(\frac{y}{1+y}\Big) - \mathrm{Li}_2\Big(\frac1{1+y}\Big).$

Hereafter, the dilogarithm function (see Lewin [, §1.1]) is defined for all complex values of $z$ by

$\mathrm{Li}_2(z) = \int_0^z\frac{\ln(1-\tau)}{-\tau} \quad\text{and}\quad \mathrm{Li}_2(z) = \sum\limits_{n = 1}^{\infty}\frac{z^n}{n^2} \quad\text{for}\quad |z|\le1.$

By making substitution and then simplifying the terms concerning logarithms and dilogarithms in the resulting expression, we find the explicit formula as in the following theorem:

Theorem 3 ( $\rho < y < 1/2$ ).

$\begin{align*} \sum\limits_{n = 1}^{\infty} \binom{3n}{n}\frac{ \mathbf{H}_{{n}}\; y^n}{n(1+y)^{3n}} & = \int_0^y\frac{-3\; dT}{1+T}\ln\Big(1-\frac{T(1+y)^3}{y(1+T)^3}\Big)\\ & = {\frac{\pi^2}2+3\Theta(y)+3\ln(1+y)\ln\Big(\frac{y}{\sqrt{(1+y)^3}}\Big)}, \end{align*}$

where, for simplicity, $\Theta(y)$ denotes the sum of five dilogarithmic values

$\begin{align*} \Theta(y)& = \mathrm{Li}_2\Big(\frac{y-\sqrt{y^2+4y}}2\Big) + \mathrm{Li}_2\Big(\frac{y+\sqrt{y^2+4y}}2\Big)\\ &- \mathrm{Li}_2\Big(\frac{y-\sqrt{y^2+4y}}{2+2y}\Big) - \mathrm{Li}_2\Big(\frac{y+\sqrt{y^2+4y}}{2+2y}\Big) - \mathrm{Li}_2\Big(\frac{1}{1+y}\Big). \end{align*}$

3. Series with $\binom{3n}{n}$ and harmonic-like numbers

According to Bailey ^[3], the classical hypergeometric series reads as

${_{1+p}F_p}\left[ \begin{array}{{r}}{a_0, a_1, a_2, \cdots, a_p}\\{b_1, b_2, \cdots, b_p} \end{array}{\Big|{z}}\right] = \sum\limits_{k = 0}^{\infty}\frac{z^k}{k!} \frac{(a_0)_k(a_1)_k(a_2)_k\cdots(a_p)_k}{(b_1)_k(b_2)_k\cdots(b_p)_k}.$

We shall confine ourselves only to examining the series with $|z| < 1$ and none of the parameters in the numerator and denominator being a non-positive integer, which imply that the corresponding series is well-defined and convergent.

There are many transformation and summation formulae for hypergeometric series (see Bailey ^[3] and Brychkov [4, Chapter 8]) in the literature. In 1928, Bailey ^[2] discovered an important cubic transformation (see the next section). Half a century later, Gessel and Stanton [15, Eq (5.4)] found the following counterpart

$_3{F_2}\left[ {\begin{array}{*{20}{c}} {\frac{{1 + a}}{3},\frac{{2 + a}}{3},\frac{{3 + a}}{3}}\\ {\frac{{1 + a + c}}{2},\;\frac{{2 + a - c}}{2}} \end{array}|\;\frac{{27y}}{{4{{(1 + y)}^3}}}} \right]\; = \frac{{{{(1 + y)}^{1 + a}}}}{{1 - 2y}}{ \times _4}{F_3}\left[ {\begin{array}{*{20}{c}} {1 + \frac{a}{3},a,\;c,\;1 - c}\\ {\frac{a}{3},\frac{{1 + a + c}}{2},\;\frac{{2 + a - c}}{2}} \end{array}|\frac{{ - y}}{4}} \right].$

(3.1)

Under the parameter replacements " $a\to 6ax, \; c\to2cx$ ", this transformation can be reformulated as

$\begin{equation} \begin{aligned}[{c}]&\sum\limits_{n = 0}^{\infty} \Big(\frac{y}{4(1+y)^3}\Big)^n \frac{(1+6ax)_{3n}}{n!(1+3ax-cx)_n(\frac12+3ax+cx)_n}\\ &\; = \frac{(1+y)^{1+6ax}}{1-2y} \Bigg\{1+6cx\sum\limits_{k = 1}^{\infty}\Big(\frac{-y}4\Big)^k \frac{(k+2ax)(1+6ax)_{k-1}(1+2cx)_{k-1}(1-2cx)_k} {k!(1+3ax-cx)_k(\frac12+3ax+cx)_k}\Bigg\}, \end{aligned} \end{equation}$

(3.2)

where " $y\in(\rho, \tfrac12)$ ", or equivalently, $\dfrac{27y}{4(1+y)^3}\in(-1, 1)$ so that the series are convergent.

The above transformation is an analytic equality with respect to $x$ in the neighborhood of $x = 0$ . Equating the constant terms on both sides yields the algebraic identity below:

$\sum\limits_{n = 0}^{\infty} \binom{3n}{n}\frac{y^n}{(1+y)^{3n}} = \frac{1+y}{1-2y}.$

As illustrated in the introduction, harmonic-like numbers can be expressed in terms of $x$ -coefficients from quotients of shifted factorials. By examining the coefficients of $x$ across the above equation, we find that the corresponding equation involving harmonic numbers can be written as

$\begin{align*} &\sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n} \Big\{6a \mathbf{H}_{{3n}}-(3a-c) \mathbf{H}_{{n}}-2(3a+c) \mathbf{O}_{{n}}\Big\}\\ &\; = 6a\frac{1+y}{1-2y}\ln(1+y) +6c\frac{1+y}{1-2y}\sum\limits_{k = 1}^{\infty} \frac{(k-1)!}{(\frac12)_k}\Big(\frac{-y}4\Big)^k. \end{align*}$

The rightmost series admits the closed form expression, as in the lemma below:

Lemma 4 ( $-4 < y < 4$ ).

$\mathrm{W}(y) = \sum\limits_{k = 1}^{\infty} \frac{(k-1)!}{(\frac12)_k}\Big(\frac{-y}4\Big)^k = -\sqrt{\frac{4y}{4+y}} \mathrm{ArcTanh}\sqrt{\frac{y}{4+y}} = \sqrt{\frac{4y}{4+y}}\ln\frac{\sqrt{4+y}-\sqrt{y}}{2}.$

Proof. Recall the beta function (see [21, §16])

$B(a, b) = \int_0^1x^{a-1}(1-x)^{b-1}dx = \frac{ \Gamma(a) \Gamma(b)}{ \Gamma(a+b)}, \quad\text{where}\quad \Re(a) > 0 \quad\text{and}\quad \Re(b) > 0.$

We can express $\mathrm{W}(y)$ as a definite integral

$\begin{align*} \mathrm{W}(y)& = \sum\limits_{k = 1}^{\infty} \frac{(k-1)!}{(\frac12)_k}\Big(\frac{-y}4\Big)^k = \sum\limits_{k = 1}^{\infty}B\Big(k, \frac12\Big) \Big(\frac{-y}4\Big)^k\\ & = \sum\limits_{k = 1}^{\infty} \Big(\frac{-y}4\Big)^k \int_0^1\frac{x^{k-1}dx}{\sqrt{1-x}} = \int_0^1\frac{dx}{x\sqrt{1-x}} \sum\limits_{k = 1}^{\infty}\Big(\frac{-xy}4\Big)^k\\ & = \int_0^1\frac{-y\; dx}{\sqrt{1-x}(4+xy)} = \int_0^1\frac{-2y}{4+y-T^2y}dT, \end{align*}$

where the last line is justified by the change of variable $\sqrt{1-x}\to T$ . Then the formula in the lemma follows by evaluating the last integral explicitly. □

Therefore, we find the general formula below:

By specifying two parameters, $a$ and $c$ , we deduce four independent series

$\begin{align} &a = 0: &&\sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n}\bar{ \mathbf{H}}_{{2n}} = 3\frac{1+y}{1-2y} \mathrm{W}(y), \\ &c = 0: &&\sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n}\Big\{ \mathbf{H}_{{3n}}- \mathbf{H}_{{2n}}\Big\} = \frac{1+y}{1-2y}\ln(1+y), \\ &c = 3a: &&\sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n}\Big\{ \mathbf{H}_{{3n}}-2 \mathbf{O}_{{n}}\Big\} = \frac{1+y}{1-2y}\ln(1+y) -3\frac{1+y}{1-2y} \mathrm{W}(y), \\ &c = -3a: &\quad&\sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n} \Big\{ \mathbf{H}_{{3n}}- \mathbf{H}_{{n}}\Big\} = \frac{1+y}{1-2y}\ln(1+y) +3\frac{1+y}{1-2y} \mathrm{W}(y). \end{align}$

Combining the above formulae with Theorem 1 and then keeping in mind the equalities

$\begin{align*} \mathbf{H}_{{2n}}& = \bar{ \mathbf{H}}_{{2n}}+ \mathbf{H}_{{n}}, \\ \mathbf{O}_{{n}}\; & = \bar{ \mathbf{H}}_{{2n}}+\frac{ \mathbf{H}_{{n}}}2, \\ \mathbf{H}_{{3n}}& = ( \mathbf{H}_{{3n}}- \mathbf{H}_{{n}})+ \mathbf{H}_{{n}}, \end{align*}$

we establish further closed formulae for four series with a free variable " $y$ ".

Theorem 5 ( $\rho < y < 1/2$ ).

$\begin{align} &\text{(a)} &&\sum\limits_{n = 1}^{\infty} \binom{3n}{n} \frac{\bar{ \mathbf{H}}_{{2n}}\; y^n}{(1+y)^{3n}} && = \frac{3+3y}{1-2y}\sqrt{\frac{4y}{4+y}} \ln\frac{\sqrt{4+y}+\sqrt{y}}{2}, \\ &\text{(b)} &&\sum\limits_{n = 1}^{\infty} \binom{3n}{n} \frac{ \mathbf{H}_{{2n}}\; y^n}{(1+y)^{3n}} && = \frac{3+3y}{2-4y}\sqrt{\frac{y}{4+y}} \ln\Big(\frac{2(1+y+\sqrt{y^2+4y})} {2-y(2+y+\sqrt{y^2+4y})}\Big) -\frac{(3+3y)\ln(1-2y)}{2(1-2y)}, \\ &\text{(c)} &&\sum\limits_{n = 1}^{\infty} \binom{3n}{n} \frac{ \mathbf{O}_{{n}}\; y^n}{(1+y)^{3n}} && = \frac{3+3y}{2-4y}\sqrt{\frac{y}{4+y}} \ln\Big(\frac{1+y+\sqrt{y^2+4y}} {\sqrt{1-2y}}\Big) -\frac{(3+3y)\ln(1-2y)}{4(1-2y)}, \\ &\text{(d)} &\quad&\sum\limits_{n = 1}^{\infty} \binom{3n}{n} \frac{ \mathbf{H}_{{3n}}\; y^n}{(1+y)^{3n}} && = \frac{1+y}{2-4y}\ln\Big(\frac{(1+y)^2}{(1-2y)^3}\Big) +\frac{3+3y}{2-4y}\sqrt{\frac{y}{4+y}} \ln\Big(\frac{2(1+y+\sqrt{y^2+4y})} {2-y(2+y+\sqrt{y^2+4y})}\Big). \end{align}$

Numerous infinite series identities with different convergence rates can be derived by assigning particular values for $y$ in Theorems 1 and 5. Five groups of representative examples are recorded as corollaries. In particular, the second formula (b) and the last formula (e) displayed in Corollary 6 confirm (1.1) and (1.2) conjectured by Sun [23, Eqs (2.12) and (2.13)].

Corollary 6 ( $y\to\tfrac{3\sqrt3-5}{2}$ in Theorems 1 and 5).

$\begin{align} &\text{(a)} &\quad&\sum\limits_{n = 1}^{\infty}\Big(\frac2{27}\Big)^n \binom{3n}{n}\bar{ \mathbf{H}}_{{2n}} && = \frac{\sqrt3-3}2\ln\big(\sqrt3-1\big), \\ &\text{(b)} &&\sum\limits_{n = 1}^{\infty}\Big(\frac{2}{27}\Big)^n \binom{3n}{n} \mathbf{H}_{{2n}} && = \frac{1+\sqrt3}2\Big\{\sqrt3\ln2-3\ln(3-\sqrt3)\Big\}, \\ &\text{(c)} &&\sum\limits_{n = 1}^{\infty}\Big(\frac{2}{27}\Big)^n \binom{3n}{n} \mathbf{H}_{{n}} && = \frac{3-\sqrt3}{4}\ln\Big(\frac43\Big) +\sqrt3\ln\Big(\frac{2+\sqrt3}3\Big), \\ &\text{(d)} &&\sum\limits_{n = 1}^{\infty}\Big(\frac2{27}\Big)^n \binom{3n}{n} \mathbf{O}_{{n}} && = \frac{3+\sqrt3}2\bigg\{\ln\Big(\frac{1+\sqrt3}{\sqrt2}\Big) -\frac{\sqrt3}4\ln3\bigg\}, \\ &\text{(e)} &&\sum\limits_{n = 1}^{\infty}\Big(\frac2{27}\Big)^n \binom{3n}{n} \mathbf{H}_{{3n}} && = \frac{1+\sqrt3}4\ln\Big(\frac{4(7+4\sqrt3)}3\Big)-\sqrt3\ln2. \end{align}$

Corollary 7 ( $y\to\sqrt5-2$ in Theorems 1 and 5).

$\begin{align} &\text{(a)} &\quad&\sum\limits_{n = 1}^{\infty}\Big(\frac18\Big)^n \binom{3n}{n}\bar{ \mathbf{H}}_{{2n}} && = \frac{3 (1-\sqrt{5}) }{\sqrt5}\ln \bigg(\frac{\sqrt{5}-1}{2} \bigg), \\ &\text{(b)} &&\sum\limits_{n = 1}^{\infty}\Big(\frac18\Big)^n \binom{3n}{n} \mathbf{H}_{{2n}} && = \frac{3 (2+\sqrt{5}) }{\sqrt5} \ln \bigg(\frac{3+\sqrt5}{2}\bigg)-\frac{3 \ln5}{\sqrt5}, \\ &\text{(c)} &&\sum\limits_{n = 1}^{\infty}\Big(\frac18\Big)^n \binom{3n}{n} \mathbf{H}_{{n}} && = \frac32 \ln\bigg(\frac{3+\sqrt5}{2} \bigg)+\frac{3 }{\sqrt5}\ln \bigg(\frac{11+5 \sqrt5}{10} \bigg), \\ &\text{(d)} &&\sum\limits_{n = 1}^{\infty}\Big(\frac18\Big)^n \binom{3n}{n} \mathbf{O}_{{n}} && = \frac{3(1+\sqrt5)}{2 \sqrt5} \ln(2+\sqrt5)-\frac{3 \ln5}{2 \sqrt5}, \\ &\text{(e)} &&\sum\limits_{n = 1}^{\infty}\Big(\frac18\Big)^n \binom{3n}{n} \mathbf{H}_{{3n}} && = \ln (11+5 \sqrt5)+\frac{3 }{\sqrt5}\ln \bigg(\frac{4+2\sqrt5}{5} \bigg). \end{align}$

Corollary 8 ( $y\to\tfrac{\sqrt5-2}{\sqrt5}$ in Theorems 1 and 5).

$\begin{align} &\text{(a)} &\quad&\sum\limits_{n = 1}^{\infty}\Big(\frac5{64}\Big)^n \binom{3n}{n}\bar{ \mathbf{H}}_{{2n}} && = \frac{6 (\sqrt{5}-1) }{11}\ln \bigg(\frac{5-\sqrt{5}}{2} \bigg), \\ &\text{(b)} &&\sum\limits_{n = 1}^{\infty}\Big(\frac5{64}\Big)^n \binom{3n}{n} \mathbf{H}_{{2n}} && = \frac{6\sqrt5}{11} \ln\bigg(\frac{25+9 \sqrt5}{22} \bigg) -\frac{3}{11} \ln \bigg(\frac{23-3 \sqrt5}{10} \bigg), \\ &\text{(c)} &&\sum\limits_{n = 1}^{\infty}\Big(\frac5{64}\Big)^n \binom{3n}{n} \mathbf{H}_{{n}} && = -\frac{3\sqrt5}{11} \ln\bigg(\frac{267-119\sqrt5}{2} \bigg) -\frac{6}{11} \ln \bigg(\frac{7+ \sqrt5}{10} \bigg), \\ &\text{(d)} &&\sum\limits_{n = 1}^{\infty}\Big(\frac5{64}\Big)^n \binom{3n}{n} \mathbf{O}_{{n}} && = \frac{3(1+\sqrt5)}{11} \ln\bigg(\frac{4+ \sqrt5}{11} \bigg) +\frac{3\sqrt5 }{11 }\ln5, \\ &\text{(e)} &&\sum\limits_{n = 1}^{\infty}\Big(\frac5{64}\Big)^n \binom{3n}{n} \mathbf{H}_{{3n}} && = \frac{6\sqrt5}{11} \ln\bigg(\frac{4(4+ \sqrt5)}{11} \bigg) +\frac1{11}\ln\bigg(\frac{200(2667-73\sqrt5)}{11^6}\bigg). \end{align}$

Corollary 9 ( $y\to-\tfrac1{10}$ in Theorems 1 and 5).

$\begin{align} &\text{(a)} &\quad&\sum\limits_{n = 1}^{\infty}\Big(-\frac{100}{729}\Big)^n \binom{3n}{n}\bar{ \mathbf{H}}_{{2n}} && = -\frac38\sqrt{\frac{3}{13}} \mathrm{ArcCos}\Big(\frac{161}{200}\Big), \\ &\text{(b)} &&\sum\limits_{n = 1}^{\infty}\Big(-\frac{100}{729}\Big)^n \binom{3n}{n} \mathbf{H}_{{2n}} && = -\frac98\ln\Big(\frac65\Big) -\frac38\sqrt{\frac{3}{13}} \mathrm{ArcCos}\Big(\frac{67}{80}\Big), \\ &\text{(c)} &&\sum\limits_{n = 1}^{\infty} \Big(-\frac{100}{729}\Big)^n \binom{3n}{n} \mathbf{H}_{{n}} && = -\frac98\ln\Big(\frac65\Big) +\frac38\sqrt{\frac{3}{13}} \mathrm{ArcCos}\Big(\frac{7987}{8000}\Big), \\ &\text{(d)} &&\sum\limits_{n = 1}^{\infty}\Big(-\frac{100}{729}\Big)^n \binom{3n}{n} \mathbf{O}_{{n}} && = -\frac9{16}\ln\Big(\frac65\Big) -\frac3{16}\sqrt{\frac{3}{13}} \mathrm{ArcCos}\Big(\frac{7}{20}\Big), \\ &\text{(e)} &&\sum\limits_{n = 1}^{\infty}\Big(-\frac{100}{729}\Big)^n \binom{3n}{n} \mathbf{H}_{{3n}} && = -\frac3{8}\ln\Big(\frac{32}{15}\Big) -\frac3{8}\sqrt{\frac{3}{13}} \mathrm{ArcCos}\Big(\frac{67}{80}\Big). \end{align}$

Corollary 10 ( $y\to\tfrac{7-3\sqrt5}{2}$ in Theorems 1 and 5).

$\begin{align} &\text{(a)} &\; &\sum\limits_{n = 1}^{\infty}\Big(\frac{3+\sqrt5}{54}\Big)^n \binom{3n}{n}\bar{ \mathbf{H}}_{{2n}} && = \sqrt{\frac{6\sqrt5-6}{\sqrt5}} \ln\Big(\frac{1+\sqrt{15+6\sqrt5}}{3+\sqrt5}\Big), \\ &\text{(b)} &&\sum\limits_{n = 1}^{\infty}\Big(\frac{3+\sqrt5}{54}\Big)^n \binom{3n}{n} \mathbf{H}_{{2n}} && = \frac{3(1+\sqrt5)}4\ln\Big(\frac{2+\sqrt5}{3}\Big) +\frac{\sqrt{3/2}}{\sqrt{5+\sqrt5}} \ln\Big(\frac{3\sqrt5-1+\sqrt{30-6\sqrt5}}{4}\Big), \\ &\text{(c)} &&\sum\limits_{n = 1}^{\infty} \Big(\frac{3+\sqrt5}{54}\Big)^n \binom{3n}{n} \mathbf{H}_{{n}} && = \frac{3(1+\sqrt5)}4\ln\Big(\frac{2+\sqrt5}{3}\Big) +{\frac{\sqrt{3/2}}{\sqrt{5+\sqrt5}} \ln\Big(\frac{11\sqrt3+\sqrt{185-82\sqrt5}} {11\sqrt3-\sqrt{185-82\sqrt5}}\Big)}, \\ &\text{(d)} &&\sum\limits_{n = 1}^{\infty}\Big(\frac{3+\sqrt5}{54}\Big)^n \binom{3n}{n} \mathbf{O}_{{n}} && = \frac{3(1+\sqrt5)}8\ln\Big(\frac{2+\sqrt5}{3}\Big) +\frac{\sqrt{3/8}}{\sqrt{5+\sqrt5}} \ln\Big(3\sqrt5-4+2\sqrt{15-6\sqrt5}\Big), \\ &\text{(e)} &&\sum\limits_{n = 1}^{\infty}\Big(\frac{3+\sqrt5}{54}\Big)^n \binom{3n}{n} \mathbf{H}_{{3n}} && = \frac{1+\sqrt5}4\ln\Big(\frac{11+5\sqrt5}{6}\Big) +\frac{\sqrt{3/2}}{\sqrt{5+\sqrt5}} \ln\Big(\frac{3\sqrt5-1+\sqrt{30-6\sqrt5}}{4}\Big). \end{align}$

We remark that the second formula (b) and the last series (e) in this corollary refine identity (1.3) conjectured by Sun [23, Eq (2.15)]:

$\sum\limits_{n = 1}^{\infty}\Big(\frac{3+\sqrt5}{54}\Big)^n \binom{3n}{n}\Big\{ \mathbf{H}_{{3n}}- \mathbf{H}_{{2n}}\Big\} = \frac{1+\sqrt5}{2}\ln\Big(\frac{9-3\sqrt5}{2}\Big).$

4. Series with $\binom{3n}{n}\Big/(2n+1)$ and harmonic-like numbers

Now, we turn to examine Bailey's cubic transformation (Bailey [2, Eq (4.05)])

$\begin{equation} {_3{F_2}}\left[ \begin{array}{{c}}{\frac{a}3, \frac{1+a}3, \frac{2+a}3}\\{\\[-3.5mm]\frac{1+a+c}2, \; \frac{2+a-c}2} \end{array}{\Big|{\; \frac{27y}{4(1+y)^3}\; }}\right] = (1+y)^{a}\times{_3{F_2}}\left[ \begin{array}{{c}}{a, \; c, \; 1-c}\\{\frac{1+a+c}2, \; \frac{2+a-c}2} \end{array}{\Big|{\; \frac{-y}4\; }}\right]. \end{equation}$

(4.1)

Under the parameter replacements " $a\to 1+6ax, \; c\to2cx$ ", it can be reformulated as

$\begin{equation} \begin{aligned}[{c}]&\sum\limits_{n = 0}^{\infty} \Big(\frac{y}{4(1+y)^3}\Big)^n \frac{(1+6ax)_{3n}}{n!(1+3ax+cx)_n(\frac32+3ax-cx)_n}\\ &\; = (1+y)^{1+6ax}\times \Bigg\{1+2cx\sum\limits_{k = 1}^{\infty}\Big(\frac{-y}4\Big)^k \frac{(1+6ax)_{k}(1+2cx)_{k-1}(1-2cx)_k} {k!(1+3ax+cx)_k(\frac32+3ax-cx)_k}\Bigg\}, \end{aligned} \end{equation}$

(4.2)

where, as before, both series converge for " $y\in(\rho, \tfrac12)$ ", or equivalently, $\dfrac{27y}{4(1+y)^3}\in(-1, 1)$ .

When $x = 0$ , the above series becomes the following simpler one:

$\sum\limits_{n = 0}^{\infty} \binom{3n}{n}\frac{y^n}{(1+2n)(1+y)^{3n}} = 1+y.$

Instead, the coefficients of $x$ across (4.2) result in the general formula as below:

$\begin{align*} &\sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n} \bigg\{\frac{6a \mathbf{H}_{{3n}}-(3a+c) \mathbf{H}_{{n}}-2(3a-c) \mathbf{O}_{{n}}} {1+2n}+\frac{4(3a-c)n}{(1+2n)^2}\bigg\}\\ &\; = 6a(1+y)\ln(1+y) +2c(1+y)\sum\limits_{k = 1}^{\infty} \frac{(k-1)!}{(\frac32)_k}\Big(\frac{-y}4\Big)^k. \end{align*}$

For brevity, we introduce the notations $\mathrm{U}(y)$ and $\mathrm{V}(y)$ for the two series:

$\begin{align*} \mathrm{U}(y)& = \sum\limits_{k = 1}^{\infty} \frac{(k-1)!}{(\frac32)_k}\Big(\frac{-y}4\Big)^k\\ \text{and} \\ \mathrm{V}(y)& = \sum\limits_{n = 0}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n}\frac{n}{(1+2n)^2}. \end{align*}$

They are evaluated in closed form by the following lemma:

Lemma 11.

$\begin{align} \text{(a)}\qquad \mathrm{U}(y)& = 2+\frac{2\sqrt{4+y}}{\sqrt y} \ln\frac{\sqrt{4+y}-\sqrt{y}}{2}, &\quad&-4 < y < 4;\\ \text{(b)}\qquad \mathrm{V}(y)& = \frac{\sqrt{(1+y)^3}}{\sqrt y} \ln\big(\sqrt y+\sqrt{1+y}\big)-(1+y), &&\rho < y < 1/2. \end{align}$

Proof. The first series can be evaluated by computing the integral as follows:

$\begin{align*} \mathrm{U}(y)& = \sum\limits_{k = 1}^{\infty} \frac{(k-1)!}{(\frac32)_k}\Big(\frac{-y}4\Big)^k = \sum\limits_{k = 1}^{\infty}B\Big(k, \frac32\Big) \Big(\frac{-y}4\Big)^k\\ & = \sum\limits_{k = 1}^{\infty} \Big(\frac{-y}4\Big)^k \int_0^1x^{k-1}\sqrt{1-x}dx = \int_0^1\frac{\sqrt{1-x}\; dx}{x} \sum\limits_{k = 1}^{\infty}\Big(\frac{-xy}4\Big)^k\\ & = \int_0^1\frac{-y\sqrt{1-x}\; dx}{4+xy} = \int_0^1\frac{-2yT^2}{4+y-T^2y}dT \qquad\boxed{\sqrt{1-x}\to T}\\ & = 2+\frac{2\sqrt{4+y}}{\sqrt y}\ln\frac{\sqrt{4+y}-\sqrt{y}}{2}. \end{align*}$

For the second series $\mathrm{V}(y)$ , by shifting forward the summation index $n\to n+1$ , we can rewrite it as

$\mathrm{V}(y) = \sum\limits_{n = 0}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n}\frac{n}{(1+2n)^2} = \frac{y}{(1+y)^{3}}\sum\limits_{n = 0}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n+3}{n}\frac{1}{3+2n}.$

Then, for " $\alpha = \beta = 3$ ", perform the replacement $\tau\to{xy}/{(1+y)^3}$ in Lambert's second series

$\sum\limits_{n = 0}^{\infty}\frac{x^ny^n}{(1+y)^{3n}} \binom{3n+3}{n} = \frac{(1+T)^4}{1-2T}.$

Now multiplying this equation by $x^\frac12$ and then integrating across over $x\in[0, 1]$ , we can proceed with

$\begin{align*} \sum\limits_{n = 0}^{\infty} \frac{y^n}{(1+y)^{3n}} \binom{3n+3}{n}\frac2{3+2n} & = \int_0^1x^{\frac12}\frac{(1+T)^4}{1-2T}dx \qquad\boxed{x = \frac{T(1+y)^3}{y(1+T)^3}}\\ & = \frac{\sqrt{(1+y)^9}}{y\sqrt y} \int_0^y\frac{\sqrt T}{\sqrt{(1+T)^3}}dT. \end{align*}$

By means of integration by parts, we can further compute

$\begin{align*} \int_0^y\frac{\sqrt T}{\sqrt{(1+T)^3}}dT & = \frac{-2\sqrt y}{\sqrt{1+y}} +\int_0^y\frac{1}{\sqrt{T(1+T)}}dT\\ & = 2\ln(\sqrt{y}+\sqrt{1+y})-\frac{2\sqrt y}{\sqrt{1+y}}. \end{align*}$

After substitution, the second formula in Lemma 11 follows accordingly. □

Therefore, we obtain the following reduced transformation:

$\begin{equation} \begin{aligned}[{c}]&\sum\limits_{n = 0}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n} \bigg\{\frac{6a \mathbf{H}_{{3n}}-(3a+c) \mathbf{H}_{{n}}-2(3a-c) \mathbf{O}_{{n}}} {1+2n}\bigg\}\\ &\; = 6a(1+y)\ln(1+y)+2c(1+y) \mathrm{U}(y)-4(3a-c) \mathrm{V}(y). \end{aligned} \end{equation}$

(4.3)

By specifying two parameters, $a$ and $c$ in the above equality, we deduce four independent series:

$\begin{align} &a = 0:&\quad &\sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n} \bigg\{\frac{\bar{ \mathbf{H}}_{{2n}}}{1+2n}\bigg\} = (1+y) \mathrm{U}(y)+2 \mathrm{V}(y), \\ &c = 0:&\quad &\sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n} \bigg\{\frac{ \mathbf{H}_{{3n}}- \mathbf{H}_{{2n}}}{1+2n}\bigg\} = (1+y)\ln(1+y)-2 \mathrm{V}(y), \\ &c = 3a:&\quad &\sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n} \bigg\{\frac{ \mathbf{H}_{{3n}}- \mathbf{H}_{{n}}}{1+2n}\bigg\} = (1+y)\ln(1+y)+(1+y) \mathrm{U}(y), \\ &c = -3a:&\quad &\sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n} \bigg\{\frac{ \mathbf{H}_{{3n}}-2 \mathbf{O}_{{n}}}{1+2n}\bigg\} = (1+y)\ln(1+y)-(1+y) \mathrm{U}(y)-4 \mathrm{V}(y). \end{align}$

For the above four equations, first replacing $\mathrm{U}(y)$ and $\mathrm{V}(y)$ by their algebraic values in Lemma 11, and then putting the resulting equations in conjunction with that in Theorem 2, we establish, after some routine simplifications, four closed formulae as in the following theorem:

Theorem 12 ( $\rho < y < 1/2$ ).

$\begin{align} &\text{(a)} &&\sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n}\frac{\bar{ \mathbf{H}}_{{2n}}}{1+2n} && = \frac{2\sqrt{(1+y)^3}}{\sqrt y}\ln\big(\sqrt y+\sqrt{1+y}\big)\\ &&&&& \nonumber +\frac{2(1+y)\sqrt{4+y}}{\sqrt y}\ln\bigg(\frac{\sqrt{4+y}-\sqrt{y}}{2}\bigg), \\ &\text{(b)} &&\sum\limits_{n = 1}^{\infty}\frac{\; y^n}{(1+y)^{3n}} \binom{3n}{n}\frac{ \mathbf{H}_{{2n}}}{1+2n} && = \frac{2\sqrt{(1+y)^3}}{\sqrt y}\ln\big(\sqrt y+\sqrt{1+y}\big) -\frac{3+3y}{2}\ln(1-2y)\\ &&&&& \nonumber +\frac{(1+y)\sqrt{4+y}}{2\sqrt{y}} \ln\bigg(\frac{2+2y-2\sqrt{y^2+4y}} {2-y(2+y-\sqrt{y^2+4y})}\bigg), \\ &\text{(c)} &&\sum\limits_{n = 1}^{\infty}\frac{\; y^n}{(1+y)^{3n}} \binom{3n}{n}\frac{ \mathbf{O}_{{n}}}{1+2n} && = \frac{2\sqrt{(1+y)^3}}{\sqrt y}\ln\big(\sqrt y+\sqrt{1+y}\big) -\frac{3+3y}{4}\ln(1-2y)\\ &&&&& \nonumber +\frac{(1+y)\sqrt{4+y}}{2\sqrt{y}} \ln\Big(\frac{1+y-\sqrt{y^2+4y}} {\sqrt{1-2y}}\Big), \\ &\text{(d)} &\quad&\sum\limits_{n = 1}^{\infty}\frac{\; y^n}{(1+y)^{3n}} \binom{3n}{n}\frac{ \mathbf{H}_{{3n}}}{1+2n} && = 2(1+y)+\frac{1+y}{2}\ln\bigg(\frac{(1+y)^2}{(1-2y)^3}\bigg)\\ &&&&& \nonumber+\frac{(1+y)\sqrt{4+y}}{2\sqrt{y}} \ln\bigg(\frac{2+2y-2\sqrt{y^2+4y}} {2-y(2+y-\sqrt{y^2+4y})}\bigg). \end{align}$

By assigning particular values for $y$ in Theorems 2 and 12, we can deduce numerous infinite series identities with different convergence rates. In order to show the variety of implications, we record five representative classes of infinite series values.

Corollary 13 ( $y\to\sqrt5-2$ in Theorems 2 and 12).

$\begin{align} &\text{(a)}\quad \sum\limits_{n = 1}^{\infty}\Big(\frac18\Big)^n \binom{3n}{n}\frac{\bar{ \mathbf{H}}_{{2n}}}{1+2n} && = 2 \sqrt2\ln \bigg(\frac{3+\sqrt{10}}{1+\sqrt2}\bigg) -(3+\sqrt5) \ln \bigg(\frac{1+\sqrt5}2\bigg), \\ &\text{(b)}\quad \sum\limits_{n = 1}^{\infty}\Big(\frac18\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{2n}}}{1+2n} && = 2 \sqrt2\ln \bigg(\frac{3+\sqrt{10}}{1+\sqrt2}\bigg) -(6-4\sqrt5) \ln \bigg(\frac{1+\sqrt5}2\bigg)-\sqrt5\ln5, \\ &\text{(c)}\quad \sum\limits_{n = 1}^{\infty}\Big(\frac18\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{n}}}{1+ 2n} && = (5\sqrt5-3) \ln \bigg(\frac{1+\sqrt5}2\bigg)-\sqrt5\ln5, \\ &\text{(d)}\quad \sum\limits_{n = 1}^{\infty}\Big(\frac18\Big)^n \binom{3n}{n}\frac{ \mathbf{O}_{{n}}}{1+2n} && = 2 \sqrt2\ln \bigg(\frac{3+\sqrt{10}}{1+\sqrt2}\bigg) -\bigg(\frac92-\frac{3\sqrt5}2\bigg) \ln \bigg(\frac{1+\sqrt5}2\bigg)-\frac{\sqrt5}2\ln5 , \\ &\text{(e)}\quad \sum\limits_{n = 1}^{\infty}\Big(\frac18\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{3n}}}{1+2n} && = (\sqrt5-1)(2+\ln2)-(5-3\sqrt5) \ln \bigg(\frac{1+\sqrt5}2\bigg)-\sqrt5\ln5. \end{align}$

Corollary 14 ( $y\to\tfrac{3\sqrt3-5}2$ in Theorems 2 and 12).

$\begin{align} &\text{(a)}\quad \sum\limits_{n = 1}^{\infty}\Big(\frac2{27}\Big)^n \binom{3n}{n}\frac{\bar{ \mathbf{H}}_{{2n}}}{1+2n} && = \frac{9+3 \sqrt3}2 \ln (\sqrt{3}-1)+\frac9{\sqrt6} \ln \bigg(\frac{3+2 \sqrt2}{\sqrt3+\sqrt2}\bigg), \\ &\text{(b)}\quad \sum\limits_{n = 1}^{\infty}\Big(\frac2{27}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{2n}}}{1+2n} && = \frac92 \ln \bigg(\frac{3-\sqrt3}2\bigg) +\frac{3 \sqrt3}2 \ln \bigg(\frac{(\sqrt{3}+1)^3}{12\sqrt3}\bigg) +\frac9{\sqrt6} \ln \bigg(\frac{3+2 \sqrt2}{\sqrt3+\sqrt2}\bigg), \\ &\text{(c)}\quad \sum\limits_{n = 1}^{\infty}\Big(\frac2{27}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{n}}}{1+2n} && = \frac94\ln\Big(\frac34\Big) +\frac{3 \sqrt3}2 \ln \bigg(\frac{(\sqrt{3}+1)^4}{24\sqrt3}\bigg), \\ &\text{(d)}\quad \sum\limits_{n = 1}^{\infty}\Big(\frac2{27}\Big)^n \binom{3n}{n}\frac{ \mathbf{O}_{{n}}}{1+2n} && = \frac34(3-\sqrt3)\ln(2\sqrt3-3)-\frac{9\ln 3 }{4\sqrt3} +\frac9{\sqrt6} \ln \bigg(\frac{3+2 \sqrt2}{\sqrt3+\sqrt2}\bigg), \\ &\text{(e)}\quad \sum\limits_{n = 1}^{\infty}\Big(\frac2{27}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{3n}}}{1+2n} && = 3\sqrt3-3+\frac32\ln\Big(\frac{2\sqrt3-3}2\Big) +\frac{3 \sqrt3}2 \ln \bigg(\frac{2+\sqrt3}{2\sqrt3}\bigg). \end{align}$

Corollary 15 ( $y\to3\sqrt2-4$ in Theorems 2 and 12).

$\begin{align} &\text{(a)}\quad \sum\limits_{n = 1}^{\infty}\Big(\frac{2+\sqrt2}{27}\Big)^n \binom{3n}{n}\frac{\bar{ \mathbf{H}}_{{2n}}}{1+2n} && = 3 \sqrt3\ln \bigg(\frac{\sqrt6-2}{\sqrt3-1}\bigg) +\frac{9\sqrt{2-\sqrt2}}{\sqrt6} \ln \Big(6\sqrt2-7+\sqrt{6 (2-\sqrt2)^3}\Big), \\ &\text{(b)}\quad \sum\limits_{n = 1}^{\infty}\Big(\frac{2+\sqrt2}{27}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{2n}}}{1+2n} && = 3 \sqrt3\ln \bigg(\frac{\sqrt3-1}{\sqrt2}\bigg) -(9-9\sqrt2)\ln\bigg(\frac{1+\sqrt2}{\sqrt3}\bigg) \\ &&&+\frac{9\sqrt{2-\sqrt2}}{\sqrt6} \ln \Big(6\sqrt2-7+\sqrt{6 (2-\sqrt2)^3}\Big), \\ &\text{(c)}\quad \sum\limits_{n = 1}^{\infty}\Big(\frac{2+\sqrt2}{27}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{n}}}{1+2n} && = 3 \sqrt3\ln \bigg(\frac{\sqrt2+\sqrt3}{2+\sqrt3}\bigg) -(9-9\sqrt2)\ln\bigg(\frac{1+\sqrt2}{\sqrt3}\bigg) , \\ &\text{(d)}\quad \sum\limits_{n = 1}^{\infty}\Big(\frac{2+\sqrt2}{27}\Big)^n \binom{3n}{n}\frac{ \mathbf{O}_{{n}}}{1+2n} && = \Big(\frac92-\frac9{\sqrt2}\Big)\ln\big(\sqrt6-\sqrt3\big) -\frac{3 \sqrt3}2\ln \big(\sqrt2+\sqrt3\big) \\ &&&+\frac{9\sqrt{2-\sqrt2}}{\sqrt6} \ln \Big(6\sqrt2-7+\sqrt{6 (2-\sqrt2)^3}\Big), \\ &\text{(e)}\quad \sum\limits_{n = 1}^{\infty}\Big(\frac{2+\sqrt2}{27}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{3n}}}{1+2n} && = 6\sqrt2-6+3 \sqrt3\ln \bigg(\frac{\sqrt3-1}{\sqrt2}\bigg) -(3-3\sqrt2)\ln\bigg(\frac{3+2\sqrt2}{\sqrt3}\bigg). \end{align}$

Corollary 16 ( $y\to\tfrac{1}5$ in Theorems 2 and 12).

$\begin{align} &\text{(a)}\quad \sum\limits_{n = 1}^{\infty}\Big(\frac{25}{216}\Big)^n \binom{3n}{n}\frac{\bar{ \mathbf{H}}_{{2n}}}{1+2n} && = \frac{12\sqrt6}{5} \ln \bigg(\frac{1+ \sqrt6}{\sqrt5}\bigg) -\frac{12\sqrt{21}}{5} \ln \bigg(\frac{1+ \sqrt{21}}{2\sqrt5}\bigg), \\ &\text{(b)}\quad \sum\limits_{n = 1}^{\infty}\Big(\frac{25}{216}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{2n}}}{1+2n} && = \frac95\ln\Big(\frac53\Big)+\frac{12\sqrt6}{5} \ln \bigg(\frac{1+ \sqrt6}{\sqrt5}\bigg) +\frac{6\sqrt{21}}{5} \ln \bigg(\frac{9- \sqrt{21}}{2\sqrt{15}}\bigg), \\ &\text{(c)}\quad \sum\limits_{n = 1}^{\infty}\Big(\frac{25}{216}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{n}}}{1+2n} && = \frac95\ln\Big(\frac53\Big) +\frac{3 \sqrt{21}}5 \ln \bigg(\frac{257-13\sqrt{21}}{250}\bigg), \\ &\text{(d)}\quad \sum\limits_{n = 1}^{\infty}\Big(\frac{25}{216}\Big)^n \binom{3n}{n}\frac{ \mathbf{O}_{{n}}}{1+2n} && = \frac9{10}\ln\Big(\frac53\Big)+\frac{12\sqrt6}{5} \ln \bigg(\frac{1+ \sqrt6}{\sqrt5}\bigg) +\frac{3\sqrt{21}}{5} \ln \bigg(\frac{6- \sqrt{21}}{\sqrt{15}}\bigg), \\ &\text{(e)}\quad \sum\limits_{n = 1}^{\infty}\Big(\frac{25}{216}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{3n}}}{1+2n} && = \frac{12}5+\frac35\ln\Big(\frac{20}3\Big) +\frac{6\sqrt{21}}{5} \ln \bigg(\frac{9- \sqrt{21}}{2\sqrt{15}}\bigg). \end{align}$

Corollary 17 ( $y\to-\tfrac1{10}$ in Theorems 2 and 12).

$\begin{align} &\text{(a)}\quad \sum\limits_{n = 1}^{\infty}\Big(-\frac{100}{729}\Big)^n \binom{3n}{n}\frac{\bar{ \mathbf{H}}_{{2n}}}{1+2n} && = \frac{27}5 \mathrm{arccot}3-\frac{9\sqrt{39}}5 \mathrm{arccot}\sqrt{39}, \\ &\text{(b)}\quad \sum\limits_{n = 1}^{\infty}\Big(-\frac{100}{729}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{2n}}}{1+2n} && = \frac{27}5 \mathrm{arccot}3 -\frac{9\sqrt{39}}{20} \mathrm{arccos}\Big(\frac{67}{80}\Big) -\frac{27}{20}\ln\Big(\frac65\Big), \\ &\text{(c)}\quad \sum\limits_{n = 1}^{\infty}\Big(-\frac{100}{729}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{n}}}{1+2n} && = \frac{9\sqrt{39}}{20} \mathrm{arccos}\Big(\frac{7987}{8000}\Big) -\frac{27}{20}\ln\Big(\frac65\Big), \\ &\text{(d)}\quad \sum\limits_{n = 1}^{\infty}\Big(-\frac{100}{729}\Big)^n \binom{3n}{n}\frac{ \mathbf{O}_{{n}}}{1+2n} && = \frac{27}5 \mathrm{arccot}3-\frac{9\sqrt{39}}{20} \mathrm{arccos}\bigg(\sqrt{\frac{27}{40}}\bigg) -\frac{27}{40}\ln\Big(\frac65\Big) , \\ &\text{(e)}\quad \sum\limits_{n = 1}^{\infty}\Big(-\frac{100}{729}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{3n}}}{1+2n} && = \frac95-\frac{9\sqrt{39}}{20} \mathrm{arccos}\Big(\frac{67}{80}\Big) -\frac{9}{20}\ln\Big(\frac{32}{15}\Big). \end{align}$

5. Series with $\binom{3n}{n}\Big/n$ and harmonic-like numbers

Under the parameter replacements " $a\to 6ax, \; c\to2cx$ ", Bailey's cubic transformation (4.1) can be reformulated alternatively as

$\begin{align*} &\sum\limits_{n = 0}^{\infty}\Big(\frac{y}{4(1+y)^3}\Big)^n \frac{(6ax)_{3n}}{n!(1+3ax-cx)_n(\frac12+3ax+cx)_n}\\ &\; = (1+y)^{6ax}\times \Bigg\{1+12acx^2\sum\limits_{k = 1}^{\infty}\Big(\frac{-y}4\Big)^k \frac{(1+6ax)_{k-1}(1+2cx)_{k-1}(1-2cx)_k} {k!(1+3ax-cx)_k(\frac12+3ax+cx)_k}\Bigg\}. \end{align*}$

Both series in the above equation converge subject to the same condition " $y\in(\rho, \tfrac12)$ ", or equivalently, $\dfrac{27y}{4(1+y)^3}\in(-1, 1)$ (as before).

Equating the constant terms from both sides yields the algebraic identity below

$\sum\limits_{n = 1}^{\infty} \binom{3n}{n}\frac{y^n}{n(1+y)^{3n}} = 3\ln(1+y).$

Instead, the coefficients of $x$ result in the following infinite series identity involving a free variable " $y$ " and two linear parameters $a$ and $c$ :

$\begin{align*} &\sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n} \bigg\{\frac{6a \mathbf{H}_{{3n}}-(3a-c) \mathbf{H}_{{n}}-2(3a+c) \mathbf{O}_{{n}}}{n}-\frac{2a}{n^2}\bigg\}\\ &\; = 9a\ln^2(1+y) +6c\sum\limits_{k = 1}^{\infty} \frac{(k-1)!}{k(\frac12)_k}\Big(\frac{-y}4\Big)^k. \end{align*}$

Let $\mathcal{{U}}(y)$ and $\mathcal{{V}}(y)$ stand for the two series

$\mathcal{{U}}(y) = \sum\limits_{k = 1}^{\infty} \frac{(k-1)!}{k(\frac12)_k}\Big(\frac{-y}4\Big)^k \quad\text{and}\quad \mathcal{{V}}(y) = \sum\limits_{n = 1}^{\infty}\frac{y^n}{n^2(1+y)^{3n}} \binom{3n}{n}.$

They admit closed formulae as in the lemma below:

Lemma 18.

$\begin{align} &\text{(a)}\quad\mathcal{{U}}(y) = -2\ln^2\bigg(\frac{2}{\sqrt y+\sqrt{4+y}}\bigg), \quad -4 < y < 4;\\ &\text{(b)}\quad\mathcal{{V}}(y) = 3\bigg\{\frac{{\pi}^2}{6}+\ln y\ln(1+y)-2\ln^2(1+y) - \mathrm{Li}_2\Big(\frac1{1+y}\Big)\bigg\}, \quad \rho < y < 1/2. \end{align}$

Proof. The first series $\mathcal{{U}}(y)$ can be reformulated as

$\begin{align*} \mathcal{{U}}(y)& = \sum\limits_{k = 1}^{\infty} \frac{(k-1)!}{k(\frac12)_k}\Big(\frac{-y}4\Big)^k = \sum\limits_{k = 1}^{\infty}B\Big(k, \frac12\Big) \Big(\frac{-y}4\Big)^k\\ & = \sum\limits_{k = 1}^{\infty} \frac{\big(-y/4)^k}k \int_0^1\frac{x^{k-1}dx}{\sqrt{1-x}} = \int_0^1\frac{dx}{x\sqrt{1-x}} \sum\limits_{k = 1}^{\infty}\frac{\big(-xy/4)^k}k\\ & = -\int_0^1\frac{\ln\big(1+\frac{xy}4\big)dx}{x\sqrt{1-x}} = -\int_0^1\frac{dT}{T}\ln\Big(\frac{1+(2+y)T+T^2}{(1+T)^2}\Big) \qquad\boxed{x\to \frac{4T}{(1+T)^2}}\\ & = \int_0^1\frac{dT}{T}\bigg\{2\ln(1+T) -\ln\Big(1+\frac{4T}{(\sqrt y+\sqrt{4+y})^2}\Big) -\ln\Big(1+\frac{4T}{(\sqrt y-\sqrt{4+y})^2}\Big)\bigg\}\\ & = \frac{{\pi}^2}6 + \mathrm{Li}_2\bigg(\frac{-4}{(\sqrt y+\sqrt{4+y})^2}\bigg) + \mathrm{Li}_2\bigg(\frac{-4}{(\sqrt y-\sqrt{4+y})^2}\bigg). \end{align*}$

Then the closed formula for $\mathcal{{U}}(y)$ in the lemma follows by applying the dilogarithmic equation

$\mathrm{Li}_2(-x)+ \mathrm{Li}_2(-x^{-1}) = -\frac{\pi^2}{6}-\frac{\ln^2x}2, \quad\text{where}\quad x > 0.$

For the second series, we can rewrite it by reindexing $n\to n+1$ as

$\mathcal{{V}}(y) = \sum\limits_{n = 1}^{\infty} \frac{y^n}{n^2(1+y)^{3n}} \binom{3n}{n} = \frac{3y}{(1+y)^{3}}\sum\limits_{n = 0}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n+2}{n}\frac{1}{(1+n)^2}.$

Now specializing Lambert's second series by " $\alpha = 2$ and $\beta = 3$ " and $\tau\to{xy}/{(1+y)^3}$

$\sum\limits_{n = 0}^{\infty}\frac{x^ny^n}{(1+y)^{3n}} \binom{3n+2}{n} = \frac{(1+T)^3}{1-2T}, \quad\text{where}\quad {x = \frac{T(1+y)^3}{y(1+T)^3}},$

then integrating twice across this equation, we can proceed with

$\begin{align*} &\sum\limits_{n = 0}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n+2}{n}\frac{1}{(1+n)^2}\\ &\; = \int_0^1\frac{dt}t\int_0^t\frac{(1+T)^3}{1-2T}dx = -\int_0^1\frac{(1+T)^3\ln x}{1-2T}dx\\ &\; = -\frac{(1+y)^3}{y} \int_0^y\frac{dT}{1+T}\ln\Big(\frac{ T(1+y)^3}{y(1+T)^3}\Big) \qquad\boxed{x = \frac{T(1+y)^3}{y(1+T)^3}}\\ &\; = -\frac{(1+y)^3}{y} \bigg\{\int_0^y\frac{dT}{1+T}\ln\Big(\frac{ T(1+y)^3}{y}\Big) -3\int_0^y\frac{\ln(1+T)}{1+T}dT\bigg\}\\ &\; = \frac{(1+y)^3}{y}\bigg\{\frac{{\pi}^2}{6}+\ln y\ln(1+y) -2\ln^2(1+y)- \mathrm{Li}_2\Big(\frac1{1+y}\Big)\bigg\}. \end{align*}$

By substitution, the closed formula for $\mathcal{{V}}(y)$ in the lemma follows consequently. □

Thus, we have the following reduced transformation:

$\begin{equation} {\sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n} \bigg\{\frac{6a \mathbf{H}_{{3n}}-(3a-c) \mathbf{H}_{{n}}-2(3a+c) \mathbf{O}_{{n}}}{n}\bigg\} = 9a\ln^2(1+y)+6c\mathcal{{U}}(y)+2a\mathcal{{V}}(y).} \end{equation}$

(5.1)

By specifying two parameters, $a$ and $c$ , we deduce four independent series

$\begin{align} &a = 0: &&\sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n}\frac{\bar{ \mathbf{H}}_{{2n}}}{n} = -3\mathcal{{U}}(y), \\ &c = 0: &&\sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n}\bigg\{\frac{ \mathbf{H}_{{3n}}- \mathbf{H}_{{2n}}}{n}\bigg\} = \frac32\ln^2(1+y)+\frac{\mathcal{{V}}(y)}3, \\ &c = 3a: &&\begin{aligned}[{t}]&\sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n}\bigg\{\frac{ \mathbf{H}_{{3n}}-2 \mathbf{O}_{{n}}}{n}\bigg\} = \frac32\ln^2(1+y)+3\mathcal{{U}}(y)+\frac{\mathcal{{V}}(y)}3, \end{aligned}\\ &c = -3a: &\quad&\begin{aligned}[{t}]&\sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n}\bigg\{\frac{ \mathbf{H}_{{3n}}- \mathbf{H}_{{n}}}{n}\bigg\} = \frac32\ln^2(1+y)-3\mathcal{{U}}(y)+\frac{\mathcal{{V}}(y)}3. \end{aligned} \end{align}$

For these four equations, replacing $\mathcal{{U}}(y)$ and $\mathcal{{V}}(y)$ first by their algebraic expressions and then relating the resulting equations to that in Theorem 3, we find, after some simplifications, four infinite series identities as below:

Theorem 19 ( $\rho < y < 1/2$ ). Letting $\Theta(y)$ be defined as in Theorem 3, the following algebraic identities hold:

$\begin{align} &\text{(a)} &&\sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n}\frac{\bar{ \mathbf{H}}_{{2n}}}{n} && = 6\ln^2\bigg(\frac{2}{\sqrt y+\sqrt{4+y}}\bigg), \\ &\text{(b)} &&\sum\limits_{n = 1}^{\infty}\frac{\; y^n}{(1+y)^{3n}} \binom{3n}{n}\frac{ \mathbf{H}_{{2n}}}{n} && = \frac{{\pi}^2}2+3\Theta(y)+3\ln y\ln(1+y)-\frac92\ln^2(1+y)+6\ln^2\bigg(\frac{2}{\sqrt y+\sqrt{4+y}}\bigg), \\ &\text{(c)} &&\sum\limits_{n = 1}^{\infty}\frac{\; y^n}{(1+y)^{3n}} \binom{3n}{n}\frac{ \mathbf{O}_{{n}}}{n} && = \frac{{\pi}^2}4+\frac32\Theta(y)+\frac32\ln y\ln(1+y)-\frac94\ln^2(1+y) +6\ln^2\bigg(\frac{2}{\sqrt y+\sqrt{4+y}}\bigg), \\ &\text{(d)} &\; &\sum\limits_{n = 1}^{\infty}\frac{\; y^n}{(1+y)^{3n}} \binom{3n}{n}\frac{ \mathbf{H}_{{3n}}}{n} && = \frac{2{\pi}^2}3+ 3\Theta(y)+4\ln y\ln(1+y)-5\ln^2(1+y) +6\ln^2\Big(\frac{2}{\sqrt y+\sqrt{4+y}}\Big)- \mathrm{Li}_2\Big(\frac{1}{1+y}\Big) . \end{align}$

By assigning five particular values for $y$ in Theorems 3 and 19, five classes of infinite series identities are displayed in the following corollaries:

Corollary 20 ( $y\to\sqrt5-2$ in Theorems 3 and 19).

$\begin{align} &\text{(a) }\; \sum\limits_{n = 1}^{\infty}\Big(\frac18\Big)^n \binom{3n}{n}\frac{\bar{ \mathbf{H}}_{{2n}}}n && = \frac32 \ln ^2\bigg(\frac{1+\sqrt5}2\bigg), \\ &\text{(b) }\; \sum\limits_{n = 1}^{\infty}\Big(\frac18\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{2n}}}n && = \frac{\pi ^2}2+\frac32 \ln ^22+3\Theta(\sqrt5-2) +6 \ln(1+\sqrt5) \ln\bigg(\frac{1+\sqrt5}4\bigg) , \\ &\text{(c) }\; \sum\limits_{n = 1}^{\infty}\Big(\frac18\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{n}}}n && = \frac{\pi ^2}2+3\Theta(\sqrt5-2) +\frac92 \ln(1+\sqrt5) \ln\bigg(\frac{1+\sqrt5}4\bigg) , \\ &\text{(d) }\; \sum\limits_{n = 1}^{\infty}\Big(\frac18\Big)^n \binom{3n}{n}\frac{ \mathbf{O}_{{n}}}n && = \frac{\pi ^2}4+\frac32 \ln ^22+\frac32\Theta(\sqrt5-2) +\frac{15}4 \ln(1+\sqrt5) \ln\bigg(\frac{1+\sqrt5}4\bigg), \\ &\text{(e)} \; \sum\limits_{n = 1}^{\infty}\Big(\frac18\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{3n}}}n && = \frac{2\pi ^2}3+\frac{11}2 \ln ^22+3\Theta(\sqrt5-2) +\frac{\ln(1+\sqrt5)}2 \ln\bigg(\frac{(1+\sqrt5)^{17}}{2^{38}}\bigg) - \mathrm{Li}_2\bigg(\frac{1+\sqrt5}4\bigg). \end{align}$

Corollary 21 ( $y\to\tfrac{\sqrt5-2}{\sqrt5}$ in Theorems 3 and 19).

$\begin{align} &\text{(a) }\; \sum\limits_{n = 1}^{\infty}\Big(\frac{5}{64}\Big)^n \binom{3n}{n}\frac{\bar{ \mathbf{H}}_{{2n}}}n && = \frac32 \ln ^2\bigg(\frac{5+\sqrt5}{10}\bigg), \\ &\text{(b) }\; \sum\limits_{n = 1}^{\infty}\Big(\frac{5}{64}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{2n}}}n && = \frac{\pi ^2}2+3\Theta\bigg(\frac{\sqrt5-2}{\sqrt5}\bigg) +\frac32 \ln\bigg(\frac{5+\sqrt5}8\bigg) \ln\bigg(64+\frac{128}{\sqrt5}\bigg) +\frac32 \ln ^2\bigg(\frac{5+\sqrt5}{10}\bigg), \\ &\text{(c) }\; \sum\limits_{n = 1}^{\infty}\Big(\frac{5}{64}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{n}}}n && = \frac{\pi ^2}2+3\Theta\bigg(\frac{\sqrt5-2}{\sqrt5}\bigg) +\frac32 \ln\bigg(\frac{5+\sqrt5}8\bigg) \ln\bigg(64+\frac{128}{\sqrt5}\bigg), \\ &\text{(d) }\; \sum\limits_{n = 1}^{\infty}\Big(\frac{5}{64}\Big)^n \binom{3n}{n}\frac{ \mathbf{O}_{{n}}}n && = \frac{\pi ^2}4+\frac32\Theta\bigg(\frac{\sqrt5-2}{\sqrt5}\bigg) +\frac34 \ln\bigg(\frac{5+\sqrt5}8\bigg) \ln\bigg(64+\frac{128}{\sqrt5}\bigg) +\frac32 \ln ^2\bigg(\frac{5+\sqrt5}{10}\bigg) , \\ &\text{(e) }\; \sum\limits_{n = 1}^{\infty}\Big(\frac{5}{64}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{3n}}}n && = \frac{2\pi ^2}3+3\Theta\bigg(\frac{\sqrt5-2}{\sqrt5}\bigg) - \mathrm{Li}_2\bigg(\frac{5+\sqrt5}8\bigg)+\frac32 \ln ^2\bigg(\frac{5+\sqrt5}{10}\bigg)\\ &&&+ \ln\bigg(\frac{5+\sqrt5}8\bigg) \ln\bigg(\frac{8(5+\sqrt5)^7}{5^4}\bigg). \end{align}$

Corollary 22 ( $y\to\tfrac13$ in Theorems 3 and 19).

$\begin{align} &\text{(a)} \; \sum\limits_{n = 1}^{\infty}\Big(\frac{9}{64}\Big)^n \binom{3n}{n}\frac{\bar{ \mathbf{H}}_{{2n}}}n && = \frac32 \ln ^2\bigg(\frac{7+\sqrt{13}}6\bigg), \\ &\text{(b) }\; \sum\limits_{n = 1}^{\infty}\Big(\frac{9}{64}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{2n}}}n && = \frac{\pi ^2}2-\frac32 \ln\Big(\frac43\Big) \ln\Big(\frac{64}3\Big) +3\Theta\Big(\frac13\Big)+\frac32 \ln ^2\bigg(\frac{7+\sqrt{13}}6\bigg), \\ &\text{(c) } \;\sum\limits_{n = 1}^{\infty}\Big(\frac{9}{64}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{n}}}n && = \frac{\pi ^2}2-\frac32 \ln\Big(\frac43\Big) \ln\Big(\frac{64}3\Big) +3\Theta\Big(\frac13\Big), \\ &\text{(d) } \;\sum\limits_{n = 1}^{\infty}\Big(\frac{9}{64}\Big)^n \binom{3n}{n}\frac{ \mathbf{O}_{{n}}}n && = \frac{\pi ^2}4-\frac34\ln\Big(\frac43\Big) \ln\Big(\frac{64}3\Big) +\frac32\Theta\Big(\frac13\Big)+\frac32 \ln ^2\bigg(\frac{7+\sqrt{13}}6\bigg) , \\ &\text{(e) }\; \sum\limits_{n = 1}^{\infty}\Big(\frac{9}{64}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{3n}}}n && = \frac{2\pi ^2}3- \ln\Big(\frac43\Big) \ln\Big(\frac{1024}3\Big) +3\Theta\Big(\frac13\Big)- \mathrm{Li}_2\Big(\frac34\Big) +\frac32 \ln ^2\bigg(\frac{7+\sqrt{13}}6\bigg). \end{align}$

Corollary 23 ( $y\to\tfrac{3 \sqrt{3}-5}2$ in Theorems 3 and 19).

$\begin{align} &\text{(a) }\; \sum\limits_{n = 1}^{\infty}\Big(\frac{2}{27}\Big)^n \binom{3n}{n}\frac{\bar{ \mathbf{H}}_{{2n}}}n && = \frac32 \ln ^2(\sqrt{3}-1), \\ &\text{(b) }\; \sum\limits_{n = 1}^{\infty}\Big(\frac{2}{27}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{2n}}}n && = \frac{\pi ^2}2+\frac32 \ln (54)\ln\Big(\frac23\Big) +6 \ln \big(1+\sqrt{3}\big) \ln \bigg(\frac{1+\sqrt3}2\bigg) +3\Theta\bigg(\frac{3\sqrt3-5}3\bigg), \\ &\text{(c) }\; \sum\limits_{n = 1}^{\infty}\Big(\frac{2}{27}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{n}}}n && = \frac{\pi ^2}2 +\frac32\ln \bigg(\frac{1+\sqrt3}3\bigg)\ln \bigg(\frac{27(1+\sqrt3)^3}4\bigg) +3\Theta\bigg(\frac{3\sqrt3-5}3\bigg) , \\ &\text{(d) }\; \sum\limits_{n = 1}^{\infty}\Big(\frac{2}{27}\Big)^n \binom{3n}{n}\frac{ \mathbf{O}_{{n}}}n && = \frac{\pi^2}4 +\frac34 \ln \bigg(\frac{1+ \sqrt3}2\bigg)\ln \bigg(\frac{(1+ \sqrt3)^5}2\bigg) +\frac34\ln (54)\ln\Big(\frac23\Big) +\frac32\Theta\bigg(\frac{3\sqrt3-5}3\bigg) , \\ &\text{(e) }\; \sum\limits_{n = 1}^{\infty}\Big(\frac{2}{27}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{3n}}}n && = \frac{2\pi^2}3+\frac32 \ln^2\bigg(\frac{1+\sqrt{3}}2 \bigg) - \mathrm{Li}_2\bigg(\frac{1+\sqrt3}3\bigg)+3\Theta\bigg(\frac{3\sqrt3-5}3\bigg)\\ &&& +\ln \bigg(\frac{1+\sqrt{3}}3 \bigg) \ln\bigg(\frac{3^5(1+\sqrt3)^7}{2^4}\bigg). \end{align}$

Corollary 24 ( $y\to 3\sqrt2-4$ in Theorems 3 and 19).

$\begin{align} &\text{(a) }\; \sum\limits_{n = 1}^{\infty}\Big(\frac{2+\sqrt2}{27}\Big)^n \binom{3n}{n}\frac{\bar{ \mathbf{H}}_{{2n}}}n && = \frac32 \ln ^2\bigg(\frac{1+\sqrt3}{2+\sqrt6}\bigg), \\ &\text{(b) }\; \sum\limits_{n = 1}^{\infty}\Big(\frac{2+\sqrt2}{27}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{2n}}}n && = \frac{\pi ^2}2+3\Theta\Big(3\sqrt2-4\Big) +\frac32 \ln ^2\bigg(\frac{1+\sqrt3}{2+\sqrt6}\bigg)\\ &&&+\frac32 \ln\bigg(\frac{1+\sqrt2}3\bigg) \ln\bigg(\frac{27(1+\sqrt2)}2\bigg), \\ &\text{(c) }\; \sum\limits_{n = 1}^{\infty}\Big(\frac{2+\sqrt2}{27}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{n}}}n && = \frac{\pi ^2}2+3\Theta\Big(3\sqrt2-4\Big) +\frac32 \ln\bigg(\frac{1+\sqrt2}3\bigg) \ln\bigg(\frac{27(1+\sqrt2)}2\bigg), \\ &\text{(d) }\; \sum\limits_{n = 1}^{\infty}\Big(\frac{2+\sqrt2}{27}\Big)^n \binom{3n}{n}\frac{ \mathbf{O}_{{n}}}n && = \frac{\pi ^2}4+\frac32\Theta\Big(3\sqrt2-4\Big) +\frac32 \ln ^2\bigg(\frac{1+\sqrt3}{2+\sqrt6}\bigg)\\ &&&+\frac34 \ln\bigg(\frac{1+\sqrt2}3\bigg) \ln\bigg(\frac{27(1+\sqrt2)}2\bigg), \\ &\text{(e) }\; \sum\limits_{n = 1}^{\infty}\Big(\frac{2+\sqrt2}{27}\Big)^n \binom{3n}{n}\frac{ \mathbf{H}_{{3n}}}n && = \frac{2\pi ^2}3+3\Theta\Big(3\sqrt2-4\Big) +\frac32 \ln ^2\bigg(\frac{1+\sqrt3}{2+\sqrt6}\bigg)\\ &&&- \mathrm{Li}_2\bigg(\frac{1+\sqrt2}3\bigg) +\ln\bigg(\frac{1+\sqrt2}3\bigg) \ln\bigg(\frac{3^5(1+\sqrt2)^3}4\bigg). \end{align}$

6. Conclusions and further comments

By integrating Lambert's series and manipulating the cubic transformations for $_3{F_2}$ -series through the "coefficient extraction method", we have shown several algebraic formulae for infinite series containing binomial coefficient $\binom{3n}{n}$ , harmonic-like numbers, and in particular, a free variable "y". As showcases, three classes of infinite series were examined in detail from Sections 3–5. More variants of these series exist and can be treated analogously. For instance, the transformation formula (3.1) under the replacements " $a\to1+6ax, \; c\to2cx$ " becomes

$\begin{aligned} & { }_3 F_2\left[\left.\begin{array}{c} 1+2 a x, \frac{2}{3}+2 a x, \frac{4}{3}+2 a x \\ 1+3 a x+c x, \frac{3}{2}+3 a x-c x \end{array} \right\rvert\, \frac{27 y}{4(1+y)^3}\right] \\ & =\frac{(1+y)^{2+6 a x}}{1-2 y} \times{ }_4 F_3\left[\left.\begin{array}{l} \frac{4}{3}+2 a x, \quad 1+6 a x, 2 c x, 1-2 c x \\ \frac{1}{3}+2 a x, 1+3 a x+c x, \frac{3}{2}+3 a x-c x \end{array} \right\rvert\, \frac{-y}{4}\right] . \\ & \end{aligned}$

Then by carrying out the same procedure as in the preceding sections, we can establish the algebraic identities for five infinite series, as in the theorem below:

Theorem 25 ( $\rho < y < 1/2$ ). Letting $\Lambda(y)$ be the function defined by

$\Lambda(y) = 3\ln(1-2y)+\frac{2-y}{\sqrt{y^2+4y}} \ln\Big(\frac{2-y(2+y+\sqrt{y^2+4y})}{2-y(2+y-\sqrt{y^2+4y})}\Big),$

we have the closed algebraic formulae for five infinite series:

$\begin{align*} &\text{(a) }\; \sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n+1}{n} \bar{ \mathbf{H}}_{{2n+1}} = \frac{2(y-2)(1+y)^2}{(1-2y)\sqrt{y^2+4y}}\ln\frac{\sqrt{4+y}-\sqrt{y}}{2}, \\ &\text{(b) }\; \sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n+1}{n} \mathbf{H}_{{2n+1}} = \frac{(1+y)^2}{1-2y}\Bigg\{\frac{2y-4}{\sqrt{y^2+4y}} \ln\frac{\sqrt{4+y}-\sqrt{y}}{2}-\frac{ \Lambda(y)}{2}\Bigg\}, \\ &\text{(c) }\; \sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n+1}{n} \mathbf{H}_{{n}} = \frac{(1+y)^2}{4y-2} \Lambda(y), \\ &\text{(d) }\; \sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n+1}{n} \mathbf{O}_{{n+1}} = \frac{(1+y)^2}{1-2y}\Bigg\{\frac{2y-4}{\sqrt{y^2+4y}}\ln\frac{\sqrt{4+y}-\sqrt{y}}{2}-\frac{ \Lambda(y)}4\Bigg\}, \\ &\text{(e) }\; \sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n+1}{n} \mathbf{H}_{{3n+1}} = \frac{(1+y)^2}{1-2y}\Bigg\{\frac{2y-4}{\sqrt{y^2+4y}}\ln\frac{\sqrt{4+y}-\sqrt{y}}{2}-\frac{ \Lambda(y)}2+\ln(1+y)\Bigg\}. \end{align*}$

From the algebraic identities exhibited for the series with harmonic-like numbers of the first order, it is natural to ask: what would happen for the corresponding series involving harmonic-like numbers of the second order? Our computations manifested that the resulting expressions are very complex and convoluted, which can be exemplified by the following attempt made by the authors. Extracting the coefficient of $x^2$ across (3.2), we can derive, after a number of reductions and simplifications, the following formulae:

$\begin{align*} \text{(a)}\quad &\sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n} \Big\{\bar{ \mathbf{H}}_{{2n}}^{{2}}+ \mathbf{H}_{{2n}}^{ \langle{{{2}}}\rangle}\Big\} = \frac{6y+6}{2y-1}\Big\{\mathcal{{U}}(y)+\mathcal{{W}}(y)\Big\}, \\ \text{(b)}\quad &\sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n} \Big\{\big( \mathbf{H}_{{3n}}- \mathbf{H}_{{2n}}\big)^2 + \mathbf{H}_{{2n}}^{ \langle{{{2}}}\rangle}- \mathbf{H}_{{3n}}^{ \langle{{{2}}}\rangle}\Big\} = \frac{1+y}{1-2y}\ln^2(1+y), \\ \text{(c)}\quad &\sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n} \Big\{\big(2 \mathbf{O}_{{n}}- \mathbf{H}_{{3n}}\big)^2 +4 \mathbf{O}_{{n}}^{ \langle{{{2}}}\rangle}- \mathbf{H}_{{3n}}^{ \langle{{{2}}}\rangle}\Big\}\\ &\; = \frac{1+y}{1-2y} \Big\{\ln^2(1+y)+6\ln(1+y) \mathrm{W}(y) -10\mathcal{{U}}(y)-12\mathcal{{W}}(y)\Big\}, \\ \text{(d)}\quad &\sum\limits_{n = 1}^{\infty}\frac{y^n}{(1+y)^{3n}} \binom{3n}{n} \Big\{\big( \mathbf{H}_{{n}}- \mathbf{H}_{{3n}}\big)^2 + \mathbf{H}_{{n}}^{ \langle{{{2}}}\rangle}- \mathbf{H}_{{3n}}^{ \langle{{{2}}}\rangle}\Big\}\\ &\; = \frac{1+y}{1-2y} \Big\{\ln^2(1+y)-6\ln(1+y) \mathrm{W}(y)-2\mathcal{{U}}(y)\Big\}; \end{align*}$

where $\mathrm{W}(y)\;\; and\;\; \mathcal{{U}}(y)$ are as in Lemmas 4 and 18, respectively, and $\mathcal{{W}}(y)$ is given by

$\begin{aligned} \mathscr{W}(y) & =\sum_{k=1}^{\infty} \frac{(k-1)!}{\left(\frac{1}{2}\right)_k}\left(\frac{-y}{4}\right)^k \overline{\mathbf{H}}_{2 k} \\ & =\left(1-\sqrt{\frac{y}{4+y}}\right) \ln ^2\left(\frac{\sqrt{4+y}+\sqrt{y}}{2}\right)-\sqrt{\frac{y}{4+y}} \operatorname{Li}_2\left(\frac{2 \sqrt{y}}{\sqrt{y}+\sqrt{4+y}}\right). \end{aligned}$

There are altogether nine harmonic-like numbers appearing in the above four series

$\Big\{ \mathbf{H}_{{n}}^{{2}}, \mathbf{O}_{{n}}^{{2}}, \mathbf{H}_{{3n}}^{{2}};\; \mathbf{H}_{{n}}^{ \langle{{{2}}}\rangle}, \mathbf{O}_{{n}}^{ \langle{{{2}}}\rangle}, \mathbf{H}_{{3n}}^{ \langle{{{2}}}\rangle};\; \mathbf{H}_{{n}} \mathbf{O}_{{n}}, \mathbf{H}_{{n}} \mathbf{H}_{{3n}}, \mathbf{H}_{{3n}} \mathbf{O}_{{n}}\Big\}.$

To evaluate the nine corresponding independent "nuclear series", it is necessary to find five similar independent equations. It seems quite a tough problem that the authors failed to arrive at solutions, even though we did succeed in determining the following long and ugly expression involving six terms in dilogarithms (that prevents us from making further simplifications). Perhaps the only value of this formula lies in its existence.

$\begin{align*} \sum\limits_{n = 1}^{\infty} \binom{3n}{n}\frac{ \mathbf{H}_{{n}}^{ \langle{{{2}}}\rangle}\; y^n}{(1+y)^{3n}} & = \int_0^y\frac{3(1+y)^3\ln\Big(\frac{T(1+y)^3}{y(1+T)^3}\Big)\; dT} {(1-2y)(1+T)(T^2y+Ty^2+3Ty-1)}\\ & = \frac{\pi ^2 (1+y)}{2 (1-2 y)} +\frac{3 (1+y) \ln^2(1+y)}{1-2 y} +\frac{3 (1+y) \ln y \ln(1+y)}{1-2 y}\\ &-\frac{3 (1+y)}{1-2 y} \mathrm{Li}_2\Big(\frac{1}{1+y}\Big) +\frac{9 (1+y)}{1-2 y}\sqrt{\frac{y}{4+y}} \ln(1+y) \ln\Big(\frac{2+y+ \sqrt{4y+y^2}}{2}\Big)\\ &+\frac{3 (1+y)}{2 (1-2 y)} \Big(1+\sqrt{\frac{y}{4+y}}\Big) \Bigg\{3 \mathrm{Li}_2\Big(\frac{-2\sqrt{y}}{(1+y) (\sqrt{y}-\sqrt{4+y})}\Big) -3 \mathrm{Li}_2\Big(\frac{-2 \sqrt{y}}{\sqrt{y}-\sqrt{4+y}}\Big) + \mathrm{Li}_2\Big(\frac{-2 y^{3/2}}{\sqrt{y}(3+y)-(1+y) \sqrt{4+y}}\Big)\Bigg\} \\ &+\frac{3 (1+y)}{2 (1-2 y)} \Big(1-\sqrt{\frac{y}{4+y}}\Big) \Bigg\{3 \mathrm{Li}_2\Big(\frac{-2 \sqrt{y}}{(1+y) (\sqrt{y}+\sqrt{4+y})}\Big) -3 \mathrm{Li}_2\Big(\frac{-2 \sqrt{y}}{\sqrt{y}+\sqrt{4+y}}\Big) + \mathrm{Li}_2\Big(\frac{-2 y^{3/2}}{\sqrt{y}(3+y)+(1+y) \sqrt{4+y}}\Big)\Bigg\} . \end{align*}$

Therefore, to evaluate the series containing binomial coefficient $\binom{3n}{n}$ and harmonic-like numbers of the second and/or higher orders, one needs to find out different approaches. The interested readers are enthusiastically encouraged to make further explorations.

Author contributions

C. Li: Writing, computation and review; W. Chu: Methodology and supervision. Both of authors have read and approved the final version of the manuscript for publication.

Use of AI tools declaration

The authors declare they have not used Artificial Intelligence tools in the creation of this article.

Acknowledgments

The authors express their gratitude to anonymous referees for their careful reading, critical comments and valuable suggestions that make the manuscript improved substantially during the revision.

Conflict of interest

Prof. Wenchang Chu is the Guest Editor of special issue "Combinatorial Analysis and Mathematical Constants" for AIMS Mathematics. Prof. Wenchang Chu was not involved in the editorial review and the decision to publish this article. The authors declare that there are no conflicts of interest regarding the publication of this paper.

References

[1]	Brown, T., Design thinking. Harvard Business Review, 2008, 86(6): 84–92. Available from: https://hbr.org/2008/06/design-thinking.
[2]	Bybee, R.W., The Case for STEM Education: Challenges and Opportunities, 2013, National Science Teachers Association Press.
[3]	Cobern, W.W., Worldview Theory and Science Education Research. In Everyday Thoughts about Nature: A Worldview Investigation of Important Concepts Students Use to Make Sense of Nature with Specific Attention of Science, 2000, 6–12. Springer Netherlands. https://doi.org/10.1007/978-94-011-4171-0_2
[4]	Costa, A., Kallick, B. and Zmuda, A., Students: Habits of Mind Explanation. The Institute for Habits of Mind. 2023. Available from: https://www.habitsofmindinstitute.org/wp-content/uploads/2022/02/Student-HOM-Explanation.pdf.
[5]	Dewey, J., Experience and Education, 1938, New York: Macmillan Company.
[6]	Glaser, B. G., The constant comparative method of qualitative analysis. Social Problems, 1965, 12(4): 436–445.
[7]	Morales, J., What is Inclusive Design? Principles and Examples. Adobe, 2021. Available from: https://xd.adobe.com/ideas/principles/design-systems/what-is-inclusive-design-principles-and-examples/.
[8]	Payling, J., Zein, A. and Zein, P., Disability unplugged, in Inclusive Design: Design for the Whole Population, J. Clarkson, S. Keates, R. Coleman, & C. Lebbon, Eds., 2003, 70–87. Springer London. https://doi.org/10.1007/978-1-4471-0001-0_5
[9]	Saldana, J., The Coding Manual for Qualitative Researchers, 2015, Thousand Oaks, CA: Sage.
[10]	Teo, T.W., Tan, A.L., Ong, Y.S. and Choy, B.H., Centricities of STEM curriculum frameworks: Variations of the S-T-E-M Quartet. STEM Education, 2021, 1(3): 141–156. https://doi.org/10.3934/steme.2021011 doi: 10.3934/steme.2021011
[11]	The World Bank, Disability Inclusion. The World Bank, 2024. Available from: https://www.worldbank.org/en/topic/disability.
[12]	Tommy Edison, What Are The Differences Between Born Blind & Becoming Blind? (Feat. Christine Ha). Tommy Edison, 2017. Available from: https://tommyedison.com/2017/04/05/what-are-the-differences-between-born-blind-becoming-blind-feat-christine-ha/.
[13]	United Nations, 10 Reduce inequality within and among countries. United Nations, 2024. Available from: https://sdgs.un.org/goals/goal10.

This article has been cited by:

1.	Chunli Li, Wenchang Chu, Generating Functions for Binomial Series Involving Harmonic-like Numbers, 2024, 12, 2227-7390, 2685, 10.3390/math12172685
2.	Artūras Dubickas, On the Range of Arithmetic Means of the Fractional Parts of Harmonic Numbers, 2024, 12, 2227-7390, 3731, 10.3390/math12233731

Author's biography Dr. Teo Tang Wee is an Associate Professor in the Natural Sciences and Science Education (Academic Group) in the National Institute of Education. She is also the Co-Head of the Multi-centric Education, Research and Industry STEM Centre in NIE. Her research interests include inclusivity in classrooms and gender issues in STEM education. She is an editorial board member of Asian Women and an Associate Editor of the Asia-Pacific Journal of Science Education, Journal of Curriculum Studies, and the Cultural Studies of Science Education. She is also the Founding Co-Editor-in-Chief of Research in Integrated STEM Education journal; Dr Tan Aik Ling is an Associate Professor of science education with the Natural Sciences and Science Education (Academic Group) at the National Institute of Education, and the Deputy Head for teaching and curriculum matters. She teaches biology education methods courses and courses related to integrated STEM curriculum. Her current research interests lie in the areas of students' ideas of science learning, science teacher professional development and STEM education; Dr Choy Ban Heng is an Assistant Professor in Mathematics Education at the National Institute of Education, Nanyang Technological University, Singapore. He is currently one of the co-Heads for Multi-centric Education, Research, and Industry STEM Centre in NIE. Dr. Choy's research focuses on developing mathematics teachers' expertise in noticing critical mathematical and instructional details during the planning, enactment, and review of lessons. More recently, he is involved in STEM education research and serves as an Associate Editor of Research in Integrated STEM Education journal published by Brill

Reader Comments

Your name:*

Email:*
© 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)