A splitstep method to include electron–electron collisions via Monte Carlo in multiple rate equation simulations
Abstract
A splitstep numerical method for calculating ultrafast freeelectron dynamics in dielectrics is introduced. The two split steps, independently programmed in C++11 and FORTRAN 2003, are interfaced via the presented open source wrapper. The first step solves a deterministic extended multirate equation for the ionization, electron–phonon collisions, and single photon absorption by freecarriers. The second step is stochastic and models electron–electron collisions using MonteCarlo techniques. This combination of deterministic and stochastic approaches is a unique and efficient method of calculating the nonlinear dynamics of 3D materials exposed to high intensity ultrashort pulses. Results from simulations solving the proposed model demonstrate how electron–electron scattering relaxes the nonequilibrium electron distribution on the femtosecond time scale.
 Authors:
 Department of Physics and OPTIMAS Research Center, University of Kaiserslautern (Germany)
 Department of Chemistry and OPTIMAS Research Center, University of Kaiserslautern (Germany)
 Department of Physics, Kennesaw State University, Kennesaw, GA 30144 (United States)
 Publication Date:
 OSTI Identifier:
 22572363
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Computational Physics; Journal Volume: 322; Other Information: Copyright (c) 2016 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; DIELECTRIC MATERIALS; ELECTRONELECTRON COLLISIONS; ELECTRONELECTRON INTERACTIONS; FORTRAN; IONIZATION; MONTE CARLO METHOD; NONLINEAR PROBLEMS; PHONONS; PHOTONS; PULSES; SIMULATION; STOCHASTIC PROCESSES
Citation Formats
Huthmacher, Klaus, Molberg, Andreas K., Rethfeld, Bärbel, and Gulley, Jeremy R., Email: jgulley@kennesaw.edu. A splitstep method to include electron–electron collisions via Monte Carlo in multiple rate equation simulations. United States: N. p., 2016.
Web. doi:10.1016/J.JCP.2016.06.043.
Huthmacher, Klaus, Molberg, Andreas K., Rethfeld, Bärbel, & Gulley, Jeremy R., Email: jgulley@kennesaw.edu. A splitstep method to include electron–electron collisions via Monte Carlo in multiple rate equation simulations. United States. doi:10.1016/J.JCP.2016.06.043.
Huthmacher, Klaus, Molberg, Andreas K., Rethfeld, Bärbel, and Gulley, Jeremy R., Email: jgulley@kennesaw.edu. 2016.
"A splitstep method to include electron–electron collisions via Monte Carlo in multiple rate equation simulations". United States.
doi:10.1016/J.JCP.2016.06.043.
@article{osti_22572363,
title = {A splitstep method to include electron–electron collisions via Monte Carlo in multiple rate equation simulations},
author = {Huthmacher, Klaus and Molberg, Andreas K. and Rethfeld, Bärbel and Gulley, Jeremy R., Email: jgulley@kennesaw.edu},
abstractNote = {A splitstep numerical method for calculating ultrafast freeelectron dynamics in dielectrics is introduced. The two split steps, independently programmed in C++11 and FORTRAN 2003, are interfaced via the presented open source wrapper. The first step solves a deterministic extended multirate equation for the ionization, electron–phonon collisions, and single photon absorption by freecarriers. The second step is stochastic and models electron–electron collisions using MonteCarlo techniques. This combination of deterministic and stochastic approaches is a unique and efficient method of calculating the nonlinear dynamics of 3D materials exposed to high intensity ultrashort pulses. Results from simulations solving the proposed model demonstrate how electron–electron scattering relaxes the nonequilibrium electron distribution on the femtosecond time scale.},
doi = {10.1016/J.JCP.2016.06.043},
journal = {Journal of Computational Physics},
number = ,
volume = 322,
place = {United States},
year = 2016,
month =
}

An extension to the multiplehistogram method (sometimes referred to as the FerrenbergSwendsen method) for use in quantum Monte Carlo simulations is presented. This method is shown to work well for the twodimensional repulsive Hubbard model, allowing measurements to be taken over a continuous region of parameters. The method also reduces the error bars over the range of parameter values due to the overlapping of multiple histograms. A continuous sweep of parameters and reduced error bars allow one to make more difficult measurements, such as Maxwell constructions used to study phase separation. Possibilities also exist for this method to be usedmore »

WEAB20901: A Monte CarloBased Method to Include Random Errors in Robust Optimization
Purpose: To develop an efficient method to implement random setup errors and organ motion in robust optimization for proton therapy treatment planning. Methods: The plans were created with an inhouse robust optimizer, coupled with a superfast Monte Carlo (MC) engine (less than 1 minute for final dose). MC simulates random errors by shifting randomly the starting point of each particle, according to their probability distribution. Such strategy assumes a sufficient number of treatment fractions. Two strategies are presented: 1) Full robust optimization with beamlets that already include the effect of random errors and 2) Mixed robust optimization, where the nominalmore » 
Accuracy of the classical trajectory Monte Carlo method for electron capture in Li{sup 3+} and Ne{sup 10+}+H(1s) collisions
The accuracy of classical trajectory Monte Carlo treatments of electron capture is studied, by focusing on collisions on H(1s) targets by Li{sup 3+} and Ne{sup 10+} projectiles, treated in a separate paper. We examine how the choice of the initial distribution, and the partition of phase space usually employed to calculate partial cross sections, influence the accuracy of the method. With respect to the former, an improvement over the singlemicrocanonical choice is advisable, but further refinements based on the electron density are not worthwhile. Regarding the latter, we illustrate the accuracy of the 'binning' method for n>2. We show thatmore »