Source author record

M. Bhuyan

M. Bhuyan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

31works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2026arXiv

Anomalous Behavior of Giant Monopole Resonance Energy with Nuclear Matter Incompressibility in the framework of Relativistic Mean Field Formalism and Coherent Density Fluctuation Model

The finite nucleus incompressibility $K^A$ is evaluated using the coherent density fluctuation model with the extended relativistic mean field density. The relativistic energy density functional for nuclear matter is replaced by the local density approximation for finite nuclei. The equation is used to calculate the finite nuclear incompressibility, which is further utilized to evaluate the isoscalar giant monopole excitation (ISGMR) energy $E_M$. This excitation energy is compared with other theoretical calculations and experimental data, wherever available. The results are comparable to the data. In contrast to the general understanding, the $E_M$ of finite nucleus is found to be maximum for the lowest nuclear matter incompressibility $K_{\infty}$, whereas it is minimum for the maximum $K_{\infty}$. These reverse results may be due to the self- and cross-interactions of the vector mesons in the nuclear potential.

preprint2022arXiv

Exploring the ground state bulk and decay properties of the nuclei in superheavy island

The $α$-decay half-lives of 204 superheavy nuclei covering the range $114 \leq Z \leq 126$ have been investigated using the relativistic mean-field model (RMF) for NL3$^*$ parameter set. The ground state bulk properties such as binding energy, quadrupole deformation parameter ($β_2$ ), and root-mean-square charge radii for these nuclei are analyzed. Four different semi-empirical formulae, namely, the universal decay law (UDL), the Viola-Seaborg (VSS) formula, the modified universal decay law (MUDL), and the modified Brown formula (MBrown), are used to obtain the $α$-decay half-lives for the considered nuclei. To examine the applicability of relativistic mean-field model within NL3$^*$ parametrization, the $α$-decay energies, and the half-lives of a few known superheavy nuclei within the range 102 $\leq$ Z $\leq$ 118 are calculated and the results are compared with the experimental data along with the theoretical predictions. The $α$-decay energies ($Q$-values) are estimated from the binding energies of the parent, and daughter from the RMF (NL3$^*$) parameter set. The calculated results are compared with macroscopic-microscopic Finite-Range-Droplet-Model (FRDM), Global Nuclear Mass Model (WS3, WS4), Weizsacker-Skyrme mass model (WS$^*$) predictions, and the experimental data, wherever available. The possible standard deviations are also estimated for experimental and various theoretical predictions. We find a good consistency for the experimental-to-UDL, FRDM-to-UDL, and WS4-to-UDL estimates of the decay energy and corresponding half-life. The present analysis provides the theoretical predictions within the microscopic model for the upcoming experiments on the superheavy region.

preprint2022arXiv

Isospin dependent properties of the isotopic chain of Scandium and Titanium nuclei within the relativistic mean-field formalism

The density-dependent nuclear symmetry energy is directly related to the isospin asymmetry for finite and infinite nuclear systems. It is critical to determine the coefficients of symmetry energy and its related observables as it holds great importance in different areas of nuclear physics, such as analyzing the structure of ground state exotic nuclei and neutron star study. The ground state bulk properties such as nuclear binding energy, quadrupole deformation, two-neutron separation energy, the differential variation of two-neutron separation energy, and root-mean-square charge radius for Scandium (Z = 21), and Titanium (Z = 22) nuclei are calculated. The coherent density fluctuation model is used to estimate the isospin-dependent properties of finite nuclei such as symmetry energy and its surface and volume components from its corresponding value in infinite nuclear matter system. Finally, we correlate the neutron-skin thickness with the coefficient of symmetry energy and its related observables corresponding to these isotopic chains of nuclei. The relativistic mean-field formalism with non-linear NL3 and Relativistic-Hartree-Bogoliubov with density-dependent DD-ME2 interaction parameters are employed in the present analysis. A shape transition is observed from spherical to prolate near N $\geq$ 44 and N $\geq$ 40 for Sc- and Ti- isotopic chain, respectively. Notable signatures of shell and/or sub-shell closures have been found for the magic neutron numbers at N = 20 and 28 for both the isotopic chain using the nuclear bulk and isospin quantities. In addition to these, a few signatures of shell/sub-shell closure are observed near the drip-line region, at N = 34 and 50 by following the surface/isospin dependent observables, namely symmetry energy and its component for both the isotopic chain of odd-A Sc- and even-even Ti- nuclei.

preprint2021arXiv

Appearance of peak in symmetry energy at N = 126 for Pb isotopic chain within relativistic energy density functional

The newly derived relativistic energy density functional [\textcolor{blue}{ Phys. Rev. C \textbf{103}, 024305 (2021)}], which stems from the effective field theory motivated relativistic mean-field (E-RMF) is employed to establish the appearance of peak/kink in the symmetry energy over the isotopic chain of Pb-nuclei. The coherent density fluctuation model parametrization procedure for finite nuclei is adopted here to obtain the relativistic energy density functional at local density. The relativistic energy density functional from E-RMFT takes precedence over the Brückner energy density functional as it accurately predicts the empirical saturation density and binding energy per nucleon $E/A$, so-called 'Coester Band Problem'. Interestingly, using the relativistic energy density functional, it is possible to predict the peak at $N=126$ for recently developed G3 and widely used NL3 parameter sets, which is not observed for Brückener's functional in-spite of using the E-RMF density. From the present analysis, the newly fitted energy density functional is found to be minutely sensitive to the choice of the parameter sets employed.

preprint2021arXiv

Application of the coherent density fluctuation model to study the nuclear matter properties of finite nuclei within the relativistic mean-field formalism

We obtained a density-dependent analytical expression of binding energy per nucleon for different neutron-proton asymmetry of the nuclear matter (NM) with a polynomial fitting, which manifests the results of effective field theory motivated relativistic mean-field (E-RMF) model. This expression has the edge over the Br$\ddot{u}$ckner energy density functional [Phys. Rev. {\bf 171}, 1188 (1968)] since it resolves the Coster-Band problem. The NM parameters like incompressibility, neutron pressure, symmetry energy, and its derivatives are calculated using the acquired expression of energy per nucleon. Further, the weight function calculated by E-RMF densities are folded with calculated NM parameters within coherent density fluctuation model to find the properties of closed/semi-closed-shell even-even $^{16}$O, $^{40}$Ca, $^{48}$Ca, $^{56}$Ni, $^{90}$Zr, $^{116}$Sn, and $^{208}$Pb nuclei. The values obtained for the neutron pressure $P^{A}$, symmetry energy $S^{A}$ and its derivative $L_{sym}^A$ known as slope parameter, lie within a narrow domain whereas there is a large variation in isoscalar incompressibility $K^{A}$ and surface incompressibility $K_{sym}^{A}$ while moving from light to heavy nuclei. The sizable variation in $K^{A}$ and $K_{sym}^{A}$ for light and heavy nuclei depicts their structural dependence due to the peculiar density distribution of each nucleus. A comparison of surface quantities calculated in the present work has also been made with ones obtained via Br$\ddot{u}$ckner energy density functional.

preprint2021arXiv

Isotopic shift and search of magic number in the superheavy region

The ground state bulk properties such as binding energy, root-mean-square radius, pairing energy, nuclear density distributions, and single-particle energies are calculated for the isotopic chain of Ca, Sn, Pb, and Z = 120 nuclei. The relativistic mean-field with recently developed G3, IOPB-1, and Relativistic-Hartree-Bogoliubov with density-dependent DD-ME1 and DD-ME2 parameter sets are used in the present analysis. The respective shifts over the isotopic chain for the structural observables and surface property like symmetry energy are also estimated using a three-point method, which is crucial for the systematic analysis of the shell/sub-shell closure. The calculated results are compared with the available experimental data for various bulk properties, wherever available. A multiple isotopic shifts leads to the shell/sub-shell closure at N = (20 \& 28), (50 \& 82), and 126 for Ca, Sn, and Pb isotopes, respectively, are observed. The analysis also supports the neutron magic at N = 40 and 184 for highly neutron-rich $^{60}$Ca, and $^{304}$120, predicted to be the next double magic beyond $^{208}$Pb, respectively. Observing the occupancy number, we notice the higher neutron orbitals are mostly occupied before the lower one, which causes the kinks at neutron magic with an amalgam in the isotopic chain trend above nuclei. We also notice the correlation between the occupation probabilities and the magicity of a nucleus and vice-versa.

preprint2020arXiv

Consistent Skyrme parametrizations constrained by GW170817

The high-density behavior of the stellar matter composed of nucleons and leptons under $β$~equilibrium and charge neutrality conditions is studied with the Skyrme parametrizations shown to be consistent (CSkP) with the nuclear matter, pure neutron matter, symmetry energy and its derivatives in a set of $11$ constraints [Dutra {\it et al.}, Phys. Rev. C 85, 035201 (2012)]. The predictions of these parametrizations on the tidal deformabilities related to the GW170817 event are also examined. The CSkP that produce massive neutron stars give a range of $11.86~\mbox{km} \leqslant R_{1.4} \leqslant 12.55~\mbox{km}$ for the canonical star radius, in agreement with other theoretical predictions. It is shown that the CSkP are compatible with the region of masses and radii obtained from the analysis of recent data from LIGO and Virgo Collaboration (LVC). A correlation between dimensionless tidal deformability and radius of the canonical star is found, namely, $Λ_{1.4} \approx 3.16\times10^{-6}R_{1.4}^{7.35}$, with results for the CSkP compatible with the recent range of $Λ_{1.4}=190_{-120}^{+390}$ from LVC. An analysis of the $Λ_1\timesΛ_2$ graph shows that all the CSkP are compatible with the recent bounds obtained by LVC. Finally, the universal correlation between the moment of inertia and the deformability of a neutron star, named as the \mbox{$I$-Love} relation, is verified for the CSkP, that are also shown to be consistent with the prediction for the moment of inertia of the \mbox{PSR J0737-3039} primary component pulsar.

preprint2020arXiv

GW170817 constraints analyzed with Gogny forces and momentum-dependent interactions

A set of equations of state obtained from finite-range Gogny forces and momentum-dependent interactions is used to investigate the recent observation of gravitational waves from the binary neutron star merger GW170817 event. For this set of interactions, we have calculated the neutron star tidal deformabilities (related to the second Love number), the mass-radius diagram, and the moment of inertia~($I$). The $I$-Love relation has been verified. We also have found strong correlations among the tidal deformability of the canonical neutron star, its radius, and the derivatives of the nuclear symmetry energy at the saturation density. Most of the obtained results are located within the constraints of the tidal deformabilities extracted from the GW170817 detection.

preprint2019arXiv

Density-dependent van der Waals model under the GW170817 constraint

We propose a density-dependent function for the attractive interaction in the original van der Waals model to correctly describe the flow constraint at the high-density regime of the symmetric nuclear matter. After a generalization to asymmetric nuclear matter, it was also possible to study the stellar matter regime from this new model. The mass-radius relation for neutron stars under $β$-equilibrium is found to agree with recent X-ray observations. The neutron star masses supported against gravity, obtained from some parametrizations of the model, are in the range of $(1.97-2.07)M_{\odot}$, compatible with observational data from the PSR J0348+0432 pulsar. Furthermore, we verify the reliability of the model in predicting tidal deformabilities of the binary system related to the GW170817 neutron star merger event and find a full agreement with the new bounds obtained by the LIGO/Virgo collaboration.

preprint2019arXiv

Effective surface properties of light, heavy, and super-heavy nuclei

Starting from light to superheavy nuclei, we have calculated the effective surface properties such as the symmetry energy, neutron pressure, and symmetry energy curvature using the coherent density fluctuation model. The isotopic chains of O, Ca, Ni, Zr, Sn, Pb, and Z = 120 are considered in the present analysis, which cover nuclei over the whole nuclear chart. The matter density distributions of these nuclei along with the ground state bulk properties are calculated within the spherically symmetric effective field theory motivated relativistic mean field model by using the recently developed IOPB-I, FSUGarnet, and G3 parameter sets. The calculated results are compared with the predictions of the widely used NL3 parameter set and found in good agreement. We observe a few signatures of shell and/or sub-shell structure in the isotopic chains of nuclei. The present investigations are quite relevant for the synthesis of exotic nuclei with high isospin asymmetry including superheavy and also to constrain an equation of state of nuclear matter.

preprint2015arXiv

The structural evolution in transitional nuclei of mass 80 $\leq$ A $\leq$ 132

In this theoretical study, we report an investigation on the behavior of two neutron separation energy, differential variation of the separation energy and the abnormality in nuclear charge radius along the isotopic and isotonic chains of transition nuclei. We have used relativistic mean field formalism with NL3 and NL3$^*$ forces for this present analysis. The study refers to {\it even-even} nuclei such as Zr, Mo, Ru and Pd with $N$ = 40$-$ 86, where a rich collective phenomena such as proton radioactivity, cluster or nucleus radioactivity, exotic shapes, {\it Island of Inversion} and etc. are observed. These non-monotonic aspects over the isotopic chain are mainly correlated with the structural properties like shell/sub-shell closures, shape transition, clustering and magicity etc. In addition to these, we have shown the internal configuration of these nuclei to get a further insight into the reason for these discrepancies.

preprint2014arXiv

Search of double shell closure in the superheavy nuclei using a simple effective interaction

This paper refers to an another attempt to search for spherical double shell closure nu- clei beyond Z=82, N=126. All calculations and results are based on a newly developed approach entitled as simple effective interaction. Our results predict the combination of magic nucleus occurs at N=182 (Z=114,120,126). All possible evidences for the oc- currence of magic nuclei are discussed systematically. And, the obtained results for all observables compared with the relativistic mean field theory for NL3 parameter.

preprint2014arXiv

Softness of Sn isotopes in relativistic semi-classical approximation

Within the frame-work of relativistic Thomas-Fermi and relativistic extended Thomas-Fermi approximations, we calculate the giant monopole resonance (GMR) excitation energies for Sn and related nuclei. A large number of non-linear relativistic force parameters are used in this calculations. We find that a parameter set is capable to reproduce the experimental monopole energy of Sn isotopes, when its nuclear matter compressibility lies within $210-230$ MeV, however fails to reproduce the GMR energy of other related nuclei. That means, simultaneously a parameter set can not reproduce the GMR values of Sn and other nuclei.

preprint2014arXiv

The effect of self interacting isoscalar-vector meson on finite nuclei and infinite nuclear matter

A detailed study is made for the nucleon-nucleon interaction based on relativistic mean field theory in which the potential is explicitly expressed in terms of mass and the coupling constant of the meson fields. A unified treatment for self-coupling of isoscalar-scalar $σ-$, isoscalar-vector $ω$-mesons and their coupling constant are given with a complete analytic form. The present investigation is focused on the effect of self-interacting higher order $σ$ and $ω$ field on nuclear properties. An attempt is made to explain the collapsing stage of nucleon by higher order $ω$-field. Both infinite nuclear matter and the finite nuclear properties are included in the present study to observe the behaviour or sensitivity of this self interacting terms.

preprint2013arXiv

Effect of isospin asymmetry in nuclear system

The effect of $δ-$ and $ω-ρ-$meson cross couplings on asymmetry nuclear systems are analyzed in the frame-work of an effective Field theory motivated relativistic mean field formalism. The calculations are done on top of the G2 parameter set, where these contributions are absent. We calculate the root mean square radius, binding energy, single particle energy (for the $1^{st}$ and last occupied orbits), density and spin-orbit interaction potential for some selected nuclei and evaluate the $L_{sym}-$ and $E_{sym}-$ coefficients for nuclear matter as function of $δ-$ and $ω-ρ-$meson coupling strengths. As expected, the influence of these effects are negligible for symmetry nuclear system and these effects are very important for systems with large isospin asymmetry.

preprint2013arXiv

Gravitational wave from rotating neutron star

Using the nuclear equation of states for a large variety of relativistic and non-relativistic force parameters, we calculate the static and rotating masses and radii of neutron stars. From these equation of states, we also evaluate the properties of rotating neutron stars, such as rotational and gravitational frequencies, moment of inertia, quadrupole deformation parameter, rotational ellipcity and gravitational wave strain amplitude. The estimated gravitational wave strain amplitude of the star is found to be $\sim 10^{-23}$.

preprint2013arXiv

Simple effective interaction: Infinite nuclear matter and finite nuclei

The mean field properties and equation of state for asymmetric nuclear matter are studied by using a simple effective interaction which has a single finite range Gaussian term. The study of finite nuclei with this effective interaction is done by means of constructing a quasilocal energy density functional for which the single particle equations take the form of Skryme-Hartree-Fock equations. The predictions of binding energies and charge radii of spherical nuclei are found to be compatible with the results of standard models as well as experimental data.

preprint2013arXiv

The evaporation residue in the fission state of Barium nuclei within relativistic mean-field theory

The evaporation residue of Barium isotopes are investigated in a microscopic study using relativistic mean field theory. The investigation includes the isotopes of Barium from the valley of stability to exotic proton-rich region. The ground as well as neck configurations for these nuclei are generated from their total nucleonic density distributions of the corresponding state. We have estimated the constituents (number of nucleons) in the elongated neck region of the fission state. We found the $α$-particle as the constituent of neck of Ba-isotopes, referred to as the evaporated residue in heavy-ion reaction studies. A strong correlation between the neutron and proton is observed throughout the isotopic chain.

preprint2013arXiv

The Oxygen core inside the Magnesium isotopes

We have studied the ground state bulk properties of magnesium isotopes using axially symmetric relativistic mean field formalism. The BCS pairing approach is employed to take care of the pairing correlation for the open shell nuclei. The contour plot of the nucleons distribution are analyzed at various parts of the nucleus, where clusters are located. The presence of an $^{16}$O core along bubble like $α$-particle(s) and few {\it nucleons} are found in the Mg isotopes.

preprint2012arXiv

An effective Nuclear Model: from Nuclear Matter to Finite Nuclei

The momentum and density dependence of mean fields in symmetric and asymmetric nuclear matter are analysed using the simple density dependent finite range effective interaction containing a single Gaussian term alongwith the zero-range terms. Within the formalism developed, it is possible to reproduce the various diverging predictions on the momentum and density dependence of isovector part of the mean field in asymmetric matter. The finite nucleus calculation is formulated for the simple Gaussian interaction in the framework of quasilocal density functional theory. The prediction of energies and charge radii of the interaction for the spherical nuclei compares well with the results of other effective theories.

preprint2012arXiv

Properties of Z=120 nuclei and the α-decay chains of the (292,304)120 isotopes using relativistic and non-relativistic formalisms

The ground state and first intrinsic excited state of superheavy nuclei with Z=120 and N=160-204 are investigated using both non-relativistic Skyrme-Hartree-Fock and the axially deformed Relativistic Mean Field formalisms. We employ a simple BCS pairing approach for calculating the energy contribution from pairing interaction. The results for isotopic chain of binding energy, quadrupole deformation parameter, two neutron separation energies and some other observables are compared with the FRDM and some recent macroscopic-microscopic calculations. We predict superdeformed ground state solutions for almost all the isotopes. Considering the possibility of magic neutron number, two different mode of α-decay chains (292)120 and (304)120 are also studied within these frameworks. The Q_α-values and the half-life T^α_{1/2} for these two different mode of decay chains are compared with FRDM and recent macroscopic-microscopic calculations. The calculation is extended for the α-decay chains of 292120 and 304120 from their exited state configuration to respective configuration, which predicts long half-life T^α_{1/2}(sec.).

preprint2012arXiv

The effect of isoscalar-isovector coupling in infinite nuclear matter

Working on the framework of Relativistic Mean Field theory, we exposed the effect of nonlinear isoscalar-isovector coupling on G2 parameter set on the density dependence of nuclear symmetry energy in infinite nuclear matter. The observables like symmetric energy and few related coefficients are studied systematically. We presented the results of stiff symmetry energy at sub-saturation densities and a soft variation at normal densities. Correlation between the symmetric energy and the isoscalar-isovector coupling parameter fully demonstrated for wide range of density. The work further extended to the octet system and showed the effect of coupling over the equation of state.

preprint2011arXiv

A new microscopic nucleon-nucleon interaction derived from relativistic mean field theory

A new microscopic nucleon-nucleon (NN) interaction has been derived for the first time from the popular relativistic mean field theory (RMFT) Lagrangian. The NN interaction so obtained remarkably relate to the inbuilt fundamental parameters of RMFT. Furthermore, by folding it with the RMFT-densities of cluster and daughter nuclei to obtain the optical potential, it's application is also examined to study the exotic cluster radioactive decays, and results obtained found comparable with the successfully used M3Y phenomenological effective NN interactions. The presently derived NN-interaction can also be used to calculate a number of other nuclear observables.

preprint2011arXiv

Magic nuclei in superheavy valley

An extensive theoretical search for the proton magic number in the superheavy valley beyond $Z=$82 and corresponding neutron magic number after $N=$126 is carried out. For this we scanned a wide range of elements $Z=112-130$ and their isotopes. The well established non-relativistic Skryme-Hartree-Fock and Relativistic Mean Field formalisms with various force parameters are used. Based on the calculated systematics of pairing gap, two neutron separation energy and the shell correction energy for these nuclei, we find $Z=$120 as the next proton magic and N=172, 182/184, 208 and 258 the subsequent neutron magic numbers.

preprint2011arXiv

The α-Decay Chains of the $^{287, 288}115$ Isotopes using Relativistic Mean Field Theory

We study the binding energy, root-mean-square radius and quadrupole deformation parameter for the synthesized superheavy element Z = 115, within the formalism of relativistic mean field theory. The calculation is dones for various isotopes of Z = 115 element, starting from A = 272 to A = 292. A systematic comparison between the binding energies and experimental data is made.The calculated binding energies are in good agreement with experimental result. The results show the prolate deformation for the ground state of these nuclei. The most stable isotope is found to be 282115 nucleus (N = 167) in the isotopic chain. We have also studied Qα and Tα for the α-decay chains of $^{287, 288}$115.

preprint2010arXiv

Application of relativistic mean field and effective field theory density to calculate scattering observables for Ca isotopes

In the frame work of relativistic mean field (RMF) theory, we calculate the density distribution of protons and neutrons for $^{40,42, 44,48}Ca$ with NL3 and G2 parameter sets. The microscopic proton-nucleus optical potential for $p+^{40,42,44,48}Ca$ system is evaluated from Dirac NN-scattering amplitude and the density of the target nucleus using Relativistic-Love-Franey and McNeil-Ray-Wallace parametrizations. Then we estimate the scattering observables, such as elastic differential scattering cross-section, analyzing power and the spin observables with relativistic impulse approximation. We compare the results with the experimental data for some selective cases and found that the use of density as well as the scattering matrix parametrization is crucial for the theoretical prediction.

preprint2010arXiv

Effects of density and parametrization on scattering observables

We calculate the density distribution of protons and neutrons for $^{40,42, 44,48}Ca$ in the frame-work of relativistic mean field (RMF) theory with NL3 and G2 parameter sets. The microscopic proton-nucleus optical potential for $p+^{40}Ca$ system is evaluted from Dirac NN-scattering amplitude and the density of the target nucleus using Relativistic-Love-Franey and McNeil-Ray-Wallace parametrizations. Then we estimate the scattering observables, such as elastic differential scattering cross-section, analysing power and the spin observables with relativistic impulse approximation. We compare the results with the experimental data for some selective cases and found that the use of density as well as the scattering matrix parametrization is crucial for the theoretical prediction.

preprint2010arXiv

Relativistic mean field study of the properties of Z=117 nucleus and the decay chains of $^{293,294}$117 isotopes

We have calculated the binding energy, root-mean-square radius and quadrupole deformation parameter for the recently synthesized superheavy element Z=117, using the axially deformed relativistic mean field (RMF) model. The calculation is extended to various isotopes of Z=117 element, strarting from A=286 till A=310. We predict almost spherical structures in the ground state for almost all the isotopes. A shape transition appears at about A=292 from prolate to a oblate shape structures of Z=117 nucleus in our mean field approach. The most stable isotope (largest binding energy per nucleon) is found to be the $^{288}$117 nucleus. Also, the Q-value of $α$-decay $Q_α$ and the half-lives $T_α$ are calculated for the $α$-decay chains of $^{293}$117 and $^{294}$117, supporting the magic numbers at N=172 and/ or 184.

preprint2009arXiv

Internal stucture of clusters in $^{112-122}$Ba nuclei within relativistic mean field theory

We study the clustering structure and the internal or sub-structure of clusters in $^{112-122}$Ba nuclei within the framework of relativistic mean field theory in an axially deformed cylindrical co-ordinate. We calculate the total density distribution, and the individual neutrons and protons density distributions. From the analysis of the clustering confugurations of the density distributions of various shapes, we find different sub-structures inside the Ba nuclei considered here. The important step, carried out for the first time, is the counting of number of protons and neutrons present in the clustering region(s). $^{12}$C is shown to consitute the cluster configuration of Ba nuclei in most cases, with $^{2,3}$H and $^4$He constituting the neck between two fissioning symmetrical fragments.