Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAn $O(N)$ $Ab~initio$ Calculation Scheme for Large-Scale Moir{e} Structures
76
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present a two-step method specifically tailored for band structure calculation of the small-angle moir{e}-pattern materials which contain tens of thousands of atoms in a unit cell. In the first step, the self-consistent field calculation for ground state is performed with $O(N)$ Krylov subspace method implemented in OpenMX. Secondly, the crystal momentum dependent Bloch Hamiltonian and overlap matrix are constructed from the results obtained in the first step and only a small number of eigenvalues near the Fermi energy are solved with shift-invert and Lanczos techniques. By systematically tuning two key parameters, the cutoff radius for electron hopping interaction and the dimension of Krylov subspace, we obtained the band structures for both rigid and corrugated twisted bilayer graphene structures at the first magic angle ($theta=1.08^circ$) and other three larger ones with satisfied accuracy on affordable costs. The band structures are in good agreement with those from tight binding models, continuum models, plane-wave pseudo-potential based $ab~initio$ calculations, and the experimental observations. This efficient two-step method is to play a crucial role in other twisted two-dimensional materials, where the band structures are much more complex than graphene and the effective model is hard to be constructed.
rate research
Read More
Using an Environmentally Friendly Renormalization Group we derive an ab initio universal scaling form for the equation of state for the O(N) model, y=f(x), that exhibits all required analyticity properties in the limits $xto 0$, $xtoinfty$ and $xto -1$. Unlike current methodologies based on a phenomenological scaling ansatz the scaling function is derived solely from the underlying Landau-Ginzburg-Wilson Hamiltonian and depends only on the three Wilson functions $gamma_lambda$, $gamma_phi$ and $gamma_{phi^2}$ which exhibit a non-trivial crossover between the Wilson-Fisher fixed point and the strong coupling fixed point associated with the Goldstone modes on the coexistence curve. We give explicit results for N=2, 3 and 4 to one-loop order and compare with known results.
The s manifold energy levels for phosphorus donors in silicon are important input parameters for the design and modelling of electronic devices on the nanoscale. In this paper we calculate these energy levels from first principles using density functional theory. The wavefunction of the donor electrons ground state is found to have a form that is similar to an atomic s orbital, with an effective Bohr radius of 1.8 nm. The corresponding binding energy of this state is found to be 41 meV, which is in good agreement with the currently accepted value of 45.59 meV. We also calculate the energies of the excited 1s(T) and 1s(E) states, finding them to be 32 and 31 meV respectively. These results constitute the first ab initio confirmation of the s manifold energy levels for phosphorus donors in silicon.
The two-body knock-out reaction 4He(e,ed)d is calculated at various momentum transfers. The full four-nucleon dynamics is taken into account microscopically both in the initial and the final states. As NN interaction the central MT-I/III potential is used. The calculation shows a strong reduction of the coincidence cross section due to the final state interaction. Nonetheless the theoretical results exhibit a considerable overestimation of the experimental cross section at lower momentum transfer. Comparisons with other, less complete, calculations suggest that consideration of a more realistic ground state might not be sufficient for a good agreement with experiment, rather a more realistic final state interaction could play an essential role.
The broken inversion symmetry at the surface of a metallic film (or, more generally, at the interface between a metallic film and a different metallic or insulating material) greatly amplifies the influence of the spin-orbit interaction on the surface properties. The best known manifestation of this effect is the momentum-dependent splitting of the surface state energies (Rashba effect). Here we show that the same interaction also generates a spin-polarization of the bulk states when an electric current is driven through the bulk of the film. For a semi-infinite jellium model, which is representative of metals with a closed Fermi surface, we prove as a theorem that, regardless of the shape of the confinement potential, the induced surface spin density at each surface is given by ${bf S} =-gamma hbar {bf hat z}times {bf j}$, where ${bf j}$ is the particle current density in the bulk, ${bf hat z}$ the unit vector normal to the surface, and $gamma=frac{hbar}{4mc^2}$ contains only fundamental constants. For a general metallic solid $gamma$ becomes a material-specific parameter that controls the strength of the interfacial spin-orbit coupling. Our theorem, combined with an {it ab initio} calculation of the spin polarization of the current-carrying film, enables a determination of $gamma$, which should be useful in modeling the spin-dependent scattering of quasiparticles at the interface.
We report on the computational characteristics of ab initio nuclear structure calculations in a symmetry-adapted no-core shell model (SA-NCSM) framework. We examine the computational complexity of the current implementation of the SA-NCSM approach, dubbed LSU3shell, by analyzing ab initio results for 6Li and 12C in large harmonic oscillator model spaces and SU(3)-selected subspaces. We demonstrate LSU3shells strong-scaling properties achieved with highly-parallel methods for computing the many-body matrix elements. Results compare favorably with complete model space calculations and significant memory savings are achieved in physically important applications. In particular, a well-chosen symmetry-adapted basis affords memory savings in calculations of states with a fixed total angular momentum in large model spaces while exactly preserving translational invariance.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
