Researcher profile

P. D. Stevenson

P. D. Stevenson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

14 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

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

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

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.

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

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.