A dielectric continuum molecular dynamics method

We introduce a novel method to simulate hydrated macromolecules with a dielectric continuum representation of the surrounding solvent. In our approach, the interaction between the solvent and the molecular degrees of freedom is described by means of a polarization density free energy functional which is minimum at electrostatic equilibrium. After a pseudospectral expansion of the polarization and a discretization of the functional, we construct the equations of motion for the system based on a Car-Parrinello technique. In the limit of the adiabatic evolution of the polarization field variables, our method provides the solution of the dielectric continuum problem "on the fly," while the molecular coordinates are propagated. In this first study, we show how our dielectric continuum molecular dynamics method can be successfully applied to hydrated biomolecules, with low cost compared to free energy simulations with explicit solvent. To our knowledge, this is the first time that stable and conservative molecular dynamic simulations of solutes can be performed for a dielectric continuum model of the solvent. (C) 2001 American Institute of Physics.

Dynamics of electric fields driving the laser acceleration of multi-MeV protons

The acceleration of multi-MeV protons from the rear surface of thin solid foils irradiated by an intense (similar to 10(18) W/cm(2)) and short (similar to 1.5 ps) laser pulse has been investigated using transverse proton probing. The structure of the electric field driving the expansion of the proton beam has been resolved with high spatial and temporal resolution. The main features of the experimental observations, namely, an initial intense sheath field and a late time field peaking at the beam front, are consistent with the results from particle-in-cell and fluid simulations of thin plasma expansion into a vacuum.

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.

Structure and dynamics of a confined ionic liquid. topics of relevance to dye-sensitized solar cells

The behavior of a model ionic liquid (IL) confined between two flat parallel walls was studied at various interwall distances using computer simulations. The results focus both on structural and dynamical properties. Mass and charge density along the confinement axis reveal a structure of layers parallel to the walls that leads to an oscillatory profile in the electrostatic potential. Orientational correlation functions indicate that cations at the interface orient tilted with respect to the surface and that any other orientational order is lost thereafter. The diffusion coefficients of the ions exhibit a maximum as a function of the confinement distance, a behavior that results from a combination of the structure of the liquid as a whole and a faster molecular motion in the vicinity of the walls. We discuss the relevance of the present results and elaborate on topics that need further attention regarding the effects of ILs in the functioning of IL-based dye-sensitized solar cells.

Instability and dynamics of two nonlinearly coupled laser beams in a plasma

The nonlinear interaction between two laser beams in a plasma is investigated in the weakly nonlinear and relativistic regime. The evolution of the laser beams is governed by two nonlinear Schrodinger equations that are coupled with the slow plasma density response. A nonlinear dispersion relation is derived and used to study the growth rates of the Raman forward and backward scattering instabilities as well of the Brillouin and self-focusing/modulational instabilities. The nonlinear evolution of the instabilities is investigated by means of direct simulations of the time-dependent system of nonlinear equations. (c) 2006 American Institute of Physics.

Intense-field dissociation dynamics of D2+ molecular ions using ultrafast laser pulses

The dynamics of dissociation of pre-ionized D2+ molecules using intense (10^12–10^15 W cm-2), ultrashort (50 fs), infrared (? = 790 nm) laser pulses are examined. Use of an intensity selective scan technique has allowed the deuterium energy spectrum to be measured over a broad range of intensity. It is found that the dominant emission shifts to lower energies as intensity is increased, in good agreement with corresponding wavepacket simulations. The results are consistent with an interpretation in terms of bond softening, which at high intensity (approximately >3 × 10^14 W cm-2) becomes dominated by dissociative ionization. Angular distribution measurements reveal the presence of slow molecular dissociation, an indication that vibrational trapping mechanisms occur in this molecule.

Effector dynamics of rhythmic wrist activity and its implications for (modeling) bimanual coordination

To examine the role of the effector dynamics of the wrist in the production of rhythmic motor activity, we estimated the phase shifts between the EMG and the task-related output for a rhythmic isometric torque production task and an oscillatory movement, and found a substantial difference (45-52degrees) between the two. For both tasks, the relation between EMG and task-related output (torque or displacement) was adequately reproduced with a physiologically motivated musculoskeletal model. The model simulations demonstrated the importance of the contribution of passive structures to the overall dynamics and provided an account for the observed phase shifts in the dynamic task. Additional simulations of the musculoskeletal model with added load suggested that particular changes in the phase relation between EMG and movement may follow largely from the intrinsic muscle dynamics, rather than being the result of adaptations in the neural control of joint stiffness. The implications of these results are discussed in relation to (models of) interlimb coordination in rhythmic tasks. (C) 2004 Elsevier B.V. All rights reserved.

The dynamics of radio-frequency driven atmospheric pressure plasma jets

The complex dynamics of radio-frequency driven atmospheric pressure plasma jets is investigated using various optical diagnostic techniques and numerical simulations. Absolute number densities of ground state atomic oxygen radicals in the plasma effluent are measured by two-photon absorption laser induced fluorescence spectroscopy (TALIF). Spatial profiles are compared with (vacuum) ultra-violet radiation from excited states of atomic oxygen and molecular oxygen, respectively. The excitation and ionization dynamics in the plasma core are dominated by electron impact and observed by space and phase resolved optical emission spectroscopy (PROES). The electron dynamics is governed through the motion of the plasma boundary sheaths in front of the electrodes as illustrated in numerical simulations using a hybrid code based on fluid equations and kinetic treatment of electrons.

Inelastic quantum transport in nanostructures: The self-consistent Born approximation and correlated electron-ion dynamics

A dynamical method for inelastic transport simulations in nanostructures is compared to a steady-state method based on nonequilibrium Green's functions. A simplified form of the dynamical method produces, in the steady state in the weak-coupling limit, effective self-energies analogous to those in the Born approximation due to electron-phonon coupling. The two methods are then compared numerically on a resonant system consisting of a linear trimer weakly embedded between metal electrodes. This system exhibits an enhanced heating at high biases and long phonon equilibration times. Despite the differences in their formulation, the static and dynamical methods capture local current-induced heating and inelastic corrections to the current with good agreement over a wide range of conditions, except in the limit of very high vibrational excitations where differences begin to emerge.

The challenge of revealing and tailoring the dynamics of radio-frequency plasmas

The complex dynamics of ionization and excitation mechanisms in capacitively coupled radio-frequency plasmas is discussed for single- and dual-frequency operations in low-pressure and atmospheric pressure plasmas. Electrons are energized through the dynamics of electric fields in the vicinity of the plasma boundary sheaths. Distinctly different power dissipation mechanisms can either co-exist or initiate mode transitions exhibiting characteristic spatio-temporal ionization structures. Phase resolved optical emission spectroscopy, in combination with adequate modelling of the population dynamics of excited states, and numerical simulations reveal dissipation associated with sheath expansion, sheath collapse, transient electron avalanches and wave–particle interactions. In dual-frequency systems the relative phase between the two frequency components provides additional strategies to tailor the plasma dynamics.

A Geostatistical Investigation into Changing Permeability of Sandstones During Weathering Simulations

The ability to predict the behavior of masonry materials is crucial to conserve building stone. Natural stone, such as sandstone, is not immune from the processes of weathering in the built environment and suffers from decay by granular disintegration, contour scaling, and multiple flaking. Spatial variation of rock properties is a major contributing factor to inconsistent responses to weathering. This has implications for moisture movement and salt input and output and storage, and results in unpredictability in the decay dynamics of masonry materials. This article explores the use of variography and kriging to investigate the spatial interactions between the trigger factors of stone decay, in particular, permeability and its effect on salt penetration. Sandstone blocks were used to represent fresh building stones from a weathering perspective and gave baseline characteristics for the interpretation of subsequent deterioration and decay pathways. Simulated weathering trials involved preloading a sandstone block with salt and subjecting a separate block to 20 cycles of a weathering trial designed to simulate a temperate weathering regime. Geostatistical analysis indicated differences in the spatial variation of permeability of the fresh rock and that subjected to the weathering regimes. Spatial prediction and visualization showed differences in the spatial continuity of permeability in a horizontal and vertical direction through the preloaded block after salt weathering. Continual wetting with salt and alternate heating increased permeability in a vertical direction, enabling the ingress and movement of salt and moisture more effectively through the stone.

Dynamics of Gompertzian tumour growth under environmental fluctuations

We investigate the effect of correlated additive and multiplicative Gaussian white noise oil the Gompertzian growth of tumours. Our results are obtained by Solving numerically the time-dependent Fokker-Planck equation (FPE) associated with the stochastic dynamics. In Our numerical approach we have adopted B-spline functions as a truncated basis to expand the approximated eigenfunctions. The eigenfunctions and eigenvalues obtained using this method are used to derive approximate solutions of the dynamics under Study. We perform simulations to analyze various aspects, of the probability distribution. of the tumour cell populations in the transient- and steady-state regimes. More precisely, we are concerned mainly with the behaviour of the relaxation time (tau) to the steady-state distribution as a function of (i) of the correlation strength (lambda) between the additive noise and the multiplicative noise and (ii) as a function of the multiplicative noise intensity (D) and additive noise intensity (alpha). It is observed that both the correlation strength and the intensities of additive and multiplicative noise, affect the relaxation time.

Plasma surface dynamics and smoothing in the relativistic few-cycle regime

Efficient production of coherent harmonic radiation from solid targets relies critically on the formation of smooth, short density scalelength plasmas. Recent experimental results (Dromey et al 2009 Nat. Phys. 5 146) suggest, however, that the target roughness on the scale of the emitted harmonic wavelength does not result in diffuse reflection-in apparent contradiction to the Rayleigh criterion for coherent reflection. In this paper we show, for the first time, using analytic theory and 2D PIC simulations, that the interaction of relativistically strong laser pulses with corrugated target surfaces results in a highly effective smoothing of the interaction surface and consequently the generation of highly collimated and temporally confined XUV pulses from rough targets, in excellent agreement with experimental observations.

Dynamics of charge-displacement channeling in intense laser-plasma interactions

The dynamics of transient electric fields generated by the interaction of high intensity laser pulses with underdense plasmas has been studied experimentally with the proton projection imaging technique. The formation of a charged channel, the propagation of its front edge and the late electric field evolution have been characterized with high temporal and spatial resolution. Particle-in-cell simulations and an electrostatic, ponderomotive model reproduce the experimental features and trace them back to the ponderomotive expulsion of electrons and the subsequent ion acceleration.