BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Nicolas Schunck

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

ResearcherUnclaimed source record
nucl-thphysics.data-anphysics.comp-ph

Catalog footprint

What is connected

10works
3topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

10 published item(s)

preprint2022arXiv

Controlling extrapolations of nuclear properties with feature selection

Predictions of nuclear properties far from measured data are inherently imprecise because of uncertainties in our knowledge of nuclear forces and in our treatment of quantum many-body effects in strongly-interacting systems. While the model bias can be directly calculated when experimental data is available, only an estimate can be made in the absence of such measurements. Current approaches to compute the estimated bias quickly lose predictive power when input variables such as proton or neutron number are extrapolated, resulting in uncontrolled uncertainties in applications such as nucleosynthesis simulations. In this letter, we present a novel technique to identify the input variables of machine learning algorithms that can provide robust estimates of model bias. Our process is based on selecting input variables, or features, based on their probability distribution functions across the entire nuclear chart. We illustrate our approach on the problem of quantifying the model bias in nuclear binding energies calculated with Density Functional Theory (DFT). We show that feature selection can systematically improve theoretical predictions without increasing uncertainties.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Theory of Nuclear Fission

Atomic nuclei are quantum many-body systems of protons and neutrons held together by strong nuclear forces. Under the proper conditions, nuclei can break into two (sometimes three) fragments which will subsequently decay by emitting particles. This phenomenon is called nuclear fission. Since different fission events may produce different fragmentations, the end-products of all fissions that occurred in a small chemical sample of matter comprise hundreds of different isotopes, including $α$ particles, together with a large number of emitted neutrons, photons, electrons and antineutrinos. The extraordinary complexity of this process, which happens at length scales of the order of a femtometer, mostly takes less than a femtosecond but is not completely over until all the lingering $β$ decays have completed - which can take years - is a fascinating window into the physics of atomic nuclei. While fission may be more naturally known in the context of its technological applications, it also plays a pivotal role in the synthesis of heavy elements in astrophysical environments. In both cases, experimental measurements are not sufficient to provide complete data. Simulations are needed, yet at levels of accuracy and precision that pose formidable challenges to nuclear theory. The goal of this article is to provide a comprehensive overview of the theoretical methods employed in the description of nuclear fission.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fission of $^{240}$Pu with symmetry-restored density functional theory

Nuclear fission plays an important role in fundamental and applied science, from astrophysics to nuclear engineering, yet it remains a major challenge to nuclear theory. Theoretical methods used so far to compute fission observables rely on symmetry-breaking schemes where basic information on the number of particles, angular momentum, and parity of the fissioning nucleus is lost. In this work, we analyze the impact of restoring broken symmetries in the benchmark case of $^{240}$Pu.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Number of Particles in Fission Fragments

In current simulations of fission, the number of protons and neutrons in a given fission fragment is almost always obtained by integrating the total density of particles in the sector of space that contains the fragment. Because of the antisymmetry of the many-body wave function of the whole nucleus, this procedure systematically gives noninteger numbers of particles in the fragments. We introduce a novel sampling method to estimate rigorously the probability of finding $Z$ protons and $N$ neutrons in a fission fragment without resorting to projectors, which can sometimes give unwieldy results. When applied on standard Hartree-Fock-Bogoliubov many-body states, we show that our approach reproduces indeed the results of full particle-number projection. We then estimate the charge and mass number dispersion of several scission configurations in $^{240}$Pu with and without pairing correlations included. We show that odd-even effects in the charge probability naturally occur within our approach, which could explain the well-known odd-even staggering of charge distributions. Our method is applicable either in static calculations of scission configurations such as, e.g., in the macroscopic-microscopic approach or energy density functional theory, but also in explicitly time-dependent density functional theory simulations of fission.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Microscopic modeling of mass and charge distributions in the spontaneous fission of 240Pu

In this letter, we outline a methodology to calculate microscopically mass and charge distributions of spontaneous fission yields. We combine the multi-dimensional minimization of collective action for fission with stochastic Langevin dynamics to track the relevant fission paths from the ground-state configuration up to scission. The nuclear potential energy and collective inertia governing the tunneling motion are obtained with nuclear density functional theory in the collective space of shape deformations and pairing. We obtain a quantitative agreement with experimental data and find that both the charge and mass distributions in the spontaneous fission of 240Pu are sensitive both to the dissipation in collective motion and to adiabatic characteristics.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Bayesian Approach for Parameter Estimation and Prediction using a Computationally Intensive Model

Bayesian methods have been very successful in quantifying uncertainty in physics-based problems in parameter estimation and prediction. In these cases, physical measurements y are modeled as the best fit of a physics-based model $η(θ)$ where $θ$ denotes the uncertain, best input setting. Hence the statistical model is of the form $y = η(θ) + ε$, where $ε$ accounts for measurement, and possibly other error sources. When non-linearity is present in $η(\cdot)$, the resulting posterior distribution for the unknown parameters in the Bayesian formulation is typically complex and non-standard, requiring computationally demanding computational approaches such as Markov chain Monte Carlo (MCMC) to produce multivariate draws from the posterior. While quite generally applicable, MCMC requires thousands, or even millions of evaluations of the physics model $η(\cdot)$. This is problematic if the model takes hours or days to evaluate. To overcome this computational bottleneck, we present an approach adapted from Bayesian model calibration. This approach combines output from an ensemble of computational model runs with physical measurements, within a statistical formulation, to carry out inference. A key component of this approach is a statistical response surface, or emulator, estimated from the ensemble of model runs. We demonstrate this approach with a case study in estimating parameters for a density functional theory (DFT) model, using experimental mass/binding energy measurements from a collection of atomic nuclei. We also demonstrate how this approach produces uncertainties in predictions for recent mass measurements obtained at Argonne National Laboratory (ANL).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Derivative-free optimization for parameter estimation in computational nuclear physics

We consider optimization problems that arise when estimating a set of unknown parameters from experimental data, particularly in the context of nuclear density functional theory. We examine the cost of not having derivatives of these functionals with respect to the parameters. We show that the POUNDERS code for local derivative-free optimization obtains consistent solutions on a variety of computationally expensive energy density functional calibration problems. We also provide a primer on the operation of the POUNDERS software in the Toolkit for Advanced Optimization.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Error Analysis in Nuclear Density Functional Theory

Nuclear density functional theory (DFT) is the only microscopic, global approach to the structure of atomic nuclei. It is used in numerous applications, from determining the limits of stability to gaining a deep understanding of the formation of elements in the universe or the mechanisms that power stars and reactors. The predictive power of the theory depends on the amount of physics embedded in the energy density functional as well as on efficient ways to determine a small number of free parameters and solve the DFT equations. In this article, we discuss the various sources of uncertainties and errors encountered in DFT and possible methods to quantify these uncertainties in a rigorous manner.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Neutrino-nucleus coherent scattering as a probe of neutron density distributions

Neutrino-nucleus coherent elastic scattering provides a theoretically appealing way to measure the neutron part of nuclear form factors. Using an expansion of form factors into moments, we show that neutrinos from stopped pions can probe not only the second moment of the form factor (the neutron radius) but also the fourth moment. Using simple Monte Carlo techniques for argon, germanium, and xenon detectors of 3.5 tonnes, 1.5 tonnes, and 300 kg, respectively, we show that the neutron radii can be found with an uncertainty of a few percent when near a neutrino flux of $3\times10^{7}$ neutrinos/cm$^{2}$/s. If the normalization of the neutrino flux is known independently, one can determine the moments accurately enough to discriminate among the predictions of various nuclear energy functionals.

0 citations0 reviewsNo code linkedOpen access