A discrete time-dependent method for metastable atoms and molecules in intense fields

The full-dimensional time-dependent Schrodinger equation for the electronic dynamics of single-electron systems in intense external fields is solved directly using a discrete method. Our approach combines the finite-difference and Lagrange mesh methods. The method is applied to calculate the quasienergies and ionization probabilities of atomic and molecular systems in intense static and dynamic electric fields. The gauge invariance and accuracy of the method is established. Applications to multiphoton ionization of positronium, the hydrogen atom and the hydrogen molecular ion are presented. At very high laser intensity, above the saturation threshold, we extend the method using a scaling technique to estimate the quasienergies of metastable states of the hydrogen molecular ion. The results are in good agreement with recent experiments. (C) 2004 American Institute of Physics.

Thin Ta films: growth, stability, and diffusion studied by molecular-dynamics simulations

We report results of classical molecular-dynamics simulations of bcc and beta-Ta thin films. Thermal PVD film growth, surface roughness, argon ion bombardment, phase stability and transformation, vacancy and adatom diffusion, and thermal relaxation kinetics are discussed. Distinct differences between the two structures are observed, including a complex vacancy diffusion mechanism in beta-Ta. Embedded atom method potentials, which were fitted to bcc properties, have been used to model the Ta-Ta interactions. In order to verify the application of these potentials to the more complex beta-Ta structure, we have also performed density functional theory calculations. Results and implications of these calculations are discussed.

Ionization dynamics of laser-driven H-2(+)

We set out aspects of a numerical algorithm used in solving the full-dimensionality time-dependent Schrodinger equation describing the electronic motion of the hydrogen molecular ion driven by an intense, linearly polarized laser pulse aligned along the molecular axis. This algorithm has been implemented within the fixed inter-nuclear separation approximation in a parallel computer code, a brief summary of which is given. Ionization rates are calculated and compared with results from other methods, notably the time-independent Floquet method. Our results compare very favourably with the precise predictions of the Floquet method, although there is some disagreement with other wavepacket calculations. Visualizations of the electron dynamics are also presented in which electron rescattering is observed.

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.

A Van der Pol-Mathieu equation for the dynamics of dust grain charge in dusty plasmas

The chaotic profile of dust grain dynamics associated with dust-acoustic oscillations in a dusty plasma is considered. The collective behaviour of the dust plasma component is described via a multi-fluid model, comprising Boltzmann distributed electrons and ions, as well as an equation of continuity possessing a source term for the dust grains, the dust momentum and Poisson's equations. A Van der Pol–Mathieu-type nonlinear ordinary differential equation for the dust grain density dynamics is derived. The dynamical system is cast into an autonomous form by employing an averaging method. Critical stability boundaries for a particular trivial solution of the governing equation with varying parameters are specified. The equation is analysed to determine the resonance region, and finally numerically solved by using a fourth-order Runge–Kutta method. The presence of chaotic limit cycles is pointed out.

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.

Double-ionization dynamics of laser-driven helium

We report a new method which allows sequential and non-sequential double-ionization rates in laser-driven helium to be distinguished and calculated separately. The method is applied to calculate such rates for two laser pulses, one of 0.236 au frequency and 8.0 × 1015 W cm-2 peak intensity, the other of 1.0 au frequency and also of 8.0 × 1015 W cm-2 peak intensity.

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.

Modelling non-adiabatic processes using correlated electron-ion dynamics

Here we survey the theory and applications of a family of methods (correlated electron-ion dynamics, or CEID) that can be applied to a diverse range of problems involving the non-adiabatic exchange of energy between electrons and nuclei. The simplest method, which is a paradigm for the others, is Ehrenfest Dynamics. This is applied to radiation damage in metals and the evolution of excited states in conjugated polymers. It is unable to reproduce the correct heating of nuclei by current carrying electrons, so we introduce a moment expansion that allows us to restore the spontaneous emission of phonons. Because of the widespread use of Non-Equilibrium Green's Functions for computing electric currents in nanoscale systems, we present a comparison of this formalism with that of CEID with open boundaries. When there is strong coupling between electrons and nuclei, the moment expansion does not converge. We thus conclude with a reworking of the CEID formalism that converges systematically and in a stable manner.

Alternative method to evaluate discharge coefficients. Part 2: case study

The purpose of this article is to apply an alternative method whereby discharge coefficients can be estimated for the flow through a poppet valve at various lifts. Presented is the development of an operational quasi-steady flow rig. An engine cylinder head poppet valve was used as the case study. The requirement to directly measure mass flowrates using a standard conventional steady flow apparatus has been eliminated. Transient mass flowrates, pressures and temperatures of air during an inflow test for a poppet valve at various lifts were measured. Mass flowrates were also calculated from measured cylinder gas pressures and corrected for heat transfer. Using both methods to determine the mass flowrates, isentropic discharge coefficients were calculated and shown to compare within +/- 4.0 per cent of steady flow data. A computational fluid dynamics (CFD) validation of the quasi-steady flow rig is also presented.

An effective method for comparing the turbulence intensity from LDA measurements and CFD predictions within a ship propeller jet

The mean velocity and turbulence intensity are the two main inputs to investigate the ship propeller induced seabed scouring resulting from a vessel is manoeuvring within a port where the underkeel clearances are low. More accurate data including the turbulence intensity is now available by using the laser doppler anemometry (LDA) measurement system and computational fluid dynamics (CFD) approach. Turbulence intensity has a loose definition, which is the velocity fluctuation as the root mean square (RMS) referenced to a mean flow velocity. However, the velocity fluctuation and mean velocity can be the overall value includingx, y and z directions or the value of a single component. LDA and CFD results were obtained from two different acquisition systems (Dantec LDA system and Fluent CFD package) and therefore the outputs cannot be compared directly. An effective method is proposed for comparing the turbulence intensity between the experimental measurements and the computational predictions within a ship propeller jet. The flow patterns of turbulence intensity within a ship propeller jet are presented by using the LDA measurements and CFD results from turbulence models of standard k-e, RNG k-e, realizable k–e, standard k–?, SST k–?and Reynolds stresses.

A robotic apparatus that dictates torque fields around joints without affecting inherent joint dynamics

This manuscript describes how motor behaviour researchers who are not at the same time expert roboticists may implement an experimental apparatus, which has the ability to dictate torque fields around a single joint on one limb or single joints on multiple limbs without otherwise interfering with the inherent dynamics of those joints. Such an apparatus expands the exploratory potential of the researcher wherever experimental distinction of factors may necessitate independent control of torque fields around multiple limbs, or the shaping of torque fields of a given joint independently of its plane of motion, or its directional phase within that plane. The apparatus utilizes torque motors. The challenge with torque motors is that they impose added inertia on limbs and thus attenuate joint dynamics. We eliminated this attenuation by establishing an accurate mathematical model of the robotic device using the Box-Jenkins method, and cancelling out its dynamics by employing the inverse of the model as a compensating controller. A direct measure of the remnant inertial torque as experienced by the hand during a 50 s period of wrist oscillations that increased gradually in frequency from 1.0 to 3.8 Hz confirmed that the removal of the inertial effect of the motor was effectively complete.

First-Order Reversal Curve Probing of Spatially Resolved Polarization Switching Dynamics in Ferroelectric Nanocapacitors

Spatially resolved polarization switching In ferroelectric nanocapacitors was studied on the sub-25 nm scale using the first-order reversal curve (FORC) method. The chosen capacitor geometry allows both high-veracity observation of the domain structure and mapping of polarization switching in a uniform field, synergistically combining microstructural observations and probing of uniform-field polarization responses as relevant to device operation. A classical Kolmogorov-Avrami-Ishibashi model has been adapted to the voltage domain, and the individual switching dynamics of the FORC response curves are well approximated by the adapted model. The comparison with microstructures suggests a strong spatial variability of the switching dynamics inside the nanocapacitors.