Source author record

P. D. Stevenson

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

24 published item(s)

preprint2019arXiv

Role of the Surface Energy in Heavy-Ion Collisions

The surface energy is one of the fundamental properties nuclei, appearing in the simplest form of the semi-empirical mass formula. The surface enery has an influence on e.g. the shape of a nucleus and its ability to deform. This in turn could be expected to have an effect in fusion reactions around the Coulomb barrier where dynamical effects such as the formation of a neck is part of the fusion process. Frozen Hartree-Fock and Time-Dependent Hartree-Fock calculations are made for a series of effective interactions in which the surface energy is systematically varied, using $^{40}$Ca + $^{48}$Ca as a test case. The dynamical lowering of the barrier is greatest for the largest surface energy, contrary to naive expectations, and we speculate that this may be due to the variation in other nuclear matter properties for these effective interactions

preprint2016arXiv

Fission dynamics within time-dependent Hartree-Fock: boost-induced fission

Background: Nuclear fission is a complex large-amplitude collective decay mode in heavy nuclei. Microscopic density functional studies of fission have previously concentrated on adiabatic approaches based on constrained static calculations ignoring dynamical excitations of the fissioning nucleus, and the daughter products. Purpose: To explore the ability of dynamic mean-field methods to describe induced fission processes, using quadrupole boosts in the nuclide $^{240}$Pu as an example. Methods: Quadrupole constrained Hartree-Fock calculations are used to create a potential energy surface. An isomeric state and a state beyond the second barrier peak are excited by means of instantaneous as well as temporally extended gauge boosts with quadrupole shapes. The subsequent deexcitation is studied in a time-dependent Hartree-Fock simulation, with emphasis on fissioned final states. The corresponding fission fragment mass numbers are studied. Results: In general, the energy deposited by the quadrupole boost is quickly absorbed by the nucleus. In instantaneous boosts, this leads to fast shape rearrangements and violent dynamics that can ultimately lead to fission. This is a qualitatively different process than the deformation-induced fission. Boosts induced within a finite time window excite the system in a relatively gentler way, and do induce fission but with a smaller energy deposition. Conclusions: The fission products obtained using boost-induced fission in time-dependent Hartree-Fock are more asymmetric than the fragments obtained in deformation-induced fission, or the corresponding adiabatic approaches.

preprint2016arXiv

Sensitivity of the fusion cross section to the density dependence of the symmetry energy

It is the aim of this paper to discuss the impact of nuclear fusion on the EOS. This is a timely subject given the expected availability of increasingly exotic beams at rare isotope facilities\,\cite{balantekin2014}. In practice, we focus on $^{48}$Ca+$^{48}$Ca fusion. We employ three different approaches to calculate fusion cross-sections for a set of energy density functionals with systematically varying nuclear matter properties. Fusion calculations are performed using frozen densities, using a dynamic microscopic method based on density-constrained time-dependent Hartree-Fock (DC-TDHF) approach, as well as direct TDHF study of above barrier cross-sections. For these studies, we employ a family of Skyrme parametrizations with systematically varied nuclear matter properties. We find a slight preference for forces which deliver a slope of symmetry energy of $L\approx 50$\,MeV that corresponds to a neutron-skin thickness of $^{48}$Ca of $R_\mathrm{skin}\!=\!(0.180\!-\!0.210)$\,fm.

preprint2015arXiv

Fission dynamics within time-dependent Hartree-Fock: deformation-induced fission

Background: Nuclear fission is a complex large-amplitude collective decay mode in heavy nuclei. Microscopic density functional studies of fission have previously concentrated on adiabatic approaches based on constrained static calculations ignoring dynamical excitations of the fissioning nucleus, and the daughter products. Purpose: To explore the ability of dynamic mean-field methods to describe fast fission processes beyond the fission barrier, using the nuclide $^{240}$Pu as an example. Methods: Time-dependent Hartree-Fock calculations based on the Skyrme interaction are used to calculate non-adiabatic fission paths, beginning from static constrained Hartree-Fock calculations. The properties of the dynamic states are interpreted in terms of the nature of their collective motion. Fission product properties are compared to data. Results: Parent nuclei constrained to begin dynamic evolution with a deformation less than the fission barrier exhibit giant-resonance-type behaviour. Those beginning just beyond the barrier explore large amplitude motion but do not fission, whereas those beginning beyond the two-fragment pathway crossing fission to final states which differ according to the exact initial deformation. Conclusions: Time-dependent Hartree-Fock is able to give a good qualitative and quantitative description of fast fission, provided one begins from a sufficiently deformed state.

preprint2015arXiv

Resonances and Reactions from Mean-Field Dynamics

The time-dependent version of nuclear density functional theory, using functionals derived from Skyrme interactions, is able to approximately describe nuclear dynamics. We present time-dependent results of calculations of dipole resonances, concentrating on excitations of valence neutrons against a proton plus neutron core in the neutron-rich doubly-magic $^{132}$Sn nucleus, and results of collision dynamics, highlighting potential routes to ternary fusion, with the example of a collision of $^{48}$Ca+$^{48}$Ca+$^{208}$Pb resulting in a compound nucleus of element 120 stable against immiedate fission

preprint2015arXiv

Shapes and Dynamics from the Time-Dependent Mean Field

Explaining observed properties in terms of underlying shape degrees of freedom is a well--established prism with which to understand atomic nuclei. Self--consistent mean--field models provide one tool to understand nuclear shapes, and their link to other nuclear properties and observables. We present examples of how the time--dependent extension of the mean--field approach can be used in particular to shed light on nuclear shape properties, particularly looking at the giant resonances built on deformed nuclear ground states, and at dynamics in highly-deformed fission isomers. Example calculations are shown of $^{28}$Si in the first case, and $^{240}$Pu in the latter case.

preprint2015arXiv

The Skyrme Tensor Force in Heavy Ion Collisions

Background: [elided for arXiv character count] Purpose: With the advancement of computational power it is now possible to undertake TDHF calculations without any symmetry assumptions and incorporate the major strides made by the nuclear structure community in improving the energy density functionals used in these calculations. In particular, time-odd and tensor terms in these functionals are naturally present during the dynamical evolution, while being absent or minimally important for most static calculations. The parameters of these terms are determined by the requirement of Galilean invariance or local gauge invariance but their significance for the reaction dynamics have not been fully studied. This work addresses this question with emphasis on the tensor force. Method: The full version of the Skyrme force, including terms arising only from the Skyrme tensor force, is applied to the study of collisions within a completely symmetry-unrestricted TDHF implementation. Results: We examine the effect on fusion thresholds with and without the tensor force terms and find an effect on the fusion threshold energy of the order several MeV. Details of the distribution of the energy within terms in the energy density functional is also discussed. Conclusions: Terms in the energy density functional linked to the tensor force can play a non-negligible role in dynamic processes in nuclei, including those in which the terms do not affect ground state properties.

preprint2015arXiv

The unrestricted Skyrme-tensor time-dependent Hartree-Fock and its application to the nuclear response from spherical to triaxial nuclei

The nuclear time-dependent Hartree-Fock model formulated in the three-dimensional space,based on the full Skyrme energy density functional and complemented with the tensor force,is presented for the first time. Full self-consistency is achieved by the model. The application to the isovector giant dipole resonance is discussed in the linear limit, ranging from spherical nuclei (16O, 120Sn) to systems displaying axial or triaxial deformation (24Mg, 28Si, 178Os, 190W, 238U). Particular attention is paid to the spin-dependent terms from the central sector of the functional, recently included together with the tensor. They turn out to be capable of producing a qualitative change on the strength distribution in this channel. The effect on the deformation properties is also discussed. The quantitative effects on the linear response are small and, overall, the giant dipole energy remains unaffected. Calculations are compared to predictions from the (quasi)-particle random phase approximation and experimental data where available, finding good agreement.

preprint2014arXiv

Nuclear Dynamics with the Sky3D code

A description is presented of how to use the Sky3D time-dependent Hartree-Fock code to calculate giant monopole resonances. This requires modification to the code, and a step-by-step guide of how to make the necessary modification is given. An example of how to analyse the output of the code to obtain quantities of physics interest is included. Together, the modifications and the post-processing are intended to serve as a typical example of how the code, which was designed to be extendable to particular users' needs, can be extended.

preprint2013arXiv

Cause of the charge radius isotope shift at the \emph{N}=126 shell gap

We discuss the mechanism causing the `kink' in the charge radius isotope shift at the N=126 shell closure. The occupation of the 1$i_{11/2}$ neutron orbital is the decisive factor for reproducing the experimentally observed kink. We investigate whether this orbital is occupied or not by different Skyrme effective interactions as neutrons are added above the shell closure. Our results demonstrate that several factors can cause an appreciable occupation of the 1$i_{11/2}$ neutron orbital, including the magnitude of the spin-orbit field, and the isoscalar effective mass of the Skyrme interaction. The symmetry energy of the effective interaction has little influence upon its ability to reproduce the kink.

preprint2013arXiv

Continuum TDHF calculation of Isoscalar and Isovector Giant Monopole Resonances

We motivate and summarise some recent results in the application of formally exact boundary conditions in nuclear time-dependent Hartree-Fock calculations, making use of Laplace transformations to calculate the values of the wave functions at the boundaries. We have realised the method in the case of giant monopole resonances of spherically-symmetric nuclei, and present strength functions of O-16 and Ca-40 using a simplified version of the Skyrme force, showing that no artefacts from discretisation occur as contaminants

preprint2013arXiv

Extension of continuum time-dependent Hartree-Fock method to proton states

This paper deals with the solution of the spherically symmetric time-dependent Hartree-Fock approximation applied to nuclear giant monopole resonances in the small amplitude regime. The problem is spatially unbounded as the resonance state is in the continuum. The practical requirement to perform the calculation in a finite-sized spatial region yields an artificial boundary, which is not present physically. The question of how to ensure the boundary does not interfere with the internal solution, while keeping the overall calculation time low is studied. Here we propose an absorbing boundary condition scheme to handle the conflict. The derivation, via a Laplace transform method, and implementation is described. An inverse Laplace transform required by the absorbing boundaries is calculated using a method of non-linear least squares. The accuracy and efficiency of the scheme is tested and results presented to support the case that they are a effective way of handling the artificial boundary.

preprint2013arXiv

Quantum measurement and thermally assisted proton tunnelling

Despite compelling evidence to the contrary in recent years, the view still persists that quantum effects cannot survive very long within a warm, noisy and complex environment that washes out quantum effects at timescales far too short for any chemically or biochemically interesting processes. It is also assumed that as the temperature of the surrounding environment increases, so the efficiency of processes such as quantum tunnelling drops. One way of viewing this has been to invoke the quantum Zeno effect: that the watched pot never boils. In this work we show that the opposite is true. For a quite general open quantum system, a proton in an asymmetric double-well potential, the action of the environment is to enhance the tunnelling rate (an anti-Zeno effect). We compare two simple mathematical models to show that, over a specific temperature range, thermally enhanced quantum tunnelling is equivalent to increasing the frequency of a von Neumann-type measurement by the environment on the system.

preprint2013arXiv

The electric dipole response of $^{76}$Se above 4 MeV

The dipole response of $^{76}_{34}$Se in the energy range 4 to 9 MeV has been analyzed using a $(\vecγ,γ')$ polarized photon scattering technique, performed at the High Intensity $γ$-Ray Source facility, to complement previous work performed using unpolarized photons. The results of this work offer both an enhanced sensitivity scan of the dipole response and an unambiguous determination of the parities of the observed J=1 states. The dipole response is found to be dominated by $E1$ excitations, and can reasonably be attributed to a pygmy dipole resonance. Evidence is presented to suggest that a significant amount of directly unobserved excitation strength is present in the region, due to unobserved branching transitions in the decays of resonantly excited states. The dipole response of the region is underestimated when considering only ground state decay branches. We investigate the electric dipole response theoretically, performing calculations in a 3D cartesian-basis time-dependent Skyrme-Hartree-Fock framework.

preprint2013arXiv

The Lipkin-Meshkov-Glick model: 'quasi-local' quantum criticality in nuclear physics

Motivated by recent work on local quantum criticality in condensed matter systems, we study the Lipkin-Meshkov-Glick (LMG) model of nuclear physics as a simple model of a kind of 'quasi-local' quantum criticality. We identify a new crossover temperature, T*(V,W), between linear and nonlinear dynamics, which is analogous to the crossover between the renormalized classical and quantum critical regimes in the condensed-matter case. This temperature T* typically vanishes logarithmically as the quantum phase transition is approached, except near the quantum tricritical point where it becomes linear. We also note a further analogy with condensed-matter quantum criticality: the LMG model exhibits quantum order-by-disorder phenomena, of the type often associated with phase reconstruction near quantum critical points.

preprint2013arXiv

The TDHF Code Sky3D

The nuclear mean-field model based on Skyrme forces or related density functionals has found wide-spread application to the description of nuclear ground states, collective vibrational excitations, and heavy-ion collisions. The code Sky3D solves the static or dynamic equations in a three-dimensional Cartesian mesh with isolated or periodic boundary conditions and no further symmetry assumptions. Pairing can be included in the BCS approximation. The code is implemented with a view to allow easy modifications for including additional physics or special analysis of the results.

preprint2012arXiv

Continuum time-dependent Hartree-Fock for giant resonances in spherical nuclei

This paper deals with the solution of the spherically symmetric time-dependent Hartree-Fock approximation applied in the case of nuclear giant monopole resonances. The problem is spatially unbounded as the resonance state is in the continuum. The practical requirement to perform the calculation in a finite-sized spatial region results in a difficulty with the spatial boundary conditions. Here we propose a absorbing boundary condition scheme to handle the conflict. The derivation, via a Laplace transform method, and implementation is described. The accuracy and efficiency of the scheme is tested and the results presented to support the case that they are a effective way of handling the artificial boundary.

preprint2012arXiv

Do Skyrme forces that fit nuclear matter work well in finite nuclei?

A shortlist of Skyrme force parameterizations, recently found to have passed a series of constraints relating to nuclear matter properties is analyzed for their ability to reproduce data in finite nuclei. We analyse binding energies, isotope shifts and fission barriers. We find that the subset of forces have no common ability to reproduce (or otherwise) properties of finite nuclei, despite passing the extensive range of nuclear matter constraints.

preprint2012arXiv

Phase space structure and dynamics within the time-dependent Hartree-Fock approach

We study the equilibration and relaxation processes within the time-dependent Hartree-Fock approach using the Wigner distribution function. On the technical side we present a geometrically unrestricted framework which allows us to calculate the full six-dimensional Wigner distribution function. With the removal of geometrical constraints, we are now able to extend our previous phase-space analysis of heavy-ion collisions in the reaction plane to unrestricted mean-field simulations of nuclear matter on a three-dimensional Cartesian lattice. From the physical point of view we provide a quantitative analysis on the stopping power in TDHF. This is linked to the effect of transparency. For the medium-heavy $^{40}$Ca+$^{40}$Ca system we examine the impact of different parametrizations of the Skyrme force, energy-dependence, and the significance of extra time-odd terms in the Skyrme functional. For the first time, transparency in TDHF is observed for a heavy system, $^{24}$Mg+$^{208}$Pb.

preprint2012arXiv

Single-particle dissipation in TDHF studied from a phase-space perspective

We study dissipation and relaxation processes within the time-dependent Hartree-Fock approach using the Wigner distribution function. On the technical side we present a geometrically unrestricted framework which allows us to calculate the full six-dimensional Wigner distribution function. With the removal of geometrical constraints, we are now able to extend our previous phase-space analysis of heavy-ion collisions in the reaction plane to unrestricted mean-field simulations of nuclear matter on a three-dimensional Cartesian lattice. From the physical point of view we provide a quantitative analysis on the stopping power in TDHF. This is linked to the effect of transparency. For the medium-heavy $^{40}$Ca+$^{40}$Ca system we examine the impact of different parametrizations of the Skyrme force, energy-dependence, and the significance of extra time-odd terms in the Skyrme functional.

preprint2012arXiv

Skyrme Interaction and Nuclear Matter Constraints

This paper presents a detailed assessment of the ability of the 240 Skyrme interaction parameter sets in the literature to satisfy a series of criteria derived from macroscopic properties of nuclear matter in the vicinity of nuclear saturation density at zero temperature and their density dependence, derived by the liquid drop model, experiments with giant resonances and heavy-ion collisions. The objective is to identify those parameterizations which best satisfy the current understanding of the physics of nuclear matter over a wide range of applications. Out of the 240 models, only 16 are shown to satisfy all these constraints. Additional, more microscopic, constraints on density dependence of the neutron and proton effective mass beta-equilibrium matter, Landau parameters of symmetric and pure neutron nuclear matter, and observational data on high- and low-mass cold neutron stars further reduce this number to 5, a very small group of recommended Skyrme parameterizations to be used in future applications of the Skyrme interaction of nuclear matter related observables. Full information on partial fulfillment of individual constraints by all Skyrme models considered is given. The results are discussed in terms of the physical interpretation of the Skyrme interaction and the validity of its use in mean-field models. Future work on application of the Skyrme forces, selected on the basis of variables of nuclear matter, in Hartree-Fock calculation of properties of finite nuclei, is outlined.

preprint2012arXiv

Why is lead so kinky?

We revisit the problem of the kink in the charge radius shift of neutron-rich even lead isotopes. We show that the ability of a Skyrme force to reproduce the isotope shift is determined by the occupation of the neutron 1i11/2 orbital beyond N=126 and the corresponding change it causes to deeply-bound protons orbitals with a principal quantum number of 1. Given the observed position of the single-particle energies, one must either ensure occupation is allowed through correlations, or not demand that the single-particle energies agree with experimental values at the mean-field level.

preprint2011arXiv

Spectroscopic calculations of the low-lying structure in exotic Os and W isotopes

Structural evolution in neutron-rich Os and W isotopes is investigated in terms of the Interacting Boson Model (IBM) Hamiltonian determined by (constrained) Hartree-Fock-Bogoliubov (HFB) calculations with the Gogny-D1S Energy Density Functional (EDF). The interaction strengths of the IBM Hamiltonian are produced by mapping the potential energy surface (PES) of the Gogny-EDF with quadrupole degrees of freedom onto the corresponding PES of the IBM system. We examine the prolate-to-oblate shape/phase transition which is predicted to take place in this region as a function of neutron number $N$ within the considered Os and W isotopic chains. The onset of this transition is found to be more rapid compared to the neighboring Pt isotopes. The calculations also allow prediction of spectroscopic variables (excited state energies and reduced transition probabilities) which are presented for the neutron-rich $^{192,194,196}$W nuclei, for which there is only very limited experimental data available to date.

preprint2010arXiv

Influence of the symmetry energy on the giant monopole resonance of neutron-rich nuclei

We analyze the influence of the density dependence of the symmetry energy on the average excitation energy of the isoscalar giant monopole resonance (GMR) in stable and exotic neutron-rich nuclei by applying the relativistic extended Thomas-Fermi method in scaling and constrained calculations. For the effective nuclear interaction, we employ the relativistic mean field model supplemented by an isoscalar-isovector meson coupling that allows one to modify the density dependence of the symmetry energy without compromising the success of the model for binding energies and charge radii. The semiclassical estimates of the average energy of the GMR are known to be in good agreement with the results obtained in full RPA calculations. The present analysis is performed along the Pb and Zr isotopic chains. In the scaling calculations, the excitation energy is larger when the symmetry energy is softer. The same happens in the constrained calculations for nuclei with small and moderate neutron excess. However, for nuclei of large isospin the constrained excitation energy becomes smaller in models having a soft symmetry energy. This effect is mainly due to the presence of loosely-bound outer neutrons in these isotopes. A sharp increase of the estimated width of the resonance is found in largely neutron-rich isotopes, even for heavy nuclei, which is enhanced when the symmetry energy of the model is soft. The results indicate that at large neutron numbers the structure of the low-energy region of the GMR strength distribution changes considerably with the density dependence of the nuclear symmetry energy, which may be worthy of further characterization in RPA calculations of the response function.