Static, dynamic, and electronic properties of liquid gallium studied by first-principles simulation

Dynamics on next-neighbour distances in liquid and undercooled gallium

Density fluctuations of liquid and 20 K undercooled gallium have been studied by neutron spectroscopy. The decay of density fluctuations has been recorded at the structure factor maximum over a wide temperature range up to twice the melting temperature. The amplitude of the scattering function falls off with rising temperature in a nonlinear way with a changing slope around [Formula: see text]. The derived generalized longitudinal viscosity shows an upturn with decreasing temperature in the same temperature range. This increase in viscosity can be understood that liquid gallium transforms from a more fluid liquid metal to a more viscous liquid metal in that temperature range upon cooling. The change in the amplitude shows a remarkable agreement with results from liquid aluminium, lead and rubidium. This study suggests a universal crossover in dynamics of liquid monatomic metals, despite the many peculiar properties of gallium.

Real-space grids and the Octopus code as tools for the development of new simulation approaches for electronic systems

Real-space grids are a powerful alternative for the simulation of electronic systems. One of the main advantages of the approach is the flexibility and simplicity of working directly in real space where the different fields are discretized on a grid, combined with competitive numerical performance and great potential for parallelization. These properties constitute a great advantage at the time of implementing and testing new physical models. Based on our experience with the Octopus code, in this article we discuss how the real-space approach has allowed for the recent development of new ideas for the simulation of electronic systems. Among these applications are approaches to calculate response properties, modeling of photoemission, optimal control of quantum systems, simulation of plasmonic systems, and the exact solution of the Schrödinger equation for low-dimensionality systems.

Dynamical and transport properties of liquid gallium at high pressures

Quantum molecular dynamics (QMD) simulations are used to calculate the equation of state, structure, and transport properties of liquid gallium along the principal shock Hugoniot. The calculated Hugoniot is in very good agreement with experimental data up to a pressure of 150 GPa as well as with our earlier classical molecular dynamics calculations using a modified embedded atom method (MEAM) potential. The self-diffusion and viscosity calculated using QMD agree with experimental measurements better than the MEAM results, which we attribute to capturing the complexity of the electronic structure at elevated temperatures. Calculations of the DC conductivity were performed around the Hugoniot. Above a density of 7.5 g/cm(3), the temperature increases rapidly along the Hugoniot, and the optical conductivity decreases, indicating simple liquid metal behavior.

Transverse acoustic excitations in liquid Ga

The transverse acoustic excitation modes were detected by inelastic x-ray scattering in liquid Ga in the Q range above 9 nm(-1) although liquid Ga is mostly described by a hard-sphere liquid. An ab initio molecular dynamics simulation clearly supports this finding. From the detailed analysis for the S(Q,omega) spectra with a good statistic quality, the lifetime of 0.5 ps and the propagating length of 0.4-0.5 nm can be estimated for the transverse acoustic phonon modes, which may correspond to the lifetime and size of cages formed instantaneously in liquid Ga.

The possibility of transverse excitation modes in liquid Ga

The dynamic structure factor S(Q,ω) of liquid Ga was measured at 100 °C using a high resolution inelastic x-ray scattering (IXS) spectrometer at 3-ID-C/APS. The spectra obtained clearly demonstrate the existence of longitudinal propagating modes at small Q values, like a previous IXS result at 42 °C obtained by Scopigno et al and an inelastic neutron scattering (INS) one at 47 °C obtained by Bove et al, but unlike an INS study at 57 °C by Bermejo et al. The dispersion relation of the excitations deviates positively from the hydrodynamic prediction by about 13%. There are two new findings from this experiment. Firstly, an additional lower energy excitation is necessary to reproduce S(Q,ω) spectra in the Q range beyond 10 nm(-1), in agreement with the result of a first-principles molecular dynamic simulation, which may indicate a transverse acoustic mode in this peculiar liquid metal. Secondly, the quasielastic line comprises a Gaussian contribution at Q values near the first maximum in S(Q), which may indicate the existence of short-lived covalent correlation in liquid Ga with a lifetime of 0.39 ps.

Ab initio molecular dynamics simulation of liquid Al88Si12 alloys

First-principles molecular dynamics simulations are carried out to study the structures, dynamics, and electronic properties of liquid Al88Si12 in the temperature ranging from 898 to 1298 K. The temperature dependence of static structure factors, pair correlation functions, and electronic density-of-states are investigated. The structural properties obtained from the simulations are in good agreement with the x-ray diffraction experimental results.

A finite temperature linear tetrahedron method for electronic structure calculations of periodic systems

A finite-temperature linear tetrahedron method for electronic structure calculations of periodic systems is developed. When compared to widely used simple temperature broadening, the number of k points necessary for accurate integration at finite temperatures is reduced. The utility of the method is demonstrated with benchmark calculations on 1D, 2D, and 3D systems.

Electrical conductivity of hot expanded aluminum: experimental measurements and ab initio calculations

Experimental measurements and theoretical calculations of the electrical conductivity of aluminum are presented in the strongly coupled partially degenerate regime (rho=0.3 g/cm(3), 5000<T<15 000 K). The experiments were performed in an isochoric plasma closed vessel designed to confine electrical plasma discharges up to 1.5 GPa. Aluminum properties were determined theoretically by ab initio molecular dynamics simulations in the local density approximation, from which the conductivity was computed using the Kubo-Greenwood formula. The theoretical results were validated in the dense coupled regime against previously published experimental results and then applied to our experimental low density regime, showing that the theoretical results overestimate the experimental conductivities.

Predicting the liquid-vapor critical point from the crystal anharmonicity

A "universal" dependence is predicted of the reduced critical parameters, k(B)T(c) / E0(gamma), V(c) / V0(gamma), and P(c)V(c)/k(B)T(c) = Z(c)(gamma), on the crystal anharmonicity gamma (closely related to the Gruneisen parameter gamma(G)). It is based on a simplified embedded-atom type approach which enables one to utilize the universal zero-temperature equation of state in a version of fluid perturbation theory. This model's critical temperature and density agree with the experimental results for both the heavy rare gases ( gamma approximately 2.85) and heavy alkali metals ( gamma approximately 1.35). Predicted critical parameters for many other liquid metals are consistent with previous estimates, but the model is not applicable when directional bonding is important.