Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAutomatic numerical evaluation of vacancy-mediated transport for arbitrary crystals: Onsager coefficients in the dilute limit using a Green function approach
301
0
0.0
(
0
)
Authors
Dallas R. Trinkle
Ask ChatGPT about the research
No Arabic abstract
A general solution for vacancy-mediated diffusion in the dilute-vacancy/dilute-solute limit for arbitrary crystal structures is derived from the master equation. A general numerical approach to the vacancy lattice Green function reduces to the sum of a few analytic functions and numerical integration of a smooth function over the Brillouin zone for arbitrary crystals. The Dyson equation solves for the Green function in the presence of a solute with arbitrary but finite interaction range to compute the transport coefficients accurately, efficiently and automatically, including cases with very large differences in solute-vacancy exchange rates. The methodology takes advantage of the space group symmetry of a crystal to reduce the complexity of the matrix inversion in the Dyson equation. An open-source implementation of the algorithm is available, and numerical results are presented for the convergence of the integration error of the bare vacancy Green function, and tracer correlation factors for a variety of crystals including wurtzite (hexagonal diamond) and garnet.
rate research
Read More
A local-orbital based ab initio approach to obtain the Green function for large heterogeneous systems is developed. First a Green function formalism is introduced based on exact diagonalization. Then the self energy is constructed from an incremental scheme, rendering the procedure feasible, while at the same time physical insight into different local correlation contributions is obtained. Subsequently the Green function is used in the frame of the Landauer theory and the wide band approximation to calculate the electronic transmission coefficient across molecular junctions. The theory is applied to meta- and para-ditholbenzene linked to gold electrodes and various correlation contributions are analyzed.
Efficient computation of lattice defect geometries such as point defects, dislocations, disconnections, grain boundaries, interfaces and free surfaces requires accurate coupling of displacements near the defect to the long-range elastic strain. Flexible boundary condition methods embedded a defect in infinite harmonic bulk through the lattice Green function. We demonstrate an efficient and accurate calculation of the lattice Green function from the force-constant matrix for general crystals with an arbitrary basis by extending a method for Bravais lattices. New terms appear due to the presence of optical modes and the possible loss of inversion symmetry. By separately treating poles and discontinuities in reciprocal space, numerical accuracy is controlled at all distances. We compute the lattice Green function for a two-dimensional model with broken symmetry to elucidate the role of different coupling terms. The algorithm is generally applicable in two and three dimensions, to crystals with arbitrary number of atoms in the unit cell, symmetry, and interactions.
The objective of this study is to assess the impact of a dense-phase treatment on the hydrodynamic description of granular, binary mixtures relative to a previous dilute-phase treatment. Two theories were considered for this purpose. The first, proposed by Garzo and Dufty (GD) [Phys. Fluids {bf 14}, 146 (2002)], is based on the Boltzmann equation which does not incorporate finite-volume effects, thereby limiting its use to dilute flows. The second, proposed by Garzo, Hrenya and Dufty (GHD) [Phys. Rev. E {bf 76}, 31303 and 031304 (2007)], is derived from the Enskog equation which does account for finite-volume effects; accordingly this theory can be applied to moderately dense systems as well. To demonstrate the significance of the dense-phase treatment relative to its dilute counterpart, the ratio of dense (GHD) to dilute (GD) predictions of all relevant transport coefficients and equations of state are plotted over a range of physical parameters (volume fraction, coefficients of restitution, material density ratio, diameter ratio, and mixture composition). These plots reveal the deviation between the two treatments, which can become quite large ($>$100%) even at moderate values of the physical parameters. Such information will be useful when choosing which theory is most applicable to a given situation, since the dilute theory offers relative simplicity and the dense theory offers improved accuracy. It is also important to note that several corrections to original GHD expressions are presented here in the form of a complete, self-contained set of relevant equations.
The disordering of an initially phase segregated system of finite size, induced by the presence of highly mobile vacancies, is shown to exhibit dynamic scaling in its late stages. A set of characteristic exponents is introduced and computed analytically, in excellent agreement with Monte Carlo data. In particular, the characteristic time scale, controlling the crossover between increasing disorder and saturation, is found to depend on the exponent scaling the number of vacancies in the sample.
In a recent publication (Phys. Rev E 77, 056705 (2008)),we have presented the stochastic Green function (SGF) algorithm, which has the properties of being general and easy to apply to any lattice Hamiltonian of the form H=V-T, where V is diagonal in the chosen occupation number basis and T has only positive matrix elements. We propose here a modified version of the update scheme that keeps the simplicity and generality of the original SGF algorithm, and enhances significantly its efficiency.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
