Fixed-node diffusion Monte Carlo computations for closed-shell jellium spheres

Fixed-node diffusion Monte Carlo computations are used to determine the ground state energy and electron density for jellium spheres with up to N = 106 electrons and background densities corresponding to the electron gas parameter 1 less than or equal to r(s)less than or equal to5.62. We analyze the density and size dependence of the surface energy, and we extrapolate our data to the thermodynamic limit. The results agree well with the predictions of density functional computations using the local density approximation. In the case of N = 20, we extend our computation to higher densities and identify a transition between atomic- and jelliumlike nodal structures occurring at the background density corresponding to r(s)=0.13. In this case the local density approximation is unable to reproduce the changes in the correlation energy due to the discontinuous transition in the ground state nodal structure. We discuss the relevance of our results for nonlocal approximations to density functional theory.

Independent particle descriptions of tunneling using the many-body quantum transport approach

Currents across thin insulators are commonly taken as single electrons moving across classically forbidden regions; this independent particle picture is well-known to describe most tunneling phenomena. Examining quantum transport from a different perspective, i.e., by explicit treatment of electron-electron interactions, we evaluate different single particle approximations with specific application to tunneling in metal-molecule-metal junctions. We find maximizing the overlap of a Slater determinant composed of single-particle states to the many-body current-carrying state is more important than energy minimization for defining single-particle approximations in a system with open boundary conditions. Thus the most suitable single particle effective potential is not one commonly in use by electronic structure methods, such as the Hartree-Fock or Kohn-Sham approximations.

Correlation Effects in the Compton Profile of Silicon

Ab initio nonlocal pseudopotential variational quantum Monte Carlo techniques are used to compute the correlation effects on the valence momentum density and Compton profile of silicon. Our results for this case are in excellent agreement with the Lam-Platzman correction computed within the local density approximation. Within the approximations used, we rule out valence electron correlations as the dominant source of discrepancies between calculated and measured Compton profiles of silicon.

Numerical ‘health check’ for scientific codes: the CADNA approach

Scientific computation has unavoidable approximations built into its very fabric. One important source of error that is difficult to detect and control is round-off error propagation which originates from the use of finite precision arithmetic. We propose that there is a need to perform regular numerical `health checks' on scientific codes in order to detect the cancerous effect of round-off error propagation. This is particularly important in scientific codes that are built on legacy software. We advocate the use of the CADNA library as a suitable numerical screening tool. We present a case study to illustrate the practical use of CADNA in scientific codes that are of interest to the Computer Physics Communications readership. In doing so we hope to stimulate a greater awareness of round-off error propagation and present a practical means by which it can be analyzed and managed.

Local and semilocal density functional computations for crystals of 1-alkyl-3-methyl-imidazolium salts

The accuracy and reliability of popular density functional approximations for the compounds giving origin to room temperature ionic liquids have been assessed by computing the T=0 K crystal structure of several 1-alkyl-3-methyl-imidazolium salts. Two prototypical exchange-correlation approximations have been considered, i.e., the local density approximation (LDA) and one gradient corrected scheme [PBE-GGA, Phys. Rev. Lett. 77, 3865 (1996)]. Comparison with low-temperature x-ray diffraction data shows that the equilibrium volume predicted by either approximations is affected by large errors, nearly equal in magnitude (~10%), and of opposite sign. In both cases the error can be traced to a poor description of the intermolecular interactions, while the intramolecular structure is fairly well reproduced by LDA and PBE-GGA. The PBE-GGA optimization of atomic positions within the experimental unit cell provides results in good agreement with the x-ray structure. The correct system volume can also be restored by supplementing PBE-GGA with empirical dispersion terms reproducing the r-6 attractive tail of the van der Waals interactions.

Quantum revivals in ultrashort intense field dissociation of molecular ions

Ultrashort (<15 fs) high intensity (1014-1016 W cm-2) laser pulses have provided novel methods for investigation of the dynamics of simple molecular ions such as H2+ and D2+. In this paper we report on simulations carried out for the D2+ molecular ion, within the Born- Oppenheimer and two-state approximations. These simulations allow one to investigate the dissociation dynamics of the D2+ molecular ion when subjected to such ultrashort, intense laser pulses. In particular, these simulations are compared to the results from recent pump-probe experiments, in which, the nuclear vibrational motion of D2+ has been imaged. Simulations suggest that the nature of the dissociation process, be it 1- or 2-photon, may be influenced by the tuning of the pump-probe delay time.

The mechanism of N2O formation via the (NO)(2) dimer: A density functional theory study

Catalytic formation of N2O via a (NO)(2) intermediate was studied employing density functional theory with generalized gradient approximations. Dimer formation was not favored on Pt(111), in agreement with previous reports. On Pt(211) a variety of dimer structures were studied, including trans-(NO)(2) and cis-(NO)(2) configurations. A possible pathway involving (NO)(2) formation at the terrace near to a Pt step is identified as the possible mechanism for low-temperature N2O formation. The dimer is stabilized by bond formation between one O atom of the dimer and two Pt step atoms. The overall mechanism has a low barrier of approximately 0.32 eV. The mechanism is also put into the context of the overall NO+H-2 reaction. A consideration of the step-wise hydrogenation of O-(ads) from the step is also presented. Removal of O-(ads) from the step is significantly different from O-(ads) hydrogenation on Pt(111). The energetically favored structure at the transition state for OH(ads) formation has an activation energy of 0.63 eV. Further hydrogenation of OH(ads) has an activation energy of 0.80 eV. (C) 2004 American Institute of Physics.

N2O and NO2 formation on Pt(111): A density functional theory study

Catalytic formation of N2O and NO2 were studied employing density functional theory with generalized gradient approximations, in order to investigate the microscopic reaction pathways of these catalytic processes on a Pt(111) surface. Transition states and reaction barriers for the addition of chemisorbed N or chemisorbed O to NO(ads) producing N2O and NO2, respectively, were calculated. The N2O transition state involves bond formation across the hcp hollow site with an associated reaction barrier of 1.78 eV. NO2 formation favors a fcc hollow site transition state with a barrier of 1.52 eV. The mechanisms for both reactions are compared to CO oxidation on the same surface. The activation of the chemisorbed NO and the chemisorbed N or O from the energetically stable initial state to the transition state are both significant contributors to the overall reaction barrier E-a, in contrast to CO oxidation in which the activation of the O-(ads) is much greater than CO(ads) activation. (C) 2002 American Institute of Physics.

Effect of quantization of vibrations on the structural properties of crystals

We study the structural effects produced by the quantization of vibrational degrees of freedom in periodic crystals at zero temperature. To this end we introduce a methodology based on mapping a suitable subspace of the vibrational manifold and solving the Schrödinger equation in it. A number of increasingly accurate approximations ranging from the quasiharmonic approximation (QHA) to the vibrational self-consistent field (VSCF) method and the exact solution are described. A thorough analysis of the approximations is presented for model monatomic and hydrogen-bonded chains, and results are presented for a linear H-F chain where the potential-energy surface is obtained via first-principles electronic structure calculations. We focus on quantum nuclear effects on the lattice constant and show that the VSCF is an excellent approximation, meaning that correlation between modes is not extremely important. The QHA is excellent for covalently bonded mildly anharmonic systems, but it fails for hydrogen-bonded ones. In the latter, the zero-point energy exhibits a nonanalytic behavior at the lattice constant where the H atoms center, which leads to a spurious secondary minimum in the quantum-corrected energy curve. An inexpensive anharmonic approximation of noninteracting modes appears to produce rather good results for hydrogen-bonded chains for small system sizes. However, it converges to the incorrect QHA results for increasing size. Isotope effects are studied for the first-principles H-F chain. We show how the lattice constant and the H-F distance increase with decreasing mass and how the QHA proves to be insufficient to reproduce this behavior.

Absolute photoionization cross sections for Xe4+, Xe5+, and Xe6+ near 13.5 nm: Experiment and theory

Absolute photoionization cross-section measurements for a mixture of ground and metastable states of Xe4+, Xe5+, and Xe6+ are reported in the photon energy range of 4d -> nf transitions, which occur within or adjacent to the 13.5 nm window for extreme ultraviolet lithography light source development. The reported values allow the quantification of opacity effects in xenon plasmas due to these 4d -> nf autoionizing states. The oscillator strengths for the 4d -> 4f and 4d -> 5f transitions in Xeq+ (q=1-6) ions are calculated using nonrelativistic Hartree-Fock and random phase approximations. These are compared with published experimental values for Xe+ to Xe3+ and with the values obtained from the present experimental cross-section measurements for Xe4+ to Xe6+. The calculations assisted in the determination of the metastable content in the ion beams for Xe5+ and Xe6+. The experiments were performed by merging a synchrotron photon beam generated by an undulator beamline of the Advanced Light Source with an ion beam produced by an electron cyclotron resonance ion source.

Discretization of highly persistent correlated AR(1) shocks

The finite state Markov-chain approximation methods developed by Tauchen (1986) and Tauchen and Hussey (1991) are widely used in economics, finance and econometrics to solve functional equations in which state variables follow autoregressive processes. For highly persistent processes, the methods require a large number of discrete values for the state variables to produce close approximations which leads to an undesirable reduction in computational speed, especially in a multivariate case. This paper proposes an alternative method of discretizing multivariate autoregressive processes. This method can be treated as an extension of Rouwenhorst's (1995) method which, according to our finding, outperforms the existing methods in the scalar case for highly persistent processes. The new method works well as an approximation that is much more robust to the number of discrete values for a wide range of the parameter space.

Three-dimensional pattern formation, multiple homogeneous soft modes, and nonlinear dielectric electroconvection

Patterns forming spontaneously in extended, three-dimensional, dissipative systems are likely to excite several homogeneous soft modes (approximate to hydrodynamic modes) of the underlying physical system, much more than quasi-one- (1D) and two-dimensional (2D) patterns are. The reason is the lack of damping boundaries. This paper compares two analytic techniques to derive the pattern dynamics from hydrodynamics, which are usually equivalent but lead to different results when applied to multiple homogeneous soft modes. Dielectric electroconvection in nematic liquid crystals is introduced as a model for 3D pattern formation. The 3D pattern dynamics including soft modes are derived. For slabs of large but finite thickness the description is reduced further to a 2D one. It is argued that the range of validity of 2D descriptions is limited to a very small region above threshold. The transition from 2D to 3D pattern dynamics is discussed. Experimentally testable predictions for the stable range of ideal patterns and the electric Nusselt numbers are made. For most results analytic approximations in terms of material parameters are given. [S1063-651X(00)09512-X].

A comparison of nonstaggered compact FDTD schemes for the 3D wave equation

This paper aims at providing a better insight into the 3D approximations of the wave equation using compact finite-difference time-domain (FDTD) schemes in the context of room acoustic simulations. A general family of 3D compact explicit and implicit schemes based on a nonstaggered rectilinear grid is analyzed in terms of stability, numerical error, and accuracy. Various special cases are compared and the most accurate explicit and implicit schemes are identified. Further considerations presented in the paper include the direct relationship with other numerical approaches found in the literature on room acoustic modeling such as the 3D digital waveguide mesh and Yee's staggered grid technique.

Chemical processes in protoplanetary disks. II. On the importance of photochemistry and X-ray ionization

We investigate the impact of photochemistry and X-ray ionization on the molecular composition of, and ionization fraction in, a protoplanetary disk surrounding a typical T Tauri star. We use a sophisticated physical model, which includes a robust treatment of the radiative transfer of UV and X-ray radiation, and calculate the time-dependent chemical structure using a comprehensive chemical network. In previous work, we approximated the photochemistry and X-ray ionization; here, we recalculate the photoreaction rates using the explicit UV wavelength spectrum and wavelength-dependent reaction cross sections. We recalculate the X-ray ionization rate using our explicit elemental composition and X-ray energy spectrum. We find that photochemistry has a larger influence on the molecular composition than X-ray ionization. Observable molecules sensitive to the photorates include OH, HCO+, N2H+, H2O, CO2, and CH3OH. The only molecule significantly affected by the X-ray ionization is N2H+, indicating that it is safe to adopt existing approximations of the X-ray ionization rate in typical T Tauri star-disk systems. The recalculation of the photorates increases the abundances of neutral molecules in the outer disk, highlighting the importance of taking into account the shape of the UV spectrum in protoplanetary disks. A recalculation of the photoreaction rates also affects the gas-phase chemistry due to the adjustment of the H/H2 and C+/C ratios. The disk ionization fraction is not significantly affected by the methods adopted to calculate the photochemistry and X-ray ionization. We determine that there is a probable "dead zone" where accretion is suppressed, present in a layer, Z/R lsim 0.1-0.2, in the disk midplane, within R ˜ 200 AU.

Dispersion interactions in room-temperature ionic liquids: Results from a non-empirical density functional

The role of dispersion or van de Waals (VDW) interactions in imidazolium-based room-temperature ionic liquids is studied within the framework of density functional theory, using a recently developed non-empirical functional [M. Dion, H. Rydberg, E. Schroder, D. C. Langreth, and B. I. Lundqvist, Phys. Rev. Lett. 92, 246401 (2004)], as efficiently implemented in the SIESTA code [G. Roman-Perez and J. M. Soler, Phys. Rev. Lett. 103, 096102 (2009)]. We present results for the equilibrium structure and lattice parameters of several crystalline phases, finding a general improvement with respect to both the local density (LDA) and the generalized gradient approximations (GGA). Similar to other systems characterized by VDW bonding, such as rare gas and benzene dimers as well as solid argon, equilibrium distances and volumes are consistently overestimated by approximate to 7%, compared to -11% within LDA and 11% within GGA. The intramolecular geometries are retained, while the intermolecular distances and orientations are significantly improved relative to LDA and GGA. The quality is superior to that achieved with tailor-made empirical VDW corrections ad hoc [M. G. Del Popolo, C. Pinilla, and P. Ballone, J. Chem. Phys. 126, 144705 (2007)]. We also analyse the performance of an optimized version of this non-empirical functional, where the screening properties of the exchange have been tuned to reproduce high-level quantum chemical calculations [J. Klimes, D. Bowler, and A. Michaelides, J. Phys.: Condens. Matter 22, 074203 (2010)]. The results for solids are even better with volumes and geometries reproduced within 2% of experimental data. We provide some insight into the issue of polymorphism of [bmim][Cl] crystals, and we present results for the geometry and energetics of [bmim][Tf] and [mmim][Cl] neutral and charged clusters, which validate the use of empirical force fields. (C) 2011 American Institute of Physics. [doi:10.1063/1.3652897]