Finite difference Coulombic potential

For the case cited above, the ponderomotive energy is approximately 1 eV. For typical short pulse experiments today, this energy can easily be hundreds of electron volts. Therefore the wave function of a photoelectron in an intense laser field does not resemble that of the normal field-free Coulomb state, but is dressed by the field, becoming, in the absence of a binding potential, a Volkov state [5], This complex motion of the photoelectrons in the continuum is very difficult to reproduce in terms of the field-free atomic basis functions, so that we have chosen to define our electron wave functions on a finite difference grid. These numerical wave functions have the flexibility to represent both the bound and continuum states in the laser field accurately. 

In classical molecular dynamics simulations, atoms are generally considered to be points which interact with other atoms by some predehned potential form. The forms of the potential can be, for example, Lennard-Jones potentials or Coulomb potentials. The atoms are given velocities in random directions with magnitudes selected from a Maxwell-Boltzman distribution, and then they are allowed to propagate via Newton s equations of motion according to a finite-difference approximation. See the following references for much more detailed discussions Allen and Tildesley (1987) and Frenkel... 

When (r — c)/a > 1, p(r) 0 and V)y(r) -> —Z/r, and the simple Coulomb potential may be used for r > lOR. Care is required in selecting an integration grid which is sufficiently dense over the nuclear skin region where the charge density is changing rapidly. An analytic Fourier-Bessel expansion which may be used in finite-difference calculations is presented in [75]. However, the complexity of (109) makes the calculation of the potential matrix required in (99)... 

The liquid structure factor of CCI4 and its derivatives with respect to temperature at fixed pressure or fixed volume, needed by eq. (2), were evaluated by Molecular Dynamics (MD) simulations. We have used the OPLS model for tetrachloromethane [9] In this model, the CCI4 molecules are described as rigid tetrahedra (dc-ci = 1 -769 A) and the intermolecular potentials are atom centered 6-12 Lennard-Jones potentials plus the coulombic interaction with partial charges on C and Cl. We performed NVT simulations with 512 molecules for about 1 ns each. The different x-ray structure factors were obtained from the accumulated partial radial distribution functions [10], using the atomic form-factors from the DABAX database [11]. In order to estimate the partial derivatives of the structure factor, we have used finite differences we considered two different temperatures, Ti = 300 K and T2 = 328 K, and two molar volumes, Vi = 97.3 cm mol and V2 = 100.65 cm mol which are the molar volumes along the liquid-vapor coexistence line for the two temperatures Tj and Tz respectively [12]. Three simulations were then run for the temperature and molar volume conditions (TiiVi), T2,V )... 

The arrows above and the symbols below the interfaces indicate the transfer of the charge at each interface when the concentration of NaF in the sample is abruptly increased. It is possible to estimate the actual number of ions that are required to establish the potential difference at the interfaces. A typical value for the doublelayer capacitor is 10 5 F cm 2. If a potential difference of n = 100 mV is established at this interface, the double-layer capacitor must be charged by the charge Q = nCdi = 10 6 coulombs. From Faraday s law (6.3), we see that it corresponds to approximately 10 11 mol cm 2 or 1012 ions cm 2 of the electrode surface area. Thus, a finite amount of the potential determining ions is removed from the sample but this charge is replenished through the liquid junction, in order to maintain electroneutrality. 

---