Molecular dynamics simulation of ionic liquids: The effect of electronic polarizability

An electronically polarizable model has been developed for the ionic liquid 1-ethyl-3-methylimidazolium nitrate (EMIM+/NO3-), Molecular dynamics simulation studies were then performed on both the polarizable and nonpolarizable versions of the model. Comparisons of shear viscosity and diffusion constants at 400 K show that the effects of polarizability are quite substantial and the polarizable model results are in better agreement with the experimental values.

Ion-Atom/Argon - Calculation of ionization cross sections by fast ion impact for neutral target atoms ranging from hydrogen to argon

A FORTRAN 90 program is presented which calculates the total cross sections, and the electron energy spectra of the singly and doubly differential cross sections for the single target ionization of neutral atoms ranging from hydrogen up to and including argon. The code is applicable for the case of both high and low Z projectile impact in fast ion-atom collisions. The theoretical models provided for the program user are based on two quantum mechanical approximations which have proved to be very successful in the study of ionization in ion-atom collisions. These are the continuum-distorted-wave (CDW) and continuum-distorted-wave eikonal-initial-state (CDW-EIS) approximations. The codes presented here extend previously published. codes for single ionization of. target hydrogen [Crothers and McCartney, Comput. Phys. Commun. 72 (1992) 288], target helium [Nesbitt, O'Rourke and Crothers, Comput. Phys. Commun. 114 (1998) 385] and target atoms ranging from lithium to neon [O'Rourke, McSherry and Crothers, Comput. Phys. Commun. 131 (2000) 129]. Cross sections for all of these target atoms may be obtained as limiting cases from the present code. Title of program: ARGON Catalogue identifier: ADSE Program summary URL: http://cpc.cs.qub.ac.uk/cpc/summaries/ADSE Program obtainable from: CPC Program Library Queen's University of Belfast, N. Ireland Licensing provisions: none Computer for which the program is designed and others on which it is operable: Computers: Four by 200 MHz Pro Pentium Linux server, DEC Alpha 21164; Four by 400 MHz Pentium 2 Xeon 450 Linux server, IBM SP2 and SUN Enterprise 3500 Installations: Queen's University, Belfast Operating systems under which the program has been tested: Red-hat Linux 5.2, Digital UNIX Version 4.0d, AIX, Solaris SunOS 5.7 Compilers: PGI workstations, DEC CAMPUS Programming language used: FORTRAN 90 with MPI directives No. of bits in a word: 64, except on Linux servers 32 Number of processors used: any number Has the code been vectorized or parallelized? Parallelized using MPI No. of bytes in distributed program, including test data, etc.: 32 189 Distribution format: tar gzip file Keywords: Single ionization, cross sections, continuum-distorted-wave model, continuum- distorted-wave eikonal-initial-state model, target atoms, wave treatment Nature of physical problem: The code calculates total, and differential cross sections for the single ionization of target atoms ranging from hydrogen up to and including argon by both light and heavy ion impact. Method of solution: ARGON allows the user to calculate the cross sections using either the CDW or CDW-EIS [J. Phys. B 16 (1983) 3229] models within the wave treatment. Restrictions on the complexity of the program: Both the CDW and CDW-EIS models are two-state perturbative approximations. Typical running time: Times vary according to input data and number of processors. For one processor the test input data for double differential cross sections (40 points) took less than one second, whereas the test input for total cross sections (20 points) took 32 minutes. Unusual features of the program: none (C) 2003 Elsevier B.V All rights reserved.

Determining the electronic structure and chemical potentials of molecules in solution

A simulation scheme is proposed for determining the excess chemical potential of a substance in solution. First, a Monte Carlo simulation is performed with classical models for solute and solvent molecules. A representative sample of these configurations is then used in a hybrid quantum/classical (QM/MM) calculation, where the solute is treated quantum-mechanically, and the average electronic structure is used to construct an improved classical model. This procedure is iterated to self-consistency in the classical model, which in practice is attained in one or two steps, depending on the quality of the initial guess. The excess free energy of the molecule within the QM/MM approach is determined relative to the classical model using thermodynamic perturbation theory with a cumulant expansion. The procedure provides a method of constructing classical point charge models appropriate for the solution and gives a measure of the importance of solvent fluctuations.

Quantum electronic transport in a configuration interaction basis

An overview of a many-body approach to calculation of electronic transport in molecular systems is given. The physics required to describe electronic transport through a molecule at the many-body level, without relying on commonly made assumptions such as the Landauer formalism or linear response theory, is discussed. Physically, our method relies on the incorporation of scattering boundary conditions into a many-body wavefunction and application of the maximum entropy principle to the transport region. Mathematically, this simple physical model translates into a constrained nonlinear optimization problem. A strategy for solving the constrained optimization problem is given. (C) 2004 Wiley Periodicals, Inc.

Optimal basis set for electronic structure calculations in periodic systems

An efficient method for calculating the electronic structure of systems that need a very fine sampling of the Brillouin zone is presented. The method is based on the variational optimization of a single (i.e., common to all points in the Brillouin zone) basis set for the expansion of the electronic orbitals. Considerations from k.p-approximation theory help to understand the efficiency of the method. The accuracy and the convergence properties of the method as a function of the optimal basis set size are analyzed for a test calculation on a 16-atom Na supercell.

A simple model of atomic interactions in noble metals based explicitly on electronic structure

A total energy tight-binding model with a basis of just one s state per atom is introduced. It is argued that this simplest of all tight-binding models provides a surprisingly good description of the structural stability and elastic constants of noble metals. By assuming inverse power scaling laws for the hopping integrals and the repulsive pair potential, it is shown that the density matrix in a perfect primitive crystal is independent of volume, and structural energy differences and equations of state are then derived analytically. The model is most likely to be of use when one wishes to consider explicitly and self-consistently the electronic and atomic structures of a generic metallic system, with the minium of computation expense. The relationship to the free-electron jellium model is described. The applicability of the model to other metals is also considered briefly.