PubMed Health⌕ Search

Biomedical subjects

B Militzer

Publications and source records attributed to B Militzer.

6 recordsLinked to original sources

First principles calculations of shock compressed fluid helium.

The properties of hot dense helium at megabar pressures are studied with two first principles computer simulation techniques: path integral Monte Carlo simulation and density functional molecular dynamics. The simulations predict that the compressibility of helium is substantially increased by electronic excitations that are present in the hot fluid at thermodynamic equilibrium. A maximum compression ratio of 5.24(4)-fold the initial density was predicted for 360 GPa and 150,000 K. This result distinguishes helium from deuterium, for which simulations predicted a maximum compression ratio of 4.3(1). Hugoniot curves for statically precompressed samples are also discussed.

Journal Article↗

Dense plasma effects on nuclear reaction rates.

Plasma density effects can cause an exponential change in charged particle nuclear reaction rates important in stellar evolution. Reaction rates in dense plasma, with emphasis on quantum aspects, are examined here using path integral Monte Carlo calculations. Quantum mechanics causes a reduction in the many body enhancement of the reaction rate compared to the value for a classical system. This can be attributed to the "quantum smearing" of the short range Coulomb interaction resulting in reduced repulsion between the reacting pair and surrounding particles. Electron screening and ion exchange effects are also examined, with screening reducing and exchange slightly increasing the many body enhancement.

Journal Article↗

Calculation of a deuterium double shock Hugoniot from ab initio simulations.

We calculate the equation of state of dense deuterium with two ab initio simulation techniques, path integral Monte Carlo and density functional theory molecular dynamics, in the density range of 0.67 < or = rho < or = 1.60 g cm(-3). We derive the double shock Hugoniot and compare with the recent laser-driven double shock wave experiments by Mostovych et al. [Phys. Rev. Lett. 85, 3870 (2000)]. We find excellent agreement between the two types of microscopic simulations, but a significant discrepancy with the laser-driven shock measurements.

Journal Article↗

Path integral Monte Carlo simulation of the low-density hydrogen plasma.

Restricted path integral Monte Carlo simulations are used to calculate the equilibrium properties of hydrogen in the density and temperature range of 9.83 x 10(-4)</=rho</=0.153 g cm(-3) and 5000</=T</=250 000 K. We test the accuracy of the pair density matrix and analyze the dependence on the system size, on the time step of the path integral, and on the type of nodal surface. We calculate the equation of state and compare with other models for hydrogen valid in this regime. Further, we characterize the state of hydrogen and describe the changes from a plasma to an atomic and molecular liquid by analyzing the pair correlation functions and estimating the number of atoms and molecules present.

Journal Article↗

Path integral monte carlo calculation of the deuterium hugoniot

Restricted path integral Monte Carlo simulations have been used to calculate the equilibrium properties of deuterium for two densities: 0.674 and 0.838 g cm(-3) ( r(s) = 2.00 and 1.86) in the temperature range of 10(5)</=T</=10(6) K. We carefully assess size effects and dependence on the time step of the path integral. Further, we compare the results obtained with a free particle nodal restriction with those from a self-consistent variational principle, which includes interactions and bound states. By using the calculated internal energies and pressures, we determine the shock Hugoniot and compare with recent laser shock wave experiments as well as other theories.

Journal Article↗

Variational density matrix method for warm, condensed matter: application to dense hydrogen

A variational principle for optimizing thermal density matrices is introduced. As a first application, the variational many-body density matrix is written as a determinant of one-body density matrices, which are approximated by Gaussians with the mean, width, and amplitude as variational parameters. The method is illustrated for the particle in an external field problem, the hydrogen molecule and dense hydrogen where the molecular, the dissociated, and the plasma regime are described. Structural and thermodynamic properties (energy, equation of state, and shock Hugoniot) are presented.

Journal Article↗