-
Minimal Basis Iterative Stockholder: Atoms in Molecules for Force-Field Development
Authors:
Toon Verstraelen,
Steven Vandenbrande,
Farnaz Heidar-Zadeh,
Louis Vanduyfhuys,
Veronique Van Speybroeck,
Michel Waroquier,
Paul W. Ayers
Abstract:
Atomic partial charges appear in the Coulomb term of many force-field models and can be derived from electronic structure calculations with a myriad of atoms-in-molecules (AIM) methods. More advanced models have also been proposed, using the distributed nature of the electron cloud and atomic multipoles. In this work, an electrostatic force field is defined through a concise approximation of the e…
▽ More
Atomic partial charges appear in the Coulomb term of many force-field models and can be derived from electronic structure calculations with a myriad of atoms-in-molecules (AIM) methods. More advanced models have also been proposed, using the distributed nature of the electron cloud and atomic multipoles. In this work, an electrostatic force field is defined through a concise approximation of the electron density, for which the Coulomb interaction is trivially evaluated. This approximate "pro-density" is expanded in a minimal basis of atom-centered s-type Slater density functions, whose parameters are optimized by minimizing the Kullback-Leibler divergence of the pro-density from a reference electron density, e.g. obtained from an electronic structure calculation. The proposed method, Minimal Basis Iterative Stockholder (MBIS), is a variant of the Hirshfeld AIM method but it can also be used as a density-fitting technique. An iterative algorithm to refine the pro-density is easily implemented with a linear-scaling computational cost, enabling applications to supramolecular systems. The benefits of the MBIS method are demonstrated with systematic applications to molecular databases and extended models of condensed phases. A comparison to 14 other AIM methods shows its effectiveness when modeling electrostatic interactions. MBIS is also suitable for rescaling atomic polarizabilities in the Tkatchenko-Sheffler scheme for dispersion interactions.
△ Less
Submitted 19 August, 2016;
originally announced August 2016.
-
Fast density-matrix based partitioning of the energy over the atoms in a molecule consistent with the Hirshfeld-I partitioning of the electron density
Authors:
Diederik Vanfleteren,
Dieter Ghillemijn,
Dimitri Van Neck,
Patrick Bultinck,
Michel Waroquier,
Paul W. Ayers
Abstract:
For the Hirshfeld-I atom-in-molecule model, associated single-atom energies and interaction energies at the Hartree-Fock level are determined efficiently in one-electron Hilbert space. In contrast to most other approaches, the energy terms are fully consistent with the partitioning of the underlying one-electron density matrix. Starting from the Hirshfeld-I atom-in-molecule model for the electron…
▽ More
For the Hirshfeld-I atom-in-molecule model, associated single-atom energies and interaction energies at the Hartree-Fock level are determined efficiently in one-electron Hilbert space. In contrast to most other approaches, the energy terms are fully consistent with the partitioning of the underlying one-electron density matrix. Starting from the Hirshfeld-I atom-in-molecule model for the electron density, the molecular one-electron density matrix is partitioned with a previously introduced double-atom scheme [Vanfleteren D. et al., J Chem Phys 2010, 132, 164111]. Single-atom density matrices are constructed from the atomic and bond contributions of the double-atom scheme. Since the Hartree-Fock energy can be expressed solely in terms of the one-electron density matrix, the partitioning of the latter over the atoms in the molecule leads naturally to a corresponding partitioning of the Hartree-Fock energy. When the size of the molecule or the molecular basis set does not grow too large, the method shows considerable computational advantages compared to other approaches that require cumbersome numerical integration of the molecular energy integrals weighted by atomic weight functions.
△ Less
Submitted 3 November, 2011;
originally announced November 2011.
-
Hilbert-space partitioning of the molecular one-electron density matrix with orthogonal projectors
Authors:
Diederik Vanfleteren,
Dimitri Van Neck,
Patrick Bultinck,
Paul W. Ayers,
Michel Waroquier
Abstract:
A double-atom partitioning of the molecular one-electron density matrix is used to describe atoms and bonds. All calculations are performed in Hilbert space. The concept of atomic weight functions (familiar from Hirshfeld analysis of the electron density) is extended to atomic weight matrices. These are constructed to be orthogonal projection operators on atomic subspaces, which has significant ad…
▽ More
A double-atom partitioning of the molecular one-electron density matrix is used to describe atoms and bonds. All calculations are performed in Hilbert space. The concept of atomic weight functions (familiar from Hirshfeld analysis of the electron density) is extended to atomic weight matrices. These are constructed to be orthogonal projection operators on atomic subspaces, which has significant advantages in the interpretation of the bond contributions. In close analogy to the iterative Hirshfeld procedure, self-consistency is built in at the level of atomic charges and occupancies. The method is applied to a test set of about 67 molecules, representing various types of chemical binding. A close correlation is observed between the atomic charges and the Hirshfeld-I atomic charges.
△ Less
Submitted 3 November, 2011;
originally announced November 2011.
-
Partitioning of the molecular density matrix over atoms and bonds
Authors:
Diederik Vanfleteren,
Dimitri Van Neck,
Patrick Bultinck,
Paul W. Ayers,
Michel Waroquier
Abstract:
A double-index atomic partitioning of the molecular first-order density matrix is proposed. Contributions diagonal in the atomic indices correspond to atomic density matrices, whereas off-diagonal contributions carry information about the bonds. The resulting matrices have good localization properties, in contrast to single-index atomic partitioning schemes of the molecular density matrix. It is s…
▽ More
A double-index atomic partitioning of the molecular first-order density matrix is proposed. Contributions diagonal in the atomic indices correspond to atomic density matrices, whereas off-diagonal contributions carry information about the bonds. The resulting matrices have good localization properties, in contrast to single-index atomic partitioning schemes of the molecular density matrix. It is shown that the electron density assigned to individual atoms, when derived from the density matrix partitioning, can be made con- sistent with well-known partitions of the electron density over AIM basins, either with sharp or with fuzzy boundaries. The method is applied to a test set of about 50 molecules, representative for various types of chemical binding. A close correlation is observed between the trace of the bond matrices and the SEDI (shared electron density index) bond index.
△ Less
Submitted 3 November, 2011;
originally announced November 2011.
-
Quasiparticle properties in a density functional framework
Authors:
D. Van Neck,
S. Verdonck,
G. Bonny,
P. W. Ayers,
M. Waroquier
Abstract:
We propose a framework to construct the ground-state energy and density matrix of an N-electron system by solving selfconsistently a set of single-particle equations. The method can be viewed as a non-trivial extension of the Kohn-Sham scheme (which is embedded as a special case). It is based on separating the Green's function into a quasi-particle part and a background part, and expressing only…
▽ More
We propose a framework to construct the ground-state energy and density matrix of an N-electron system by solving selfconsistently a set of single-particle equations. The method can be viewed as a non-trivial extension of the Kohn-Sham scheme (which is embedded as a special case). It is based on separating the Green's function into a quasi-particle part and a background part, and expressing only the background part as a functional of the density matrix. The calculated single-particle energies and wave functions have a clear physical interpretation as quasiparticle energies and orbitals.
△ Less
Submitted 5 September, 2006;
originally announced September 2006.
-
On the nuclear symmetry energy and the neutron skin in neutron-rich nuclei
Authors:
A. E. L. Dieperink,
Y. Dewulf,
D. Van Neck,
M. Waroquier,
V. Rodin
Abstract:
The symmetry energy for nuclear matter and its relation to the neutron skin in finite nuclei is discussed. The symmetry energy as a function of density obtained in a self-consistent Green function approach is presented and compared to the results of other recent theoretical approaches. A partial explanation of the linear relation between the symmetry energy and the neutron skin is proposed. The…
▽ More
The symmetry energy for nuclear matter and its relation to the neutron skin in finite nuclei is discussed. The symmetry energy as a function of density obtained in a self-consistent Green function approach is presented and compared to the results of other recent theoretical approaches. A partial explanation of the linear relation between the symmetry energy and the neutron skin is proposed. The potential of several experimental methods to extract the neutron skin is examined.
△ Less
Submitted 10 December, 2003;
originally announced December 2003.
-
Saturation of nuclear matter and short-range correlations
Authors:
Y. Dewulf,
W. H. Dickhoff,
D. Van Neck,
E. R. Stoddard,
M. Waroquier
Abstract:
A fully self-consistent treatment of short-range correlations in nuclear matter is presented. Different implementations of the determination of the nucleon spectral functions for different interactions are shown to be consistent with each other. The resulting saturation densities are closer to the empirical result when compared with (continuous-choice) Brueckner-Hartree-Fock values. Arguments fo…
▽ More
A fully self-consistent treatment of short-range correlations in nuclear matter is presented. Different implementations of the determination of the nucleon spectral functions for different interactions are shown to be consistent with each other. The resulting saturation densities are closer to the empirical result when compared with (continuous-choice) Brueckner-Hartree-Fock values. Arguments for the dominance of short-range correlations in determining the nuclear-matter saturation density are presented. A further survey of the role of long-range correlations suggests that the inclusion of pionic contributions to ring diagrams in nuclear matter leads to higher saturation densities than empirically observed. A possible resolution of the nuclear-matter saturation problem is suggested.
△ Less
Submitted 19 March, 2003;
originally announced March 2003.
-
Electromagnetic interaction in chiral quantum hadrodynamics and decay of vector and axial-vector mesons
Authors:
A. Yu. Korchin,
D. Van Neck,
M. Waroquier
Abstract:
The chiral invariant QHD-III model of Serot and Walecka is applied in the calculation of some meson properties. The electromagnetic interaction is included by extending the symmetry of the model to the local U(1) \times SU(2)_{R} \times SU(2)_{L} group. The minimal and nonminimal contributions to the electromagnetic Lagrangian are obtained in a new representation of QHD-III. Strong decays of the…
▽ More
The chiral invariant QHD-III model of Serot and Walecka is applied in the calculation of some meson properties. The electromagnetic interaction is included by extending the symmetry of the model to the local U(1) \times SU(2)_{R} \times SU(2)_{L} group. The minimal and nonminimal contributions to the electromagnetic Lagrangian are obtained in a new representation of QHD-III. Strong decays of the axial-vector meson, a_{1} \to πρ, a_{1} \to πσ, and the electromagnetic decays ρ\to ππγ, a_{1} \to πγand ρ\to πγare calculated. The low-energy parameters for the π-πscattering are calculated in the tree-level approximation. The effect of the auxiliary Higgs bosons, introduced in QHD-III in order to generate masses of the vector and axial-vector mesons via the Higgs mechanism, is studied as well. This is done on the tree level for π-πscattering and on the level of one-loop diagrams for the a_{1} \to πγdecay. It is demonstrated that the model successfully describes some features of meson phenomenology in the non-strange sector.
△ Less
Submitted 17 February, 2003;
originally announced February 2003.
-
Coherent Compton scattering on light nuclei in the delta resonance region
Authors:
L. Van Daele,
A. Yu. Korchin,
D. Van Neck,
O. Scholten,
M. Waroquier
Abstract:
Coherent Compton scattering on light nuclei in the delta resonance region is studied in the impulse approximation and is shown to be a sensitive probe of the in-medium properties of the delta resonance. The elementary amplitude on a single nucleon is calculated from the unitary K-matrix approach developed previously. Modifications of the properties of the delta resonance due to the nuclear mediu…
▽ More
Coherent Compton scattering on light nuclei in the delta resonance region is studied in the impulse approximation and is shown to be a sensitive probe of the in-medium properties of the delta resonance. The elementary amplitude on a single nucleon is calculated from the unitary K-matrix approach developed previously. Modifications of the properties of the delta resonance due to the nuclear medium are accounted for through the self-energy operator of the delta, calculated from the one-pion loop. The dominant medium effects such as the Pauli blocking, mean-field modification of the nucleon and delta masses, and particle-hole excitations in the pion propagator are consistently included in nuclear matter.
△ Less
Submitted 23 October, 2001;
originally announced October 2001.
-
Improved lower bounds for the ground-state energy of many-body systems
Authors:
D. Van Neck,
Y. Dewulf,
M. Waroquier
Abstract:
New lower bounds for the binding energy of a quantum-mechanical system of interacting particles are presented. The new bounds are expressed in terms of two-particle quantities and improve the conventional bounds of the Hall-Post type. They are constructed by considering not only the energy in the two-particle system, but also the structure of the pair wave function. We apply the formal results t…
▽ More
New lower bounds for the binding energy of a quantum-mechanical system of interacting particles are presented. The new bounds are expressed in terms of two-particle quantities and improve the conventional bounds of the Hall-Post type. They are constructed by considering not only the energy in the two-particle system, but also the structure of the pair wave function. We apply the formal results to various numerical examples, and show that in some cases dramatic improvement over the existing bounds is reached.
△ Less
Submitted 28 February, 2001;
originally announced February 2001.
-
Short-range correlations in nuclear matter using Green's functions within a discrete pole approximation
Authors:
Y. Dewulf,
D. Van Neck,
M. Waroquier
Abstract:
We treat short-range correlations in nuclear matter, induced by the repulsive core of the nucleon-nucleon potential, within the framework of a self-consistent Green's function theory. The effective in-medium interaction sums the ladder diagrams of both the particle-particle and hole-hole type. The demand of self-consistency results in a set of nonlinear equations which must be solved by iteratio…
▽ More
We treat short-range correlations in nuclear matter, induced by the repulsive core of the nucleon-nucleon potential, within the framework of a self-consistent Green's function theory. The effective in-medium interaction sums the ladder diagrams of both the particle-particle and hole-hole type. The demand of self-consistency results in a set of nonlinear equations which must be solved by iteration. We explore the possibility of approximating the single-particle Green's function by a limited number of poles and residues.
△ Less
Submitted 24 April, 2001; v1 submitted 6 December, 2000;
originally announced December 2000.
-
Production of e+e- pairs in proton-deuteron capture to 3He
Authors:
A. Yu. Korchin,
D. Van Neck,
M. Waroquier,
O. Scholten,
A. E. L. Dieperink
Abstract:
The process p+d \leftrightarrow 3He + γ* at intermediate energies is described using a covariant and gauge-invariant model, and a realistic pd3He vertex. Both photodisintegration of 3He and proton-deuteron capture with production of e+e- pairs are studied, and results for cross sections and response functions are presented. The effect of time-like formfactors on the dilepton cross sections is in…
▽ More
The process p+d \leftrightarrow 3He + γ* at intermediate energies is described using a covariant and gauge-invariant model, and a realistic pd3He vertex. Both photodisintegration of 3He and proton-deuteron capture with production of e+e- pairs are studied, and results for cross sections and response functions are presented. The effect of time-like formfactors on the dilepton cross sections is investigated as well.
△ Less
Submitted 1 October, 1998;
originally announced October 1998.
-
Center-of-mass effects on the quasi-hole spectroscopic factors in the 16O(e,e'p) reaction
Authors:
D. Van Neck,
M. Waroquier,
A. E. L. Dieperink,
Steven C. Pieper,
V. R. Pandharipande
Abstract:
The spectroscopic factors for the low-lying quasi-hole states observed in the 16O(e,e'p)15N reaction are reinvestigated with a variational Monte Carlo calculation for the structure of the initial and final nucleus. A computational error in a previous report is rectified. It is shown that a proper treatment of center-of-mass motion does not lead to a reduction of the spectroscopic factor for $p$-…
▽ More
The spectroscopic factors for the low-lying quasi-hole states observed in the 16O(e,e'p)15N reaction are reinvestigated with a variational Monte Carlo calculation for the structure of the initial and final nucleus. A computational error in a previous report is rectified. It is shown that a proper treatment of center-of-mass motion does not lead to a reduction of the spectroscopic factor for $p$-shell quasi-hole states, but rather to a 7% enhancement. This is in agreement with analytical results obtained in the harmonic oscillator model. The center-of-mass effect worsens the discrepancy between present theoretical models and the experimentally observed single-particle strength. We discuss the present status of this problem, including some other mechanisms that may be relevant in this respect.
△ Less
Submitted 28 April, 1998;
originally announced April 1998.
-
Electroinduced two-nucleon knockout and correlations in nuclei
Authors:
J. Ryckebusch,
V. Van der Sluys,
K. Heyde,
H. Holvoet,
W. Van Nespen,
M. Waroquier,
M. Vanderhaeghen
Abstract:
We present a model to calculate cross sections for electroinduced two-nucleon emission from finite nuclei. Short-range correlations in the wave functions and meson-exchange contributions to the photoabsorption process are implemented. Effects of the short-range correlations are studied with the aid of a perturbation expansion method with various choices of the Jastrow correlation function. The m…
▽ More
We present a model to calculate cross sections for electroinduced two-nucleon emission from finite nuclei. Short-range correlations in the wave functions and meson-exchange contributions to the photoabsorption process are implemented. Effects of the short-range correlations are studied with the aid of a perturbation expansion method with various choices of the Jastrow correlation function. The model is used to investigate the relative importance of the different reaction mechanisms contributing to the A(e,e$'$pn) and A(e,e$'$pp) process. Representative examples for the target nuclei $^{12}$C and $^{16}$O and for kinematical conditions accessible with contemporary high-duty cycle electron accelerators are presented. A procedure is outlined to calculate the two-nucleon knockout contribution to the semi-exclusive (e,e$'$p) cross section. Using this technique we investigate in how far semi-exclusive (e,e$'$p) reactions can be used to detect high-momentum components in the nuclear spectral function.
△ Less
Submitted 24 February, 1997;
originally announced February 1997.
-
Systematic study of Coulomb distortion effects in exclusive (e,e'p) reactions
Authors:
V. Van der Sluys,
K. Heyde,
J. Ryckebusch,
M. Waroquier
Abstract:
A technique to deal with Coulomb electron distortions in the analysis of (e,e'p) reactions is presented. Thereby, no approximations are made. The suggested technique relies on a partial-wave expansion of the electron wave functions and a multipole decomposition of the electron and nuclear current in momentum space. In that way, we succeed in keeping the computational times within reasonable limi…
▽ More
A technique to deal with Coulomb electron distortions in the analysis of (e,e'p) reactions is presented. Thereby, no approximations are made. The suggested technique relies on a partial-wave expansion of the electron wave functions and a multipole decomposition of the electron and nuclear current in momentum space. In that way, we succeed in keeping the computational times within reasonable limits. This theoretical framework is used to calculate the quasielastic (e,e'p) reduced cross sections for proton knockout from the valence shells in $^{16}$O, $^{40}$Ca, $^{90}$Zr and $^{208}$Pb. The final-state interaction of the ejected proton with the residual nucleus is treated within an optical potential model. The role of electron distortion on the extracted spectroscopic factors is discussed.
△ Less
Submitted 31 January, 1997;
originally announced February 1997.
-
Long-range correlations in finite nuclei: comparison of two self-consistent treatments
Authors:
Y. Dewulf,
D. Van Neck,
L. Van Daele,
M. Waroquier
Abstract:
Long-range correlations, which are partially responsible for the observed fragmentation and depletion of low-lying single-particle strength, are studied in the Green's function formalism. The self-energy is expanded up to second order in the residual interaction. We compare two methods of implementing self-consistency in the solution of the Dyson equation beyond Hartree-Fock, for the case of the…
▽ More
Long-range correlations, which are partially responsible for the observed fragmentation and depletion of low-lying single-particle strength, are studied in the Green's function formalism. The self-energy is expanded up to second order in the residual interaction. We compare two methods of implementing self-consistency in the solution of the Dyson equation beyond Hartree-Fock, for the case of the 16O nucleus. It is found that the energy-bin method and the BAGEL method lead to globally equivalent results. In both methods the final single-particle strength functions are characterized by exponential tails at energies far from the Fermi level.
△ Less
Submitted 29 January, 1997;
originally announced January 1997.
-
The exclusive (e,e$'$p) reaction at high missing momenta
Authors:
V. Van der Sluys,
J. Ryckebusch,
M. Waroquier
Abstract:
The reduced (e,e$'$p) cross section is calculated for kinematics that probe high missing momenta. The final-state interaction is handled within a non-relativistic many-body framework. One- and two-body nuclear currents are included. Electron distortion effects are treated in an exact distorted wave calculation. It is shown that at high missing momenta the calculated (e,e$'$p) cross sections exhi…
▽ More
The reduced (e,e$'$p) cross section is calculated for kinematics that probe high missing momenta. The final-state interaction is handled within a non-relativistic many-body framework. One- and two-body nuclear currents are included. Electron distortion effects are treated in an exact distorted wave calculation. It is shown that at high missing momenta the calculated (e,e$'$p) cross sections exhibit a pronounced sensitivity to ground-state correlations of the RPA type and two-body currents. The role of these mechanisms is found to be relatively small at low missing momenta.
△ Less
Submitted 12 June, 1996;
originally announced June 1996.
-
Two-body currents in inclusive electron scattering
Authors:
Veerle Van der Sluys,
Jan Ryckebusch,
Michel Waroquier
Abstract:
The longitudinal and transverse structure functions are calculated for inclusive electron scattering from ${}^{12}$C and ${}^{40}$Ca in the quasielastic and dip region. The microscopic model presented here incorporates two-body currents derived from one-pion exchange and intermediate $Δ$ excitation. The reaction mechanism involves both one-nucleon and two-nucleon knockout processes. It is demonstr…
▽ More
The longitudinal and transverse structure functions are calculated for inclusive electron scattering from ${}^{12}$C and ${}^{40}$Ca in the quasielastic and dip region. The microscopic model presented here incorporates two-body currents derived from one-pion exchange and intermediate $Δ$ excitation. The reaction mechanism involves both one-nucleon and two-nucleon knockout processes. It is demonstrated that, even for quasielastic kinematics, two-body currents give a substantial contribution to the transverse structure function. Furthermore, the observed excess of transverse strength in the dip region between the quasielastic and delta peak can be partially ascribed to direct two-nucleon knockout following photoabsorption on a two-body current.
△ Less
Submitted 9 March, 1995;
originally announced March 1995.
-
Short-range correlations in (e,e$'$p) and (e,e$'$pp) reactions on complex nuclei
Authors:
Jan Ryckebusch,
Marc Vanderhaeghen,
Kris Heyde,
Michel Waroquier
Abstract:
The influence of short-range correlations (SRC) on the triple-coincidence (e,e$'$pp) reactions is studied. The non-relativistic model uses a mean-field potential to account for the distortions that the escaping particles undergo. Apart from the SRC, that are implemented through a Jastrow ansatz with a realistic correlation function, we incorporate the contribution from pion exchange and intermedia…
▽ More
The influence of short-range correlations (SRC) on the triple-coincidence (e,e$'$pp) reactions is studied. The non-relativistic model uses a mean-field potential to account for the distortions that the escaping particles undergo. Apart from the SRC, that are implemented through a Jastrow ansatz with a realistic correlation function, we incorporate the contribution from pion exchange and intermediate $Δ_{33}$ currents. The (e,e$'$pp) cross sections are predicted to exhibit a sizeable sensitivity to the SRC. The contribution from the two-nucleon breakup channel to the semi-exclusive $^{12}$C(e,e$'$p) cross section is calculated in the kinematics of a recent NIKHEF-K experiment. In the semi-exclusive channel, a selective sensitivity in terms of the missing energy and momentum to the SRC is found.
△ Less
Submitted 9 March, 1995;
originally announced March 1995.
-
Two-nucleon knockout contributions to the $^{12}$C$(e,e'p)$ reaction in the dip and {$Δ$}(1232) regions
Authors:
J. Ryckebusch,
V. Van der Sluys,
M. Waroquier,
L. J. H. M. Kester,
W. H. A. Hesselink,
E. Jans,
A. Zondervan
Abstract:
The contributions from $^{12}$C$(e,e'pn)$ and $^{12}$C$(e,e'pp)$ to the semi-exclusive $^{12}$C$(e,e'p)$ cross section have been calculated in an unfactorized model for two-nucleon emission. We assume direct two-nucleon knockout after virtual photon coupling with the two-body pion-exchange currents in the target nucleus. Results are presented at several kinematical conditions in the dip and $Δ$(…
▽ More
The contributions from $^{12}$C$(e,e'pn)$ and $^{12}$C$(e,e'pp)$ to the semi-exclusive $^{12}$C$(e,e'p)$ cross section have been calculated in an unfactorized model for two-nucleon emission. We assume direct two-nucleon knockout after virtual photon coupling with the two-body pion-exchange currents in the target nucleus. Results are presented at several kinematical conditions in the dip and $Δ$(1232) regions. The calculated two-nucleon knockout strength is observed to account for a large fraction of the measured $(e,e'p)$ strength above the two-nucleon emission threshold.
△ Less
Submitted 13 June, 1994;
originally announced June 1994.
-
Multinucleon Mechanisms in ($γ$,N) and ($γ$,NN) Reactions
Authors:
J. Ryckebusch,
L. Machenil,
M. Vanderhaeghen,
V. Van der Sluys,
M. Waroquier
Abstract:
The similarities in the experimental indications for multinucleon mechanisms in $(γ,p)$ and $(e,e'p)$ processes are pointed out. For both types of reactions, the substantial role of two-nucleon emission processes for transitions to high excitation energies in the residual nucleus is stressed. A microscopic model for the calculation of the two-body knockout contributions to the inclusive $(γ,N)$…
▽ More
The similarities in the experimental indications for multinucleon mechanisms in $(γ,p)$ and $(e,e'p)$ processes are pointed out. For both types of reactions, the substantial role of two-nucleon emission processes for transitions to high excitation energies in the residual nucleus is stressed. A microscopic model for the calculation of the two-body knockout contributions to the inclusive $(γ,N)$ reaction is presented. It is based on an unfactorized formalism for the calculation of electromagnetically induced two-nucleon emission cross sections. The model is shown to yield a reasonable description of the overall behaviour of the $^{12}$C$(γ,p)$ and $^{12}$C$(γ,n)$ data at high excitation energies in the residual nucleus. In the calculations effects from non-resonant and resonant pion exchange currents are included. Photoabsorption on these currents are predicted to produce the major contributions to the exclusive $^{16}$O$(γ,n_0)^{15}$O process at photonenergies above the pion threshold. Double differential cross sections for photon induced $pp$ and $pn$ emission from $^{16}$O are calculated and compared with the data.
△ Less
Submitted 21 October, 1993;
originally announced October 1993.
-
Effects of meson-exchange currents on the $(\protect\vec{e},e'p)$ structure functions
Authors:
V. Van der Sluys,
J. Ryckebusch,
M. Waroquier
Abstract:
The response functions for the unpolarized $(e,e'p)$ and polarized $(\vec{e},e'p)$ reaction are calculated for medium-heavy nuclei under quasifree conditions. The formalism presented here incorporates two-body currents related to meson-exchange and the $Δ(1232)$ excitation. The final-state interaction of the outgoing nucleon with the residual nucleus is handled in a HF-RPA formalism. The sensiti…
▽ More
The response functions for the unpolarized $(e,e'p)$ and polarized $(\vec{e},e'p)$ reaction are calculated for medium-heavy nuclei under quasifree conditions. The formalism presented here incorporates two-body currents related to meson-exchange and the $Δ(1232)$ excitation. The final-state interaction of the outgoing nucleon with the residual nucleus is handled in a HF-RPA formalism. The sensitivity of the results to the two-body currents is discussed for the five structure functions in quasielastic $(\vec{e},e'p)$ scattering off the target nuclei $^{16}O$ and $^{40}Ca$. A selective sensitivity to the two-body currents is obtained in the longitudinal-transverse interference term $W_{LT}$ where two-body currents can explain part of the discrepancy between the impulse-approximation calculations and the data.
△ Less
Submitted 13 July, 1993;
originally announced July 1993.
-
Effects of the final-state interaction in ($γ$,pn) and ($γ$,pp) processes
Authors:
J. Ryckebusch,
M. Vanderhaeghen,
L. Machenil,
M. Waroquier
Abstract:
A model is presented to describe electromagnetically induced two-nucleon emission processes in a shell-model picture. Distortions in the outgoing nucleon waves are accounted for by performing a partial-wave expansion in a real mean-field potential. The antisymmetry condition for the A-body wavefunctions is shown to be naturally preserved.
The model is used to calculate ($γ$,pn) and ($γ$,pp) cr…
▽ More
A model is presented to describe electromagnetically induced two-nucleon emission processes in a shell-model picture. Distortions in the outgoing nucleon waves are accounted for by performing a partial-wave expansion in a real mean-field potential. The antisymmetry condition for the A-body wavefunctions is shown to be naturally preserved.
The model is used to calculate ($γ$,pn) and ($γ$,pp) cross sections off the target nuclei $^{16}$O and $^{12}$C for photon energies ranging from 50 MeV up to the $\bigtriangleup$(1232) isobar threshold. Effects due to the pionic currents and intermediate $\bigtriangleup$ creation are implemented. The impact of the distortions due to the interaction of the outgoing nucleon waves with the (A-2) core is examined. Hadronic form factors are introduced to regularize the $π$NN vertices and the sensitivity of the cross section to the pion cut-off mass is examined. The relative contribution of the ($γ$,pp) and the ($γ$,pn) channel to the total photoabsorption strength is discussed. Further, the photon energy dependence of the ($γ$,pp)/($γ$,pn) ratio is investigated.
△ Less
Submitted 5 July, 1993;
originally announced July 1993.
-
Absorption mechanisms in photon induced two-body knockout
Authors:
L. Machenil,
M. Vanderhaeghen,
J. Ryckebusch,
M. Waroquier
Abstract:
Calculations have been performed for the $^{16}$O($γ$,pn) and the $^{16}$O($γ$,pp) reaction in the photon-energy range $E_γ$ = 60-300 MeV. Besides the contribution from the more common photoabsorption on the pionic degrees of freedom, we have investigated the influence of heavier meson exchange ($ρ, σ, ω$) and intermediate $Δ$ creation with $π$ and $ρ$ exchange. Whereas the $π$ meson is found to…
▽ More
Calculations have been performed for the $^{16}$O($γ$,pn) and the $^{16}$O($γ$,pp) reaction in the photon-energy range $E_γ$ = 60-300 MeV. Besides the contribution from the more common photoabsorption on the pionic degrees of freedom, we have investigated the influence of heavier meson exchange ($ρ, σ, ω$) and intermediate $Δ$ creation with $π$ and $ρ$ exchange. Whereas the $π$ meson is found to set the main trends, the $ρ$ meson is found not to be discardable in a theoretical description of the ($γ$,pn) reaction. The incorporation of an energy dependence and a decay width in the $Δ$ propagator is observed to be essential in order to arrive at a more realistic description of ($γ$,NN) reactions at higher photon energies.
△ Less
Submitted 2 July, 1993; v1 submitted 2 July, 1993;
originally announced July 1993.