Source author record

Kenneth J. Roche

Kenneth J. Roche 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

13works
6topics
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

13 published item(s)

preprint2016arXiv

Induced Fission of $^{240}$Pu within a Real-Time Microscopic Framework

We describe the fissioning dynamics of Pu240 from a configuration in the proximity of the outer fission barrier to full scission and the formation of the fragments within an implementation of density functional theory extended to superfluid systems and real-time dynamics. The fission fragments emerge with properties similar to those determined experimentally, while the fission dynamics appears to be quite complex, with many excited shape and pairing modes. The evolution is found to be much slower than previously expected, and the ultimate role of the collective inertia is found to be negligible in this fully nonadiabatic treatment of nuclear dynamics, where all collective degrees of freedom (CDOF) are included (unlike adiabatic treatments with a small number of CDOF).

preprint2016arXiv

Object Kinetic Monte Carlo Simulations of Radiation Damage in Bulk Tungsten Part-II: With a PKA Spectrum Corresponding to 14-MeV Neutrons

Object kinetic Monte Carlo was employed to study the effect of dose rate on the evolution of vacancy microstructure in polycrystalline tungsten under neutron bombardment. The evolution was followed up to 1.0 displacement per atom (dpa) with point defects generated in accordance with a primary knock-on atom (PKA) spectrum corresponding to 14-MeV neutrons. The present study includes the effect of grain size (2.0 and 4.0 $μ$m) but excludes the impact of transmutation or pre-existing defects beyond grain boundary sinks. Vacancy cluster density increases with dose rate, while the density of vacancies decreases. Consequently, the average vacancy cluster size and the fraction of vacancies in visible clusters decrease with increasing dose rate. The density of vacancies and vacancy clusters decrease with grain size such that the average size of the clusters remains similar. However, the average size is larger for larger grains at dose rates < 4.5 x 10-7 dpa/s. The trend of vacancy accumulation as a function of dose, dose rate, and grain size is similar to that obtained with the High Flux Isotope Reactor (HFIR) PKA spectrum. However, the amount of vacancy accumulation and the vacancy microstructure are quite different. Compared to the HFIR case, we find that even though the dose rates are 2.5 times higher, the density of vacancies and the average vacancy cluster sizes are lower. In addition, a void lattice forms only for the lowest two dose rates (4.5 x 10-8 and 4.5 x 10-9 dpa/s). In contrast, a void lattice formed at all dose rates studied using the HFIR PKA spectrum. We discuss in detail the factors that lead to these different microstructures.

preprint2016arXiv

Object Kinetic Monte Carlo Simulations of Radiation Damage in Neutron-Irradiated Tungsten Part-I: Neutron Flux with a PKA Spectrum Corresponding to the High-flux Isotope Reactor

Object kinetic Monte Carlo simulations were performed to study the impact of varying dose rate and grain size up to a dose of 1.0 dpa in pure, polycrystalline tungsten, subjected to a neutron irradiation having a PKA spectrum corresponding to the High Flux Isotope Reactor. The present study models defect cluster accumulation in tungsten, but does not consider the impact of transmutation or pre-existing defects beyond the grain boundary sinks, with varying grain size. With increasing dose rate, the vacancy cluster density increases, while the number density of vacancies decreases. Accordingly, the average vacancy cluster size and the fraction of vacancies that are part of visible clusters decreases with increasing dose rate. With increasing grain size, both the number densities of vacancies and vacancy clusters decrease, while both the fraction of vacancies in visible clusters and the average vacancy cluster size increase. This is caused by the pseudo-ripening of the vacancy clusters due to the longer-lived self-interstitial clusters in larger grains. The spatial ordering of vacancy clusters along {110} planes was observed for both grain sizes and all dose rates studied. Interplanar spacing increases with grain size; however, no clear dependence on dose or dose rate was observed. The results of this study show that 1D diffusion of self-interstitial clusters, while necessary, is not sufficient to form a void lattice, and that the diffusion of vacancies is also required. A methodology is suggested for choosing the simulation box dimensions so as to represent more faithfully the effects of one-dimensional migrating self-interstitial-atom clusters.

preprint2014arXiv

Displacement cascades and defects annealing in tungsten, Part I: defect database from molecular dynamics simulations

Molecular dynamics simulations have been used to generate a comprehensive database of surviving defects due to displacement cascades in bulk tungsten. Twenty one data points of primary knock-on atom (PKA) energies ranging from 100 eV (sub-threshold energy) to 100 keV ($\sim$780$\times E_d$, where $E_d$ is the average displacement threshold energy) have been completed at 300 K, 1025 K and 2050 K. Within this range of PKA energies, two regimes of power-law energy-dependence of the defect production are observed. A distinct power-law exponent characterizes the number of Frenkel pairs produced within each regime. The two regimes intersect at a transition energy which occurs at approximately 250$\times E_d$. The transition energy also marks the onset of the formation of large self-interstitial atom (SIAs) clusters (size 14 or more). The observed defect clustering behavior is asymmetric, with SIA clustering increasing with temperature, while the vacancy clustering decreases. This asymmetry increases with temperature such that at 2050 K ($\sim 0.5 T_m$) practically no large vacancy clusters are formed, meanwhile large SIA clusters appear in all simulations. The implication of such asymmetry on the long-term defect survival and damage accumulation is discussed. In addition, rare $<$100$>${110} SIA loops are observed.

preprint2014arXiv

Displacement cascades and defects annealing in tungsten, Part II: Object kinetic Monte Carlo Simulation of Tungsten Cascade Aging

We describe the results of object kinetic Monte Carlo (OKMC) simulations of the annealing of primary cascade damage in bulk tungsten using a comprehensive database of cascades obtained from molecular dynamics [1] as a function of primary knock-on atom (PKA) energy and direction, and temperatures of 300, 1025 and 2050 K. An increase in SIA clustering but decrease in vacancy clustering with temperature combined with disparate mobilities of SIAs versus vacancies causes an interesting temperature effect on cascade annealing, which is quite different from what one would expect. The annealing efficiency (ratio of number of defects after and before annealing) exhibits an inverse U-shape curve as a function of temperature. In addition, we will also describe the capabilities of our newly developed OKMC code; KSOME (kinetic simulations of microstructure evolution) used to carryout these simulations

preprint2014arXiv

Life Cycle of Superfluid Vortices and Quantum Turbulence in the Unitary Fermi Gas

The unitary Fermi gas (UFG) offers an unique opportunity to study quantum turbulence both experimentally and theoretically in a strongly interacting fermionic superfluid. It yields to accurate and controlled experiments, and admits the only dynamical microscopic description via time-dependent density functional theory (DFT) - apart from dilute bosonic gases - of the crossing and reconnection of superfluid vortex lines conjectured by Feynman in 1955 to be at the origin of quantum turbulence in superfluids at zero temperature. We demonstrate how various vortex configurations can be generated by using well established experimental techniques: laser stirring and phase imprinting. New imagining techniques demonstrated by the MIT group [Ku et al. arXiv:1402.7052] should be able to directly visualize these crossings and reconnections in greater detail than performed so far in liquid helium. We demonstrate the critical role played by the geometry of the trap in the formation and dynamics of a vortex in the UFG and how laser stirring and phase imprint can be used to create vortex tangles with clear signatures of the onset of quantum turbulence.

preprint2014arXiv

Quantized Superfluid Vortex Rings in the Unitary Fermi Gas

In a recent article, Yefsah et al. [Nature 499, 426 (2013)] report the observation of an unusual excitation in an elongated harmonically trapped unitary Fermi gas. After phase imprinting a domain wall, they observe oscillations almost an order of magnitude slower than predicted by any theory of domain walls which they interpret as a "heavy soliton" of inertial mass some 200 times larger than the free fermion mass or 50 times larger than expected for a domain wall. We present compelling evidence that this "soliton" is instead a quantized vortex ring by showing that the main aspects of the experiment can be naturally explained within the framework of time-dependent superfluid DFT.

preprint2013arXiv

Computational Nuclear Quantum Many-Body Problem: The UNEDF Project

The UNEDF project was a large-scale collaborative effort that applied high-performance computing to the nuclear quantum many-body problem. UNEDF demonstrated that close associations among nuclear physicists, mathematicians, and computer scientists can lead to novel physics outcomes built on algorithmic innovations and computational developments. This review showcases a wide range of UNEDF science results to illustrate this interplay.

preprint2013arXiv

Cooper pairing above the critical temperature in a unitary Fermi gas

We present an ab initio determination of spin responses of the unitary Fermi gas, based on finite temperature quantum Monte Carlo calculations and the Kubo linear-response formalism. We determine the temperature dependence of the spin susceptibility and the spin conductivity. We show that both quantities exhibit suppression above the critical temperature of the superfluid to normal phase transition due to presence of the Cooper pairs. The spin diffusion transport coefficient does not display the existence of a minimum in the vicinity of the critical temperature and it drops to very low values D_s approx 0.8hbar/m in the superfluid phase. All these spin observables show a smooth and monotonic behavior with temperature when crossing the critical temperature T_c, until the Fermi liquid regime is attained at the temperature T*, where the pseudogap regime disappears.

preprint2013arXiv

Quantum Friction: Cooling Quantum Systems with Unitary Time Evolution

We introduce a type of quantum dissipation -- local quantum friction -- by adding to the Hamiltonian a local potential that breaks time-reversal invariance so as to cool the system. Unlike the Kossakowski-Lindblad master equation, local quantum friction directly effects unitary evolution of the wavefunctions rather than the density matrix: it may thus be used to cool fermionic many-body systems with thousands of wavefunctions that must remain orthogonal. In addition to providing an efficient way to simulate quantum dissipation and non-equilibrium dynamics, local quantum friction coupled with adiabatic state preparation significantly speeds up many-body simulations, making the solution of the time-dependent Schrödinger equation significantly simpler than the solution of its stationary counterpart.

preprint2013arXiv

The temperature evolution of the shear viscosity in a unitary Fermi gas

We present an ab initio determination of the shear viscosity for the unitary Fermi gas based on finite temperature quantum Monte Carlo (QMC) calculations and the Kubo linear-response formalism. The results are confronted with the bound for the shear viscosity originating from hydrodynamic fluctuations. Assuming smoothness of the frequency dependent shear viscosity eta(omega), we show that the bound is violated in the low temperature regime and the violation occurs simultaneously with the onset of the Cooper paring in the system. In order to preserve the hydrodynamic bound in QMC eta(omega) has to possess a sharp structure located in the vicinity of zero frequency which is not resolved by an analytic continuation procedure.

preprint2012arXiv

Quantum Shock Waves and Domain Walls in the Real-Time Dynamics of a Superfluid Unitary Fermi Gas

We show that in the collision of two superfluid fermionic atomic clouds one observes the formation of quantum shock waves as discontinuities in the number density and collective flow velocity. Domain walls, which are topological excitations of the superfluid order parameter, are also generated and exhibit abrupt phase changes by $π$ and slower motion than the shock waves. The domain walls are distinct from the gray soliton train or number density ripples formed in the wake of the shock waves and observed in the collisions of superfluid bosonic atomic clouds. Domain walls with opposite phase jumps appear to collide elastically.

preprint2010arXiv

Generation and Dynamics of Quantized Vortices in a Unitary Fermi Superfluid

Superfluidity and superconductivity are remarkable manifestations of quantum coherence at a macroscopic scale. The dynamics of superfluids has dominated the study of these systems for decades now, but a comprehensive theoretical framework is still lacking. We introduce a local extension of the time-dependent density functional theory to describe the dynamics of fermionic superfluids. Within this approach one can correctly represent vortex quantization, generation, and dynamics, the transition from a superfluid to a normal phase and a number of other large amplitude collective modes which are beyond the scope of two-fluid hydrodynamics, Ginzburg-Landau and/or Gross-Pitaevskii approaches. We illustrate the power of this approach by studying the generation of quantized vortices, vortex rings, vortex reconnection, and transition from a superfluid to a normal state in real time for a unitary Fermi gas. We predict the emergence of a new qualitative phenomenon in superfluid dynamics of gases, the existence of stable superfluidity when the systems are stirred with velocities significantly exceeding the nominal Landau critical velocity in these systems.