Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSinc-based method for an efficient solution in the direct space of quantum wave equations with periodic boundary conditions
131
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The solution of differential problems, and in particular of quantum wave equations, can in general be performed both in the direct and in the reciprocal space. However, to achieve the same accuracy, direct-space finite-difference approaches usually involve handling larger algebraic problems with respect to the approaches based on the Fourier transform in reciprocal space. This is the result of the errors that direct-space discretization formulas introduce into the treatment of derivatives. Here, we propose an approach, relying on a set of sinc-based functions, that allows us to achieve an exact representation of the derivatives in the direct space and that is equivalent to the solution in the reciprocal space. We apply this method to the numerical solution of the Dirac equation in an armchair graphene nanoribbon with a potential varying only in the transverse direction.
rate research
Read More
We describe a method and its implementation for calculating electronic structure and electron transport without approximating the structure using periodic super-cells. This effectively removes spurious periodic images and interference effects. Our method is based on already established methods readily available in the non-equilibrium Green function formalism and allows for non-equilibrium transport. We present examples of a N defect in graphene, finite voltage bias transport in a point-contact to graphene, and a graphene-nanoribbon junction. This method is less costly, in terms of CPU-hours, than the super-cell approximation.
We developed a micromagnetic method for modeling magnetic systems with periodic boundary conditions along an arbitrary number of dimensions. The main feature is an adaptation of the Ewald summation technique for evaluation of long-range dipolar interactions. The method was applied to investigate the hysteresis process in hard-soft magnetic nanocomposites with various geometries. The dependence of the results on different micromagnetic parameters was studied. We found that for layered structures with an out-of-plane hard phase easy axis the hysteretic properties are very sensitive to the strength of the interlayer exchange coupling, as long as the spontaneous magnetization for the hard phase is significantly smaller than for the soft phase. The origin of this behavior was discussed. Additionally, we investigated the soft phase size optimizing the energy product of hard-soft nanocomposites.
We address periodic-image errors arising from the use of periodic boundary conditions to describe systems that do not exhibit full three-dimensional periodicity. The difference between the periodic potential, as straightforwardly obtained from a Fourier transform, and the potential satisfying any other boundary conditions can be characterized analytically. In light of this observation, we present an efficient real-space method to correct periodic-image errors, based on a multigrid solver for the potential difference, and demonstrate that exponential convergence of the energy with respect to cell size can be achieved in practical calculations. Additionally, we derive rapidly convergent expansions for determining the Madelung constants of point-charge assemblies in one, two, and three dimensions.
The existence of an inertial manifold for the 3D Cahn-Hilliard equation with periodic boundary conditions is verified using the proper extension of the so-called spatial averaging principle introduced by G. Sell and J. Mallet-Paret. Moreover, the extra regularity of this manifold is also obtained.
An approximate method based on adiabatic time dependent density functional theory (TDDFT) is presented, that allows for the description of the electron dynamics in nanoscale junctions under arbitrary time dependent external potentials. In this scheme, the density matrix of the device region is propagated according to the Liouville-von Neumann equation. The semi-infinite leads give rise to dissipative terms in the equation of motion which are calculated from first principles in the wide band limit. In contrast to earlier {em ab-initio} implementations of this formalism, the Hamiltonian is here approximated by a density expansion in the spirit of the density functional based tight-binding (DFTB) method without introducing empirical parameters. Results are presented for two prototypical molecular devices and compared to calculations at the full TDDFT level. The issue of non-existence of a steady state under certain conditions is also briefly touched on.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
