## Brazilian Journal of Physics

*Print version* ISSN 0103-9733

### Braz. J. Phys. vol.34 no.3a São Paulo Sept. 2004

#### http://dx.doi.org/10.1590/S0103-97332004000500018

**Coulomb nuclear interference with deuterons in even palladium isotopes **

**M. R. D. Rodrigues ^{I}; C. L. Rodrigues^{I}; T. Borello-Lewin^{I}; L. B. Horodynski-Matsushigue^{I}; J. L. M. Duarte^{I}; G. M. Ukita^{I, II} **

^{I}Instituto de Física, Universidade de São Paulo, Caixa Postal 66318, CEP 05389-970, São Paulo, SP, Brazil

^{II}Faculdade de Psicologia, Universidade de Santo Amaro, Rua Professor Enéas da Siqueira Neto, 340, CEP 04829-300, São Paulo, SP, Brazil

**ABSTRACT**

Angular distributions for the inelastic scattering of 13.0 MeV deuterons on ^{104;106;108;110}Pd were measured with the Sao Paulo Pelletron-Enge-Spectrograph facility in the range of 12^{0} __<__ q_{lab}__<__ 64^{0}. A Coulomb-Nuclear Interference analysis, employing the Distorted Wave Born Approximation with the Deformed Optical Model as transition potential, under well established global optical parameters, was applied to the first quadrupolar excitations. The values of *C* = / , the ratio of charge to isoscalar deformation lengths, and of ()^{2} were extracted through the comparison of experimental and predicted cross section angular distributions. The ratios of reduced charge to isoscalar transition probabilities, B(EL) to B(ISL) respectively, are related to the square of the parameter C and were thus obtained with the advantage of scale uncertainties cancellation. For 104Pd, and preliminary for ^{108}Pd, the respective values of C = 1.18(3) and C = 1.13(4) reveal an enhanced contribution of the protons relative to the neutrons to the excitation, while a smaller effect is found for ^{106}Pd, C = 1.06(3) and for ^{110}Pd, C = 1.07(3), in comparison with the value C = 1.00 expected for homogenous collective excitations.

**1 Introduction**

Inelastic scattering is known to be a powerful tool for studying collective states. In particular, the Coulomb-Nuclear Interference (CNI), using a probe with isoscalar character, allows for the simultaneous extraction of the isoscalar and electric reduced transition probabilities, B(ISL) and B(EL), respectively. In the *A* ~ 100 transitional mass region, due to the relevant role played by the neutrons, an a-priori assumption of the dominance of simple homogenous collective effects (relative contributions of protons and neutrons of about Z/N) is not indicated and a direct access to the reduced isoscalar transition probability B(ISL) (mass), besides B(EL), is demanded. The São Paulo Nuclear Spectroscopy with Light Ions Group developed a research line that uses CNI in the inelastic scattering of isoscalar probes, in particular deuterons, by nuclei, with the aim of following the evolution of the collective behavior throughout isotopic chains [1, 2]. The adopted procedure applies the deformed optical potential model (DOMP) with global optical parameters as the nuclear transition potential in the analysis of the inelastic scattering. It is to be noted that, for the first quadrupolar excitation, the majority of the microscopic form factors do not differ substantially from the macroscopic ones in the important tail region. Through this macroscopic CNI analysis in the DWBA approach, ()^{2} , the square of mass deformation length, is extracted as a scale factor from the fit of the predicted cross sections to the experimental data of the inelastic scattering reaction and, analyzing the characteristic changes in the angular distribution shape, the value of the ratio between charge () and mass () deformation lengths, C, is also obtained. These quantities can be put in correspondence with the value of B(ISL) and of the ratio B(EL)/B(ISL), for which, therefore, a scale uncertainty cancellation occurs, favoring more accurate results. Furthermore, for deuterons there are global optical potential parameter sets, which were extensively tested, allowing for a good control of the free parameters in the analyses.

The present work applies the CNI method to the inelastic scattering of deuterons on even palladium isotopes. Data were obtained for the ^{104,106,108,110}Pd(*d, d*')Pd () reactions at 13.0 MeV in the Pelletron Laboratory, using nuclear emulsion plates on the focal plane of the Enge Split Pole Spectrograph. For these excitations, no results of B(IS2) had been formerly reported in the literature.

**2 Experimental Procedure**

Data for the inelastic scattering of deuterons of 13.0 MeV by ^{104,106,108,110}Pd were obtained, using the Pelletron-Enge-Spectrograph system. The deuteron beam was extracted from a ''MC-SNICS" ion source, provided with a solid titanium deuteride cathode. The beam was focused on the target after passing defining slits of 1.2 mm width and 3.0 mm height (spectrograph object). An adequate profile of the beam is important for a good quality of the data.

Targets of isotopically enriched materials of ^{104,106,108}Pd (98.98%, 98.53%,99%), with thicknesses between 20 and 35 mg/cm^{2}, were produced in São Paulo by the electron bombardment method. Only rather small amounts of the usual contaminations by C, N, O, Si, S, K e W were diagnosed in these targets and also in the ^{110}Pd target (64mg/cm^{2}) produced formerly in Germany by a sputtering procedure.

The scattered deuterons were momentum analyzed by the Enge Magnetic Spectrograph and detected in nuclear emulsion plates (Fuji G6B, 2G and 7D, 50 mm thick), which covered 25 cm along the focal plane of the spectrograph. The greatest advantage of nuclear emulsion is that this detector does not respond to the abundant background, mostly X and g rays from (*n*, g) reactions in the spectrograph iron core after deuteron break-up. Spectra were taken at judiciously chosen scattering angles in a range of 12^{0} __<__q_{lab}__<__ 64^{0}, in order to characterize the CNI in the angular distribution corresponding to the first quadrupolar excitation. After processing, the exposed plates were scanned in strips of 200 mm across the plates. An energy resolution of 7 keV was achieved. For instance, Fig. 1 displays portions of deuteron spectra showing the peaks associated with the state for ^{106}Pd and ^{110}Pd each at two scattering angles. It can be seen the strong influence of elastic tail at 15^{0} as well as the low background at greater angles. At 33^{0} the elastic peaks of the contaminants C, N, O are near to the peak of interest.

The total number of projectiles corresponding to each spectrum was obtained with the help of a calibrated current integrator connected to an aligned Faraday cup with electron suppression. The beam direction was continuously monitored. The relative normalization of the spectra was thus obtained, while the absolute normalization of the cross sections was referred to optical model predictions for the elastic scattering of deuterons on the same target, measured under similar conditions. The elastic spectra were measured at laboratory scattering angles in the range of 30^{0} __<__ q_{lab}__<__ 78 ^{0}.

**3 Analysis**

The excitation of the first quadrupolar state in ^{104,106,108,110}Pd was interpreted in terms of the distorted wave Born approximation (DWBA) - deformed optical potential model (DOMP). The parameters of Perey & Perey [3], were chosen, both, for the absolute normalization of the elastic cross section, and for the generation of the distorted incident and outgoing waves in the reaction analysis and also for the transition potential to the collective state. These parameters are being systematically used by the S. Paulo Group in the CNI studies with deuterons, in this mass region. The DWBA-DOMP calculations were performed with the code DWUCK4 [4] with the usual corrections to account for non-locality effects.

The DWBA-DOMP description of the experimental cross section is:

s* _{exp}*(q) = as

_{DWBA}_{–}

*(q).*

_{DOMP}In the fit of each experimental angular distribution, through the factor a, the square of the mass deformation length is obtained, while the characteristic shape is related to the parameter C.

In Fig. 2, the curves associated with the values of C_{fitted} are results of linked parameters fitting, since C and are correlated. Due to the cross section not being a linear function of C, the procedure chosen for the c^{2} minimization was the iterative method of Gauss [5], that presents a rapid convergence. From the values of C, it is possible to calculate the ratio:

where, *r _{C}* and

*r*are, respectively, the characteristic reduced radii of the charge and the mass distributions of the nucleus. In this work, the equivalent sharp cutoff radii

_{m}*r*= 1.22 fm and

_{C}*r*= 1.16 fm [7], respectively, were used.

_{m}The reduced isoscalar transition probability was therein related to , using the Bernstein [6] definition, through:

with *R _{m}* =

*r*

_{m}A^{1/3}.

Under further assumptions [1], the modulus of the ratio between the quadrupolar moments of the neutron and proton distributions, | *M _{n}*/

*M*|, an indicator of the isospin character of the nuclear transition, can be obtained through the relation:

_{p} The comparison between | *M _{n}*/

*M*| and

_{p}*N/Z*(expected to be approximately equal for a homogeneous collective excitation) is another manner to characterize the role played by the neutrons relative to the protons in the excitations.

**4 Some Results and Discussion**

In the analysis of ^{104,106,110}Pd, the experimental data for the inelastic scattering to the state are discriminative for C and when compared with the predicted angular distributions, as shown in Fig. 2. In order to illustrate the sensibility of the method, two curves associated with values of *C* @ *C _{fitted}* ± 0.2, are also represented in Fig. 2, with the respective value of c

^{2}.

In order to evaluate the adequacy of the method employed to extract the two correlated parameters, C and , and respective statistical uncertainties obtained, a direct statistical test was performed through a Monte Carlo simulation of 5,000 new data sets. This procedure provides equivalent results, for Gaussian data distributions. Each new data set was generated starting from the measured data set, by randomly choosing values around each point from Gaussian distributions with the given standard deviation. So the was determined and then the respective C and values were obtained. Two different views of the c^{2} surface associated with the ^{106}Pd in the space determined by the (C, ) plane is seen in Fig. 3. The c^{2} contour lines obtained are almost perfect ellipses, in very good agreement with outcomes of the Gauss approximation.

The values of C and extracted through the inelastic scattering of deuterons of 13.0 MeV on ^{104}Pd, ^{106}Pd and ^{110}Pd, for the states, are presented in Table 1. The preliminary results for ^{108}Pd are also included. The uncertainties shown are only the statistical ones. The values of the ratios *B*(*E*2)/*B*(*IS*2)e^{2}, of *B*(*IS*2)(*b*^{2}) and of |*M _{n}/M_{p}*| obtained are also shown in Table 1. The uncertainties indicated were calculated by simple propagation.

For ^{104}Pd and preliminary for ^{108}Pd, the extracted C values clearly reveal a non-homogeneous behavior, with an enhanced contribution of the protons relative to the neutrons to the first quadrupolar excitation, while a less pronounced effect is found for ^{106,110}Pd, although the values of C are also slightly higher than 1.00 for both.

Figure 4 shows the values of the *B*(*E*2)/*B*(*IS*2)*e*^{2} ratios, extracted in this work for the four even isotopes of Pd analyzed, in comparison with the results obtained by the São Paulo group in the *A* ~ 100 region, with the CNI methodology applied to the inelastic scattering of isoscalar projectiles, for further consideration. The isotopes ^{94}Mo and ^{98}Mo were studied using deuterons as projectiles [1], while ^{100,102,104}Ru were investigated through inelastic scattering of alpha particles [2]. Comparing the nuclei studied in the present work with their isotones in the Ru chain, formerly investigated, it is to be noted that the nucleus ^{102}Ru [2], in the same manner as the isotone ^{104}Pd does, reveals a more important contribution of protons than of neutrons to the first quadrupolar excitation. Following the Ru chain, although a rather intense non-homogeneity also in favor of the protons is presented by ^{104}Ru, this behavior is in strong contrast to its isotone ^{106}Pd, for which the non-homogeneity is diluted. Following the evolution of *B*(*E*2)/*B*(*IS*2)*e*^{2} ratios along the isotopic chains, it is seen that for the even Ru isotopes there is a growth as N increases, while for the even Pd isotopes a different behavior emerges.

Figure 5 shows the evolution of the B(E2)(*e*^{2}*b*^{2}) values, taken from the comprehensive survey by Raman [8]. The nuclei studied recently by the São Paulo group are circled in the figure for easier identification. It is seen that, between ^{94}Mo, on one side, and ^{104}Ru and ^{110}Pd, on the other, the B(E2)(*e*^{2}*b*^{2}) values increase by more than a factor of 4, demonstrating the well known tendency to deformation, verified in this mass region for *N* ~ 60. As shown by Fig. 4, the values of the *B*(E2)/*B*(*IS*2)*e*^{2} ratio vary from around 1.2 e^{2}(^{106,110}Pd and ^{100}Ru) to approximately 1.6 *e*^{2}(^{104}Pd and ^{104}Ru), therefore the neutron contribution to the first quadrupolar excitations generally lags behind the proton one in the three isotopic chains investigated for *N* ~ 60. No monotonic behavior is evident for this ratio, pointing out that, in this long established region of shape transitions, neutrons may be responsible for some interesting effects outside the simple collective interpretation.

**Acknowledgements**

This work was partially supported by Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP).

**References**

[1] G. M. Ukita, T. Borello-Lewin, L. B. Horodynski-Matsushigue, J. L. M. Duarte, and L. C. Gomes, Phys. Rev. C **64**, 014316 (2001). [ Links ]

[2] L. C. Gomes *et al.*, Phys. Rev. C **54**, 2296 (1996). [ Links ]

[3] C. M. Perey and F. G. Perey, At. Data and Nucl. Data Tables **17**, 1 (1976). [ Links ]

[4] P. D. Kunz, computer code DWUCK4, University of Colorado, 1974. [ Links ]

[5] P. R. Bevington, Data Reduction and Error Analysis for the Physical Sciences (McGraw Hill, New York, 1969). [ Links ]

[6] A. M. Bernstein, V. R. Brown, and V. A. Madsen, Phys. Lett. **71B**, 48 (1977); **103B**, 255 (1981). [ Links ]

[7] H. S. Chung and W. D. Myers, Nucl. Phys. A **513**, 283 (1990). [ Links ]

[8] S. Raman *et al.*, At. Data and Nucl. Data Tables, **78**, 1 (2001). [ Links ]

Received on 11 September, 2003