In the context of an active topical research in laser-related physics,1‒5 the problem of charge emission from bound-states under the action of the electromagnetic radiation is receiving an increasing interest. Some investigations focus especially on the effect the optical-laser radiation may have on the spontaneous alpha-particle decay of the atomic nuclei,6‒11 or nuclear proton emission,12,13 but the area may be extended to atom ionization or molecular or atomic clusters fragmentation.14‒17 The aim of the present paper is to estimate the effect of the adiabatically-applied electromagnetic radiation upon the rate of spontaneous nuclear alpha decay and proton emission. Specifically, the paper is motivated by the interest in computing the rate of tunneling through a Coulomb potential barrier in the presence of electric fields. It is claimed that the rate of alpha decay is practically not affected by electric fields,8 nor it is greatly enhanced by strong electric fields.18 On the other side, the atomic electron cloud may screen appreciably the external electric fields, such that the atomic nucleus may experience, in fact, rather low electric fields. It is this point, related to low electric fields, which may raise technical difficulties in estimating the small effect of these external fields upon the alpha decay.
We adopt a nuclear model with Z protons and A − Z neutrons, where A is the mass number of the nucleus, moving in the nuclear mean field. The experiments proceed usually by placing a collection of heavy atoms in the focal region of a laser beam, and focusing radiation pulses upon that collection of atoms. We consider an optical-laser radiation with a typical frequency of the order 1015s−1 (corresponding to a period and a wavelength ). We assume that the motion of the charges under the action of the electromagnetic radiation remains non-relativistic, i.e. where q is the particle charge, m is the particle mass and A0 is the amplitude of the vector potential (c denotes the speed of light in vacuum). For protons in atomic nuclei this condition yields a very high electric field (1011 electro-static units), which corresponds to a maximum intensity of the laser beam in the focal region of the order Typically, the duration of the laser pulse is of the order of tens of radiation period (or longer), such that we may consider the action of the electromagnetic radiation much longer than the period of the radiation. The repetition rate of the laser pulses is usually much longer than the pulse duration. For simplification we consider linearly-polarized radiation plane waves; the calculations can be extended to a general polarization. The laser-beam shape or multi-mode operation has little relevance upon the results presented here.
The electric fields are appreciably screened by the electronic cloud of the heavy atoms. The screening effects on the thermonuclear reactions, alpha decay and lifetimes have been considered previously.10,11,19,20 A convenient means of treating the electron cloud in heavy atoms is the linearized Thomas-Fermi model.21 According to this model, the radial electron distribution is concentrated at distance (screening distance), where is the Bohr radius and Z is the atomic number (Z ≫ 1); q and m denote the electron charge and mass, respectively. The atomic binding energy depends on R, and the atom exhibits an eigenmode related to the change in R (a breathing-type mode), with an eigenfrequency We recognize in the plasma frequency of a mean electron density It corresponds to the atomic giant-dipole oscillations discussed in Ref.21 In the presence of an external electric field E oriented along the z-direction the electrons are displaced by u(with fixed nucleus), a displacement which produces an energy change By integrating over z, we get a factor in the eigenfrequency , as expected. It follows that the displacement u obeys the equation of motion where the internal field is (polarization and the dipole moment p=Zqu). For the solution of this equation is and the internal field is the total electric field acting upon the atomic nucleus is
(1)
since we may use the approximation (where ); we can see that the total field acting upon the nucleus is appreciably reduced by the electron screening. For Z = 50 this reduction factor is the maximum field 3 × 1013V /cm is reduced to 107V /cm (electrostatic units). It follows that we may limit ourselves here to low fields acting upon the atomic nuclei. The cases of strong fields have been analyzed in Refs.8,9,18,22,23 At the same time, an induced electric field generated by the dipolar eigenmodes occurs inside the atom, which oscillates with the higher eigenfrequency
If the field is low, the bound-state charge oscillates, emits higher-order harmonics of electromagnetic radiation and tunneling may appear; in this latter case, the charge accommodates itself in the field, in a long time, which amounts to an adiabatically-introduced interaction; this regime allows the usual, standard application of the tunneling approach. As we shall see below, the threshold field which separates the two regimes (low-field regime from high-field regime) can be obtained from where a is a distance of the order of the bound-state dimension (a = 10−13cm) (for protons, the threshold field is (102 electrostatic units)).
Originally, the charge emission from bound states, like atom ionization, has been treated by using adiabatic hypothesis, either by time-dependent perturbation theory, or by imaginary-time tunneling, or other equivalent approaches.24‒32 Quasi-classical tunneling through the potential barrier generated by the field has been applied in classical works to static fields and the hydrogen atom (in parabollic coordinates).33‒35 For alpha-particle decay or proton emission the situation is different. First, in spontaneous decay, the alpha particle (and, in general, the ejected charge) is preformed and, second, the tunneling through the Coulomb potential barrier must be included.36‒40 We analyze below the spontaneous charge emission, affected by the presence of an adiabatically-introduced electromagnetic radiation, in the presence of a Coulomb barrier; the problem may exhibit relevance for studies of alpha-particle decay or proton emission.
The standard model of spontaneous alpha decay is based on Bohr’s concept of compound nuclei.41 In an alpha-unstable nucleus the pre-formed alpha particle acquires a kinetic energy and penetrates (tunnels through) the Coulomb potential barrier. Consequently, the alpha-unstable nucleus is in fact in a "metastable state". In this simple model, the spontaneous alpha-particle decay and proton emission proceed by tunneling through the Coulomb potential barrier, as a result of many "attempts" the charge makes to penetrate the barrier. The (high) frequency of this process is of the order 1/ta, where ta corresponds, approximately, to the energy level spacing; in atomic nuclei this spacing, for the relevant energy levels, is of the order which gives41 also, the broadening of the charge energy levels introduces an energy uncertainty (we leave aside the so-called tunneling through the internal potential barrier and the pre-formation factor of the alpha particle). The order of magnitude of the energy of the charge is a few M eV , which ensures a quasi-classical tunneling. The effect of the electromagnetic radiation upon the initial preparation of the charge for tunneling may be neglected. Similarly, we consider a sufficiently low electromagnetic radiation, such that we may neglect its effects upon the mean-field potential. We limit ourselves to the effect of the electromagnetic interaction on the tunneling rate.
Let us consider a charge q >0 with mass m in the potential barrier V(r) in the presence of an electromagnetic radiation with the vector potential where A0 is the amplitude of the vector potential, is the radiation frequency and k is the radiation wavevector the electromagnetic field is transverse, i.e. kA = 0. Since the phase velocity of the non-relativistic charge is much smaller than the speed of light c in vacuum, we may neglect the spatial phase kr in comparison with the temporal phase Kr consequently, the vector potential may be approximated by This approximation amounts to neglecting the effects of the magnetic field. It is assumed that this potential is introduced adiabatically. The charge is immersed in the radiation field, such that we may start with the standard non-relativistic hamiltonian
(2)
where the momentum p includes the electromagnetic momentum qA/c beside the mechanical momentum mv, where v is the velocity of the particle. We consider the Schrodinger equation
(3)
since the interaction is time-dependent we need the time evolution of the wavefunction. Consequently, in equation(3) we perform the well-known Kramers-Henneberger transform42‒45 (with a vanishing electromagnetic interaction for t → −∞)
(4)
the Schrodinger equation becomes
(5)
it is convenient to introduce the electric field we get
(6)
and (7)
We can see that for high-intensity fields the potential (including the mean- field potential) is rapidly vanishing along the field direction. Here we assume that the field intensity is low; specifically we assume where a is the dimension of the region the charge moves in (the atomic nucleus); for protons this inequality means (102 electrostatic units), as stated above. The preformed alpha particle (or emitted proton) may tunnel through the potential barrier given by equation (7); the "attempt" frequency to penetrate the barrier and the energy uncertainty are practically not affected by the low-intensity field.
We adopt a model of nuclear decay by assuming a Coulumb potential barrier (with the center-of-mass of the original nucleus placed at the origin); in the absence of the field the tunneling proceeds from r1 = a to where is the radial energy of the charge; it is convenient to introduce the parameter which includes the effect of the field. In the presence of the field these limits become
(8)
and where a = ar/r. We expand in powers of and get
(9)
where θ is the angle the radius vector r makes with the electric field E0.
To continue, we assume that the free charge attempting to penetrate the potential barrier has momentum pn and kinetic energy where n is a generic notation for its state; we may leave aside the orbital motion and denote by prn the radial momentum and by the radial energy. Let pr and be the highest radial momentum and, respectively, the highest radial energy; they correspond to the total momentum p and, respectively, total energy (in general, a degeneration may exist). This charge may tunnel through the potential barrier V(r) fromto The relevant factors in the wave function given by equation (4) are
(10)
Where it is easy to see that p2 = 0. It follows that the tunneling probability (transmission coefficient) is given by where
(11)
and (the condition ensures the existence of the bound state). We expand the coefficient A in powers of and take the average with respect to time; we get
(12)
the same procedure applied to the coefficient B leads to
(13)
where corresponds to the absence of the radiation; finally, we get
(14)
We can see that the effect of the radiation is to increase the rate of charge emission by a factor proportional to the square of the electric field and to introduce a slight anisotropy. It is worth noting that the radiation field contributes not only to the tunneling factor, as expressed by the coefficient B, but it is present also in the coefficient A, via the time-dependence of the wave function provided by the Kramers-Henneberger transform.
We can define a total disintegration probability
(15)
by integrating over angle θ, where The disintegration rate per unit time is where is related to the time ta estimated above and the time introduced by the energy uncertainty.41
The exponent corresponding to the absence of the radiation, is
(16)
since (for protons q2/a = 2.5M eV ) , we may use the approximate formulae
(17)
and (18)
As it is well know the interplay between the very large values of and the very small values of makes the disintegration rate to be very sensitive to the energy values, and to vary over a wide range.41 The result can be cast in the form of the Geiger-Nuttall law, which, in the absence of the radiation, can be written as and being well-known constants;41 the only effect of the radiation is to modify the constant into The correction to can also be written as for The maximum value of this correction is of the order of the unity; it follows that the decay rate is enhanced by the radiation by a factor of the order
After the emission of the charge, the mean-field potential suffers a reconfiguration (re-arrangement) process and the potential V (r) is modified; this is the well-known process of "core shake-up" (or "core excitation"); a new bound state is formed and a new transformation process may begin for the modified potential V(r). The tunneling probability w given above is a transmission coefficient (we can check that w<1); with probability 1-w the charge is reflected from the potential barrier; in these conditions the bound state is "shaken-up" and the charge resumes its motion, or its pre-formation process, until it tunnels, or is rescattered back to the core; the latter is the well-known recollision process.8,46‒50
The case of a static field requires a special discussion. Within the present formalism a static electric field E can be obtained from a vector potential A=-cEt; the position vector in the mean-field potential is shifted to where the special discussion is necessary because the parameter is unbounded in time. The distance a is covered in time for proton, a is of the order a = 10−13cm and the threshold field is E=E0=3×104 V/cm (102 electrostatic units) given above; we get This is a very long duration, in comparison with the relevant nuclear times, in particular the attempt time ( estimated above). In general, the condition of adiabatic interaction reads where is the mean separation of the energy levels; it implies which allows for high static fields. In these conditions the protons accommodate themselves to the electric field, which is absorbed into slightly modified energy levels; this change, which can be estimated by perturbation theory, is irrelevant for our discussion, since the field strength is small. However, it has an important consequence in that the electric field, once taken in the energy levels, is not available anymore for the Kramers-Henneberger transform given by equation (4); therefore, the present time-dependent formalism cannot be applied. Instead of using the hamiltonian given by equation (2), we start with the (equivalent) dipole hamiltonian which includes the interaction term −qEr. Consequently, the potential barrier is changed into
(19)
We compute the tunneling rate by using this potential barrier. In view of the small value of the correction parameter proportional to E in equation (19), we may expand the momentum in powers of this parameter and replace the powers of r2 by their mean values over the tunneling range from r1=a to since r2≫a, we get the small parameter in equation (19). For Z=100 and this parameter is which is much smaller than unity for any usual static field. Integrating over angles and assuming where is given by equation (17), we get finally
(20)
We can see that the correction brought by a static electric field to the decay rate is extremely small, as expected.
Finally, it is worth discussing the case of intermediate fields, i.e. field strengths which satisfy the inequality (in our case, fields from 3×104 V/cm to 107 V/cm).51 In this case the adiabatic hypothesis cannot be used anymore, and the initial conditions for introducing the interaction are important. The corresponding Kramers-Henneberger transform diminishes appreciably the potential barrier and the charge is set free in a short time, which is the reciprocal of the decay rate; this rate may exhibit oscillations as a function of the field strength.