Development of complex classical force fields through force matching to ab initio data: Application to a room-temperature ionic liquid

Recent experimental neutron diffraction data and ab initio molecular dynamics simulation of the ionic liquid dimethylimidazolium chloride ([dmim]Cl) have provided a structural description of the system at the molecular level. However, partial radial distribution functions calculated from the latter, when compared to previous classical simulation results, highlight some limitations in the structural description offered by force fieldbased simulations. With the availability of ab initio data it is possible to improve the classical description of [dmim]Cl by using the force matching approach, and the strategy for fitting complex force fields in their original functional form is discussed. A self-consistent optimization method for the generation of classical potentials of general functional form is presented and applied, and a force field that better reproduces the observed first principles forces is obtained. When used in simulation, it predicts structural data which reproduces more faithfully that observed in the ab initio studies. Some possible refinements to the technique, its application, and the general suitability of common potential energy functions used within many ionic liquid force fields are discussed.

Current-driven magnetic rearrangements in spin-polarized point contacts

A method for investigating the dynamics of atomic magnetic moments in current-carrying magnetic point contacts under bias is presented. This combines the nonequilibrium Green's function (NEGF) method for evaluating the current and the charge density with a description of the dynamics of the magnetization in terms of quasistatic thermally activated transitions between stationary configurations. This method is then implemented in a tight-binding (TB) model with parameters chosen to simulate the main features of the electronic structures of magnetic transition metals. We investigate the domain wall (DW) migration in magnetic monoatomic chains sandwiched between magnetic leads, and for realistic parameters find that collinear arrangement of the magnetic moments of the chain is always favorable. Several stationary magnetic configurations are identified, corresponding to a different number of Bloch walls in the chain and to a different current. The relative stability of these configurations depends on the geometrical details of the junction and on the bias; however, we predict transitions between different configurations with activation barriers of the order of a few tens of meV. Since different magnetic configurations are associated with different resistances, this suggests an intrinsic random telegraph noise at microwave frequencies in the I-V curves of magnetic atomic point contacts at room temperature. Finally, we investigate whether or not current-induced torques are conservative.

Correlated electron-ion dynamics with open boundaries: formalism

We extend a new formalism, which allows correlated electron-ion dynamics to be applied to the problem of open boundary conditions. We implement this at the first moment level (allowing heating of ions by electrons) and observe the expected cooling in the classical part of the ionic kinetic energy and current-induced heating in the quantum contribution. The formalism for open boundaries should be easily extended to higher moments of the correlated electron-ion fluctuations.

Correlated electron-ion dynamics: the excitation of atomic motion by energetic electrons

Correlated electron-ion dynamics (CEID) is an extension of molecular dynamics that allows us to introduce in a correct manner the exchange of energy between electrons and ions. The formalism is based on a systematic approximation: small amplitude moment expansion. This formalism is extended here to include the explicit quantum spread of the ions and a generalization of the Hartree-Fock approximation for incoherent sums of Slater determinants. We demonstrate that the resultant dynamical equations reproduce analytically the selection rules for inelastic electron-phonon scattering from perturbation theory, which control the mutually driven excitations of the two interacting subsystems. We then use CEID to make direct numerical simulations of inelastic current-voltage spectroscopy in atomic wires, and to exhibit the crossover from ionic cooling to heating as a function of the relative degree of excitation of the electronic and ionic subsystems.

Beyond Ehrenfest: correlated non-adiabatic molecular dynamics

A method for introducing correlations between electrons and ions that is computationally affordable is described. The central assumption is that the ionic wavefunctions are narrow, which makes possible a moment expansion for the full density matrix. To make the problem tractable we reduce the remaining many-electron problem to a single-electron problem by performing a trace over all electronic degrees of freedom except one. This introduces both one- and two-electron quantities into the equations of motion. Quantities depending on more than one electron are removed by making a Hartree-Fock approximation. Using the first-moment approximation, we perform a number of tight binding simulations of the effect of an electric current on a mobile atom. The classical contribution to the ionic kinetic energy exhibits cooling and is independent of the bias. The quantum contribution exhibits strong heating, with the heating rate proportional to the bias. However, increased scattering of electrons with increasing ionic kinetic energy is not observed. This effect requires the introduction of the second moment.

Time delay between singly and doubly ionizing wavepackets in laser-driven helium

We present calculations of the time delay between single and double ionization of helium, obtained from full-dimensionality numerical integrations of the helium-laser Schroedinger equation. The notion of a quantum mechanical time delay is defined in terms of the interval between correlated bursts of single and double ionization. Calculations are performed at 390 and 780 nm in laser intensities that range from 2 X 10^14 to 14 X 10^14 W /cm^2. We find results consistent with the rescattering model of double ionization but supporting its classical interpretation only at 780 nm.

Impulsive electric fields driven by high-intensity interactions

The interaction of high-intensity laser pulses with matter releases instantaneously ultra-large currents of highly energetic electrons, leading to the generation of highly-transient, large-amplitude electric and magnetic fields. We report results of recent experiments in which such charge dynamics have been studied by using proton probing techniques able to provide maps of the electrostatic fields with high spatial and temporal resolution. The dynamics of ponderomotive channeling in underdense plasmas have been studied in this way, as also the processes of Debye sheath formation and MeV ion front expansion at the rear of laser-irradiated thin metallic foils. Laser-driven impulsive fields at the surface of solid targets can be applied for energy-selective ion beam focusing.

Sub-fossil Insects and ecosystem dynamics in wetlands; implications for biodiversity and conservation.

We review the uses of fossil insects, particularly Coleoptera (beetles) and Chironomidae (non-biting midges) from ancient deposits to inform the study of wetland ecosystems and their ecological and restoration processes. In particular, we focus on two contrasting ecosystems, drawing upon research undertaken by us on British raised mire peats and shallow lake systems, one an essentially terrestrial ecosystem, the other aquatic, but in which wetland insects play an important and integral part. The study of raised mire peats suggests that faunal stability is a characteristic of these wetland systems, over what appear to be extensive periods of time (up to several millennia), whilst studies of shallow lake ecosystems over recent timescales indicates that faunal instability appears to be more common, usually driven by increasing eutrophication. Drawing upon a series of fossil Coleoptera records spanning several thousand years from Hatfield Moors, south Yorkshire, we reconstruct in some detail the mire’s ontogeny and fluctuations in site hydrology and vegetation cover, illustrating the intimate association between substrate, topography and peat development. A comparison between fossil and modern beetle populations indicates that the faunal characteristics of this mire and its adjacent neighbour, Thorne Moors, become established during the early phases of peat development, including its rare endemics, and that the faunal biodiversity on the sites today is dictated by complex site histories. The over-riding characteristic of these faunas is of stability over several thousand years, which has important implications for the restoration of degraded sites, especially those where refugial areas are limited. In contrast, analyses of fossil Chironomidae from shallow lakes allow researchers to track changes in limnological status and while attempts have been made to reconstruct changes in nutrient levels quantitatively, the chironomids respond indirectly to such changes, typically mediated through complex ecosystem dynamics such as changes in fish and/or macrophyte communities. These changes are illustrated via historic chironomid stratigraphies and diversity indices from a range of shallow lakes located across Britain: Slapton Ley, Frensham Great Pond, Fleet Pond, Kyre Pool and Barnes Loch. These sites have shown varying degrees of eutrophication over recent timescales which tends to be associated with a decline in chironomid diversity. While complex functional processes exist within these ecosystems, our evidence suggests that one of the key drivers in the loss of shallow lake chironomid diversity appears to be the loss of aquatic macrophytes. Overall, while chironomids do show a clear response to altered nutrient regimes, multi-proxy reconstructions are recommended for a clear interpretation of past change. We conclude that if we are to have a better understanding of biota at the ecosystem level we need to know more of the complex interactions between different insect groups as well as with other animal and plant communities. A palaeoecological approach is thus crucial in order to assess the role of insect groups in ecosystem processes, both in the recent past and over long time scales, and is essential for wetland managers and conservation organisations involved in long term management and restoration of wetland systems.

2.4 GHz Class-E power amplifier with transmissionline harmonic terminations

The design procedure, fabrication and measurement of a Class-E power amplifier with excellent second- and third-harmonic suppression levels are presented. A simplified design technique offering compact physical layout is proposed. With a 1.2 mm gate-width GaAs MESFET as a switching device, the amplifier is capable of delivering 19.2 dBm output power at 2.41 GHz, achieves peak PAE of 60% and drain efficiency of 69%, and exhibits 9 dB power gain when operated from a 3 V DC supply voltage. When compared to the classical Class-E two-harmonic termination amplifier, the Class-E amplifier employing three-harmonic terminations has more than 10% higher drain efficiency and 23 dB better third-harmonic suppression level. Experimental results are presented and good agreement with simulation is obtained. Further, to verify the practical implementation in communication systems, the Bluetooth-standard GFSK modulated signal is applied to both two- and three-harmonic amplifiers. The measured RMS FSK deviation error and RMS magnitude error were, for the three-harmonic case, 1.01 kHz and 0.122%, respectively, and, for the two-harmonic case, 1.09 kHz and 0.133%. © 2007 The Institution of Engineering and Technology.

Liquid structure of the ionic liquid, 1-methyl-4-cyanopyridinium bis{(trifluoromethyl)sulfonyl}imide determined from neutron scattering and molecular dynamics simulations

The liquid structure of 1-methyl-4-cyanopyridinium bis {(trifluoromethyl)sulfonyl}imide, a prototypical ionic liquid containing an electron-withdrawing group on the cation, has been investigated at 368 K. Experimental neutron scattering combined with empirical potential structure refinement analysis of the data and classical molecular dynamics simulations have been used to probe the liquid structure in detail. Both techniques generated highly consistent results that provide valuable validation of the force fields and refinement approaches. A significant degree of apparent charge ordering is found in the liquid structure, although the nonspherical shape of the ions results in interpenetration of cations into the first shell of adjacent cations, with much shorter closest contact distances than the averaged center-of-mass cation-cation and cation-anion separations.

Coherent X-ray production via pulse reflection from laser-driven dense electron sheets

A scheme to obtain brilliant x-ray sources by coherent reflection of a counter-propagating pulse from laser-driven dense electron sheets is theoretically and numerically investigated in a self-consistent manner. A radiation pressure acceleration model for the dynamics of the electron sheets blown out from laser-irradiated ultrathin foils is developed and verified by PIC simulations. The first multidimensional and integral demonstration of the scheme by 2D PIC simulations is presented. It is found that the reflected pulse undergoes Doppler-upshift by a factor 4?z2, where ?z = (1- vz2/c2)-1/2 is the effective Lorentz factor of the electron sheet al ong its normal direction. Meanwhile the pulse electric field is intensified by a factor depending on the electron density of the sheet in its moving frame ne/?, where ? is the full Lorentz factor.

Diagnostic based modelling of radio-frequency driven atmospheric pressure plasmas

Diagnostic based modelling (DBM) actively combines complementary advantages of numerical plasma simulations and relatively simple optical emission spectroscopy (OES). DBM is employed to determine absolute atomic oxygen ground state densities in a helium–oxygen radio-frequency driven atmospheric pressure plasma jet. A comparatively simple one-dimensional simulation yields detailed information on electron properties governing the population dynamics of excited states. Important characteristics of the electron dynamics are found to be largely insensitive to details of the chemical composition and to be in very good agreement with space and phase-resolved OES. Benchmarking the time and space resolved simulation allows us to subsequently derive effective excitation rates as the basis for DBM with simple space and time integrated OES. The population dynamics of the upper O 3p 3P (? = 844 nm) atomic oxygen state is governed by direct electron impact excitation, dissociative excitation, radiation losses and collisional induced quenching. Absolute values for atomic oxygen densities are obtained through tracer comparison with the upper Ar 2p1 (? = 750.4 nm) state. The presented results for the atomic oxygen density show excellent quantitative agreement with independent two-photon laser-induced fluorescence measurements.