Search

We present the results of the first Charged-Particle Transport Coefficient Code Comparison Workshop, which was held in Albuquerque, NM October 4-6, 2016. In this first workshop, scientists from eight institutions and four countries gathered to compare calculations of transport coefficients including thermal and electrical conduction, electron-ion coupling, inter-ion diffusion, ion viscosity, and charged particle stopping powers. Here, we give general background on Coulomb coupling and computational expense, review where some transport coefficients appear in hydrodynamic equations, and present the submitted data. Large variations are found when either the relevant Coulomb coupling parameter is large or computational expense causes difficulties. Understanding the general accuracy and uncertainty associated with such transport coefficients is important for quantifying errors in hydrodynamic simulations of inertial confinement fusion and high-energy density experiments., Comment: 45 pages, 17 figures

Using the Klein--Majda--Damodaran model of nearly parallel vortex filaments, we construct vortex knots and links on a torus involving periodic boundary conditions and analyze their stability. For a...

The high cost of density functional theory has hitherto limited the ab initio prediction of equation of state (EOS). In this article, we employ a combination of large scale computing, advanced simulation techniques, and smart data science strategies to provide an unprecedented, ab initio performance analysis of the high explosive pentaerythritol tetranitrate (PETN). Comparison to both experiment and thermochemical predictions reveals important quantitative limitations of DFT for EOS prediction, and thus the assessment of high explosives. In particular, we find DFT predicts the energy of PETN detonation products to be systematically too high relative to the unreacted neat crystalline material, resulting in an underprediction of the detonation velocity, pressure, and temperature at the Chapman-Jouguet (CJ) state. The energetic bias can be partially accounted for by high-level electronic structure calculations of the product molecules. We also demonstrate a modeling strategy for mapping chemical composition across a wide parameter space with limited numerical data, the results of which suggest additional molecular species to consider in thermochemical modeling., 16 pages, 8 figures

Here we outline the methods used in LANL's thermochemical equation of state modeling code. We lay out the basic numerical modeling used to describe equilibrium molecular fluid equation of state. We look at an example of PBX-9501 and find good agreement between experiments and other theoretical modeling. This work serves as the basis from which we will further develop modeling of reactive matter equation of state for materials such as high explosives.

A calibration of the AWSD (Arrhenius Wescott-Stewart-Davis) reactive flow model, applied to the HMX based explosive PBX 9501 is presented. Requisite equation of state parameters for the high explosive reactants and products are based on a combination of experimental data and theoretical calculations. The experimental data needed to calibrate the reaction rate model, specifically shock to detonation transition and detonation propagation at a variety of initial temperatures, are presented along with model results and planned future activities.

Here we present reactive Monte Carlo benchmarks of LANL’s new thermochemical code (Magpie). Thermochemical codes make various approximations such as representing the non-ideal portions of the equation of state with perturbation or integral equation theories. Further, such codes implement various mixing rules, and solve the chemical equilibrium problem by minimizing the free energy. The Monte Carlo (MC) calculations presented here forego the aforementioned approximations and provide an exact benchmark result for thermochemical EOSs (within certain shared approximations). We present results of Monte Carlo calculations of the CJ state, the overdriven Hugoniot, and the release isentrope of PETN, and find that the errors due to perturbation theory and the ideal-mixing approximation implemented in Magpie are quite acceptable for PETN.

Complex plasma mixtures with three or more components are often encountered in astrophysics or in inertial confinement fusion (ICF) experiments. For mixtures containing species with large differences in atomic number Z, the modeling needs to consider at the same time the kinetic theory for low-Z elements combined with the theory of strongly coupled plasma for high-Z elements, as well as all the intermediate situations that can appear in multi-component systems. For such cases, we study the pair distribution functions, self-diffusions, mutual diffusion and viscosity for ternary mixtures at extreme conditions. These quantities can be produced from first principles using orbital free molecular dynamics at the computational expense of very intensive simulations to reach good statistics. Utilizing the first-principles results as reference data, we assess the merit of a global analytic model for transport coefficients, Pseudo-Ions in Jellium (PIJ), based on an iso-electronic assumption ( iso-n e ). With a multi-component hypernetted-chain integral equation, we verify the quality of the iso-n e prescription for describing the static structure of the mixtures. This semi-analytical modeling compares well with the simulation results and allows one to consider plasma mixtures not accessible to simulations. Applications are given for the mix of materials in ICF experiments. A reduction of a multicomponent mixture to an effective binary mixture is also established in the hydrodynamic limit and compared with PIJ estimations for ICF relevant mixtures., Comment: 17 pages, 10 Figures

We present the formulation, simulations, and results for multicomponent mutual diffusion coefficients in the warm, dense matter regime. While binary mixtures have received considerable attention for mass transport, far fewer studies have addressed ternary and more complex systems. We therefore explicitly examine ternary systems utilizing the Maxwell-Stefan formulation that relates diffusion to gradients in the chemical potential. Onsager coefficients then connect the macroscopic diffusion to microscopic particle motions, evinced in trajectories characterized by positions and velocities, through various autocorrelation functions (ACFs). These trajectories are generated by molecular dynamics (MD) simulations either through the Born-Oppenheimer approximation, which treats the ions classically and the electrons quantum-mechanically by an orbital-free density-functional theory, or through a classical MD approach with Yukawa pair-potentials, whose effective ionizations and electron screening length derive from quantal considerations. We employ the reference-mean form of the ACFs and determine the center-of-mass coefficients through a simple reference-frame-dependent similarity transformation. The Onsager terms in turn determine the mutual diffusion coefficients. We examine a representative sample of ternary mixtures as a function of density and temperature from those with only light elements (D-Li-C, D-Li-Al) to those with highly asymmetric mass components (D-Li-Cu, D-Li-Ag, H-C-Ag). We also follow trends in the diffusion as a function of number concentration and evaluated the efficacy of various approximations such as the Darken approximation.

A novel set of experiments and reactive flow modeling of pentaerythritol tetranitrate (PETN) is presented. Here, the specific phenomenon of shock to detonation transition is examined, where an initial, relatively weak shock is propagated into pressed PETN powder at 1.65 g/cm 3 and the subsequent buildup to detonation is observed experimentally. These experiments, in conjunction with reactant and products’ equations of state, are utilized for building reactive flow models.

A new density functional theory, average atom based model for the electrical conductivity of dense plasmas with a mixture of ion species, containing no adjustable parameters, is presented. The model takes the temperature, mass density and relative abundances of the species as input. It takes into account partial ionization, ionic structure, and core-valence orthogonality, and uses quantum mechanical calculations of cross sections. Comparison to an existing high fidelity but computationally expensive method reveals good agreement. The new model is computationally efficient and can reach high temperatures. A new mixing rule is also presented that gives reasonably accurate conductivities for high temperature plasma mixtures.

The long-time existence of topologically nontrivial configurations of quantum vortices in the form of torus knots and links in trapped Bose-Einstein condensates is demonstrated numerically within the three-dimensional Gross-Pitaevskii equation with an external anisotropic parabolic potential. We identify parametric domains of trap anisotropy, characterized by the axial over planar frequency ratio $\ensuremath{\lambda}\ensuremath{\approx}1.5--1.6$, where the lifetime of such quasistationary rotating vortex structures is many hundreds of typical rotation times. This suggests the potential experimental observability of the structures. We quantify the relevant lifetimes as a function of the model parameters (e.g., $\ensuremath{\lambda}$) and initial condition parameters of the knot profile.

In the present work, motivated by numerous recent experimental developments we revisit the dynamics of a single vortex ring in anisotropic harmonic traps. At the theoretical level, we start from a general Lagrangian dynamically capturing the evolution of a vortex ring and not only consider its spectrum of linearized excitations, but also explore the full nonlinear dynamical evolution of the ring as a vortical filament. The theory predicts that the ring is stable for $1 \leq \lambda \leq 2$, where $\lambda=\omega_z/\omega_r$ is the ratio of the trapping frequencies along the $z$ and $r$ axes, i.e., for spherical to slightly oblate condensates. We compare this prediction with direct numerical simulations of the full 3D Gross-Pitaevskii equation (GPE) capturing the linearization spectrum of the ring for different values of the chemical potential and as a function of the anisotropy parameter $\lambda$. We identify this result as being only asymptotically valid as the chemical potential $\mu \rightarrow \infty$, revealing how the stability interval narrows and, in particular, its upper bound decreases for finite $\mu$. Finally, we compare at the dynamical level the results of the GPE with the ones effectively capturing the ring dynamics, revealing the unstable evolution for different values of $\lambda$, as well as the good agreement between the two., Comment: Corrected citation, 10 pages and many fun figures

Topological materials are of interest to both fundamental science and advanced technologies, because topological states are robust with respect to perturbations and dissipation. Experimental detection of topological invariants is thus in great demand, but it remains extremely challenging. Ultrafast laser-matter interactions, and in particular high-harmonic generation (HHG), meanwhile, were proposed several years ago as tools to explore the structural and dynamical properties of various matter targets. Here we show that the high-harmonic emission signal produced by a circularly-polarized laser contains signatures of the topological phase transition in the paradigmatic Haldane model. In addition to clear shifts of the overall emissivity and harmonic cutoff, the high-harmonic emission shows a unique circular dichroism, which exhibits clear changes in behavior at the topological phase boundary. Our findings pave the way to understand fundamental questions about the ultrafast electron-hole pair dynamics in topological materials via non-linear high-harmonic generation spectroscopy., 9 pages including references, then 16 pages of Supplementary Material

For quantum fluids, the role of quantum fluctuations may be significant in several regimes such as when the dimensionality is low, the density is high, the interactions are strong, or for low particle numbers. In this paper we propose a fundamentally different regime for enhanced quantum fluctuations without being restricted by any of the above conditions. Instead, our scheme relies on the engineering of an effective attractive interaction in a dilute, two-component Bose-Einstein condensate (BEC) consisting of thousands of atoms. In such a regime, the quantum spin fluctuations are significantly enhanced (atom bunching with respect to the noninteracting limit) since they act to reduce the interaction energy - a remarkable property given that spin fluctuations are normally suppressed (anti-bunching) at zero temperature. In contrast to the case of true attractive interactions, our approach is not vulnerable to BEC collapse. We numerically demonstrate that these quantum fluctuations are experimentally accessible by either spin or single-component Bragg spectroscopy, offering a useful platform on which to test beyond-mean-field theories. We also develop a variational model and use it to analytically predict the shift of the immiscibility critical point, finding good agreement with our numerics., 12 pages (main body ~ 7 pages), 6 figures

The validity of the widely used linear mixing approximation for the equations of state (EOS) of planetary ices is investigated at pressure-temperature conditions typical for the interior of Uranus and Neptune. The basis of this study are ab initio data ranging up to 1000 GPa and 20 000 K calculated via density functional theory molecular dynamics simulations. In particular, we calculate a new EOS for methane and EOS data for the 1:1 binary mixtures of methane, ammonia, and water, as well as their 2:1:4 ternary mixture. Additionally, the self-diffusion coefficients in the ternary mixture are calculated along three different Uranus interior profiles and compared to the values of the pure compounds. We find that deviations of the linear mixing approximation from the results of the real mixture are generally small; for the thermal EOS they amount to 4% or less. The diffusion coefficients in the mixture agree with those of the pure compounds within 20% or better. Finally, a new adiabatic model of Uranus with an inner layer of almost pure ices is developed. The model is consistent with the gravity field data and results in a rather cold interior ($\mathrm{T_{core}} \mathtt{\sim}$ 4000 K).

Transport properties of mixtures of elements in the dense plasma regime play an important role in natural astrophysical and experimental systems, e.g., inertial confinement fusion. We present a series of orbital-free molecular dynamics simulations on dense plasma mixtures with comparison to a global pseudo ion in jellium model. Hydrogen is mixed with elements of increasingly high atomic number (lithium, carbon, aluminum, copper, and silver) at a fixed temperature of 100 eV and constant pressure set by pure hydrogen at 2g/cm^{3}, namely, 370 Mbars. We compute ionic transport coefficients, such as self-diffusion, mutual diffusion, and viscosity for various concentrations. Small concentrations of the heavy atoms significantly change the density of the plasma and decrease the transport coefficients. The structure of the mixture evidences a strong Coulomb coupling between heavy ions and the appearance of a broad correlation peak at short distances between hydrogen atoms. The concept of an effective one component plasma is used to quantify the overcorrelation of the light element induced by the admixture of a heavy element.

When heated at constant volume (isochoric heating) a hot and dense plasma (10-5000 eV, 1-50 g/cm3) exhibits the same persistent microscopic structure over a wide range of temperatures as intuited long time ago by Laughlin [1]. In this steady-state regime, which depends on the chosen density and on the atomic number, the static structure is essentially independent of the temperature and results in the subtle balance between ionization and temperature leading to a constant coupling between ions. This behavior, suggested by simulations, is confirmed by an analysis in the framework of the Thomas-Fermi scaling laws and is driven by the ionization dynamics which regulates the coupling between ions and electrons. A simple fit is derived allowing for predicting the occurrence of this self-organized regime: the Γ-plateau. (© 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

We construct new equations of state (EOS) for high density and ultrahigh molecular weight polyethylene and their chemical decomposition products under shock loading. The former were built using the SESAME framework, based in part on new specific heat and thermal expansion data reported here. The products EOS was based on thermochemical modeling under the assumption of full thermodynamic and chemical equilibrium. The products are represented as the ideal mixture of bulk carbon in the form of diamond, H2, H, and CH4. In the process of building a new EOS for the products, we recalibrated our exponential-6 pair potential for methane in order to better agree with data that have appeared since its original parameterization. The polyethylene EOS were calibrated to thermal, thermomechanical, and shock data, and their performance was evaluated in hydrodynamic modeling of deep release experiments reported previously.

In the present work, we explore the existence, stability and dynamics of single and multiple vortex ring states that can arise in Bose-Einstein condensates. Earlier works have illustrated the bifurcation of such states, in the vicinity of the linear limit, for isotropic or anisotropic three-dimensional harmonic traps. Here, we extend these states to the regime of large chemical potentials, the so-called Thomas-Fermi limit, and explore their properties such as equilibrium radii and inter-ring distance, for multi-ring states, as well as their vibrational spectra and possible instabilities. In this limit, both the existence and stability characteristics can be partially traced to a particle picture that considers the rings as individual particles oscillating within the trap and interacting pairwise with one another. Finally, we examine some representative instability scenarios of the multi-ring dynamics including breakup and reconnections, as well as the transient formation of vortex lines., 10 pages, 8 figures

We study how concentration changes ionic transport properties along isobars-isotherms for a mixture of hydrogen and silver, representative of turbulent layers relevant to inertial confinement fusion and astrophysics. Hydrogen will typically be fully ionized while silver will be only partially ionized but can have a large effective charge. This will lead to very different physical conditions for the H and Ag. Large first principles orbital free molecular dynamics simulations are performed and the resulting transport properties are analyzed. Comparisons are made with transport theory in the kinetic regime and in the coupled regime. The addition of a small amount of heavy element in a light material has a dramatic effect on viscosity and diffusion of the mixture. This effect is explained through kinetic theory as a manifestation of a crossover between classical diffusion and Lorentz diffusion.

We present an alternative method for interpreting the velocity autocorrelation function (VACF) of a fluid with application to extracting the entropy in a manner similar to the methods developed by Lin et al. [J. Chem. Phys. 119, 11792 (2003)]JCPSA60021-960610.1063/1.1624057 and improved upon by Desjarlais [Phys. Rev. E 88, 062145 (2013)]PLEEE81539-375510.1103/PhysRevE.88.062145. The liquid VACF is decomposed into two components, one gas and one solid, and each contribution's entropic portion is calculated. However, we fit both the gas and solid portions of the VACF in the time domain. This approach is applied to a single-component liquid (a two-phase model of liquid Al at the melt line) and two different two-component systems: a superionic-to-superionic (bcc to fcc) phase transition in H_{2}O at high temperatures and pressures and a metastable liquid state of MgO. For all three examples, comparisons to existing results in the literature demonstrate the validity of our alternative.

We performed a series of plate impact experiments on NH3 gas initially at room temperature and at a pressure of ∼100 psi. Shocked states were determined by optical velocimetry and the temperatures by optical pyrometry, yielding compression ratios of ∼5-10 and second shock temperatures in excess of 7500 K. A first-principles statistical mechanical (thermochemical) approach that included chemical dissociation yielded reasonable agreement with experimental results on the principal Hugoniot, even with interparticle interactions neglected. Theoretical analysis of reshocked states, which predicts a significant degree of chemical dissociation, showed reasonable agreement with experimental data for higher temperature shots; however, reshock calculations required the use of interaction potentials. We rationalize the very different shock temperatures obtained, relative to previous results for argon, in terms of atomic versus molecular heat capacities.

We investigate the plasma enhancement of nuclear reactions in the intermediate coupling regime using orbital free molecular dynamics (OFMD) simulations. Mixtures of H-Cu and H-Ag serve as prototypes of simultaneous weak and strong couplings due to the charge asymmetry. Of particular importance is the partial ionization of Cu and Ag and the free electron polarization captured by OFMD simulations. By comparing a series of OFMD simulations at various concentrations and constant pressure to multi-component hyper-netted chain (MCHNC) calculations of effective binary ionic mixtures, we set a general procedure for computing enhancement factors. The MCHNC procedure allows extension to very low concentrations (5% or less) and to very high temperatures (few keV) unreachable by the simulations. Enhancement factors for nuclear reaction rates extracted from the MCHNC approach are compared with the Salpeter theory in the weak and strong coupling regimes, and a new interpolation is proposed.We investigate the plasma enhancement of nuclear reactions in the intermediate coupling regime using orbital free molecular dynamics (OFMD) simulations. Mixtures of H-Cu and H-Ag serve as prototypes of simultaneous weak and strong couplings due to the charge asymmetry. Of particular importance is the partial ionization of Cu and Ag and the free electron polarization captured by OFMD simulations. By comparing a series of OFMD simulations at various concentrations and constant pressure to multi-component hyper-netted chain (MCHNC) calculations of effective binary ionic mixtures, we set a general procedure for computing enhancement factors. The MCHNC procedure allows extension to very low concentrations (5% or less) and to very high temperatures (few keV) unreachable by the simulations. Enhancement factors for nuclear reaction rates extracted from the MCHNC approach are compared with the Salpeter theory in the weak and strong coupling regimes, and a new interpolation is proposed.

Orbital-free molecular dynamics simulations are used to benchmark two popular models for hot dense plasmas: the one component plasma (OCP) and the Yukawa model. A unified concept emerges where an effective OCP (eOCP) is constructed from the short-range structure of the plasma. An unambiguous ionization and the screening length can be defined and used for a Yukawa system, which reproduces the long range structure with finite compressibility. Similarly, the dispersion relation of longitudinal waves is consistent with the screened model at vanishing wavenumber but merges with the OCP at high wavenumber. Additionally, the eOCP reproduces the overall relaxation timescales of the correlation functions associated with ionic motion. In the hot dense regime, this unified concept of eOCP can be fruitfully applied to deduce properties such as the equation of state, ionic transport coefficients, and the ion feature in x-ray Thomson scattering experiments., 5 pages, 4 Figures and Supplemental with 1 Figure Accepted for publication in PRL, Physical Review Letters 2016

Performing a systematic Bogoliubov--de Gennes spectral analysis, we illustrate that stationary vortex lines, vortex rings, and more exotic states, such as hopfions, are robust in three-dimensional atomic Bose-Einstein condensates, for large parameter intervals. Importantly, we find that the hopfion can be stabilized in a simple parabolic trap, without the need for trap rotation or inhomogeneous interactions. We supplement our spectral analysis by studying the dynamics of such stationary states; we find them to be robust against significant perturbations of the initial state. In the unstable regimes, we not only identify the unstable mode, such as a quadrupolar or hexapolar mode, but we also observe the corresponding instability dynamics. Furthermore, deep in the Thomas-Fermi regime, we investigate the particlelike behavior of vortex rings and hopfions.

We have performed finite-temperature density functional theory molecular dynamics simulations on dense methane, ammonia, and water mixtures (CH4:NH3:H2O) for various compositions and temperatures (2000 K ≤ T ≤ 10,000 K) that span a set of possible conditions in the interiors of ice-giant exoplanets. The equation-of-state, pair distribution functions, and bond autocorrelation functions (BACF) were used to probe the structure and dynamics of these complex fluids. In particular, an improvement to the choice of the cutoff in the BACF was developed that allowed analysis refinements for density and temperature effects. We note the relative changes in the nature of these systems engendered by variations in the concentration ratios. A basic tenet emerges from all these comparisons that varying the relative amounts of the three heavy components (C,N,O) can effect considerable changes in the nature of the fluid and may in turn have ramifications for the structure and composition of various planetary layers.

Mass transport properties (shear viscosity and diffusion coefficients) for Pu/DT mixtures were calculated with Orbital Free Molecular Dynamics (OFMD). The results were fitted to simple functions of mass density (for ρ=10.4 to 62.4 g/cm3) and temperature (for T=100 up to 3,000 eV) for Pu/DT mixtures consisting of 100/0, 25/75, 50/50, and 75/25 by number.

We present simulations of a four-component mixture of HCNO with orbital free molecular dynamics (OFMD). These simulations were conducted for 5-200 eV with densities ranging between 0.184 and 36.8 g/cm3. We extract the equation of state from the simulations and compare to average atom models. We found that we only need to add a cold curve model to find excellent agreement. Additionally, we studied mass transport properties. We present fits to the self-diffusion and shear viscosity that are able to reproduce the transport properties over the parameter range studied. We compare these OFMD results to models based on the Coulomb coupling parameter and one-component plasmas.

In the present work, we investigate how single- and multi-vortex-ring states can emerge from a planar dark soliton in three-dimensional (3D) Bose-Einstein condensates (confined in isotropic or anisotropic traps) through bifurcations. We characterize such bifurcations quantitatively using a Galerkin-type approach, and find good qualitative and quantitative agreement with our Bogoliubov-de Gennes (BdG) analysis. We also systematically characterize the BdG spectrum of the dark solitons, using perturbation theory, and obtain a quantitative match with our 3D BdG numerical calculations. We then turn our attention to the emergence of single- and multi-vortex-ring states. We systematically capture these as stationary states of the system and quantify their BdG spectra numerically. We find that although the vortex ring may be unstable when bifurcating, its instabilities weaken and may even eventually disappear, for sufficiently large chemical potentials and suitable trap settings. For instance, we demonstrate the stability of the vortex ring for an isotropic trap in the large chemical potential regime., 14 pages, 9 figures

In the warm dense matter (WDM) regime, material properties like diffusion and viscosity can be obtained from lengthy quantum molecular dynamics simulations, where the quantum behavior of the electrons is represented using either Kohn-Sham or orbital-free density functional theory. To reduce the simulation duration, we fit the time dependence of the autocorrelation functions (ACFs) and then use the fit to find values of the diffusion and viscosity. This fitting procedure avoids noise in the long time behavior of the ACFs. We present a detailed analysis of the functional form used to fit the ACFs, which is always a more efficient means to obtain mass transport properties. We use the fits to estimate the statistical error of the transport properties. We apply this methodology to a dense correlated plasma of copper and a mixture of carbon and hydrogen. Both systems show structure in their ACFs and exhibit multiple time scales. The mixture contains different structural forms of the ACFs for each component in the mixture.

We present a theory for the number fluctuations of a quasi-two-dimensional (quasi-2D) dipolar Bose-Einstein condensate measured with finite resolution cells. We show that when the dipoles are tilted to have a component parallel to the plane of the trap, the number fluctuations become anisotropic, i.e. depend on the in-plane orientation of the measurement cell. We develop analytic results for the quantum and thermal fluctuations applicable to the cell sizes accessible in experiments. We show that as cell size is increased the thermodynamic fluctuation result is approached much more slowly than in condensates with short range interactions, so experiments would not require high numerical aperture imaging to observe the predicted effect., Comment: 5 pages, 4 figures

We have performed nonequilibrium classical and quantum-mechanical molecular dynamics simulations that follow the interpenetration of deuterium-tritium (DT) and carbon (C) components through an interface initially in hydrostatic and thermal equilibrium. We concentrate on the warm, dense matter regime with initial densities of 2.5--5.5 g/cm${}^{3}$ and temperatures from 10 to 100 eV. The classical treatment employs a Yukawa pair-potential with the parameters adjusted to the plasma conditions, and the quantum treatment rests on an orbital-free density functional theory at the Thomas-Fermi-Dirac level. For times greater than about a picosecond, the component concentrations evolve in accordance with Fick's law for a classically diffusing fluid with the motion, though, described by the mutual diffusion coefficient of the mixed system rather than the self-diffusion of the individual components. For shorter times, microscopic processes control the clearly non-Fickian dynamics and require a detailed representation of the electron probability density in space and time.

A recent and unexpected discrepancy between \textit{ab initio} simulations and the interpretation of a laser shock experiment on aluminum, probed by X-ray Thomson scattering (XRTS), is addressed. The ion-ion structure factor deduced from the XRTS elastic peak (ion feature) is only compatible with a strongly coupled out-of-equilibrium state. Orbital free molecular dynamics simulations with ions colder than the electrons are employed to interpret the experiment. The relevance of decoupled temperatures for ions and electrons is discussed. The possibility that it mimics a transient, or metastable, out-of-equilibrium state after melting is also suggested., Comment: 12 pages, 3 figures

We demonstrate that measurements of number fluctuations within finite cells provide a direct means to study fluctuation scaling in a trapped two-component condensate. This quantum system supports a second-order phase transition between miscible (co-spatial) and immiscible (symmetry-broken) states that is driven by a diverging susceptibility to magnetic fluctuations. As the transition is approached from the miscible side the magnetic susceptibility is found to depend strongly on the geometry and orientation of the observation cell. However, a scaling exponent consistent with that for the homogenous gas ($\gamma = 1$) can be recovered, for all cells considered, as long as the fit excludes the region in the immediate vicinity of the critical point. As the transition is approached from the immiscible side, the magnetic fluctuations exhibit a non-trivial scaling exponent $\gamma \simeq 1.30$. Experimentally, the observation cells may be formed either by considering individual imaging pixels or by combining pixels to form larger cells, and fluctuation statistics can be obtained by repeated \emph{in situ} images. Interestingly, on both sides of the transition, we find it best to extract the exponents using an observation cell that covers half of the trapped system. This implies that relatively low-resolution \emph{in situ} imaging will be adequate for the investigation of these exponents. We also investigate the gap energy and find exponents $\nu z$ = 0.505 on the miscible side and, unexpectedly, $\nu z$ = 0.60(3) for the immiscible phase., Comment: 7 pages, 3 figures

We present a study for the dispersion relation and character of the excitations of a single- and two-component Bose-Einstein condensate (BEC). We study the single-component dispersion for a finite BEC system and look at examples of quasiparticles to understand and characterize the dispersion relation. Next we present the dispersion relation for a two-component BEC in both the miscible and immiscible parameter regimes. Then we present examples of the quasiparticles for both regimes.

We consider the fluctuations in atom number that occur within finite-sized measurement cells in a trapped Bose-Einstein condensate (BEC). This approximates the fluctuation measurements made in current experiments with finite resolution in situ imaging. A numerical scheme is developed to calculate these fluctuations using the quasiparticle modes of a cylindrically symmetric three-dimensionally trapped condensate with either contact or dipole-dipole interactions (DDIs). We use this scheme to study the properties of a pancake shaped condensate using cylindrical cells. The extension of the theory to washer shaped cells with azimuthal weighting is made and used to discriminate between the low energy roton modes in a dipolar condensate according to their pro- jection of angular momentum. Our results are based on the Bogoliubov approach valid for zero and small finite temperatures., 11 pages, 9 figures

Extending the well-known Thomas-Fermi $Z$-scaling laws to the Coulomb coupling parameter, we investigate the stabilization of the ionic coupling in isochoric heating [Cl\'erouin et al., Phys. Rev. E 87, 061101 (2013)]. This stabilization is restricted to a domain in atomic number $Z$, temperature, and density, including strong limitations on high couplings, that can only be obtained for high-$Z$ elements. Contact is made with recent isochoric heating experiments. The consequences for corresponding states with respect to ionic coupling are also quantified via orbital free molecular dynamics simulations. This opens avenues for future isochoric heating experiments.