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.

Scattering theory of multilevel atoms interacting with arbitrary radiation fields

We present a generic transfer matrix approach for the description of the interaction of atoms possessing multiple ground state and excited state sublevels with light fields. This model allows us to treat multi-level atoms as classical scatterers in light fields modified by, in principle, arbitrarily complex optical components such as mirrors, resonators, dispersive or dichroic elements, or filters. We verify our formalism for two prototypical sub-Doppler cooling mechanisms and show that it agrees with the standard literature.

Optical Cooling Using the Dipole Force

The term `laser cooling' is applied to the use of optical means to cool the motional energies of either atoms and molecules, or micromirrors. In the literature, these two strands are kept largely separate; both, however suffer from severe limitations. Laser cooling of atoms and molecules largely relies on the internal level structure of the species being cooled. As a result, only a small number of elements and a tiny number of molecules can be cooled this way. In the case of micromirrors, the problem lies in the engineering of micromirrors that need to satisfy a large number of constraints---these include a high mechanical Q-factor, high reflectivity and very good optical quality, weak coupling to the substrate, etc.---in order to enable efficient cooling. During the course of this thesis, I will draw these two sides of laser cooling closer together by means of a single, generically applicable scattering theory that can be used to explain the interaction between light and matter at a very general level. I use this `transfer matrix' formalism to explore the use of the retarded dipole--dipole interaction as a means of both enhancing the efficiency of micromirror cooling systems and rendering the laser cooling of atoms and molecules less species selective. In particular, I identify the `external cavity cooling' mechanism, whereby the use of an optical memory in the form of a resonant element (such as a cavity), outside which the object to be cooled sits, can potentially lead to the construction of fully integrated optomechanical systems and even two-dimensional arrays of translationally cold atoms, molecules or even micromirrors.

Logics for belief functions on MV-algebras

In this paper we present a generalization of belief functions over fuzzy events. In particular we focus on belief functions defined in the algebraic framework of finite MV-algebras of fuzzy sets. We introduce a fuzzy modal logic to formalize reasoning with belief functions on many-valued events. We prove, among other results, that several different notions of belief functions can be characterized in a quite uniform way, just by slightly modifying the complete axiomatization of one of the modal logics involved in the definition of our formalism. © 2012 Elsevier Inc. All rights reserved.

Constraint rule-based programming of norms for electronic institutions

Norms constitute a powerful coordination mechanism among heterogeneous agents. In this paper, we propose a rule language to specify and explicitly manage the normative positions of agents (permissions, prohibitions and obligations), with which distinct deontic notions and their relationships can be captured. Our rule-based formalism includes constraints for more expressiveness and precision and allows to supplement (and implement) electronic institutions with norms. We also show how some normative aspects are given computational interpretation. © 2008 Springer Science+Business Media, LLC.

Monte Carlo radiation hydrodynamics:Methods, tests and application to Type Ia supernova ejecta

In astrophysical systems, radiation-matter interactions are important in transferring energy and momentum between the radiation field and the surrounding material. This coupling often makes it necessary to consider the role of radiation when modelling the dynamics of astrophysical fluids. During the last few years, there have been rapid developments in the use of Monte Carlo methods for numerical radiative transfer simulations. Here, we present an approach to radiation hydrodynamics that is based on coupling Monte Carlo radiative transfer techniques with finite-volume hydrodynamical methods in an operator-split manner. In particular, we adopt an indivisible packet formalism to discretize the radiation field into an ensemble of Monte Carlo packets and employ volume-based estimators to reconstruct the radiation field characteristics. In this paper the numerical tools of this method are presented and their accuracy is verified in a series of test calculations. Finally, as a practical example, we use our approach to study the influence of the radiation-matter coupling on the homologous expansion phase and the bolometric light curve of Type Ia supernova explosions. © 2012 The Authors Monthly Notices of the Royal Astronomical Society © 2012 RAS.

Using sketch-map coordinates to analyze and bias molecular dynamics simulations

When examining complex problems, such as the folding of proteins, coarse grained descriptions of the system drive our investigation and help us to rationalize the results. Oftentimes collective variables (CVs), derived through some chemical intuition about the process of interest, serve this purpose. Because finding these CVs is the most difficult part of any investigation, we recently developed a dimensionality reduction algorithm, sketch-map, that can be used to build a low-dimensional map of a phase space of high-dimensionality. In this paper we discuss how these machine-generated CVs can be used to accelerate the exploration of phase space and to reconstruct free-energy landscapes. To do so, we develop a formalism in which high-dimensional configurations are no longer represented by low-dimensional position vectors. Instead, for each configuration we calculate a probability distribution, which has a domain that encompasses the entirety of the low-dimensional space. To construct a biasing potential, we exploit an analogy with metadynamics and use the trajectory to adaptively construct a repulsive, history-dependent bias from the distributions that correspond to the previously visited configurations. This potential forces the system to explore more of phase space by making it desirable to adopt configurations whose distributions do not overlap with the bias. We apply this algorithm to a small model protein and succeed in reproducing the free-energy surface that we obtain from a parallel tempering calculation.

Molecular dynamics simulation model for the quantitative assessment of tool wear during single point diamond turning of cubic silicon carbide

Silicon carbide (SiC) is a material of great technological interest for engineering applications concerning hostile environments where silicon-based components cannot work (beyond 623 K). Single point diamond turning (SPDT) has remained a superior and viable method to harness process efficiency and freeform shapes on this harder material. However, it is extremely difficult to machine this ceramic consistently in the ductile regime due to sudden and rapid tool wear. It thus becomes non trivial to develop an accurate understanding of tool wear mechanism during SPDT of SiC in order to identify measures to suppress wear to minimize operational cost.

In this paper, molecular dynamics (MD) simulation has been deployed with a realistic analytical bond order potential (ABOP) formalism based potential energy function to understand tool wear mechanism during single point diamond turning of SiC. The most significant result was obtained using the radial distribution function which suggests graphitization of diamond tool during the machining process. This phenomenon occurs due to the abrasive processes between these two ultra hard materials. The abrasive action results in locally high temperature which compounds with the massive cutting forces leading to sp3–sp2 order–disorder transition of diamond tool. This represents the root cause of tool wear during SPDT operation of cubic SiC. Further testing led to the development of a novel method for quantitative assessment of the progression of diamond tool wear from MD simulations.