This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Quadrupole Moment of 37K

K. Minamisono minamiso@nscl.msu.edu P. F. Mantica H. L. Crawford J. S. Pinter J. B. Stoker Y. Utsuno R. R. Weerasiri National Superconducting Cyclotron Laboratory, Michigan State University, East Lansing, MI 48824, USA Department of Chemistry, Michigan State University, East Lansing, MI 48824, USA Japan Atomic Energy Agency, Tokai, Ibaraki 319-1195, Japan
Abstract

The electric quadrupole coupling constant of the ground state of 37K(IπI^{\pi} = 3/2+, T1/2T_{1/2} = 1.22 s) in a tetragonal KH2PO4 single crystal was measured to be |eqQ/h|=2.99±0.07|eqQ/h|=2.99\pm 0.07 MHz using the β\beta-ray detecting nuclear quadrupole resonance technique. The electric quadrupole moment of 37K was determined to be |Q(37K)|=10.6±0.4|Q(^{37}{\rm K})|=10.6\pm 0.4 efm2, where the known electric quadrupole coupling constant of stable 39K in the KH2PO4 crystal was used as a reference. The present experimental result is larger than that predicted by shell-model calculations in the sdsd or the sdsd and fpfp model spaces. A possible variation of effective charges was explored to explain the discrepancy.

keywords:
37K , electric coupling constant , quadrupole moment , β\beta-ray detecting nuclear magnetic resonance (β\betaNMR) , β\beta-ray detecting nuclear quadrupole resonance (β\betaNQR) , shell model
PACS:
21.10.Ky , 21.60.Cs , 27.30. +t , 27.40.+z , 29.27.Hj , 29.38.-c

, , , , , , and

thanks: Current address: Technova Corp. Lansing, MI 48906, USA

1 Introduction

The spectroscopic electric-quadrupole moment (QQ) provides a direct measure of the deviation of the charge distribution in a nucleus from spherical symmetry. Proton and neutron effective charges are used in shell-model calculations to obtain quadrupole moments. Within a given shell-model space, the effective charges represent 2ω\hbar\omega, IπI^{\pi} = 2+ excitations of the core nucleons to valence orbits [1] and reflect the virtual excitation of the isoscalar and isovector giant quadrupole resonances. Values of effective charges for sdsd-shell nuclei have been obtained from a systematic analysis of experimental EE2 matrix elements [1, 2] and give reasonable agreement between experimental and theoretical quadrupole moments. However, some nuclei show significant disagreement between experiment and theory. The neutron-rich B isotopes, for example, show a reduction of the neutron effective charge [3] because the loosely-bound valence neutrons, far removed from the core, have less probability to excite the core than well bounded neutrons. Neutron-deficient nuclei, especially those with small proton separation energies, may also be expected to show variation of effective charges. Experimental QQ of such neutron-deficient nuclei are still scarce, even in the sdsd shell. Additional systematic data of QQ of neutron-deficient nuclei are important to further improve our knowledge about the exotic structure of dripline nuclei.

The previously measured quadrupole moment of neutron-deficient 37K, QQ(37K) = 11 ±\pm 4 efm2, was obtained from the hyperfine coupling constant of the 4p24p^{2} P3/2P_{3/2} state of the 37K atom [4]. In the evaluation of Q(37Q(^{37}K), the 4p24p^{2} P3/2P_{3/2} coupling constant of the 39K [5] atom and Q(39Q(^{39}K) evaluated with a calculated electric field gradient [6] were used as reference. The OXBASH shell-model calculation code [7] in the sdsd model space with the USDA interaction [8] gives QQ(37K)theory=7.5 efm2. The experimental value is larger in magnitude, but within 1σ\sigma error of the shell-model result.

Q(37Q(^{37}K) was remeasured using a newly-developed β\beta-ray detecting nuclear quadrupole resonance (β\beta NQR) technique, to improve the precision of the previous measurement and to clarify the possible deviation from shell-model expectations. The new result, Q(37Q(^{37}K) = 10.6 ±\pm 0.4 efm2, is well away from the shell-model result calculated in the sdsd model space. An extension of the shell-model space to include both sdsd and fpfp shells also does not account for the large QQ observed for the 37K ground state. A possible variation of effective charges from the traditional values was explored to explain the discrepancy between theory and experiment.

2 Experiment

The Q(37Q(^{37}K) measurement was performed at National Superconducting Cyclotron Laboratory (NSCL) at Michigan State University. 37K ions were produced from a primary beam of 36Ar accelerated to 150 MeV/nucleon by the coupled cyclotrons and impinging on a 564 mg/cm2 9Be target. A charge pick-up reaction (36Ar + p \rightarrow 37K) was employed to produce 37K because a large nuclear polarization is expected at the peak of the momentum-yield curve of the fragment [9]. The primary beam was steered at an angle of +2+2^{\circ} relative to the normal beam axis at the production target in order to produce the polarized beam. 37K ions were separated from other reaction products by the A1900 fragment separator [10]. The full angular acceptance (±2.5\pm 2.5^{\circ}) of the A1900 was used. An achromatic wedge (1200 mg/cm2 Al) was placed at the second intermediate dispersive image to separate 37K based on relative energy loss in the wedge. The central momentum of the 37K beam was selected with a 1% momentum acceptance to optimize the polarization. A typical counting rate of 37K ions at the experimental apparatus was 50 particles/s/pnA and 40 pnA of primary beam was available. The major contamination in the secondary beam following the A1900 was stable 36Ar, which did not adversely impact the β\beta-NMR measurement.

The polarized 37K ions were delivered to the β\beta-NMR apparatus [11] and implanted into a tetragonal KH2PO4 (KDP) single crystal. The crystal was 2 mm thick and tilted by 45 relative to the β\beta-detector surface. An external magnetic field of H0H_{0} = 0.45 T was applied parallel to the direction of polarization and orthogonal to the beam direction to provide the Zeeman splitting field and to maintain the polarization in KDP. 37K decays to the daughter 37Ar by emitting β+\beta^{+} rays with a half life of 1.226 s. The branching ratio to the ground state (IπI^{\pi} = 3/2+) is 98.2%. The maximum β\beta-ray energy is 5.13 MeV. The asymmetry parameter for direct β\beta decay to the 37Ar ground state is AA = -0.57. β\beta rays from the stopped 37K ions were detected by a set of plastic scintillator telescopes placed at 0(u) and 180(d) relative to the external field direction. The counting rate between u and d counters is asymmetric for a polarized source. The angular distribution,

W(φ)1+APcosφ,W(\varphi)\sim 1+AP\cos{\varphi}, (1)

depends on the polarization PP, AA and the angle φ\varphi between the direction of momentum of decaying β\beta ray and the polarization axis.

The KDP crystal used for implantation had its c-axis set parallel to the external magnetic field. The electric field gradient, qq, at the K lattice site in KDP is axially symmetric (asymmetry parameter η=0\eta=0) and parallel to the c-axis. The magnetic sub-level energies, EmE_{m}, of the 37K implanted in KDP are given in this condition by

Em=gμNH0m+hνQ12(3cos2θ1){3m2I(I+1)},E_{m}=-g\mu_{N}H_{0}m+\frac{h\nu_{Q}}{12}(3\cos^{2}\theta-1)\{3m^{2}-I(I+1)\}, (2)

where mm is the magnetic quantum number, gg is the nuclear gg factor, νQ=3eqQ/{2I(2I1)h}\nu_{Q}=3eqQ/\{2I(2I-1)h\} is the normalized electric coupling constant, where eqQ/heqQ/h is the electric coupling constant, and θ\theta is the angle between the c-axis and external magnetic field. Eq. (2) is given to the first-order of eqQ/heqQ/h, taking the electric interaction as a perturbation to the main magnetic interaction. Since the 37K nucleus implanted into KDP has a large electric interaction, the electric interaction can not be considered as a perturbation to the magnetic interaction. Higher-order terms of the electric interaction have to be considered. However, Eq. (2) and (3) derived later may be used for the system with θ=0\theta=0 and η=0\eta=0, where the first order perturbation calculation gives the exact transition frequencies. The first term in Eq. (2) indicates the 2I+12I+1 magnetic sublevels separated by a fixed energy value due solely to the magnetic interaction (Zeeman splitting). These sublevels are further shifted by the electric interaction and the energy spacing between adjacent sublevels is no longer constant. The 2II separate transition frequencies that appear due to the electric interaction are determined as

fm1m=νLνQ4(3cos2θ1)(2m1).f_{m-1\leftrightarrow m}=\nu_{L}-\frac{\nu_{Q}}{4}(3\cos^{2}\theta-1)(2m-1). (3)

The transition frequencies correspond to the energy difference between two adjacent energy levels in Eq. (2), since the allowed transitions have Δm=±1\Delta m=\pm 1. Here νL=gμNH0/h\nu_{L}=g\mu_{N}H_{0}/h is the Larmor frequency.

The NMR signal is obtained from a double ratio defined as

R=[W(0)W(180)]off/[W(0)W(180)]on1+2AΔP,R=\left[\frac{W(0^{\circ})}{W(180^{\circ})}\right]_{\rm off}\Bigg{/}\left[\frac{W(0^{\circ})}{W(180^{\circ})}\right]_{\rm on}\sim 1+2A\Delta P, (4)

where the subscript, off(on), stands for without(with) RF applied and ΔP\Delta P is the measured change in polarization between RF on and off conditions. The NMR signal can be most effectively searched by applying the 2II transition frequencies in Eq. (3) simultaneously, as a function of the electric coupling constant. Such application ensures total destruction (ΔP=P0\Delta P=P_{0}, where P0P_{0} is the initial polarization produced in the nuclear reaction) of the polarization at the actual electric coupling constant and makes the measurement efficient. The NMR signal for a nucleus with II = 3/2 is 10 times bigger for ΔP=P0\Delta P=P_{0} than that obtained by applying a single transition frequency, where only a partial destruction of P0P_{0} can be achieved. A one day measurement with the multiple frequency technique would be at least a 100 days measurement with the single frequency technique. In this calculation, a linear distribution of populations in EmE_{m} is assumed as amam1=ϵa_{m}-a_{m-1}=\epsilon and m=IIam=1\sum^{I}_{m=-I}a_{m}=1, where ϵ\epsilon is a constant. This multiple frequency version of the β\beta-NMR technique is known as the β\beta-ray detecting nuclear quadrupole resonance (β\beta NQR) technique and is discussed in detail elsewhere [12].

A schematic of the β\beta-NQR system at NSCL is shown in Fig. 1. An RF signal generated by a function generator (labeled FG 1) is selected by a gate (labeled DBM) and sent to an RF amplifier. The amplified signal is then applied to an RF coil, which is part of an LCR resonance circuit, with an impedance matching transformer and one of five variable capacitors and one fixed (\sim 10000 pF) capacitor. After 100 ms irradiation time, the frequency from FG 2 is selected by a gate signal and sent to the same LCR resonance circuit. A different capacitor, which has been tuned to the second frequency to satisfy the LCR resonance condition, is selected by the fast switching relay system. The system ensures sufficient power for any set of three transition frequencies over the expected search region of eqQ/h=eqQ/h= 1.9 \sim 4.0 MHz, so that a high transition probability is achieved. The transition frequency search range spans from 0.465 MHz to 2.465 MHz, as deduced from Eq. (3).

The electric coupling constant of 37K in KDP was searched by applying three transition frequencies to the 37K ions. The three transition frequencies for each electric coupling constant calculated using Eq. (3) were repeatedly applied, in sequence, through the RF coil. Each signal was frequency modulated (FM = ±\pm 20 kHz) to cover a certain region of the electric coupling constant. The signal was frequency modulated at a rate of 50 Hz, and each signal was applied for 100 ms. The oscillating magnetic field strength was \sim 0.6 mT. The 37K ions were continuously implanted into KDP and data were collected for consecutive RF on and off periods of 30 s each.

3 Result

The β\beta-NQR spectrum of 37K in KDP is shown in Fig. 2, where the NMR signal, RR, is plotted as a function of the electric coupling constant. The solid circles are the experimental data and the electric coupling constant of each point corresponds to a set of three transition frequencies, as noted earlier. The horizontal bar on each point is the range of electric coupling constant covered by the frequency modulation of the applied frequencies. The shaded area is the baseline of the resonance obtained from the data outside the resonance region. The solid line is a Gaussian fit to the data. The centroid of the fit gives the electric coupling constant of 37K in KDP:

|eqQ(37K)h|=2.99±0.07(stat.)±0.02(syst.)MHz.\left|\frac{eqQ(^{37}{\rm K})}{h}\right|=2.99\pm 0.07({\rm stat.})\pm 0.02({\rm syst.}){\rm\ \ MHz}. (5)

The dashed line is a theoretical line shape taking into account the frequency modulation of each transition frequency and the distribution of the electric field gradient Δq/q\Delta q/q. The latter broadens the natural line width and may lead to a different centroid. The Δq/q\Delta q/q in KDP at the K lattice site is not known but several percent is expected in ionic crystals [13]. It is difficult to fit the theoretical line shape to data with Δq/q\Delta q/q as a fitting parameter since the data around the resonance are too sparse. Therefore, the theoretical line shape was obtained with an assumption of Δq/q\Delta q/q = 3%, giving reasonable reproduction of the data. The systematic error in the coupling constant reflects the difference between the centroids obtained with the Gaussian and theoretical line shape fits.

The electric coupling constant of 39K in KDP was measured by the conventional NMR technique as |eqQ(39K)/h|=1.69±0.01|eqQ(^{39}{\rm K})/h|=1.69\pm 0.01 MHz [14]. The ratio of the electric quadrupole coupling constants of 37K and 39K in KDP was determined to be:

[eqQ(37K)h/eqQ(39K)h]KDP=Q(37K)Q(39K)=1.77±0.04.\left[\frac{eqQ(^{37}{\rm K})}{h}\Bigg{/}\frac{eqQ(^{39}{\rm K})}{h}\right]_{\rm KDP}=\frac{Q(^{37}{\rm K})}{Q(^{39}{\rm K})}=1.77\pm 0.04. (6)

Q(37Q(^{37}K) can be extracted from Eq. (6) if Q(39Q(^{39}K) is known. Ab initio calculations were performed to obtain the electric field gradient of the 3p53d4s3p^{5}3d4s F9/24{}^{4}F_{9/2} state in the 39K atom [6]. Applying the electric field gradient to the electric coupling constant of the F9/24{}^{4}F_{9/2} state [15], Q(39K)=6.01±0.15Q(^{39}{\rm K})=6.01\pm 0.15 efm2 was deduced. This Q(39Q(^{39}K) was used to deduce Q(37Q(^{37}K) reported here. The method to extract QQ from the measured eqQ/heqQ/h with a precisely calculated qq has been successfully applied to many atomic systems [16]. Other reported values of Q(39Q(^{39}K) obtained from the 4P3/2P_{3/2}, 5P3/2P_{3/2} or 6P3/2P_{3/2} hyperfine coupling constants of the 39K atom [5] with the Sternheimer correction [17], and from molecular microwave spectra combined with calculations of electric field gradient [18] are available as well. The latter gives QQ values, which are consistent with and as precise as the ab initio calculations. The values obtained with Sternheimer correction are subject to rather large error bars due to the accuracy of the correction [5]. However, these values of QQ are in good agreement within quoted errors and do not increase the systematic error of the present result.

Q(37K)Q(^{37}{\rm K}) was extracted from Eq. (6) and Q(39Q(^{39}K) discussed above as

Q(37K)=10.6±0.4efm2.Q(^{37}{\rm K})=10.6\pm 0.4{\rm\ \ efm}^{2}. (7)

The error includes the statistical and systematic errors that resulted from the line shape of the β\beta-NQR spectrum presented in Fig. 2. A temperature dependence of the electric field gradient at the K lattice site in KDP was reported in Ref. [14]. The temperature dependence was not considered for the Q(37Q(^{37}K) reported here, since its contribution is smaller than the uncertainty in the present result.

The present value of Q(37Q(^{37}K) is consistent and more precise than the previous value, Q(37K)Q(^{37}{\rm K}) = 11 ±\pm 4 efm2, which was extracted from the hyperfine coupling constant of the 4p4p P3/22{}^{2}P_{3/2} state of 37K atom, BP3/2(37K)=5.4±1.8B_{P_{3/2}}(^{37}{\rm K})=5.4\pm 1.8 MHz reported in Ref. [4], with the 4p4p P3/22{}^{2}P_{3/2} coupling constant of 39K atom, BP3/2(39K)=2.83±0.13B_{P_{3/2}}(^{39}{\rm K})=2.83\pm 0.13 MHz [5] and Q(39Q(^{39}K) = 6.01 ±\pm 0.15 efm2 evaluated with a calculated electric field gradient [6].

4 Discussion

Theoretical calculations were performed using the OXBASH shell-model code [7] in the sdsd-model space with the USDA interaction [8] and Woods-Saxson single-particle wave functions. Theoretical quadrupole moments were calculated with Qtheory=epQ0p+enQ0nQ_{\rm theory}=e_{p}Q_{0}^{p}+e_{n}Q_{0}^{n}, where Q0p(n)Q_{0}^{p(n)} is the bare quadrupole moment of the proton (neutron) and epe_{p} and ene_{n} are the effective charges for the proton and neutron, respectively. Values of epe_{p} = 1.3ee and ene_{n} = 0.5ee were used for the shell-model results. These effective charges are typical for sdsd shell nuclei [1, 2]. The experimental and theoretical results are summarized in Table 1 and shown in Fig. 3 for the K isotopes (ZZ = 19) and their mirror nuclei (the NN = 19 isotones). Calculations with the USD [20] or the USDB [8] interactions gave results within the expected theoretical uncertainty and therefore only the results of the calculations performed with USDA interaction are discussed here. The quadrupole moments of the mirror nuclei were obtained by exchanging the bare quadrupole moments of the proton and neutron, assuming charge symmetry. The theoretical result Q(37Q(^{37}K)theory = 7.5 efm2 is significantly smaller than the present experimental value [Eq. (7)]. The theoretical ratio [Q(37Q(^{37}K)/Q(39Q(^{39}K)]theory = 1.1 also disagrees with the experimental result [Eq.(6)] as well. Finally, the experimental trend of K isotopes is poorly reproduced by the shell model results as depicted in Fig. 3.

The smaller theoretical ratio of Q(37Q(^{37}K) to Q(39Q(^{39}K) than deduced from experiment is understood as follows. An increase in Q(37Q(^{37}K)theory is expected when two neutrons are removed from 39K. This is realized from the Q0nQ_{0}^{n} value reported in Table 1. However, a decrease in Q0pQ_{0}^{p} is also calculated between 39K and 37K, which is nearly half the increase in Q0nQ_{0}^{n}. The reduction of Q0pQ_{0}^{p} is attributed to two neutrons coupled to 2+ in the d3/2d_{3/2} shell. The wave function of the ground state of 37K is approximated by

|K37(3/2+)\displaystyle\left|{}^{37}{\rm K}(3/2^{+})\right> α|[πd3/21×ν(38Ca:0+)]3/2\displaystyle\sim\alpha\left|\left[\pi d^{-1}_{3/2}\times\nu(^{\rm 38}{\rm Ca}:0^{+})\right]^{3/2}\right>
+β|[πd3/21×ν(38Ca:2+)]3/2,\displaystyle+\beta\left|\left[\pi d^{-1}_{3/2}\times\nu(^{\rm 38}{\rm Ca}:2^{+})\right]^{3/2}\right>, (8)

where Ca38:0+{}^{38}{\rm Ca}:0^{+} and Ca38:2+{}^{38}{\rm Ca}:2^{+} represent the ground and first excited states, respectively. Eq. (4) can be rewritten as |K37(3/2+)α|1+β|2\left|{}^{37}{\rm K}(3/2^{+})\right>\sim\alpha\left|1\right>+\beta\left|2\right>. The amplitudes of these states, α2\alpha^{2} and β2\beta^{2}, are calculated to be 0.76 and 0.16, respectively. The |K37(3/2+)\left|{}^{37}{\rm K}(3/2^{+})\right> covers \sim92% of all the possible wave functions in the sdsd model space. The contribution from each term to the Q0pQ_{0}^{p} can be evaluated as α21Q^p1=1×α2d3/21Q^pd3/21\alpha^{2}\left<1||\widehat{Q}_{p}||1\right>=1\times\alpha^{2}\left<d_{3/2}^{-1}||\widehat{Q}_{p}||d_{3/2}^{-1}\right> and β22Q^p2=0.6×β2d3/21Q^pd3/21\beta^{2}\left<2||\widehat{Q}_{p}||2\right>=-0.6\times\beta^{2}\left<d_{3/2}^{-1}||\widehat{Q}_{p}||d_{3/2}^{-1}\right>, where Q^p\widehat{Q}_{p} is the quadrupole operator. The factors 1 and -0.6 in the righthand side of each equation are the phase factors given by the initial and final state spins. The opposite sign of the 2Q^p2\left<2||\widehat{Q}_{p}||2\right> term to that of the 1Q^p1\left<1||\widehat{Q}_{p}||1\right> term reduces the Q0pQ^{p}_{0} of 37K compared to that of 39K, where the 1Q^p1\left<1||\widehat{Q}_{p}||1\right> term is dominant. The neutron configuration in 37K counteracts the Q0pQ_{0}^{p} and the reduction results in the small increase of the overall Q(37Q(^{37}K)theory relative to Q(39Q(^{39}K)theory. Consequently, the experimental ratio, namely the trend in QQ, would be difficult to reproduce in this framework of the shell model, even if other interactions are considered.

Cross shell (2p-2h) excitations were considered for protons and neutrons with the Monte Carlo Shell Model (MCSM) in the sdsd and fpfp model space with the SDPF-M interaction [21]. The fractions of valence protons and neutrons excited to the fpfp shell in the calculation are 0.29 and 0.28 for 37K and 0.32 and 0.42 for 39K, respectively. The resulting [Q(37Q(^{37}K)/Q(39Q(^{39}K)]theory = 1.1 and QQ(37K)theory = 7.47 efm2 are still significantly smaller than experimental values. The MCSM results are surprisingly similar to the OXBASH result ( [Q(37Q(^{37}K)/Q(39Q(^{39}K)]theory = 1.1 and QQ(37K)theory = 7.52 efm2, respectively), which were limited to the sdsd model space. Although \sim 30% of all the nucleons in the MCSM calculation are excited to the fpfp shell, the Q(37Q(^{37}K)theory is not increased, compared to that calculated in the sdsd model space. The contribution of the 2p-2h excitations to QQ(37K)theory can therefore be considered negligible.

The experimental result may be better reproduced by the shell-model calculations discussed above with an adjustment of the effective charges. A quadrupole deformation of the valence orbitals will lead to a self-consistent polarization of the core. The deformation is described as the coupling of 2ω\hbar\omega, IπI^{\pi}=2+ excitations of the core to the valence orbits [1] and is represented within the model space as effective charges for protons and neutrons. As discussed in the introduction, there is evidence of reduced effective neutron charges in the B isotopes near the neutron dripline [3]. But no precedent has been established for changes in effective charges in neutron-deficient nuclei like 37K.

The experimental ratio between Q(37Q(^{37}K) and Q(39Q(^{39}K) may be used as a guide in estimating effective charges. The ratio can be represented as:

Q(37K)Q(39K)=Q0p(37K)Q0p(39K)+enepQ0n(37K)Q0p(39K),\frac{Q(^{37}{\rm K})}{Q(^{39}{\rm K})}=\frac{Q_{0}^{p}(^{37}{\rm K})}{Q_{0}^{p}(^{39}{\rm K})}+\frac{e_{n}}{e_{p}}\frac{Q_{0}^{n}(^{37}{\rm K})}{Q_{0}^{p}(^{39}{\rm K})}, (9)

where the NN = 20 neutron core in 39K couples to 0+. The ratio of the effective charges can be extracted as en/epe_{n}/e_{p}\sim 1.14, using Eq. (6) and the bare quadrupole moments reported in Table 1. The typical effective charges ep=1.3ee_{p}=1.3e, en=0.5ee_{n}=0.5e give en/ep0.4e_{n}/e_{p}\sim 0.4. An increase in ene_{n} and/or a decrease in epe_{p} is required to reproduce the experimental Q(37Q(^{37}K)/Q(39Q(^{39}K) [Eq. (6)].

The suggestion of alternate effective charges is also seen in the experimental quadrupole moments of mirror nuclei near 37K. 39K and 39Ca have one proton hole and one neutron hole, respectively, in doubly-magic 40Ca. The simple shell-model configurations of these nuclei give the effective charge as the ratio of experimental and bare quadrupole moments for protons (39K) and neutrons (39Ca). The effective charges obtained in this approach are epe_{p} \sim 1.13ee and ene_{n} \sim 0.71ee. The increase in ene_{n} and decrease in epe_{p}, relative to the typical values, lead to [Q(37Q(^{37}K)/Q(39Q(^{39}K)]theory1.3{}_{\rm theory}\sim 1.3 and Q(37Q(^{37}K)theory = 7.9 efm2. These values are increased relative to the shell-model results with the typical values of effective charges, but still remain smaller than the experimental values. The QQ evaluated with these new effective charges derived from the AA = 39, TT = 1/2 mirror partners are summarized in Table 1.

The driving force for the variation of effective charges may lie in the interaction between the valence neutrons and core protons. The effective charges can be expanded to include the isoscalar, epol(0)e_{\rm pol}^{(0)}, and isovector, epol(1)e_{\rm pol}^{(1)}, polarization charges arising from the virtual excitation of the isoscalar and isovector giant quadrupole resonances of the core:

ep=1e+epol(0)epol(1),en=epol(0)+epol(1).e_{p}=1e+e_{\rm pol}^{(0)}-e_{\rm pol}^{(1)},\ \ \ e_{n}=e_{\rm pol}^{(0)}+e_{\rm pol}^{(1)}. (10)

The epol(0)e_{\rm pol}^{(0)} is well determined from the EE2 matrix elements between low-lying states [2], which is dominated by the isoscalar component, and can be derived from typical effective charges as epol(0)0.4ee_{\rm pol}^{(0)}\sim 0.4e. Therefore, using en/ep1.14e_{n}/e_{p}\sim 1.14 obtained by considering the ratio of experimental QQ of 37K and 39K given in Eq. (9), a value epol(1)0.6ee_{\rm pol}^{(1)}\sim 0.6e is determined, much larger than the value epol(1)0.1ee_{\rm pol}^{(1)}\sim 0.1e from the typical effective charges, ep=1.3ee_{p}=1.3e and en=0.5ee_{n}=0.5e. The large epol(1)e_{\rm pol}^{(1)} indicates stronger coupling to the isovector quadrupole giant resonance. It is also noted here that such a large epol(1)e_{\rm pol}^{(1)} has also been discussed for the AA = 51 mirror nuclei [22], which reside in the fpfp shell.

An alternative way to evaluate QQ is by a collective model approach. The ground state wave function of 37K given in Eq. (4) may also be considered as a single-hole configuration in the d3/2d_{3/2} shell interacting with the 2+ vibrational state of the 38Ca core. The influence of particle-core coupling can be calculated as a correction to the QQ by taking into account the effect induced by the collective EE2 operator. A simple relation between QQ in the coupled particle-core system and the QQ of the single particle configuration can be derived [23] as:

Q(37K;3/2+)\displaystyle Q(^{37}{\rm K};3/2^{+}) =Q(|πd3/21;3/2+)epeff\displaystyle=Q(\left|\pi d^{-1}_{3/2}\right>;3/2^{+})e_{p}^{\rm eff} (11)
=Q(|πd3/21;3/2+){1+2r24π15rdVdrB(E2;01+21+)ZeR02ω2}.\displaystyle=Q(\left|\pi d^{-1}_{3/2}\right>;3/2^{+})\left\{1+\frac{2}{r^{2}}\frac{4\pi}{15}\left<r\frac{dV}{dr}\right>\frac{B(E2;0^{+}_{1}\rightarrow 2^{+}_{1})}{ZeR_{0}^{2}\hbar\omega_{2}}\right\}.

The neutron-dependent proton effective charge, epeffe_{p}^{\rm eff}, can be derived from the experimental B(E2)B(E2) and excitation energy ω2=E21+\hbar\omega_{2}=E_{2_{1}^{+}} of the first-excited 21+2^{+}_{1} state in the 38Ca core. Eq. (11) can be evaluated using the values, Q(|πd3/21;3/2+)5.2Q(\left|\pi d^{-1}_{3/2}\right>;3/2^{+})\sim 5.2 fm2, r2d3/2\left<r^{2}\right>_{d_{3/2}} \sim 13 fm2, rdV/dr\left<rdV/dr\right>\sim 40 MeV, E21+E_{2^{+}_{1}} = 2.21 ±\pm 0.01 MeV and B(E2;01+21+)=96±21e2B(E2;0_{1}^{+}\rightarrow 2_{1}^{+})=96\pm 21e^{2}fm4 [24], to obtain Q(37K)=8.7±0.8Q(^{37}{\rm K})=8.7\pm 0.8 efm2. This collective approach using empirical data gives a larger value of QQ than the shell model. The agreement with the present experimental result is improved, however, the error is large. A B(E2)B(E2) with higher precision would be useful for further discussion.

The large experimental Q(37Q(^{37}K) and relatively small proton separation energy of 37K, Sp=1.858S_{p}=1.858 MeV, would be an indication of extended valence proton orbit or even the development of a proton halo [25, 26]. In such a proton halo system, epe_{p} should converge to 1ee due to a smaller probability of the valence protons to polarize the core. Though the experimental Q(37Q(^{37}K)/Q(39Q(^{39}K) may indicate a smaller epe_{p}, a corresponding increase in ene_{n} is essential to better reproduce Q(37Q(^{37}K) and the ratio Q(37Q(^{37}K)/Q(39Q(^{39}K) as discussed above. Also, the mean-square charge radius of 37K relative to that of 39K is known [4] and does not show evidence of an extended distribution of protons. Therefore, the large Q(37Q(^{37}K), the small epe_{p} and the small proton separation energy of 37K are not attributed to an extended proton orbit or a proton halo. The angular momentum of the valence nucleon (d3/2d_{3/2}) and the large atomic number (ZZ = 19) are expected to suppress an extended proton orbit in 37K, due to large centrifugal and Coulomb barriers, respectively.

5 Summary

A newly-developed β\beta-NQR system was used at NSCL to measure the charge distribution of the neutron-deficient nucleus 37K. The electric quadrupole coupling constant of 37K, implanted in a tetragonal KDP single crystal, was measured to be |eqQ/h||eqQ/h| = 2.99 ±\pm 0.07 MHz. Together with the known electric quadrupole coupling constant of 39K in KDP and quadrupole moment of 39K, the quadrupole moment of 37K was deduced as |Q(37K)||Q(^{37}{\rm K})| = 10.6 ±\pm 0.4 efm2. The present result is consistent with and significantly more precise than the previous value (11 ±\pm 4 efm2), evaluated from the hyperfine coupling constant measured by laser spectroscopy. Shell-model results in the sdsd or the sdsd and fpfp model spaces are not able to reproduce the large Q(37Q(^{37}K) value. A increased neutron effective charge and a decreased proton effective charge relative to the typical values in the sdsd shell are required for a better agreement between experiment and theory. Such changes in the effective charges for the QQ of neutron-deficient 37K indicate that the well bound valence neutrons strongly polarize the core protons. More substantial coupling to the isovector giant resonance is inferred, beyond that indicated by the typical effective charges for the sdsd shell nuclei.

Acknowledgement

This work was supported in part by the National Science Foundation, Grant PHY06-06007. The authors would like to express their thanks to the NSCL operations staff for providing the primary and secondary beams and the electronics group for their help in developing the RF system.

References

  • [1] H. Sagawa and B. A. Brown, Nucl. Phys. A 430 (1984) 84.
  • [2] B. A. Brown and B. H. Wildenthal, Ann. Rev. Nucl. Part. Sci. 38 (1988) 29.
  • [3] H. Ogawa et al., Phys. Rev. C 67 (2004) 064308.
  • [4] J. A. Behr et al., Phys. Rev. Lett. 79 (1997) 375.
  • [5] E. Arimondo et al., Rev. Mod. Phys. 49 (1977) 31.
  • [6] D. Sundholm et al., J. Chem. Phys. 98 (1993) 7152.
  • [7] B. A. Brown et al., MSU-NSCL Report No. 1289.
  • [8] B. A. Brown and W. A. Richter, Phys. Rev. C 74 (2006) 034315.
  • [9] D. E. Groh et al., Phys. Rev. Lett. 90 (2003) 202502.
  • [10] D. J. Morrissey, Nucl. Instrum. Methods Phys. Res. B 126 (1997) 316.
  • [11] P. F. Mantica et al., Phys. Rev. C 55 (1997) 2501.
  • [12] K. Minamisono et al., Nucl. Instrum. Methods Phys. Res. A in press.
  • [13] T. Minamisono, Y. Nojiri and K. Sugimoto, Nucl. Phys. A 236 (1984) 416.
  • [14] J. Seliger and V. Zagar, Phys. Rev. B 49 (1994) 14918.
  • [15] G. Spprott and R. Novick, Phys. Rev. Lett. 21 (1968) 336.
  • [16] M. Tokman et al., Chem. Phys. Lett. 291 (1998) 414, and references therein.
  • [17] R. M. Sternheimer and R. F. Peierls, Phys. Rev. A 3 (1971) 837.
  • [18] V. Kello and A. J. Sadlej, Chem. Phys. Lett. 292 (1998) 403.
  • [19] N. J. Stone, At. Data Nucl. Data Tables 90 (2005) 75.
  • [20] B. H. Wildenthal, Phys. Rev. C 28 (1983) 1343.
  • [21] Y. Utsuno et al., Phys. Rev. C 60 (1999) 054315.
  • [22] R. du Rietz et al., Phys. Rev. Lett. 93 (2004) 222501.
  • [23] G. Neyens et al., Nucl. Phys. A 625 (1997) 668.
  • [24] P. D. Cottle et al., Phys. Rev. C 60 (1999) 031301.
  • [25] T. Minamisono et al., Phys. Rev. Lett. 69 (1992) 2058.
  • [26] T. Sumikama et al., Phys. Rev. C 74 (2006) 024327.
Table 1: Experimental and theoretical quadrupole moments of the ZZ = 19 isotopes and NN = 19 isotones. The Q0Q_{0} values, calculated with shell model using the USDA interaction, are bare quadrupole moments without effective charges. Typical effective charges, epe_{p} = 1.3ee and ene_{n} = 0.5ee in the sdsd shell and effective charges obtained from QexpQ_{\rm exp} for the AA = 39 mirror pair, epe_{p} = 1.13ee and ene_{n} = 0.71ee, were used to calculate QtheoryQ_{\rm theory}.
nucleus Q0pQ^{p}_{0} Q0nQ^{n}_{0} QtheoryQ_{\rm theory} (efm2) Qexp.Q_{\rm exp.} (efm2) Ref.
epe_{p}, ene_{n} epe_{p}, ene_{n}
1.3ee, 0.5ee 1.13ee, 0.71ee
39K 5.32 0 6.92 6.01 6.01 ±\pm 0.15 [6]
37K 3.94 4.80 7.52 7.86 10.6 ±\pm 0.4 present
35K 4.71 2.91 7.58 7.39 -
39Ca 0 5.32 2.66 3.78 3.8 ±\pm 0.65 [19]
37Ar 4.80 3.94 8.21 8.22 7.6 ±\pm 0.9 [19]
35S 2.91 4.71 6.14 6.63 4.71 ±\pm 0.09 [19]
Refer to caption
Figure 1: Schematic illustration of the β\beta-NQR system at NSCL. One of the RF signals generated by the function generators is selected by a gate and sent to the amplifier. The amplified signal is then applied to an RF coil through an LCR resonance circuit. One of five variable capacitors and a fixed capacitor is selected by a fast relay switch to satisfy the resonance condition. A pulse pattern generator, RPV-071, controls the function generators, the timing of DBM gating, the fast switching relay system and the pulsing of the short-lived radioactive beam.
Refer to caption
Figure 2: Quadrupole resonance spectrum of 37K in KDP. The NMR signal is plotted as a function of the electric coupling constant. The solid circles are the data and the x-axis value of each point corresponds to a set of three transition frequencies. The horizontal bar of each point is the range of electric coupling constant covered by the frequency modulation of the applied frequencies. The shaded area is the base line of the resonance obtained from the data outside the resonance region. The solid line is a Gaussian fit to the data and the dashed line is the theoretical line shape.
Refer to caption
Figure 3: Experimental and theoretical quadrupole moments of the K isotopes in the sdsd shell and their mirror nuclei. The solid circles and line are the data and theory of the K isotopes, respectively. The open circles and dashed line are the data and theory of their mirror partners, respectively.