-
Experimental and theoretical confirmation of an orthorhombic phase transition in niobium at high pressure and temperature
Authors:
Daniel Errandonea,
Leonid Burakovsky,
Dean L. Preston,
Simon G. MacLeod,
David Santamaria-Perez,
Shaoping Chen,
Hyunchae Cynn,
Sergey I. Simak,
Malcolm I. McMahon,
John E. Proctor,
Mohamed Mezouar
Abstract:
Compared to other body-centered cubic (bcc) transition metals Nb has been the subject of fewer compression studies and there are still aspects of its phase diagram which are unclear. Here, we report a combined theoretical and experimental study of Nb under high pressure and temperature. We present the results of static laser-heated diamond anvil cell experiments up to 120 GPa using synchrotron-bas…
▽ More
Compared to other body-centered cubic (bcc) transition metals Nb has been the subject of fewer compression studies and there are still aspects of its phase diagram which are unclear. Here, we report a combined theoretical and experimental study of Nb under high pressure and temperature. We present the results of static laser-heated diamond anvil cell experiments up to 120 GPa using synchrotron-based fast x-ray diffraction combined with ab initio quantum molecular dynamics simulations. The melting curve of Nb is determined, and evidence for a solid-solid phase transformation in Nb with increasing temperature is found. The high-temperature phase of Nb is orthorhombic Pnma. The bcc-Pnma transition is clearly seen in the experimental data on the Nb principal Hugoniot. The bcc-Pnma coexistence observed in our experiments is explained. Agreement between the measured and calculated melting curves is very good except at 40-60 GPa where three experimental points lie below the theoretical melting curve by 250 K (or 7%); a possible explanation is given.
△ Less
Submitted 24 January, 2024;
originally announced January 2024.
-
On the temperature and density dependence of dislocation drag from phonon wind
Authors:
Daniel N. Blaschke,
Leonid Burakovsky,
Dean L. Preston
Abstract:
At extreme strain rates, where fast moving dislocations govern plastic deformation, anharmonic phonon scattering imparts a drag force on the dislocations. In this paper, we present calculations of the dislocation drag coefficients of aluminum and copper as functions of temperature and density. We discuss the sensitivity of the drag coefficients to changes in the third-order elastic constants with…
▽ More
At extreme strain rates, where fast moving dislocations govern plastic deformation, anharmonic phonon scattering imparts a drag force on the dislocations. In this paper, we present calculations of the dislocation drag coefficients of aluminum and copper as functions of temperature and density. We discuss the sensitivity of the drag coefficients to changes in the third-order elastic constants with temperature and density.
△ Less
Submitted 12 June, 2021; v1 submitted 17 April, 2021;
originally announced April 2021.
-
Analytic model of the remobilization of pinned glide dislocations: including dislocation drag from phonon wind
Authors:
Daniel N. Blaschke,
Abigail Hunter,
Dean L. Preston
Abstract:
In this paper we discuss the effect of a non-constant dislocation drag coefficient on the very high strain rate regime within an analytic model describing mobile-immobile dislocation intersections applicable to polycrystals. Based on previous work on dislocation drag, we estimate its temperature and pressure dependence and its effects on stress-strain rate relations. In the high temperature regime…
▽ More
In this paper we discuss the effect of a non-constant dislocation drag coefficient on the very high strain rate regime within an analytic model describing mobile-immobile dislocation intersections applicable to polycrystals. Based on previous work on dislocation drag, we estimate its temperature and pressure dependence and its effects on stress-strain rate relations. In the high temperature regime, we show that drag can remain the dominating effect even down to intermediate strain rates. We also discuss the consequences of having a limiting dislocation velocity, a feature which is typically predicted by analytic models of dislocation drag, but which is somewhat under debate because a number of MD simulations predict supersonic dislocations.
△ Less
Submitted 2 March, 2020; v1 submitted 18 December, 2019;
originally announced December 2019.
-
Dislocation drag from phonon wind in an isotropic crystal at large velocities
Authors:
Daniel N. Blaschke,
Emil Mottola,
Dean L. Preston
Abstract:
The anharmonic interaction and scattering of phonons by a moving dislocation, the photon wind, imparts a drag force $v B(v, T, ρ)$ on the dislocation. In early studies the drag coefficient $B$ was computed and experimentally determined only for dislocation velocities $v$ much less than transverse sound speed, $c_t$. In this paper we derive analytic expressions for the velocity dependence of $B$ up…
▽ More
The anharmonic interaction and scattering of phonons by a moving dislocation, the photon wind, imparts a drag force $v B(v, T, ρ)$ on the dislocation. In early studies the drag coefficient $B$ was computed and experimentally determined only for dislocation velocities $v$ much less than transverse sound speed, $c_t$. In this paper we derive analytic expressions for the velocity dependence of $B$ up to $c_t$ in terms of the third-order continuum elastic constants of an isotropic crystal, in the continuum Debye approximation, valid for dislocation velocities approaching the sound speed. In so doing we point out that the most general form of the third order elastic potential for such a crystal and the dislocation-phonon interaction requires two additional elastic constants involving asymmetric local rotational strains, which have been neglected previously. We compute the velocity dependence of the transverse phonon wind contribution to $B$ in the range 1%-90% $c_t$ for Al, Cu, Fe, and Nb in the isotropic Debye approximation. The drag coefficient for transverse phonons scattering from screw dislocations is finite as $v \rightarrow c_t$, whereas $B$ is divergent for transverse phonons scattering from edge dislocations in the same limit. This divergence indicates the breakdown of the Debye approximation and sensitivity of the drag coefficient at very high velocities to the microscopic crystalline lattice cutoff. We compare our results to experimental results wherever possible and identify ways to validate and further improve the theory of dislocation drag at high velocities with realistic phonon dispersion relations, inclusion of lattice cutoff effects, MD simulation data, and more accurate experimental measurements.
△ Less
Submitted 1 October, 2019; v1 submitted 28 June, 2019;
originally announced July 2019.
-
Thermoelastic-Plastic Flow Equations in General Coordinates
Authors:
Daniel N. Blaschke,
Dean L. Preston
Abstract:
The equations governing the thermoelastic-plastic flow of isotropic solids in the Prandtl-Reuss and small anisotropy approximations in Cartesian coordinates are generalized to arbitrary coordinate systems. In applications the choice of coordinates is dictated by the symmetry of the solid flow. The generally invariant equations are evaluated in spherical, cylindrical (including uniaxial), and both…
▽ More
The equations governing the thermoelastic-plastic flow of isotropic solids in the Prandtl-Reuss and small anisotropy approximations in Cartesian coordinates are generalized to arbitrary coordinate systems. In applications the choice of coordinates is dictated by the symmetry of the solid flow. The generally invariant equations are evaluated in spherical, cylindrical (including uniaxial), and both prolate and oblate spheroidal coordinates.
△ Less
Submitted 22 March, 2018; v1 submitted 12 September, 2017;
originally announced September 2017.
-
Leading Relativistic Corrections to the Kompaneets Equation
Authors:
Lowell S. Brown,
Dean L. Preston
Abstract:
We calculate the first relativistic corrections to the Kompaneets equation for the evolution of the photon frequency distribution brought about by Compton scattering. The Lorentz invariant Boltzmann equation for electron-photon scattering is first specialized to isotropic electron and photon distributions, the squared scattering amplitude and the energy-momentum conserving delta function are each…
▽ More
We calculate the first relativistic corrections to the Kompaneets equation for the evolution of the photon frequency distribution brought about by Compton scattering. The Lorentz invariant Boltzmann equation for electron-photon scattering is first specialized to isotropic electron and photon distributions, the squared scattering amplitude and the energy-momentum conserving delta function are each expanded to order v^/c^4, averages over the directions of the electron and photon momenta are then carried out, and finally an integration over the photon energy yields our Fokker- Planck equation. The Kompaneets equation, which involves only first- and second-order derivatives with respect to the photon energy, results from the order v^2/c^2 terms, while the first relativistic corrections of order v^4/c^4 introduce third- and fourth-order derivatives. We emphasize that our result holds when neither the electrons nor the photons are in thermal equilibrium; two effective temperatures characterize a general, non-thermal electron distribution. When the electrons are in thermal equilibrium our relativistic Fokker-Planck equation is in complete agreement with the most recent published results, but we both disagree with older work.
△ Less
Submitted 31 January, 2012; v1 submitted 26 January, 2012;
originally announced January 2012.
-
Charged Particle Motion in a Plasma: Electron-Ion Energy Partition
Authors:
Lowell S. Brown,
Dean L. Preston,
Robert L. Singleton Jr
Abstract:
A charged particle traversing a plasma loses its energy to both plasma electrons and ions. We compute the energy partition, the fractions $E_e/E_0$ and $E_\smI/E_0$ of the initial energy $E_0$ of this `impurity particle' that are deposited into the electrons and ions when it has slowed down into an equilibrium distribution that we shall determine. We use a well-defined Fokker-Planck equation for t…
▽ More
A charged particle traversing a plasma loses its energy to both plasma electrons and ions. We compute the energy partition, the fractions $E_e/E_0$ and $E_\smI/E_0$ of the initial energy $E_0$ of this `impurity particle' that are deposited into the electrons and ions when it has slowed down into an equilibrium distribution that we shall determine. We use a well-defined Fokker-Planck equation for the phase space distribution of the charged impurity particles in a weakly to moderately coupled plasma. The Fokker-Planck equation holds to first sub-leading order in the dimensionless plasma coupling constant, which means we compute to order $n\ln n$ (leading) and $n$ (sub-leading) in the plasma density $n$. Previously, the order $n$ terms had been estimated, not calculated. Since the charged particle does not come to rest, the energy loss obtained by an integration of a $dE/dx$ has an ambiguity of order of the plasma temperature. Our Fokker-Planck formulation provides an unambiguous, precise definition of the energy fractions. For equal electron and ion temperatures, we find that our precise results agree well with a fit obtained by Fraley, Linnebur, Mason, and Morse. The case with differing electron and ion temperatures, a case of great importance for nuclear fusion, will be investigated in detail in the present paper. The energy partitions for this general case, partitions that have not been obtained before, will be presented. We find that now the proper solution of the Fokker-Planck equation yields a quasi-static equilibrium distribution to which fast particles relax that has neither the electron nor the ion temperature. This "schizophrenic" final ensemble of slowed particles gives a new mechanism to bring the electron and ion temperatures together. The rate at which this new mechanism brings the electrons and ions in the plasma into thermal equilibrium will be computed.
△ Less
Submitted 19 September, 2011; v1 submitted 7 June, 2011;
originally announced June 2011.
-
Highly Charged Ions in a Dilute Plasma: An Exact Asymptotic Solution Involving Strong Coupling
Authors:
Lowell S. Brown,
David C. Dooling,
Dean L. Preston
Abstract:
The ion sphere model introduced long ago by Salpeter is placed in a rigorous theoretical setting. The leading corrections to this model for very highly charged but dilute ions in thermal equilibrium with a weakly coupled, one-component background plasma are explicitly computed, and the subleading corrections shown to be negligibly small. This is done using effective field theory methods advocate…
▽ More
The ion sphere model introduced long ago by Salpeter is placed in a rigorous theoretical setting. The leading corrections to this model for very highly charged but dilute ions in thermal equilibrium with a weakly coupled, one-component background plasma are explicitly computed, and the subleading corrections shown to be negligibly small. This is done using effective field theory methods advocated by Brown and Yaffe. Thus, corrections to nuclear reaction rates that such highly charged ions may undergo can be computed precisely. Moreover, their contribution to the equation of state can also be computed with precision. Such analytic results for very strong coupling are rarely available, and they can serve as benchmarks for testing computer models in this limit.
△ Less
Submitted 19 April, 2006; v1 submitted 29 March, 2006;
originally announced March 2006.
-
Highly charged ions in a weakly coupled plasma: an exact solution
Authors:
Lowell S. Brown,
David C. Dooling,
Dean L. Preston
Abstract:
The ion sphere model introduced long ago by Salpeter is placed in a rigorous theoretical setting. The leading corrections to this model for very highly charged but dilute ions in thermal equilibrium with a weakly coupled, one-component background plasma are explicitly computed, and the subleading corrections are shown to be negligibly small. Such analytic results for very strong coupling are rar…
▽ More
The ion sphere model introduced long ago by Salpeter is placed in a rigorous theoretical setting. The leading corrections to this model for very highly charged but dilute ions in thermal equilibrium with a weakly coupled, one-component background plasma are explicitly computed, and the subleading corrections are shown to be negligibly small. Such analytic results for very strong coupling are rarely available, and they can serve as benchmarks for testing computer models in this limit.
△ Less
Submitted 22 September, 2005;
originally announced September 2005.
-
Rigorous theory of nuclear fusion rates in a plasma
Authors:
Lowell S. Brown,
David C. Dooling,
Dean L. Preston
Abstract:
Real-time thermal field theory is used to reveal the structure of plasma corrections to nuclear reactions. Previous results are recovered in a fashion that clarifies their nature, and new extensions are made. Brown and Yaffe have introduced the methods of effective quantum field theory into plasma physics. They are used here to treat the interesting limiting case of dilute but very highly charge…
▽ More
Real-time thermal field theory is used to reveal the structure of plasma corrections to nuclear reactions. Previous results are recovered in a fashion that clarifies their nature, and new extensions are made. Brown and Yaffe have introduced the methods of effective quantum field theory into plasma physics. They are used here to treat the interesting limiting case of dilute but very highly charged particles reacting in a dilute, one-component plasma. The highly charged particles are very strongly coupled to this background plasma. The effective field theory proves that this mean field solution plus the one-loop term dominate; higher loop corrections are negligible even though the problem involves strong coupling. Such analytic results for very strong coupling are rarely available, and they can serve as benchmarks for testing computer models.
△ Less
Submitted 22 September, 2005;
originally announced September 2005.
-
dE/dx to Subleading Order in the Plasma Density
Authors:
Lowell S. Brown,
Dean L. Preston,
Robert L. Singleton Jr
Abstract:
Dimensional continuation is employed to compute the energy loss rate for a non-relativistic particle moving through a highly ionized plasma. No restriction is made on the charge, mass, or speed of this particle, but it is assumed that the plasma is not strongly coupled in that the dimensionless plasma coupling parameter g= e^2 κ_D/4πT is small, where κ_D is the Debye wave number. To leading orde…
▽ More
Dimensional continuation is employed to compute the energy loss rate for a non-relativistic particle moving through a highly ionized plasma. No restriction is made on the charge, mass, or speed of this particle, but it is assumed that the plasma is not strongly coupled in that the dimensionless plasma coupling parameter g= e^2 κ_D/4πT is small, where κ_D is the Debye wave number. To leading order in this coupling, dE/dx is of the generic form g^2 \ln[g^2 C]. The prefactor of the logarithm is well known. We compute the constant C under the logarithm exactly. Our result differs from approximations given in the literature, with differences in the range of about 20% for cases relevant to inertial confinement fusion experiments.
△ Less
Submitted 20 September, 2005;
originally announced September 2005.
-
Charged Particle Motion in a Highly Ionized Plasma
Authors:
Lowell S. Brown,
Dean L. Preston,
Robert L. Singleton Jr
Abstract:
A recently introduced method utilizing dimensional continuation is employed to compute the energy loss rate for a non-relativistic particle moving through a highly ionized plasma. No restriction is made on the charge, mass, or speed of this particle. It is, however, assumed that the plasma is not strongly coupled in the sense that the dimensionless plasma coupling parameter g=e^2κ_D/ 4πT is smal…
▽ More
A recently introduced method utilizing dimensional continuation is employed to compute the energy loss rate for a non-relativistic particle moving through a highly ionized plasma. No restriction is made on the charge, mass, or speed of this particle. It is, however, assumed that the plasma is not strongly coupled in the sense that the dimensionless plasma coupling parameter g=e^2κ_D/ 4πT is small, where κ_D is the Debye wave number of the plasma. To leading and next-to-leading order in this coupling, dE/dx is of the generic form g^2 \ln[C g^2]. The precise numerical coefficient out in front of the logarithm is well known. We compute the constant C under the logarithm exactly for arbitrary particle speeds. Our exact results differ from approximations given in the literature. The differences are in the range of 20% for cases relevant to inertial confinement fusion experiments. The same method is also employed to compute the rate of momentum loss for a projectile moving in a plasma, and the rate at which two plasmas at different temperatures come into thermal equilibrium. Again these calculations are done precisely to the order given above. The loss rates of energy and momentum uniquely define a Fokker-Planck equation that describes particle motion in the plasma. The coefficients determined in this way are thus well-defined, contain no arbitrary parameters or cutoffs, and are accurate to the order described. This Fokker-Planck equation describes the longitudinal straggling and the transverse diffusion of a beam of particles. It should be emphasized that our work does not involve a model, but rather it is a precisely defined evaluation of the leading terms in a well-defined perturbation theory.
△ Less
Submitted 20 March, 2007; v1 submitted 17 January, 2005;
originally announced January 2005.
-
An analytic model of the shear modulus at all densities and temperatures
Authors:
L. Burakovsky,
C. W. Greeff,
D. L. Preston
Abstract:
An analytic model of the shear modulus applicable at temperatures up to melt and at all densities is presented. It is based in part on a relation between the melting temperature and the shear modulus at melt. Experimental data on argon are shown to agree with this relation to within 1%. The model of the shear modulus involves seven parameters, all of which can be determined from zero-pressure ex…
▽ More
An analytic model of the shear modulus applicable at temperatures up to melt and at all densities is presented. It is based in part on a relation between the melting temperature and the shear modulus at melt. Experimental data on argon are shown to agree with this relation to within 1%. The model of the shear modulus involves seven parameters, all of which can be determined from zero-pressure experimental data. We obtain the values of these parameters for 11 elemental solids. Both the experimental data on the room-temperature shear modulus of argon to compressions of \sim 2.5, and theoretical calculations of the zero-temperature shear modulus of aluminum to compressions of \sim 3.5 are in good agreement with the model. Electronic structure calculations of the shear moduli of copper and gold to compressions of 2, performed by us, agree with the model to within uncertainties.
△ Less
Submitted 7 November, 2002; v1 submitted 29 August, 2002;
originally announced August 2002.
-
An analytic model of the Gruneisen parameter at all densities
Authors:
L. Burakovsky,
D. L. Preston
Abstract:
We model the density dependence of the Gruneisen parameter as gamma(rho) = 1/2 + gamma_1/rho^{1/3} + gamma_2/rho^{q}, where gamma_1, gamma_2, and q>1 are constants. This form is based on the assumption that gamma is an analytic function of V^{1/3}, and was designed to accurately represent the experimentally determined low-pressure behavior of gamma. The numerical values of the constants are obta…
▽ More
We model the density dependence of the Gruneisen parameter as gamma(rho) = 1/2 + gamma_1/rho^{1/3} + gamma_2/rho^{q}, where gamma_1, gamma_2, and q>1 are constants. This form is based on the assumption that gamma is an analytic function of V^{1/3}, and was designed to accurately represent the experimentally determined low-pressure behavior of gamma. The numerical values of the constants are obtained for 20 elemental solids. Using the Lindemann criterion with our model for gamma, we calculate the melting curves for Al, Ar, Ni, Pd, and Pt and compare them to available experimental melt data. We also determine the Z (atomic number) dependence of gamma_1. The high-compression limit of the model is shown to follow from a generalization of the Slater, Dugdale-MacDonald, and Vashchenko-Zubarev forms for the dependence of the Gruneisen parameter.
△ Less
Submitted 4 November, 2002; v1 submitted 11 June, 2002;
originally announced June 2002.
-
Dislocation-Mediated Melting: The One-Component Plasma Limit
Authors:
Leonid Burakovsky,
Dean L. Preston
Abstract:
The melting parameter $Γ_m$ of a classical one-component plasma is estimated using a relation between melting temperature, density, shear modulus, and crystal coordination number that follows from our model of dislocation-mediated melting. We obtain $Γ_m=172\pm 35,$ in good agreement with the results of numerous Monte-Carlo calculations.
The melting parameter $Γ_m$ of a classical one-component plasma is estimated using a relation between melting temperature, density, shear modulus, and crystal coordination number that follows from our model of dislocation-mediated melting. We obtain $Γ_m=172\pm 35,$ in good agreement with the results of numerous Monte-Carlo calculations.
△ Less
Submitted 19 September, 2000;
originally announced September 2000.
-
Analysis of Dislocation Mechanism for Melting of Elements: Pressure Dependence
Authors:
Leonid Burakovsky,
Dean L. Preston,
Richard R. Silbar
Abstract:
In the framework of melting as a dislocation-mediated phase transition we derive an equation for the pressure dependence of the melting temperatures of the elements valid up to pressures of order their ambient bulk moduli. Melting curves are calculated for Al, Mg, Ni, Pb, the iron group (Fe, Ru, Os), the chromium group (Cr, Mo, W), the copper group (Cu, Ag, Au), noble gases (Ne, Ar, Kr, Xe, Rn),…
▽ More
In the framework of melting as a dislocation-mediated phase transition we derive an equation for the pressure dependence of the melting temperatures of the elements valid up to pressures of order their ambient bulk moduli. Melting curves are calculated for Al, Mg, Ni, Pb, the iron group (Fe, Ru, Os), the chromium group (Cr, Mo, W), the copper group (Cu, Ag, Au), noble gases (Ne, Ar, Kr, Xe, Rn), and six actinides (Am, Cm, Np, Pa, Th, U). These calculated melting curves are in good agreement with existing data. We also discuss the apparent equivalence of our melting relation and the Lindemann criterion, and the lack of the rigorous proof of their equivalence. We show that the would-be mathematical equivalence of both formulas must manifest itself in a new relation between the Grüneisen constant, bulk and shear moduli, and the pressure derivative of the shear modulus.
△ Less
Submitted 5 May, 2000;
originally announced May 2000.
-
On the Bragg, Leibfried, and Modified Leibfried Numbers
Authors:
Leonid Burakovsky,
Dean L. Preston,
Richard R. Silbar
Abstract:
The Bragg, Leibfried, and modified Leibfried numbers are defined in the context of a theory of dislocation-mediated melting, and their values are determined from the properties of the dislocation ensemble at the melting temperature. The approximate numerical coincidence of the Bragg and modified Leibfried numbers is explained. The parameter K in the definition of the modified Leibfried number is…
▽ More
The Bragg, Leibfried, and modified Leibfried numbers are defined in the context of a theory of dislocation-mediated melting, and their values are determined from the properties of the dislocation ensemble at the melting temperature. The approximate numerical coincidence of the Bragg and modified Leibfried numbers is explained. The parameter K in the definition of the modified Leibfried number is shown to be the natural logarithm of the effective coordination number. Our analysis reveals that the Bragg number can be considered an elemental constant, in contrast to the Leibfried and modified Leibfried numbers.
△ Less
Submitted 4 April, 2000;
originally announced April 2000.
-
Melting as a String-Mediated Phase Transition
Authors:
Leonid Burakovsky,
Dean L. Preston,
Richard R. Silbar
Abstract:
We present a theory of the melting of elemental solids as a dislocation-mediated phase transition. We model dislocations near melt as non-interacting closed strings on a lattice. In this framework we derive simple expressions for the melting temperature and latent heat of fusion that depend on the dislocation density at melt. We use experimental data for more than half the elements in the Period…
▽ More
We present a theory of the melting of elemental solids as a dislocation-mediated phase transition. We model dislocations near melt as non-interacting closed strings on a lattice. In this framework we derive simple expressions for the melting temperature and latent heat of fusion that depend on the dislocation density at melt. We use experimental data for more than half the elements in the Periodic Table to determine the dislocation density from both relations. Melting temperatures yield a dislocation density of (0.61\pm 0.20) b^{-2}, in good agreement with the density obtained from latent heats, (0.66\pm 0.11) b^{-2}, where b is the length of the smallest perfect-dislocation Burgers vector. Melting corresponds to the situation where, on average, half of the atoms are within a dislocation core.
△ Less
Submitted 3 April, 2000;
originally announced April 2000.
-
Analysis of Dislocation Mechanism for Melting of Elements
Authors:
Leonid Burakovsky,
Dean L. Preston
Abstract:
The melting of elemental solids is modelled as a dislocation-mediated transition on a lattice. Statistical mechanics of linear defects is used to obtain a new relation between melting temperature, crystal structure, atomic volume, and shear modulus that is accurate to 17% for at least half of the Periodic Table.
The melting of elemental solids is modelled as a dislocation-mediated transition on a lattice. Statistical mechanics of linear defects is used to obtain a new relation between melting temperature, crystal structure, atomic volume, and shear modulus that is accurate to 17% for at least half of the Periodic Table.
△ Less
Submitted 30 March, 2000;
originally announced March 2000.
