QMCPACK is an open source quantum Monte Carlo package for ab initio electronic structure calculations. It supports calculations of metallic and insulating solids, molecules, atoms, and some model Hamiltonians. Implemented real space quantum Monte Carlo algorithms include variational, diffusion, and reptation Monte Carlo. QMCPACK uses Slater-Jastrow type trial wavefunctions in conjunction with a sophisticated optimizer capable of optimizing tens of thousands of parameters. The orbital space auxiliary-field quantum Monte Carlo method is also implemented, enabling cross validation between different highly accurate methods. The code is specifically optimized for calculations with large numbers of electrons on the latest high performance computing architectures, including multicore central processing unit and graphical processing unit systems. We detail the program's capabilities, outline its structure, and give examples of its use in current research calculations. The package is available at http://qmcpack.org.

In multiphase flow in porous media the consistent pore to Darcy scale description of two-fluid flow processes has been a long-standing challenge. Immiscible displacement processes occur at the scale of individual pores. However, the larger scale behavior is described by phenomenological relationships such as relative permeability, which typically uses only fluid saturation as a state variable. As a consequence pore scale properties such as contact angle cannot be directly related to Darcy scale flow parameters. Advanced imaging and computational technologies are closing the gap between the pore and Darcy scale, supporting the development of new theory. We utilize fast x-ray microtomography to observe pore-scale two-fluid configurations during immiscible flow and initialize lattice Boltzmann simulations that demonstrate that the mobilization of disconnected nonwetting phase clusters can account for a significant fraction of the total flux. We show that fluid topology can undergo substantial changes during flow at constant saturation, which is one of the underlying causes of hysteretic behavior. Traditional assumptions about fluid configurations are therefore an oversimplification. Our results suggest that the role of fluid connectivity cannot be ignored for multiphase flow. On the Darcy scale, fluid topology and phase connectivity are accounted for by interfacial area and Euler characteristic as parameters that are missing from our current models.

Models that describe two-fluid flow in porous media suffer from a widely-recognized problem that the constitutive relationships used to predict capillary pressure as a function of the fluid saturation are non-unique, thus requiring a hysteretic description. As an alternative to the traditional perspective, we consider a geometrical description of the capillary pressure, which relates the average mean curvature, the fluid saturation, the interfacial area between fluids, and the Euler characteristic. The state equation is formulated using notions from algebraic topology and cast in terms of measures of the macroscale state. Synchrotron-based X-ray micro-computed tomography (µCT) and highresolution pore-scale simulation is applied to examine the uniqueness of the proposed relationship for six different porous media. We show that the geometric state function is able to characterize the microscopic fluid configurations that result from a wide range of simulated flow conditions in an averaged sense. The geometric state function can serve as a closure relationship within macroscale models to effectively remove hysteretic behavior attributed to the arrangement of fluids within a porous medium. This provides a critical missing component needed to enable a new generation of higher fidelity models to describe two-fluid flow in porous media. PACS numbers: 92.40Cy, 92.40.Kf, 91.60.Tn, 91.60.Fe

Multiphase flows in porous medium systems are typically modeled at the macroscale by applying the principles of continuum mechanics to develop models that describe the behavior of averaged quantities, such as fluid pressure and saturation. These models require closure relations to produce solvable forms. One of these required closure relations is an expression relating the capillary pressure to fluid saturation and, in some cases, other topological invariants such as interfacial area and the Euler characteristic (or average Gaussian curvature). The forms that are used in traditional models, which typically consider only the relationship between capillary pressure and saturation, are hysteretic. An unresolved question is whether the inclusion of additional morphological and topological measures can lead to a nonhysteretic closure relation. Relying on the lattice Boltzmann (LB) method, we develop an approach to investigate equilibrium states for a two-fluid-phase porous medium system, which includes disconnected nonwetting phase features. A set of simulations are performed within a random close pack of 1964 spheres to produce a total of 42 908 distinct equilibrium configurations. This information is evaluated using generalized additive models to quantitatively assess the degree to which functional relationships can explain the behavior of the equilibrium data. The variance of various model estimates is computed, and we conclude that, except for the limiting behavior close to a single fluid regime, capillary pressure can be expressed as a deterministic and nonhysteretic function of fluid saturation, interfacial area between the fluid phases, and the Euler characteristic. To our knowledge, this work is unique in the methods employed, the size of the data set, the resolution in space and time, the true equilibrium nature of the data, the parametrizations investigated, and the broad set of functions examined. The conclusion of essentially nonhysteretic behavior provides support for an evolving class of two-fluid-phase flow in porous medium systems models.

No abstract

Plasma-based x-ray lasers allow single-shot nano-scale imaging and other experiments requiring a large number of photons per pulse to be conducted in compact facilities. However, compact repetitively fired gain-saturated x-ray lasers have been limited to wavelengths above λ 8.85 nm. Here we extend their range to λ 6.85 nm by transient traveling wave excitation of Ni-like Gd ions in a plasma created with an optimized pre-pulse followed by rapid heating with an intense sub-picosecond pump pulse. Isoelectronic scaling also produced strong lasing at 6.67 nm and 6.11 nm in Ni-like Tb and amplification at 6.41 nm and 5.85 nm in Ni-like Dy. This scaling to shorter wavelengths was obtained by progressively increasing the pump pulse grazing incidence angle to access increased plasma densities. We experimentally demonstrate that the optimum grazing incidence angle increases linearly with atomic number from 17 deg for Z 42 (Mo) to 43 deg for Z 66 (Dy). The results will enable applications of sub-7 nm lasers at compact facilities.

We report a comprehensive study of the emission from Gd and Tb laser produced plasmas in the 6.5-6.7 nm wavelength region for a broad range of laser irradiation parameters using a single λ = 1030 nm laser with tunable pulse duration in the 120 ps to 4 ns range. The results are of interest for beyond extreme ultraviolet (BEUV) lithography of integrated circuits. BEUV emission spectra were measured as a function of laser pulse duration, emission angle, and

We have developed a set of reduced, proxy applications (''MiniApps'') based on large-scale application codes supported at the Oak Ridge Leadership Computing Facility (OLCF). The MiniApps are designed to encapsulate the details of the most important (i.e. the most time-consuming and/or unique) facets of the applications that run in production mode on the OLCF. In each case, we have produced or plan to produce individual versions of the MiniApps using different specific programing models (e.g., OpenACC, CUDA, OpenMP). We describe some of our initial observations regarding these different implementations along with estimates of how closely the MiniApps track the actual performance characteristics (in particular, the overall scalability) of the large-scale applications from which they are derived.