Applications of the generalized Langevin equation: Towards a realistic description of the baths

The generalized Langevin equation (GLE) method, as developed previously [L. Stella et al., Phys. Rev. B 89, 134303 (2014)], is used to calculate the dissipative dynamics of systems described at the atomic level. The GLE scheme goes beyond the commonly used bilinear coupling between the central system and the bath, and permits us to have a realistic description of both the dissipative central system and its surrounding bath. We show how to obtain the vibrational properties of a realistic bath and how to convey such properties into an extended Langevin dynamics by the use of the mapping of the bath vibrational properties onto a set of auxiliary variables. Our calculations for a model of a Lennard-Jones solid show that our GLE scheme provides a stable dynamics, with the dissipative/relaxation processes properly described. The total kinetic energy of the central system always thermalizes toward the expected bath temperature, with appropriate fluctuation around the mean value. More importantly, we obtain a velocity distribution for the individual atoms in the central system which follows the expected canonical distribution at the corresponding temperature. This confirms that both our GLE scheme and our mapping procedure onto an extended Langevin dynamics provide the correct thermostat. We also examined the velocity autocorrelation functions and compare our results with more conventional Langevin dynamics.

Electrostatic Solitary Waves in Relativistic Degenerate Electron-Positron-Ion Plasma

The linear and nonlinear properties of ion acoustic excitations propagating in warm dense electron-positron-ion plasma are investigated. Electrons and positrons are assumed relativistic and degenerate, following the Fermi-Dirac statistics, whereas the warm ions are described by a set of classical fluid equations. A linear dispersion relation is derived in the linear approximation. Adopting a reductive perturbation method, the Korteweg-de Vries equation is derived, which admits a localized wave solution in the form of a small-amplitude weakly super-acoustic pulse-shaped soliton. The analysis is extended to account for arbitrary amplitude solitary waves, by deriving a pseudoenergy-balance like equation, involving a Sagdeev-type pseudopotential. It is shown that the two approaches agree exactly in the small-amplitude weakly super-acoustic limit. The range of allowed values of the pulse soliton speed (Mach number), wherein solitary waves may exist, is determined. The effects of the key plasma configuration parameters, namely, the electron relativistic degeneracy parameter, the ion (thermal)-to-the electron (Fermi) temperature ratio, and the positron-to-electron density ratio, on the soliton characteristics and existence domain, are studied in detail. Our results aim at elucidating the characteristics of ion acoustic excitations in relativistic degenerate plasmas, e.g., in dense astrophysical objects, where degenerate electrons and positrons may occur.

GPU Acceleration of Partial Differential Equation Solvers

Differential equations are often directly solvable by analytical means only in their one dimensional version. Partial differential equations are generally not solvable by analytical means in two and three dimensions, with the exception of few special cases. In all other cases, numerical approximation methods need to be utilized. One of the most popular methods is the finite element method. The main areas of focus, here, are the Poisson heat equation and the plate bending equation. The purpose of this paper is to provide a quick walkthrough of the various approaches that the authors followed in pursuit of creating optimal solvers, accelerated with the use of graphical processing units, and comparing them in terms of accuracy and time efficiency with existing or self-made non-accelerated solvers.

Millimeter Wave Amplitude Comparison Monopulse Using a Waveguide Rotman Lens

A 94 GHz waveguide Rotman lens is described which can be used to implement an amplitude comparison monopulse RADAR. In transmit mode, adjacent dual beam ports are excited with equal amplitude and phase to form a sum radiation pattern, and in receive mode, the outputs of the beam port pairs are combined using magic tees to provide a sum and a difference signal which can be used to calculate an angular error estimate for target acquisition and tracking. This approach provides an amplitude comparison monopulse system which can be scanned in azimuth and which has a low component count, with no requirement for phase shift circuitry in the array feed lines, making it suitable for mm-wave frequencies. A 12 input (beam ports), 12 output (array ports) lens is designed using CST Microwave Studio, and the predicted results are presented.

Electron-scale dissipative electrostatic solitons in multi-species plasmas

The linear and nonlinear properties of small-amplitude electron-acoustic solitary waves are investigated via the fluid dynamical approach. A three-component plasma is considered, composed of hot electrons, cold electrons, and ions (considered stationary at the scale of interest). A dissipative (wave damping) effect is assumed due to electron-neutral collisions. The background (hot) electrons are characterized by an energetic (excessively superthermal) population and are thus modeled via a κ-type nonthermal distribution. The linear characteristics of electron-acoustic excitations are discussed, for different values of the plasma parameters (superthermality index κ and cold versus hot electron population concentration β). Large wavelengths (beyond a threshold value) are shown to be overdamped. The reductive perturbation technique is used to derive a dissipative Korteweg de-Vries (KdV) equation for small-amplitude electrostatic potential disturbances. These are expressed by exact solutions in the form of dissipative solitary waves, whose dynamics is investigated analytically and numerically. Our results should be useful in elucidating the behavior of space and experimental plasmas characterized by a coexistence of electron populations at different temperatures, where electron-neutral collisions are of relevance.

Solutions of fourth-order parabolic equation modeling thin film growth

In this paper we study the well-posedness for a fourth-order parabolic equation modeling epitaxial thin film growth. Using Kato's Method [1], [2] and [3] we establish existence, uniqueness and regularity of the solution to the model, in suitable spaces, namelyC⁰([0,T];L^p(Ω)) where  with 1<α<2, n∈N and n≥2. We also show the global existence solution to the nonlinear parabolic equations for small initial data. Our main tools are L_p–L_q-estimates, regularization property of the linear part of e^−tΔ2 and successive approximations. Furthermore, we illustrate the qualitative behavior of the approximate solution through some numerical simulations. The approximate solutions exhibit some favorable absorption properties of the model, which highlight the stabilizing effect of our specific formulation of the source term associated with the upward hopping of atoms. Consequently, the solutions describe well some experimentally observed phenomena, which characterize the growth of thin film such as grain coarsening, island formation and thickness growth.

On the numerical solution of the 2d wave equation with compact fdtd schemes

This paper discusses compact-stencil finite difference time domain (FDTD) schemes for approximating the 2D wave equation in the context of digital audio. Stability, accuracy, and efficiency are investigated and new ways of viewing and interpreting the results are discussed. It is shown that if a tight accuracy constraint is applied, implicit schemes outperform explicit schemes. The paper also discusses the relevance to digital waveguide mesh modelling, and highlights the optimally efficient explicit scheme.