Combined R-matrix eigenstate basis set and finite-difference propagation method for the time-dependent Schrodinger equation: The one-electron case

In this work we present the theoretical framework for the solution of the time-dependent Schrödinger equation (TDSE) of atomic and molecular systems under strong electromagnetic fields with the configuration space of the electron’s coordinates separated over two regions; that is, regions I and II. In region I the solution of the TDSE is obtained by an R-matrix basis set representation of the time-dependent wave function. In region II a grid representation of the wave function is considered and propagation in space and time is obtained through the finite-difference method. With this, a combination of basis set and grid methods is put forward for tackling multiregion time-dependent problems. In both regions, a high-order explicit scheme is employed for the time propagation. While, in a purely hydrogenic system no approximation is involved due to this separation, in multielectron systems the validity and the usefulness of the present method relies on the basic assumption of R-matrix theory, namely, that beyond a certain distance (encompassing region I) a single ejected electron is distinguishable from the other electrons of the multielectron system and evolves there (region II) effectively as a one-electron system. The method is developed in detail for single active electron systems and applied to the exemplar case of the hydrogen atom in an intense laser field.

Electron-acoustic plasma waves: Oblique modulation and envelope solitons

Theoretical and numerical studies are presented of the amplitude modulation of electron-acoustic waves (EAWs) propagating in space plasmas whose constituents are inertial cold electrons, Boltzmann distributed hot electrons, and stationary ions. Perturbations oblique to the carrier EAW propagation direction have been considered. The stability analysis, based on a nonlinear Schrodinger equation, reveals that the EAW may become unstable; the stability criteria depend on the angle theta between the modulation and propagation directions. Different types of localized EA excitations are shown to exist.

Fokker-Planck equation for a test-particle weakly coupled to a magnetized plasma

We consider the derivation of a kinetic equation for a charged test particle weakly interacting with an electrostatic plasma in thermal equilibrium, subject to a uniform external magnetic field. The Liouville equation leads to a generalized master equation to second order in the `weak' interaction; a Fokker-Planck-type equation then follows as a `Markovian' approximation. It is shown that such an equation does not preserve the positivity of the distribution function f(x,v;t). By applying techniques developed in the theory of open systems, a correct Fokker-Planck equation is derived. Explicit expressions for the diffusion and drift coefficients, depending on the magnetic field, are obtained.

A comparison of nonstaggered compact FDTD schemes for the 3D wave equation

This paper aims at providing a better insight into the 3D approximations of the wave equation using compact finite-difference time-domain (FDTD) schemes in the context of room acoustic simulations. A general family of 3D compact explicit and implicit schemes based on a nonstaggered rectilinear grid is analyzed in terms of stability, numerical error, and accuracy. Various special cases are compared and the most accurate explicit and implicit schemes are identified. Further considerations presented in the paper include the direct relationship with other numerical approaches found in the literature on room acoustic modeling such as the 3D digital waveguide mesh and Yee's staggered grid technique.

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.

CRITICAL PROCESSES, LANGEVIN EQUATION AND UNIVERSALITY

In nature there are ubiquitous systems that can naturally approach critical states, The Langevin equation in the discrete version can be used to describe a class of critical processes, which are characterized by power-law behaviors and scaling relations. As an example, we present a simple model for a clinical thermometer, whose reading cannot fall even when its temperature decreases. The fibers bundle model and the spring-block model are also shown to belong to such a class.

Note on the single-shock solutions of the Korteweg-de Vries-Burgers equation

The well-known shock solutions of the Korteweg-de Vries-Burgers equation are revisited, together with their limitations in the context of plasma (astro)physical applications. Although available in the literature for a long time, it seems to have been forgotten in recent papers that such shocks are monotonic and unique, for a given plasma configuration, and cannot show oscillatory or bell-shaped features. This uniqueness is contrasted to solitary wave solutions of the two parent equations (Korteweg-de Vries and Burgers), which form a family of curves parameterized by the excess velocity over the linear phase speed.

Current-induced atomic dynamics, instabilities, and Raman signals: Quasiclassical Langevin equation approach

We derive and employ a semiclassical Langevin equation obtained from path integrals to describe the ionic dynamics of a molecular junction in the presence of electrical current. The electronic environment serves as an effective nonequilibrium bath. The bath results in random forces describing Joule heating, current-induced forces including the nonconservative wind force, dissipative frictional forces, and an effective Lorentz-type force due to the Berry phase of the nonequilibrium electrons. Using a generic two-level molecular model, we highlight the importance of both current-induced forces and Joule heating for the stability of the system. We compare the impact of the different forces, and the wide-band approximation for the electronic structure on our result. We examine the current-induced instabilities (excitation of runaway "waterwheel" modes) and investigate the signature of these in the Raman signals.

Real-time approach to the optical properties of solids and nanostructures: Time-dependent Bethe-Salpeter equation

Many-body effects are known to play a crucial role in the electronic and optical properties of solids and nanostructures. Nevertheless, the majority of theoretical and numerical approaches able to capture the influence of Coulomb correlations are restricted to the linear response regime. In this work, we introduce an approach based on a real-time solution of the electronic dynamics. The proposed approach reduces to the well-known Bethe-Salpeter equation in the linear limit regime and it makes it possible, at the same time, to investigate correlation effects in nonlinear phenomena. We show the flexibility and numerical stability of the proposed approach by calculating the dielectric constants and the effect of a strong pulse excitation in bulk h-BN.

Parameter estimation and equation formulation in Business Dynamics

System Dynamics enables modelling and simulation of highly non-linear feedback systems to predict future system behaviour. Parameter estimation and equation formulation are techniques in System Dynamics, used to retrieve the values of parameters or the equations for ?ows and/or variables. These techniques are crucial for the annotations and thereafter the simulation. This paper critically examines existing and well established approaches in parameter estimation and equation formulation along with their limitations, identifying performance gaps as well as providing directions for potential future research.

Ion-acoustic supersolitons in plasmas with two-temperature electrons: Boltzmann and kappa distributions

Acoustic supersolitons arise when a plasma model is able to support three consecutive local extrema of the Sagdeev pseudopotential between the undisturbed conditions and an accessible root. This leads to a characteristic electric field signature, where a simple bipolar shape is enriched by subsidiary maxima. Large-amplitude nonlinear acoustic modes are investigated, using a pseudopotential approach, for plasmas containing two-temperature electrons having Boltzmann or kappa distributions, in the presence of cold fluid ions. The existence domains for positive supersolitons are derived in a methodological way, both for structure velocities and amplitudes, in terms of plasma compositional parameters. In addition, typical pseudopotentials, soliton, and electric field profiles have been given to illustrate that positive supersolitons can be found in the whole range of electron distributions from Maxwellian to a very hard nonthermal spectrum in kappa. However, it is found that the parameter ranges that support supersolitons vary significantly over the wide range of kappa considered. VC 2013 AIP Publishing LLC. [http://dx.doi.org/10.1063/1.4818888]

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.