Tight-binding simulation of current-carrying nanostructures

The tight-binding (TB) approach to the modelling of electrical conduction in small structures is introduced. Different equivalent forms of the TB expression for the electrical current in a nanoscale junction are derived. The use of the formalism to calculate the current density and local potential is illustrated by model examples. A first-principles time-dependent TB formalism for calculating current-induced forces and the dynamical response of atoms is presented. An earlier expression for current-induced forces under steady-state conditions is generalized beyond local charge neutrality and beyond orthogonal TB. Future directions in the modelling of power dissipation and local heating in nanoscale conductors are discussed.

Current-induced forces in atomic-scale conductors

We present a self-consistent tight-binding formalism to calculate the forces on individual atoms due to the flow of electrical current in atomic-scale conductors. Simultaneously with the forces, the method yields the local current density and the local potential in the presence of current flow, allowing a direct comparison between these quantities. The method is applicable to structures of arbitrary atomic geometry and can be used to model current-induced mechanical effects in realistic nanoscale junctions and wires. The formalism is implemented within a simple Is tight-binding model and is applied to two model structures; atomic chains and a nanoscale wire containing a vacancy.

Multi-zone TAP-reactors theory and application - IV. Ideal and non-ideal boundary conditions

We study the influence of non-ideal boundary and initial conditions (BIC) of a temporal analysis of products (TAP) reactor model on the data (observed exit flux) analysis. The general theory of multi-response state-defining experiments for a multi-zone TAP reactor is extended and applied to model several alternative boundary and initial conditions proposed in the literature. The method used is based on the Laplace transform and the transfer matrix formalism for multi-response experiments. Two non-idealities are studied: (1) the inlet pulse not being narrow enough (gas pulse not entering the reactor in Dirac delta function shape) and (2) the outlet non-ideality due to imperfect vacuum. The effect of these non-idealities is analyzed to the first and second order of approximation. The corresponding corrections were obtained and discussed in detail. It was found that they are negligible. Therefore, the model with ideal boundary conditions is proven to be completely adequate to the description and interpretation of transport-reaction data obtained with TAP-2 reactors.

Linear quadratic and tumour control probability modelling in external beam radiotherapy

The standard linear-quadratic (LQ) survival model for external beam radiotherapy is reviewed with particular emphasis on studying how different schedules of radiation treatment planning may be affected by different tumour repopulation kinetics. The LQ model is further examined in the context of tumour control probability (TCP) models. The application of the Zaider and Minerbo non-Poissonian TCP model incorporating the effect of cellular repopulation is reviewed. In particular the recent development of a cell cycle model within the original Zaider and Minerbo TCP formalism is highlighted. Application of this TCP cell-cycle model in clinical treatment plans is explored and analysed.

The Combination of Multiple Classifiers Using an Evidential Reasoning Approach

In many domains when we have several competing classifiers available we want to synthesize them or some of them to get a more accurate classifier by a combination function. In this paper we propose a ‘class-indifferent’ method for combining classifier decisions represented by evidential structures called triplet and quartet, using Dempster's rule of combination. This method is unique in that it distinguishes important elements from the trivial ones in representing classifier decisions, makes use of more information than others in calculating the support for class labels and provides a practical way to apply the theoretically appealing Dempster–Shafer theory of evidence to the problem of ensemble learning. We present a formalism for modelling classifier decisions as triplet mass functions and we establish a range of formulae for combining these mass functions in order to arrive at a consensus decision. In addition we carry out a comparative study with the alternatives of simplet and dichotomous structure and also compare two combination methods, Dempster's rule and majority voting, over the UCI benchmark data, to demonstrate the advantage our approach offers. (A continuation of the work in this area that was published in IEEE Trans on KDE, and conferences)

Electrochemical Oxidation of Hydrogen Sulfide at Platinum Electrodes in Room Temperature Ionic Liquids: Evidence for Significant Accumulation of H2S at the Pt/1-Butyl-3-methylimidazolium Trifluoromethylsulfonate Interface

Electrochemical oxidation of hydrogen sulfide gas (H2S) has been studied at a platinum microelectrode (10 mu m diameter) in five room temperature ionic liquids (RTILs): [C(4)mim][OTf], [C(4)dmim][NTf2], [C(4)mim][PF6],. [C(6)mim][FAP], and [P-14,P-6,P-6,P-6][FAP] (where [C-n mim](+) = 1-alkyl-3-methylimidazolium, [C(n)dmim](+) = 1-alkyl-2,3-dimethylimidazolium, [P-14,P-6,P-6,P-6](+) = tris(p-hexyl)-tetradecylphosphonium, [OTf](-) = trifluoromethlysulfonate, [NTf2](-) = bis(trifluoromethylsulfonyl)imide, [PF6](-) = hexafluorophosphate, and [FAP](-) = trifluorotris(pentafluoroethyl)phosphate). In four of the RTILs ([C(4)dmim][NTf2], [C(4)mim][PF6], [C(6)mim][FAP], and [P-14,P-6,P-6,P-6][FAP]), no clear oxidative signal was observed. In [C(4)mim][OTf], a chemically irreversible oxidation peak was observed on the oxidative sweep with no signal seen on the reverse scan. The oxidative signal showed an adsorptive stripping peak type followed by near steady-state limiting current behavior. Potential step chronoamperometry was carried out on the reductive wave, giving a diffusion coefficient and solubility of 1.6 x 10(-11) m(2) s(-1) and 7 mM, respectively (at 25 degrees C). Using these data, we modeled the oxidation signal kinetically, assuming adsorption preceded oxidation and that adsorption was approximately Langmuirian. The oxidation step was described by an electrochemically fully irreversible Tafel law/Butler-Volmer formalism. Modeling indicated a substantial buildup of H2S in the double layer in excess of the coverage that would be expected for a monolayer of chemisorbed H2S, reflecting high solubility of the gas in [C(4)mim][OTf] and possible attractive interactions with the [OTf](-) anions accumulated at the electrode at potentials positive of the potential of zero charge. Solute enrichment of the double layer in the solution adjacent to the electrode appears a novel feature of RTIL electrochemistry.

Passing quantum correlations to qubits using any two-mode state

We draw an explicit connection between the statistical properties of an entangled two-mode continuous variable (CV) resource and the amount of entanglement that can be dynamically transferred to a pair of noninteracting two-level systems. More specifically, we rigorously reformulate entanglement-transfer process by making use of covariance matrix formalism. When the resource state is Gaussian, our method makes the approach to the transfer of quantum correlations much more flexible than in previously considered schemes and allows the straightforward inclusion of the effects of noise affecting the CV system. Moreover, the proposed method reveals that the use of de-Gaussified two-mode states is almost never advantageous for transferring entanglement with respect to the full Gaussian picture, despite the entanglement in the non-Gaussian resource can be much larger than in its Gaussian counterpart. We can thus conclude that the entanglement-transfer map overthrows the

Many-body theory calculations of positron binding to negative ions

A many-body theory approach developed by the authors [Phys. Rev. A 70, 032720 (2004)] is applied to positron bound states and annihilation rates in atomic systems. Within the formalism, full account of virtual positronium (Ps) formation is made by summing the electron-positron ladder diagram series, thus enabling the theory to include all important many-body correlation effects in the positron problem. Numerical calculations have been performed for positron bound states with the hydrogen and halogen negative ions, also known as Ps hydride and Ps halides. The Ps binding energies of 1.118, 2.718, 2.245, 1.873 and 1.393 eV and annihilation rates of 2.544, 2.482, 1.984, 1.913 and 1.809 ns^{-1}, have been obtained for PsH, PsF, PsCl, PsBr and PsI, respectively.

Modulational instability criteria for two-component Bose-Einstein condensates

The stability of colliding Bose-Einstein condensates is investigated. A set of coupled Gross-Pitaevskii equations is thus considered, and analyzed via a perturbative approach. No assumption is made on the signs ( or magnitudes) of the relevant parameters like the scattering lengths and the coupling coefficients. The formalism is therefore valid for asymmetric as well as symmetric coupled condensate wave states. A new set of explicit criteria is derived and analyzed. An extended instability region, in addition to an enhanced instability growth rate, is predicted for unstable two component bosons, as compared to the individual ( uncoupled) state.

Acoustic solitary waves in dusty and/or multi-ion plasmas with cold, adiabatic, and hot constituents

Large nonlinear acoustic waves are discussed in a four-component plasma, made up of two superhot isothermal species, and two species with lower thermal velocities, being, respectively, adiabatic and cold. First a model is considered in which the isothermal species are electrons and ions, while the cooler species are positive and/or negative dust. Using a Sagdeev pseudopotential formalism, large dust-acoustic structures have been studied in a systematic way, to delimit the compositional parameter space in which they can be found, without restrictions on the charges and masses of the dust species and their charge signs. Solitary waves can only occur for nonlinear structure velocities smaller than the adiabatic dust thermal velocity, leading to a novel dust-acoustic-like mode based on the interplay between the two dust species. If the cold and adiabatic dust are oppositely charged, only solitary waves exist, having the polarity of the cold dust, their parameter range being limited by infinite compression of the cold dust. However, when the charges of the cold and adiabatic species have the same sign, solitary structures are limited for increasing Mach numbers successively by infinite cold dust compression, by encountering the adiabatic dust sonic point, and by the occurrence of double layers. The latter have, for smaller Mach numbers, the same polarity as the charged dust, but switch at the high Mach number end to the opposite polarity. Typical Sagdeev pseudopotentials and solitary wave profiles have been presented. Finally, the analysis has nowhere used the assumption that the dust would be much more massive than the ions and hence, one or both dust species can easily be replaced by positive and/or negative ions and the conclusions will apply to that plasma model equally well. This would cover a number of different scenarios, such as, for example, very hot electrons and ions, together with a mix of adiabatic ions and dust (of either polarity) or a very hot electron-positron mix, together with a two-ion mix or together with adiabatic ions and cold dust (both of either charge sign), to name but some of the possible plasma compositions.

Relativistic laser pulse compression in plasmas with a linear axial density gradient

The self-compression of a relativistic Gaussian laser pulse propagating in a non-uniform plasma is investigated. A linear density inhomogeneity (density ramp) is assumed in the axial direction. The nonlinear Schrodinger equation is first solved within a one-dimensional geometry by using the paraxial formalism to demonstrate the occurrence of longitudinal pulse compression and the associated increase in intensity. Both longitudinal and transverse self-compression in plasma is examined for a finite extent Gaussian laser pulse. A pair of appropriate trial functions, for the beam width parameter (in space) and the pulse width parameter (in time) are defined and the corresponding equations of space and time evolution are derived. A numerical investigation shows that inhomogeneity in the plasma can further boost the compression mechanism and localize the pulse intensity, in comparison with a homogeneous plasma. A 100 fs pulse is compressed in an inhomogeneous plasma medium by more than ten times. Our findings indicate the possibility for the generation of particularly intense and short pulses, with relevance to the future development of tabletop high-power ultrashort laser pulse based particle acceleration devices and associated high harmonic generation. An extension of the model is proposed to investigate relativistic laser pulse compression in magnetized plasmas.

Electron beam-plasma interaction and ion-acoustic solitary waves in plasmas with a superthermal electron component

The study of non-Maxwellian plasmas is crucial to the understanding of space and astrophysical plasma dynamics. In this paper, we investigate the existence of arbitrary amplitude ion-acoustic solitary waves in an unmagnetized plasma consisting of ions and excess superthermal electrons (modelled by a kappa-type distribution), which is penetrated by an electron beam. A kappa (kappa-) type distribution is assumed for the background electrons. A (Sagdeev-type) pseudopotential formalism is employed to derive an energy-balance like equation. The range of allowed values of the soliton speed (Mach number), wherein solitary waves may exist, is determined. The Mach number range (allowed soliton speed values) becomes narrower under the combined effect of the electron beam and of the superthermal electrons, and may even be reduced to nil (predicting no solitary wave existence) for high enough beam density and low enough kappa (significant superthermality). For fixed values of all other parameters (Mach number, electron beam-to-ion density ratio and electron beam velocity), both soliton amplitude and (electric potential perturbation) profile steepness increase as kappa decreases. The combined occurrence of small-amplitude negative potential structures and larger amplitude positive ones is pointed out, while the dependence of either type on the plasma parameters is investigated.

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.

Shape-induced phase transition of domain patterns in ferroelectric platelets

Using a combination of experimental and computational techniques, changes in the domain structures seen infreestanding single-crystal platelets of BaTiO3 have been described in terms of a second-order phase transition.The transition is driven by the change in the length-to-width ratio of the platelet sidewalls and results in a symmetrybreaking of a complex, quadrant domain pattern. The phenomenon can be described by a Landau formalism inwhich (1) the order parameter is not the polarization but rather is the degree to which the domain pattern becomesoff-centered, and (2) the shape anisotropy of the platelet substitutes for temperature in the conventional Landauexpansion as the controlling thermodynamic variable. Bistability, in terms of the direction in which the domainpattern moves off center, coupled with the spontaneous macroscopic polarization and toroidal moment that resultfrom this off-centering, prompt the possibility of a new form of memory storage.

Excitation of the 3p^4(4s,3d,4p) Ar^+ states during Ar photoionization:Intensity, alignment, and orientation

Ar photoionization is studied using the R-matrix formalism with emphasis on the simultaneous excitation of the residual A^r+ ion. Cross sections have been obtained for excitation of the 3p^4(3d,4s,4p) states. A comparison with experiments having a resolution of 70 meV shows reasonable agreement for the position and shape of resonance structures. The relative magnitude of the resonances proves to be more elusive. The partial cross section for excitation of the 3p^4(3Pe)4p(2P_3/2^o) and (2D_3/2^o) levels is treated in more detail. A comparison of LS-coupling calculations with high-resolution experimental results shows good agreement for both the excitation cross sections and the polarization of the fluorescence. We also predict the orientation for both levels. We demonstrate that the polarization of the fluorescence originating from the (2D_3/2^o) level can be employed to study spin-orbit effects in Ar photoionization.