Modulated transverse off-plane dust-lattice wave packets in hexagonal two-dimensional dusty plasma crystals

The propagation of nonlinear dust-lattice waves in a two-dimensional hexagonal crystal is investigated. Transverse (off-plane) dust grain oscillatory motion is considered in the form of a backward propagating wave packet whose linear and nonlinear characteristics are investigated. An evolution equation is obtained for the slowly varying amplitude of the first (fundamental) harmonic by making use of a two-dimensional lattice multiple scales technique. An analysis based on the continuum approximation (spatially extended excitations compared to the lattice spacing) shows that wave packets will be modulationally stable and that dark-type envelope solitons (density holes) may occur in the long wavelength region. Evidence is provided of modulational instability and of the occurrence of bright-type envelopes (pulses) at shorter wavelengths. The role of second neighbor interactions is also investigated and is shown to be rather weak in determining the modulational stability region. The effect of dissipation, assumed negligible in the algebra throughout the article, is briefly discussed.

Fully kinetic simulation of ion acoustic and dust-ion acoustic waves

A series of numerical simulations is presented, based on a recurrence-free Vlasov kinetic model using kinetic phase point trajectories. All plasma components are modeled kinetically via a Vlasov evolution equation, then coupled through Poisson’s equation. The dynamics of ion acoustic waves in an electron-ion and in a dusty (electron-ion-dust) plasma configuration are investigated, focusing on wave decay due to Landau damping and, in particular, on the parametric dependence of the damping rate on the dust concentration and on the electron-to-ion temperature ratio. In the absence of dust, the occurrence of damping was observed, as expected, and its dependence to the relative magnitude of the electron vs ion temperature(s) was investigated. When present, the dust component influences the charge balance, enabling dust-ion acoustic waves to survive Landau damping even in the extreme regime where Te???Ti. The Landau damping rate is shown to be minimized for a strong dust concentration or/and for a high value of the electron-to-ion temperature ratio. Our results confirm earlier theoretical considerations and contribute to the interpretation of experimental observations of dust-ion acoustic wave characteristics.

Nonlinear wavepackets in pair-ion and electron-positron-ion plasmas

The occurrence of amplitude-modulated electrostatic and electromagnetic
wavepackets in pair plasmas is investigated. A static additional charged background species is considered, accounting for dust defects or for heavy ion
presence in the background. Relying on a two-fluid description, a nonlinear
Schrodinger type evolution equation is obtained and analyzed, in terms of the
slow dynamics of the wave amplitude. Exact envelope excitations are obtained,
modelling envelope pulses or holes, and their characteristics are discussed.

Ion-acoustic envelope modes in a degenerate relativistic electron-ion plasma

A self-consistent relativistic two-fluid model is proposed for one-dimensional electron-ion plasma dynamics. A multiple scales perturbation technique is employed, leading to an evolution equation for the wave envelope, in the form of a nonlinear Schrödinger type equation (NLSE). The inclusion of relativistic effects is shown to introduce density-dependent factors, not present in the non-relativistic case - in the conditions for modulational instability. The role of relativistic effects on the linear dispersion laws and on envelope soliton solutions of the NLSE is discussed.

Self-consistent non-Markovian theory of a quantum-state evolution for quantum-information processing

We study non-Markovian decoherence phenomena by employing projection-operator formalism when a quantum system (a quantum bit or a register of quantum bits) is coupled to a reservoir. By projecting out the degree of freedom of the reservoir, we derive a non-Markovian master equation for the system, which is reduced to a Lindblad master equation in Markovian limit, and obtain the operator sum representation for the time evolution. It is found that the system is decohered slower in the non- Markovian reservoir than the Markovian because the quantum information of the system is memorized in the non-Markovian reservoir. We discuss the potential importance of non-Markovian reservoirs for quantum-information processing.

Escape times for rigid Brownian rotators in a bistable potential from the time evolution of the Green function and the characteristic time of the probability evolution

The greatest relaxation time for an assembly of three- dimensional rigid rotators in an axially symmetric bistable potential is obtained exactly in terms of continued fractions as a sum of the zero frequency decay functions (averages of the Legendre polynomials) of the system. This is accomplished by studying the entire time evolution of the Green function (transition probability) by expanding the time dependent distribution as a Fourier series and proceeding to the zero frequency limit of the Laplace transform of that distribution. The procedure is entirely analogous to the calculation of the characteristic time of the probability evolution (the integral of the configuration space probability density function with respect to the position co-ordinate) for a particle undergoing translational diffusion in a potential; a concept originally used by Malakhov and Pankratov (Physica A 229 (1996) 109). This procedure allowed them to obtain exact solutions of the Kramers one-dimensional translational escape rate problem for piecewise parabolic potentials. The solution was accomplished by posing the problem in terms of the appropriate Sturm-Liouville equation which could be solved in terms of the parabolic cylinder functions. The method (as applied to rotational problems and posed in terms of recurrence relations for the decay functions, i.e., the Brinkman approach c.f. Blomberg, Physica A 86 (1977) 49, as opposed to the Sturm-Liouville one) demonstrates clearly that the greatest relaxation time unlike the integral relaxation time which is governed by a single decay function (albeit coupled to all the others in non-linear fashion via the underlying recurrence relation) is governed by a sum of decay functions. The method is easily generalized to multidimensional state spaces by matrix continued fraction methods allowing one to treat non-axially symmetric potentials, where the distribution function is governed by two state variables. (C) 2001 Elsevier Science B.V. All rights reserved.

Molecular Hydrogen Emission from Protoplanetary Disks: Effects of X-ray Irradiation and Dust Evolution

Detailed models for the density and temperature profiles of gas and dust in protoplanetary disks are constructed by taking into account X-ray and UV irradiation from a central T Tauri star, as well as dust size growth and settling toward the disk midplane. The spatial and size distributions of dust grains are numerically computed by solving the coagulation equation for settling dust particles, with the result that the mass and total surface area of dust grains per unit volume of the gas in the disks are very small, except at the midplane. The H2 level populations and line emission are calculated using the derived physical structure of the disks. X-ray irradiation is the dominant heating source of the gas in the inner disk and in the surface layer, while the UV heating dominates otherwise. If the central star has strong X-ray and weak UV radiation, the H2 level populations are controlled by X-ray pumping, and the X-rayinduced transition lines could be observable. If the UV irradiation is strong, the level populations are controlled by thermal collisions or UV pumping, depending on the dust properties. As the dust particles evolve in the disks, the gas temperature at the disk surface drops because the grain photoelectric heating becomes less efficient. This makes the level populations change from LTE to non-LTE distributions, which results in changes to the line ratios. Our results suggest that dust evolution in protoplanetary disks could be observable through the H2 line ratios. The emission lines are strong from disks irradiated by strong UV and X-rays and possessing small dust grains; such disks will be good targets in which to observe H2 emission.

Mapping the evolution of optically generated rotational wave packets in a room-temperature ensemble of D2

A coherent superposition of rotational states in D2 has been excited by nonresonant, ultrafast (12 fs), intense (2×1014 W cm-2) 800 nm laser pulses, leading to impulsive dynamic alignment. Field-free evolution of this rotational wave packet has been mapped to high temporal resolution by a time-delayed pulse, initiating rapid double ionization, which is highly sensitive to the angle of orientation of the molecular axis with respect to the polarization direction, . The detailed fractional revivals of the neutral D2 wave packet as a function of and evolution time have been observed and modeled theoretically.

Evolution of linearly polarized electromagnetic pulses in laser plasmas

An analytical and numerical investigation is presented of the behavior of a linearly polarized electromagnetic pulse as it propagates through a plasma. Considering a weakly relativistic regime, the system of one-dimensional fluid-Maxwell equations is reduced to a generalized nonlinear Schrodinger type equation, which is solved numerically using a split step Fourier method. The spatio-temporal evolution of an electromagnetic pulse is investigated. The evolution of the envelope amplitude of density harmonics is also studied. An electromagnetic pulse propagating through the plasma tends to broaden due to dispersion, while the nonlinear frequency shift is observed to slow down the pulse at a speed lower than the group velocity. Such nonlinear effects are more important for higher density plasmas. The pulse broadening factor is calculated numerically, and is shown to be related to the background plasma density. In particular, the broadening effect appears to be stronger for dense plasmas. The relation to existing results on electromagnetic pulses in laser plasmas is discussed. (c) 2008 American Institute of Physics.

Berry phase in open quantum systems: a quantum Langevin equation approach

The evolution of a two level system with a slowly varying Hamiltonian, modeled as a spin 1/2 in a slowly varying magnetic field, and interacting with a quantum environment, modeled as a bath of harmonic oscillators is analyzed using a quantum Langevin approach. This allows to easily obtain the dissipation time and the correction to the Berry phase in the case of an adiabatic cyclic evolution.

1-D constitutive model for evolution of stress-induced R-phase and localized Lüders-like stress-induced martensitic transformation of super-elastic NiTi wires

NiTi alloys have been widely used in the applications for micro-electro-mechanical-systems (MEMS), which often involve some precise and complex motion control. However, when using the NiTi alloys in MEMS application, the main problem to be considered is the degradation of functional property during cycling loading. This also stresses the importance of accurate prediction of the functional behavior of NiTi alloys. In the last two decades, a large number of constitutive models have been proposed to achieve the task. A portion of them focused on the deformation behavior of NiTi alloys under cyclic loading, which is a practical and non-negligible situation. Despite of the scale of modeling studies of the field in NiTi alloys, two experimental observations under uniaxial tension loading have not received proper attentions. First, a deviation from linearity well before the stress-induced martensitic transformation (SIMT) has not been modeled. Recent experiments confirmed that it is caused by the formation of stress-induced R phase. Second, the influence of the well-known localized Lüders-like SIMT on the macroscopic behavior of NiTi alloys, in particular the residual strain during cyclic loading, has not been addressed. In response, we develop a 1-D phenomenological constitutive model for NiTi alloys with two novel features: the formation of stress-induced R phase and the explicit modeling of the localized Lüders-like SIMT. The derived constitutive relations are simple and at the same time sufficient to describe the behavior of NiTi alloys. The accumulation of residual strain caused by R phase under different loading schemes is accurately described by the proposed model. Also, the residual strain caused by irreversible SIMT at different maximum loading strain under cyclic tension loading in individual samples can be explained by and fitted into a single equation in the proposed model. These results show that the proposed model successfully captures the behavior of R phase and the essence of localized SIMT.

Magma Evolution in the Primitive, Intra-oceanic Tonga Arc: Petrogenesis of Basaltic Andesites at Tofua Volcano

Tofua volcano is situated midway along the Tonga oceanic arc and has undergone two phases of ignimbrite-forming activity. The eruptive products are almost entirely basaltic andesites (52 center dot 5-57 wt % SiO2) with the exception of a volumetrically minor pre-caldera dacite. The suite displays a strong tholeiitic trend with K2O <1 wt %. Phenocryst assemblages typically comprise plagioclase + clinopyroxene +/- orthopyroxene with microlites of Ti-magnetite. Olivine (Fo(83-88)) is rare and believed to be dominantly antecrystic. An increase in the extent and frequency of reverse zoning in phenocrysts, sieve-textured plagioclase and the occurrence of antecrystic phases in post-caldera lavas record a shift to dynamic conditions, allowing the interaction of magma batches that were previously distinct. Pyroxene thermobarometry suggests crystallization at 950-1200 degrees C and 0 center dot 8-1 center dot 8 kbar. Volatile measurements of glassy melt inclusions indicate a maximum H2O content of 4 center dot 16 wt % H2O, and CO2-H2O saturation curves indicate that crystallization occurred at two levels, at depths of 4-5 center dot 5 km and 1 center dot 5-2 center dot 5 km. Major and trace element models suggest that the compositions of the majority of the samples represent a differentiation trend whereby the dacite was produced by 65% fractional crystallization of the most primitive basaltic andesite. Trace element models suggest that the sub-arc mantle source is the residuum of depleted Indian mid-ocean ridge basalt mantle (IDMM-1% melt), whereas radiogenic isotope data imply addition of 0 center dot 2% average Tongan sediment melt and a fluid component derived from the subducted altered Pacific oceanic crust. A horizontal array on the U-Th equiline diagram and Ra excesses of up to 500% suggest fluid addition to the mantle wedge within the last few thousand years. Time-integrated (Ra-226/Th-230) vs Sr/Th and Ba/Th fractionation models imply differentiation timescales of up to 4500 years for the dacitic magma compositions at Tofua.

Magma Evolution in the Primitive, Intra-oceanic Tonga Arc: Rapid Petrogenesis of Dacites at Fonualei Volcano

Fonualei is unusual amongst subaerial volcanoes in the Tonga arc because it has erupted dacitic vesicular lavas, tuffs and phreomagmatic deposits for the last 165 years. The total volume of dacite may approach 5 km(3) and overlies basal basaltic andesite and andesite lavas that are constrained to be less than a few millennia in age. All of the products are crystal-poor and formed from relatively low-viscosity magmas inferred to have had temperatures of 1100-1000 degrees C, 2-4 wt % H2O and oxygen fugacities 1-2 log units above the quartz-fayalite-magnetite buffer. Major and trace element data, along with Sr-Nd-Pb and U-Th-Ra isotope data, are used to assess competing models for the origin of the dacites. Positive correlations between Sc and Zr and Sr rule out evolution of the within-dacite compositional array by closed-system crystal fractionation of a single magma batch. An origin by partial melting of lower crustal amphibolites cannot reproduce these data trends or, arguably, any of the dacites either. Instead, we develop a model in which the dacites reflect mixing between two dacitic magmas, each the product of fractional crystallization of basaltic andesite magmas formed by different degrees of partial melting. Mixing was efficient because the two magmas had similar temperatures and viscosities. This is inferred to have occurred at shallow (2-6 km) depths beneath the volcano. U-Th-Ra disequilibria in the basaltic andesite and andesite indicate that the parental magmas had fluids added to their mantle source regions less than 8 kyr ago and that fractionation to the dacitic compositions took less than a few millennia. The 165 year eruption period for the dacites implies that mixing occurred on a similar timescale, possibly during ascent in conduits. The composition of the dacites renders them unsuitable candidates as contributors to average continental crust.

Generalized Langevin equation: An efficient approach to nonequilibrium molecular dynamics of open systems

The generalized Langevin equation (GLE) has been recently suggested to simulate the time evolution of classical solid and molecular systems when considering general nonequilibrium processes. In this approach, a part of the whole system (an open system), which interacts and exchanges energy with its dissipative environment, is studied. Because the GLE is derived by projecting out exactly the harmonic environment, the coupling to it is realistic, while the equations of motion are non-Markovian. Although the GLE formalism has already found promising applications, e. g., in nanotribology and as a powerful thermostat for equilibration in classical molecular dynamics simulations, efficient algorithms to solve the GLE for realistic memory kernels are highly nontrivial, especially if the memory kernels decay nonexponentially. This is due to the fact that one has to generate a colored noise and take account of the memory effects in a consistent manner. In this paper, we present a simple, yet efficient, algorithm for solving the GLE for practical memory kernels and we demonstrate its capability for the exactly solvable case of a harmonic oscillator coupled to a Debye bath.

Constitutive modeling, microstructure evolution, and processing map for a nitride-strengthened heat-resistant steel

A constitutive equation was established to describe the deformation behavior of a nitride-strengthened (NS) steel through isothermal compression simulation test. All the parameters in the constitutive equation including the constant and the activation energy were precisely calculated for the NS steel. The result also showed that from the stress-strain curves, there existed two different linear relationships between critical stress and critical strain in the NS steel due to the augmentation of auxiliary softening effect of the dynamic strain-induced transformation. In the calculation of processing maps, with the change of Zener-Hollomon value, three domains of different levels of workability were found, namely excellent workability region with equiaxed-grain microstructure, good workability region with “stripe” microstructure, and the poor workability region with martensitic-ferritic blend microstructure. With the increase of strain, the poor workability region first expanded, then shrank to barely existing, but appeared again at the strain of 0.6.