Total: 1
The Born-Oppenheimer approximation is widely used to investigate the properties of hydrogen-like systems and doubly heavy hadrons. However, the extent to which this approximation reliably captures the true features of such systems remains an open question. In this work, we adopt the results obtained with the Gaussian expansion method as a benchmark to assess the validity of the Born-Oppenheimer approximation in hadronic systems. We also investigate the dependence of the Born-Oppenheimer approximation results on the choice of trial wave functions. A comprehensive study of the Born-Oppenheimer approximation is carried out by performing calculations using Slater-type functions and Gaussian-type functions as trial wave functions, and by comparing the resulting predictions with those obtained from the Gaussian expansion method. We find that the calculations performed within the Born-Oppenheimer approximation are close to those obtained with the Gaussian expansion method when the heavy-quark mass is relatively small. However, as the heavy-quark mass increases, calculations employing Slater-type functions yield larger values than those from the Gaussian expansion method, whereas those using Gaussian-type functions lead to smaller ones. The use of Slater-type functions generally leads to an enhanced binding energy. The underestimation observed in Born-Oppenheimer approximation calculations with Gaussian-type functions primarily stems from the neglect of non-adiabatic corrections. This comparative study provides deeper insight into the structure of doubly heavy hadrons and clarifies the applicability and limitations of the Born-Oppenheimer treatment in these systems.