Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion problems”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Steady-State Thermal-Hydraulic Analysis and Bowing Reactivity Evaluation Methods Based on Neutron and Gamma Transport Calculations

With the advancement of computer technology, the variational nodal transport code VARIANT is now routinely used for sodium-cooled fast reactor (SFR) design analysis. In particular, VARIANT transport calculations are usually performed for whole core neutronics analyses and for fuel cycle analyses with REBUS-3. However, the steady-state thermal hydraulics analysis with the SE2-ANL coed and the assembly bowing reactivity calculation with the structural analysis with the NUBOW-3D code are still performed with the DIF3D finite difference diffusion theory code. Furthermore, the DIF3D diffusion calculation for these applications is performed with only six triangular meshes per hexagonal assembly, and the SE2-ANL code determines the pin power distributions approximately by assuming that the power distribution within an assembly is separable in the radial and axial directions. In order to overcome the limitations and to improve the accuracy of the existing methods, it was proposed to develop the relevant computational methods based on the transport calculations with the VARIANT and PROTEUS-SN codes. The objectives of this study are to develop a new heating calculation method based on the coupled neutron and gamma transport calculations with VARIANT, to implement the new heating calculation method into the SE2-ANL steady-state thermal-hydraulics analysis code with improved numerical algorithms and additional capabilities for automatic flow allocation calculations, and to develop a method to calculate bowing reactivity coefficients based on the VARIANT core solutions and the PROTEUS-SN assembly solutions, which can be used in NUBOW-3D calculations. A new, coupled neutron and gamma heating calculation procedure has been developed based on VARIANT transport calculations. In the new heating calculation procedure, the neutron flux distribution is first determined by solving an eigenvalue transport problem using the VARIANT transport code. Using the calculated neutron flux, the intra-nodal gamma source distributions are calculated to be consistent with the neutron flux distribution with the GAMSOR code. For this, the GAMSOR code has been updated to generate the intra-nodal gamma source distributions in the form of the VARSRC dataset of VARIANT. With this gamma source distribution, the gamma flux distribution is determined by solving a fixed source gamma transport problem using VARIANT again. With the neutron and gamma flux distributions, the heat generation rates are evaluated at individual fuel pins and assembly duct walls, without assuming separable axial and radial profiles for the intra-assembly power distribution.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Addressing Critical Problems in Materials Science Through Multiscale and Multimode Characterization (Project 1); Characterization and Optimization of Novel Triple-Conducting Oxide Materials for Energy Applications (Project 2) (CRADA Final Report)

PROJECT 1: Address critical problems in materials science and simultaneously advance the state-of-the-art in multiscale and multimode characterization using the combined advanced analytical capabilities and expertise of Colorado School of Mines (CSM) and the National Renewable Energy Laboratory (NREL). The primary effort of the Phase I of this CRADA is to establish the International Center for Multiscale Characterization using shared resources at both NREL and CSM. Phase II will focus on capability development and marketing, choosing candidate materials science issues in the areas of structure imaging, chemical composition mapping, and correlating properties and performance of materials for impact in energy-related, environmental and critical materials areas. The CRADA will be modified to include specific topics of concern in materials science to industry member partners. Advanced analytical capabilities and expertise at CSM and NREL will be used to advance materials understanding and performance through characterization of multiscale phenomena including structural imaging, chemical composition mapping, and other techniques correlating properties and performance of materials. PROJECT 2: As part of the International Center for Materials Characterization, work under Modification #1 will be led by Colorado School of Mines (CSM), working in collaboration with NREL staff to mentor and advise CSM postdoctoral researchers on set up of diffusion annealing experiments. The purpose of the modification is to provide for NREL staff to mentor and advise CSM postdoctoral researchers on set up of diffusion annealing experiments, including mentoring and advising the CSM-NREL team on proper Secondary Ion Mass Spectrometry (SIMS) data analysis as needed. SIMS measurements of 10-20 samples will be performed at NREL during the project duration.

08 HYDROGEN↗
DOE OSTI2017

Influence of Sedimentation Length on the Convective Stability of a Colloidal Suspension

The convective stability of a colloidal suspension is studied in the case when the vertical dimension of the cavity is less than or comparable with the sedimentation length of nanoparticles. The analysis is carried out within the Boussinesq approximation on the basis of a modified model that takes into account the dependence of a thermodiffusion flow on the local value of impurity concentration. A new parameter of the problem is the ratio of the sedimentation length to the vertical dimension of the cavity. For a quiescent colloidal suspension, exact and approximate (in the case of small concentrations) solutions are obtained that describe the distributions of nanoparticles. A transformation is obtained that allows one to investigate the convective stability of a colloidal suspension stratified in the gravitational field by the Galerkin method with a set of simple trial functions. Instability boundaries and the characteristics of critical perturbations are determined. It is shown that, in the case of negative thermodiffusion, a decrease in the sedimentation length leads to a decrease in the convection threshold and the frequency of neutral oscillations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

High thermal stability and sluggish crystallization kinetics of high-entropy bulk metallic glasses

Metallic glasses are metastable and their thermal stability is critical for practical applications, particularly at elevated temperatures. The conventional bulk metallic glasses (BMGs), though exhibiting high glass-forming ability (GFA), crystallize quickly when being heated to a temperature higher than their glass transition temperature. This problem may potentially be alleviated due to the recent developments of high-entropy (or multi-principle-element) bulk metallic glasses (HE-BMGs). In this work, we demonstrate that typical HE-BMGs, i.e., ZrTiHfCuNiBe and ZrTiCuNiBe, have higher kinetic stability, as compared with the benchmark glass Vitreoy1 (Zr{sub 41.2}Ti{sub 13.8}Cu{sub 12.5}Ni{sub 10}Be{sub 22.5}) with a similar chemical composition. The measured activation energy for glass transition and crystallization of the HE-BMGs is nearly twice that of Vitreloy 1. Moreover, the sluggish crystallization region ΔT{sub pl-pf}, defined as the temperature span between the last exothermic crystallization peak temperature T{sub pl} and the first crystallization exothermic peak temperature T{sub pf}, of all the HE-BMGs is much wider than that of Vitreloy 1. In addition, high-resolution transmission electron microscopy characterization of the crystallized products at different temperatures and the continuous heating transformation diagram which is proposed to estimate the lifetime at any temperature below the melting point further confirm high thermal stability of the HE-BMGs. Surprisingly, all the HE-BMGs show a small fragility value, which contradicts with their low GFA, suggesting that the underlying diffusion mechanism in the liquid and the solid of HE-BMGs is different.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Optimal white-noise stochastic forcing for linear models of turbulent channel flow

In the present study an optimisation problem is formulated to determine the forcing of an eddy-viscosity-based linearised Navier–Stokes model in channel flow at $Re_\tau \approx 5200$ ( $Re_\tau$ is the friction Reynolds number), where the forcing is white-in-time and spatially decorrelated. The objective functional is prescribed such that the forcing drives a response to best match a set of velocity spectra from direct numerical simulation (DNS), as well as remaining sufficiently smooth. Strong quantitative agreement is obtained between the velocity spectra from the linear model with optimal forcing and from DNS, but only qualitative agreement between the Reynolds shear stress co-spectra from the model and DNS. The forcing spectra exhibit a level of self-similarity, associated with the primary peak in the velocity spectra, but they also reveal a non-negligible amount of energy spent in phenomenologically mimicking the non-self-similar part of the velocity spectra associated with energy cascade. By exploiting linearity, the effect of the individual forcing components is assessed and the contributions from the Orr mechanism and the lift-up effect are also identified. Finally, the effect of the strength of the eddy viscosity on the optimisation performance is investigated. The inclusion of the eddy viscosity diffusion operator is shown to be essential in modelling of the near-wall features, while still allowing the forcing of the self-similar primary peak. In particular, reducing the strength of the eddy viscosity results in a considerable increase in the near-wall forcing of wall-parallel components.

Mechanics↗
DOE OSTI2016

Neutrino physics with JUNO

The Jiangmen Underground Neutrino Observatory (JUNO), a 20 kton multi-purpose underground liquid scintillator detector, was proposed with the determination of the neutrino mass hierarchy (MH) as a primary physics goal. The excellent energy resolution and the large fiducial volume anticipated for the JUNO detector offer exciting opportunities for addressing many important topics in neutrino and astro-particle physics. In this document, we present the physics motivations and the anticipated performance of the JUNO detector for various proposed measurements. Following an introduction summarizing the current status and open issues in neutrino physics, we discuss how the detection of antineutrinos generated by a cluster of nuclear power plants allows the determination of the neutrino MH at a 3–4σ significance with six years of running of JUNO. The measurement of antineutrino spectrum with excellent energy resolution will also lead to the precise determination of the neutrino oscillation parameters ${\mathrm{sin}}^{2}{\theta }_{12}$, ${\rm{\Delta }}{m}_{21}^{2}$, and $| {\rm{\Delta }}{m}_{{ee}}^{2}| $ to an accuracy of better than 1%, which will play a crucial role in the future unitarity test of the MNSP matrix. The JUNO detector is capable of observing not only antineutrinos from the power plants, but also neutrinos/antineutrinos from terrestrial and extra-terrestrial sources, including supernova burst neutrinos, diffuse supernova neutrino background, geoneutrinos, atmospheric neutrinos, and solar neutrinos. As a result of JUNO's large size, excellent energy resolution, and vertex reconstruction capability, interesting new data on these topics can be collected. For example, a neutrino burst from a typical core-collapse supernova at a distance of 10 kpc would lead to ~5000 inverse-beta-decay events and ~2000 all-flavor neutrino–proton ES events in JUNO, which are of crucial importance for understanding the mechanism of supernova explosion and for exploring novel phenomena such as collective neutrino oscillations. Detection of neutrinos from all past core-collapse supernova explosions in the visible universe with JUNO would further provide valuable information on the cosmic star-formation rate and the average core-collapse neutrino energy spectrum. Antineutrinos originating from the radioactive decay of uranium and thorium in the Earth can be detected in JUNO with a rate of ~400 events per year, significantly improving the statistics of existing geoneutrino event samples. Atmospheric neutrino events collected in JUNO can provide independent inputs for determining the MH and the octant of the ${\theta }_{23}$ mixing angle. Detection of the 7Be and 8B solar neutrino events at JUNO would shed new light on the solar metallicity problem and examine the transition region between the vacuum and matter dominated neutrino oscillations. Regarding light sterile neutrino topics, sterile neutrinos with ${10}^{-5}\;{{\rm{eV}}}^{2}\lt {\rm{\Delta }}{m}_{41}^{2}\lt {10}^{-2}\;{{\rm{eV}}}^{2}$ and a sufficiently large mixing angle ${\theta }_{14}$ could be identified through a precise measurement of the reactor antineutrino energy spectrum. Meanwhile, JUNO can also provide us excellent opportunities to test the eV-scale sterile neutrino hypothesis, using either the radioactive neutrino sources or a cyclotron-produced neutrino beam. The JUNO detector is also sensitive to several other beyondthe-standard-model physics. Examples include the search for proton decay via the $p\to {K}^{+}+\bar{\nu }$ decay channel, search for neutrinos resulting from dark-matter annihilation in the Sun, search for violation of Lorentz invariance via the sidereal modulation of the reactor neutrino event rate, and search for the effects of non-standard interactions. The proposed construction of the JUNO detector will provide a unique facility to address many outstanding crucial questions in particle and astrophysics in a timely and cost-effective fashion. Furthermore, it holds the great potential for further advancing our quest to understanding the fundamental properties of neutrinos, one of the building blocks of our Universe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Experimental Observation of Thin-shell Instability in a Collisionless Plasma

We report on the experimental observation of the instability of a plasma shell, which formed during the expansion of a laser-ablated plasma into a rarefied ambient medium. By means of a proton radiography technique, the evolution of the instability is temporally and spatially resolved on a timescale much shorter than the hydrodynamic one. The density of the thin shell exceeds that of the surrounding plasma, which lets electrons diffuse outward. An ambipolar electric field grows on both sides of the thin shell that is antiparallel to the density gradient. Ripples in the thin shell result in a spatially varying balance between the thermal pressure force mediated by this field and the ram pressure force that is exerted on it by the inflowing plasma. This mismatch amplifies the ripples by the same mechanism that drives the hydrodynamic nonlinear thin-shell instability (NTSI). Our results thus constitute the first experimental verification that the NTSI can develop in colliding flows.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Measurement of Kinetic Parameters by Oscillated Experiments

The availability of reliable values of the kinetic parameters (effective delayed neutron fraction {beta} and effective mean neutron generation time {Lambda}) is important for reactor kinetic analyses and for reactivity monitoring when adopting the point kinetic approach. The kinetic parameters may be obtained either by simulations (Monte Carlo or deterministic) of the system or by direct measurements. Several methods have been proposed for the direct experimental determination of the kinetic parameters, e.g. control rod drops and noise analysis approaches, e.g. In this work, the problem of the determination of kinetic parameters in subcritical systems is tackled. The system is subjected to an oscillated source and the value of the ratio {beta}/{Lambda} is related to a measurable quantity when the system response reaches an equilibrium condition. In the analysis presented, a pure frequency is considered: when analyzing more complicated source periodic behaviours, a Fourier analysis can allow to isolate the contributions of individual frequencies. The phase difference between the oscillating source and the system response is adopted as the experimental observable. The method for the interpretation of the measurements is based on the standard point kinetic model. The response of the system is obtained by localized neutron detectors and, therefore, spatial and spectral effects may also affect the results. This drawback needs some further discussion and will be considered in the future work. The phase shift between the source signal and the system response may be analytically retrieved by Laplace transforming the point kinetic equations in the presence of an oscillated source. An explicit formula relates the phase shift to the ratio {beta}/{Lambda} and to the frequency of the source. Based on this relationship, the ratio {beta}/{Lambda} can be determined. To validate the technique, the experimental data are simulated by a numerical code solving the time-dependent diffusion equations for a source-driven system. The work allows to assess the possibility of the technique and its limits of applicability. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

MODELING DUST IN THE MAGELLANIC CLOUDS

We model the extinction profiles observed in the Small and Large Magellanic clouds with a synthetic population of dust grains consisting of core-mantle particles and a collection of free-flying polycyclic aromatic hydrocarbons (PAHs). All different flavors of the extinction curves observed in the Magellanic Clouds (MCs) can be described by the present model, which has been previously (successfully) applied to a large sample of diffuse and translucent lines of sight in the Milky Way. We find that in the MCs the extinction produced by classical grains is generally larger than absorption by PAHs. Within this model, the nonlinear far-UV rise is accounted for by PAHs, whose presence in turn is always associated with a gap in the size distribution of classical particles. This hints either at a physical connection between (e.g., a common cause for) PAHs and the absence of middle-sized dust particles or the need for an additional component in the model that can account for the nonlinear far-UV rise without contributing to the UV bump at ∼217 nm such as, e.g., nanodiamonds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Fixed-node diffusion Monte Carlo description of nitrogen defects in zinc oxide

Using the fixed-node diffusion Monte Carlo (FN-DMC) method, we evaluate the formation energies and charge transition levels of substitutional nitrogen defects in the wide-band-gap semiconductor zinc oxide. The use of a direct-solution, many-body approach inherently secures a good description of electron-electron interactions, achieving high accuracy without adjustable parameters. According to FN-DMC nitrogen is a deep acceptor with a charge transition level 1.0(3) eV above the valence-band maximum when 72-atom supercells are used. This result falls on the lower end of typically reported hybrid density functional results for the same size supercells, which range from 1.0 to 1.8 eV. Further, residual finite-size effects due to charged defect image interactions in the 72-atom supercells are estimated by supercell extrapolation within hybrid density functional theory. When the finite-size correction is included, we obtain a deep acceptor at 1.6(3) eV. This result is in good agreement with recent experimental measurements. We also analyze the local compressibility of charge according to FN-DMC and common density functionals and find that the use of hybrid functionals obtains compressibilities in better agreement with the many-body theory. Our work illustrates the application of the FN-DMC method to a challenging point defect problem, demonstrating that uncertainties and approximations can be well controlled.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Self-screening of the quantum confined Stark effect by the polarization induced bulk charges in the quantum barriers

InGaN/GaN light-emitting diodes (LEDs) grown along the polar orientations significantly suffer from the quantum confined Stark effect (QCSE) caused by the strong polarization induced electric field in the quantum wells, which is a fundamental problem intrinsic to the III-nitrides. Here, we show that the QCSE is self-screened by the polarization induced bulk charges enabled by designing quantum barriers. The InN composition of the InGaN quantum barrier graded along the growth orientation opportunely generates the polarization induced bulk charges in the quantum barrier, which well compensate the polarization induced interface charges, thus avoiding the electric field in the quantum wells. Consequently, the optical output power and the external quantum efficiency are substantially improved for the LEDs. The ability to self-screen the QCSE using polarization induced bulk charges opens up new possibilities for device engineering of III-nitrides not only in LEDs but also in other optoelectronic devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Simultaneous Momentum and Position Measurement and the Instrumental Weyl-Heisenberg Group

The canonical commutation relation, [Q,P]=iℏ, stands at the foundation of quantum theory and the original Hilbert space. The interpretation of P and Q as observables has always relied on the analogies that exist between the unitary transformations of Hilbert space and the canonical (also known as contact) transformations of classical phase space. Now that the theory of quantum measurement is essentially complete (this took a while), it is possible to revisit the canonical commutation relation in a way that sets the foundation of quantum theory not on unitary transformations but on positive transformations. This paper shows how the concept of simultaneous measurement leads to a fundamental differential geometric problem whose solution shows us the following. The simultaneous P and Q measurement (SPQM) defines a universal measuring instrument, which takes the shape of a seven-dimensional manifold, a universal covering group we call the instrumental Weyl-Heisenberg (IWH) group. The group IWH connects the identity to classical phase space in unexpected ways that are significant enough that the positive-operator-valued measure (POVM) offers a complete alternative to energy quantization. Five of the dimensions define processes that can be easily recognized and understood. The other two dimensions, the normalization and phase in the center of the IWH group, are less familiar. The normalization, in particular, requires special handling in order to describe and understand the SPQM instrument.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Quantum Monte Carlo benchmarking of large noncovalent complexes in the L7 benchmark set

We have used diffusion Monte Carlo (DMC) to perform calculations on the L7 benchmark set. DMC is a stochastic numerical integration scheme in real-space and part of a larger set of quantum Monte Carlo methods. The L7 set was designed to test the ability of electronic structure methods to include dispersive interactions. While the agreement between DMC and quantum-chemical state-of-the-art methods is excellent for some of the structures, there are significant differences in others. In contrast to wavefunction-based quantum chemical methods, DMC is a first-principle many-body method with the many-body wavefunction evolving in real space. It includes explicitly all electron–electron interactions and is relatively insensitive to the size of the basis set.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Enhancing Unknown Waveform Detection by Learning Intra and Inter-domain Dependencies with Advanced Attention Fusion Mechanisms

Detection of unknown waveforms in mission-critical communications is a crucial area of interest for the Department of Energy (DoE). Traditional methods and recent deep learning-based approaches often assume that the training set includes all possible classes, which is impractical for detecting new waveforms. This limitation gives rise to the problem of open-set recognition (OSR), which involves correctly identifying known classes while detecting and rejecting unknown or unseen classes. To address this limitation, we propose a novel dual-domain complex-valued neural architecture that jointly processes time-domain and frequency-domain signal representations using transformer mechanisms. A transformer model is a deep learning architecture that uses self-attention mechanisms to process and learn relationships in sequential data. Our model employs a cosine similarity loss to extract domain-specific features and incorporates a transformer architecture in the latent space to weigh the importance of different features from the time and frequency domains. The transformer layer includes stacked self-attention and cross-attention modules to learn intra-domain and inter-domain dependencies, creating a more holistic signal representation. An attention-based fusion module intelligently combines the time and frequency-domain features using multi-head attention, enabling the network to learn the optimal feature for each domain in each input signal. Quantitative results demonstrate the impact of these architectural choices on overall performance, showing significant improvement after incorporating self and cross-attention modules and using complex attention fusion over simple weighted fusion. Our ongoing work will focus on addressing the limitations of threshold-based OSR methods by developing a novel generative framework that integrates a conditional diffusion probabilistic model (DPM). DPM is a generative framework that learns to synthesize complex data by reversing a gradual noising process using a neural network trained to denoise step-by-step. Our goal is to leverage the inherent strengths of DPMs for identifying unknown signals more robustly. One primary advantage of using a DPM is its ability to provide a more reliable anomaly score based on the model's reconstruction error, rather than relying solely on classifier confidence. Additionally, the iterative denoising process of DPMs makes this approach naturally resilient to low Signal-to-Noise Ratio (SNR) conditions, where traditional methods often fail. By implementing this generative framework, we aim to enhance the model's capability to accurately detect unknown waveforms and maintain performance in challenging environments.

99 - GENERAL AND MISCELLANEOUS↗
DOE OSTI2013

MONTE CARLO SIMULATIONS OF NONLINEAR PARTICLE ACCELERATION IN PARALLEL TRANS-RELATIVISTIC SHOCKS

We present results from a Monte Carlo simulation of a parallel collisionless shock undergoing particle acceleration. Our simulation, which contains parameterized scattering and a particular thermal leakage injection model, calculates the feedback between accelerated particles ahead of the shock, which influence the shock precursor and 'smooth' the shock, and thermal particle injection. We show that there is a transition between nonrelativistic shocks, where the acceleration efficiency can be extremely high and the nonlinear compression ratio can be substantially greater than the Rankine-Hugoniot value, and fully relativistic shocks, where diffusive shock acceleration is less efficient and the compression ratio remains at the Rankine-Hugoniot value. This transition occurs in the trans-relativistic regime and, for the particular parameters we use, occurs around a shock Lorentz factor γ{sub 0} = 1.5. We also find that nonlinear shock smoothing dramatically reduces the acceleration efficiency presumed to occur with large-angle scattering in ultra-relativistic shocks. Our ability to seamlessly treat the transition from ultra-relativistic to trans-relativistic to nonrelativistic shocks may be important for evolving relativistic systems, such as gamma-ray bursts and Type Ibc supernovae. We expect a substantial evolution of shock accelerated spectra during this transition from soft early on to much harder when the blast-wave shock becomes nonrelativistic.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Quantal diffusion approach for multinucleon transfer processes in the 58,64 Ni + 208 Pb reactions: Toward the production of unknown neutron-rich nuclei

In recent years, substantial efforts have been made for the study of multinucleon transfer reactions at energies around the Coulomb barrier both experimentally and theoretically, aiming at the production of unknown neutron-rich heavy nuclei. It is crucial to provide reliable theoretical predictions based on microscopic theories with sufficient predictive power. Purpose: This article aims to clarify the applicability of the quantal diffusion approach based on the stochastic mean-field (SMF) theory for multinucleon transfer processes. Isotope production cross sections are evaluated for the reactions of 64 Ni + 208 Pb at E c.m. = 268 MeV and 58 Ni + 208 Pb at E c.m. = 270 MeV and are compared with available experimental data. Methods: Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are carried out for a range of initial orbital angular momenta with Skyrme SLy4d functional. Quantal diffusion equations, derived based on the SMF theory, for variances and covariance of neutron and proton numbers of reaction products are solved, with microscopic drift and diffusion coefficients obtained from time evolution of occupied single-particle orbitals in TDHF. Secondary de-excitation processes, both particle evaporation and fission, are simulated by a statistical compound-nucleus de-excitation model, gemini++. Dynamics of a fast isospin equilibration process followed by a slow drift toward the mass symmetry are commonly observed, as expected. Various reaction outcomes are evaluated, including average mass and charge numbers of reaction products, total kinetic energy loss (TKEL), scattering angle, contact time, and production cross sections for primary and secondary products. By comparing with the experimental data, we find that SMF and TDHF quantitatively reproduce experimental data for few-nucleon-transfer channels around the average values. In contrast, for many-nucleon-transfer channels, we find that the SMF approach provides much better description of the experimentally measured isotopic distributions. The results underline the importance of beyond-mean-field effects, especially one-body (mean-field) fluctuations and correlations, in describing multinucleon transfer processes. Moreover, through a combined analysis of SMF with a statistical model, gemini++, we find a significant contribution of transfer-induced fission, which is consistent with the experimental observation. In some cases, the SMF approach overestimates the isotopic width, requiring further improvements of the theoretical description. Possible ways to improve the description are discussed. The SMF approach is designed to describe the quantum many-body problem according to an ensemble of mean-field trajectories, taking into account part of many-body correlations in the description. As it requires feasible computational costs comparable to the ordinary TDHF approach, together with further model improvements, it will be a promising tool in the search for optimal reaction conditions to produce yet-unknown neutron-rich heavy nuclei through the multinucleon transfer reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Stochastic Schrödinger evolution over piecewise enlarged filtrations

This paper constructs a nonlinear filtering framework that admits appearances of new information processes at random times by introducing piecewise enlargements of filtrations and proposes a new energy-based Schrodinger evolution expressed as a stochastic differential equation on a complex Hilbert space. Each information process is modeled as the sum of a random variable taking the eigenvalues of a Hamiltonian and an independent Brownian bridge noise. It is shown that under a piecewise enlarged filtration, the wave function is a jump-diffusion process until it collapses at some terminal time. In between discontinuities, the dynamics of the state vector are governed by different Wiener processes and diffusion coefficients. This motivates the introduction of an inclusive chain of Kolmogorov probability spaces or a *-isomorphic chain of commutative von Neumann probability spaces, on which the quantum system evolves differently based on the number of active information processes. The expectation of the Hamiltonian at a given state is the solution of a second-order nonlinear differential equation determined by one of the possible regimes that the quantum system belongs to. It is shown that the collapse rate is a submartingale with positive jumps and the Shannon entropy process is a supermartingale with expected negative jumps when passing to higher-order probability spaces. The framework is extended to the case when the Hamiltonian is modeled as a function of a set of commutative operators, where each operator is associated with a different piecewise enlarged filtration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗