Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userFast calculation of the electrostatic potential in ionic crystals by direct summation method
219
0
0.0
(
0
)
Added by
Marie-Bernadette Lepetit
Publication date
2008
fields
Physics
and research's language is
English
Authors
Alain Gelle
Ask ChatGPT about the research
No Arabic abstract
An efficient real space method is derived for the evaluation of the Madelungs potential of ionic crystals. The proposed method is an extension of the Evjens method. It takes advantage of a general analysis for the potential convergence in real space. Indeed, we show that the series convergence is exponential as a function of the number of annulled multipolar momenta in the unit cell. The method proposed in this work reaches such an exponential convergence rate. Its efficiency is comparable to the Ewalds method, however unlike the latter, it uses only simple algebraic functions.
rate research
Read More
We present a fast summation method for lattice sums of the type which arise when solving wave scattering problems with periodic boundary conditions. While there are a variety of effective algorithms in the literature for such calculations, the approach presented here is new and leads to a rigorous analysis of Woods anomalies. These arise when illuminating a grating at specific combinations of the angle of incidence and the frequency of the wave, for which the lattice sums diverge. They were discovered by Wood in 1902 as singularities in the spectral response. The primary tools in our approach are the Euler-Maclaurin formula and a steepest descent argument. The resulting algorithm has super-algebraic convergence and requires only milliseconds of CPU time.
In the molecular dynamics calculations for the free energy of ions and ionic molecules, we often encounter wet charged molecular systems where electrical neutrality condition is broken. This causes a problem in the evaluation of electrostatic interaction under periodic boundary condition. A standard remedy for the problem is to consider a hypothetical homogeneous background charge density to neutralize the total system. Here, we present a new expression for the evaluation of electrostatic interactions for the system including the background charge by fast multipole method (FMM). Further, an efficient scheme to evaluate solute-solvent interaction energy by FMM has been developed to reduce the computation of far-field part. We have calculated hydration free energy of ions, Mg$^{2+}$, Na$^{+}$, and Cl$^{-}$ dissolved in neutral solvent using the new expression. The calculated free energy showed a good agreement with the result using well-established particle mesh Ewald method, demonstrating the validity of the present expression in the framework of FMM. An advantage of the present scheme is in an efficient free energy calculation of a large-scale charged systems (particularly over million particles) based on highly parallel computations.
Previous theoretical studies of calamitic (i.e., rod-like) ionic liquid crystals (ILCs) based on an effective one-species model led to indications of a novel smectic-A phase with a layer spacing being much larger than the length of the mesogenic (i.e., liquid-crystal forming) ions. In order to rule out the possibility that this wide smectic-A phase is merely an artifact caused by the one-species approximation, we investigate an extension which accounts explicitly for cations and anions in ILCs. Our present findings, obtained by grand canonical Monte Carlo simulations, show that the phase transitions between the isotropic and the smectic-A phases of the cation-anion system are in qualitative agreement with the effective one-species model used in the preceding studies. In particular, for ILCs with mesogenes (i.e., liquid-crystal forming species) carrying charged sites at their tips, the wide smectic-A phase forms, at low temperatures and within an intermediate density range, in between the isotropic and a hexagonal crystal phase. We find that in the ordinary smectic-A phase the spatial distribution of the counterions of the mesogens is approximately uniform, whereas in the wide smectic-A phase the small counterions accumulate in between the smectic layers. Due to this phenomenology the wide smectic-A phase could be interesting for applications which hinge on the presence of conductivity channels for mobile ions.
We present a detailed analysis of the modulated-carrier quantum phase gate implemented with Wigner crystals of ions confined in Penning traps. We elaborate on a recent scheme, proposed by two of the authors, to engineer two-body interactions between ions in such crystals. We analyze for the first time the situation in which the cyclotron (w_c) and the crystal rotation (w_r) frequencies do not fulfill the condition w_c=2w_r. It is shown that even in the presence of the magnetic field in the rotating frame the many-body (classical) Hamiltonian describing small oscillations from the ion equilibrium positions can be recast in canonical form. As a consequence, we are able to demonstrate that fast and robust two-qubit gates are achievable within the current experimental limitations. Moreover, we describe a realization of the state-dependent sign-changing dipole forces needed to realize the investigated quantum computing scheme.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
