Quantum simulation studies of metal–ammonia solutions

Quantifying solvated electrons' delocalization

The electron delocalization range EDR(r;u_av) (left) captures the spin density (right) of an electron delocalized over u_av = 5.77 Å on the surface of an (H₂O)₂₀⁻ cluster.

The identification of a solvated electron pair in the gaseous clusters of Na(-)(H2O)n and Li(-)(H2O)n

By first principles calculations, we explore the possibility that Na(-)(H(2)O)(n) and Li(-)(H(2)O)(n) clusters, which have been measured previously by photoelectron experiments, could serve as gas-phase molecular models for the solvation of two electrons. Such models would capture the electron-electron interaction in a solution environment, which is missed in the well-known anionic water clusters (H(2)O)(n) (-). Our results show that by n = 10, the two loosely bound s electrons in Li(-)(H(2)O)(n) are indeed detached from lithium, and they could exist in either the singlet (spin-paring) or the triplet (spin-coupling) state. In contrast, the two electrons would prefer to stay on the sodium atom in Na(-)(H(2)O)(n) and on the surface of the cluster. The formation of a solvated electron pair and the variation in solvation structures make these two cluster series interesting subjects for further experimental investigation.

Dispersion forces between solvated electrons

Using the path integral centroid approach, we investigate dispersion interactions between electrons solvated in metal-ammonia solutions. We have argued that at finite metal concentrations, the behavior of the solvated electrons is controlled by these interactions. The latter result in a peculiar nonmetal-metal transition, which appears as a sharp dielectric enhancement and a mechanical instability of the system. Our results indicate also that the solvated electrons are to be considered as a two-component mixture consisting of localized and delocalized electrons beyond the critical density corresponding to this mechanical instability.

A sodium atom in a large water cluster: electron delocalization and infrared spectra

Ab initio molecular dynamics simulations modeling low-energy collisions of a sodium atom with a cluster with more than 30 water molecules are presented. We follow the dynamics of the atom-cluster interaction and the delocalization of the valence electron of sodium together with the changes in the electron binding energy. This electron tends to be shared by the nascent sodium cation and the water cluster. IR spectra of the sodium-water cluster are both computationally and experimentally obtained, with a good agreement between the two approaches.

Nature of metal-nonmetal transition in metal-ammonia solutions. II. From uniform metallic state to inhomogeneous electronic microstructure

Applying semianalytical models of nonideal plasma, we evaluate the behavior of the metallic phase in metal-ammonia solutions (MAS). This behavior is mainly controlled by the degenerate electron gas, which remains stable down to 5 MPM due to high solvent polarizability and strong dielectric screening of solvated ions. Comparing the behavior of the metallic state with those of localized solvated electrons, we have estimated the miscibility gap Delta n for various alkali metals and found Delta n(Na)>Delta n(K). It is rather narrow in Rb-NH3 and does not occur in Cs-NH3 solutions, which is in full agreement with the experiments. The case of Li is discussed separately. The difference calculated in the excess free energies of the metallic and nonmetallic phases is in the order of kBT, yielding a thermally fluctuating mixed state at intermediate metal concentrations. It results in a continuous metal-nonmetal (MNM) transition above the consolute point Tc and a phase separation below Tc. We propose a criterion for the MNM transition which may be attributed to the line of the maximum of compressibility above Tc. This line crosses the spinodal one at the critical temperature. Finally, we assert that a new electronic phase similar to microemulsion should also arise between the spinodal and the binodal lines.

Low variance energy estimators for systems of quantum Drude oscillators: Treating harmonic path integrals with large separations of time scales

In the effort to develop atomistic models capable of accurately describing nanoscale systems with complex interfaces, it has become clear that simple treatments with rigid charge distributions and dispersion coefficients selected to generate bulk properties are insufficient to predict important physical properties. The quantum Drude oscillator model, a system of one-electron pseudoatoms whose "pseudoelectrons" are harmonically bound to their respective "pseudonuclei," is capable of treating many-body polarization and dispersion interactions in molecular systems on an equal footing due to the ability of the pseudoatoms to mimic the long-range interactions that characterize real materials. Using imaginary time path integration, the Drude oscillator model can, in principle, be solved in computer operation counts that scale linearly with the number of atoms in the system. In practice, however, standard expressions for the energy and pressure, including the commonly used virial estimator, have extremely large variances that require untenably long simulation times to generate converged averages. In this paper, low-variance estimators for the internal energy are derived, in which the large zero-point energy of the oscillators does not contribute to the variance. The new estimators are applicable to any system of harmonic oscillators coupled to one another (or to the environment) via an arbitrary set of anharmonic interactions. The variance of the new estimators is found to be much smaller than standard estimators in three example problems, a one-dimensional anharmonic oscillator and quantum Drude models of the xenon dimer and solid (fcc) xenon, respectively, yielding 2-3 orders of magnitude improvement in computational efficiency.

Ionization potentials of large sodium doped ammonia clusters

In a continuous neat supersonic expansion ammonia clusters are generated and doped with sodium atoms in a pickup cell. Thus clusters of the form Na(NH(3))(n) are produced that are photoionized by a tunable dye laser system. The ions are mass analyzed in a reflectron time-of-flight mass spectrometer, and the wavelength dependent ion signals serve for the determination of the ionization potentials (IP) of the different clusters in the size range 10< or =n< or =1500. Aside from a plateau for 10< or =n< or =17 and smaller steps at n=24, 35, and 59 on the average a continuous decrease of the IP with cluster size is observed. The IPs in this size range are linear with (n+1)(-13) and extrapolate to IP(n=infinity)=1.66+/-0.01 eV. The slope is consistent with a dielectric continuum model of the solvated electron and the dielectric constant of the solid. The extrapolated IPs are compared with results obtained for negative ammonia cluster ions and metallic solutions in liquid ammonia. Differences are explained by the presence of counterions and their various distances from the solvated electron.

Scalable fine-grained parallelization of plane-wave-based ab initio molecular dynamics for large supercomputers

Many systems of great importance in material science, chemistry, solid-state physics, and biophysics require forces generated from an electronic structure calculation, as opposed to an empirically derived force law to describe their properties adequately. The use of such forces as input to Newton's equations of motion forms the basis of the ab initio molecular dynamics method, which is able to treat the dynamics of chemical bond-breaking and -forming events. However, a very large number of electronic structure calculations must be performed to compute an ab initio molecular dynamics trajectory, making the efficiency as well as the accuracy of the electronic structure representation critical issues. One efficient and accurate electronic structure method is the generalized gradient approximation to the Kohn-Sham density functional theory implemented using a plane-wave basis set and atomic pseudopotentials. The marriage of the gradient-corrected density functional approach with molecular dynamics, as pioneered by Car and Parrinello (R. Car and M. Parrinello, Phys Rev Lett 1985, 55, 2471), has been demonstrated to be capable of elucidating the atomic scale structure and dynamics underlying many complex systems at finite temperature. However, despite the relative efficiency of this approach, it has not been possible to obtain parallel scaling of the technique beyond several hundred processors on moderately sized systems using standard approaches. Consequently, the time scales that can be accessed and the degree of phase space sampling are severely limited. To take advantage of next generation computer platforms with thousands of processors such as IBM's BlueGene, a novel scalable parallelization strategy for Car-Parrinello molecular dynamics is developed using the concept of processor virtualization as embodied by the Charm++ parallel programming system. Charm++ allows the diverse elements of a Car-Parrinello molecular dynamics calculation to be interleaved with low latency such that unprecedented scaling is achieved. As a benchmark, a system of 32 water molecules, a common system size employed in the study of the aqueous solvation and chemistry of small molecules, is shown to scale on more than 1500 processors, which is impossible to achieve using standard approaches. This degree of parallel scaling is expected to open new opportunities for scientific inquiry.

Nose-Hoover chain method for nonequilibrium molecular dynamics simulation

It is demonstrated that the Nose-Hoover chain (NHC) thermostat method influences the properties of a system out of equilibrium. A simple modification of the NHC equations is proposed which enables both steady-state averages and time correlation functions to be obtained correctly for nonequilibrium states, and consistent with those computed using Gaussian isokinetic and Nose-Hoover dynamics.