Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe Moving-Grid Effect in the Harmonic Vibrational Frequency Calculations with Numeric Atom-Centered Orbitals
236
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
When using atom-centered integration grids, the portion of the grid that belongs to a certain atom also moves when this atom is displaced. In the paper, we investigate the moving-grid effect in the calculation of the harmonic vibrational frequencies when using all-electron full-potential numeric atomic-centered orbitals as the basis set. We find that, unlike the first order derivative (i.e., forces), the moving-grid effect plays an essential role for the second order derivatives (i.e., vibrational frequencies). Further analysis reveals that predominantly diagonal force constant terms are affected, which can be bypassed efficiently by invoking translational symmetry. Our approaches have been demonstrated in both finite (molecules) and extended (periodic) systems.
rate research
Read More
The Bethe-Salpeter equation (BSE) based on GW quasiparticle levels is a successful approach for calculating the optical gaps and spectra of solids and also for predicting the neutral excitations of small molecules. We here present an all-electron implementation of the GW+BSE formalism for molecules, using numeric atom-centered orbital (NAO) basis sets. We present benchmarks for low-lying excitation energies for a set of small organic molecules, denoted in the literature as Thiels set. Literature reference data based on Gaussian-type orbitals are reproduced to about one meV precision for the molecular benchmark set, when using the same GW quasiparticle energies and basis sets as the input to the BSE calculations. For valence correlation consistent NAO basis sets, as well as for standard NAO basis sets for ground state density-functional theory with extended augmentation functions, we demonstrate excellent convergence of the predicted low-lying excitations to the complete basis set limit. A simple and affordable augmented NAO basis set denoted tier2+aug2 is recommended as a particularly efficient formulation for production calculations. We finally demonstrate that the same convergence properties also apply to linear-response time-dependent density functional theory within the NAO formalism.
We present an efficient, linear-scaling implementation for building the (screened) Hartree-Fock exchange (HFX) matrix for periodic systems within the framework of numerical atomic orbital (NAO) basis functions. Our implementation is based on the localized resolution of the identity approximation by which two-electron Coulomb repulsion integrals can be obtained by only computing two-center quantities -- a feature that is highly beneficial to NAOs. By exploiting the locality of basis functions and efficient prescreening of the intermediate three- and two-index tensors, one can achieve a linear scaling of the computational cost for building the HFX matrix with respect to the system size. Our implementation is massively parallel, thanks to a MPI/OpenMP hybrid parallelization strategy for distributing the computational load and memory storage. All these factors add together to enable highly efficient hybrid functional calculations for large-scale periodic systems. In this work we describe the key algorithms and implementation details for the HFX build as implemented in the ABACUS code package. The performance and scalability of our implementation with respect to the system size and the number of CPU cores are demonstrated for selected benchmark systems up to 4096 atoms.
Although usually considered as a technique for predicting electron states in dense plasmas, atom-in-jellium calculations can be used to predict the mean displacement of the ion from its equilibrium position in colder matter, as a function of compression and temperature. The Lindemann criterion of a critical displacement for melting can then be employed to predict the melt locus, normalizing for instance to the observed melt temperature or to more direct simulations such as molecular dynamics (MD). This approach reproduces the high pressure melting behavior of Al as calculated using the Lindemann model and thermal vibrations in the solid. Applied to Fe, we find that it reproduces the limited-range melt locus of a multiphase equation of state (EOS) and the results of ab initio MD simulations, and agrees less well with a Lindemann construction using an older EOS. The resulting melt locus lies significantly above the older melt locus for pressures above 1.5,TPa, but is closer to recent ab initio MD results and extrapolations of an analytic fit to them. This study confirms the importance of core freezing in massive exoplanets, predicting that a slightly smaller range of exoplanets than previously assessed would be likely to exhibit dynamo generation of magnetic fields by convection in the liquid portion of the core.
Atom-in-jellium calculations of the Einstein frequency were used to calculate the mean displacement of an ion over a wide range of compression and temperature. Expressed as a fraction of the Wigner-Seitz radius, the displacement is a measure of the asymptotic freedom of the ion at high temperature, and thus of the change in heat capacity from 6 to 3 quadratic degrees of freedom per atom. A functional form for free energy was proposed based on the Maxwell-Boltzmann distribution as a correction to the Debye free energy, with a single free parameter representing the effective density of potential modes to be saturated. This parameter was investigated using molecular dynamics simulations, and found to be ~0.2 per atom. In this way, the ion-thermal contribution can be calculated for a wide-range equation of state (EOS) without requiring a large number of molecular dynamics simulations. Example calculations were performed for carbon, including the sensitivity of key EOS loci to ionic freedom.
Low-loss optical communication requires light sources at 1.5um wavelengths. Experiments showed without much theoretical guidance that InAs/GaAs quantum dots (QDs) may be tuned to such wavelengths by adjusting the In fraction in an InxGa1-xAs strain-reducing capping layer (SRCL). In this work systematic multimillion atom electronic structure calculations qualitatively and quantitatively explain for the first time available experimental data. The NEMO 3-D simulations treat strain in a 15 million atom system and electronic structure in a subset of ~9 million atoms using the experimentally given nominal geometries and without any further parameter adjustments the simulations match the nonlinear behavior of experimental data very closely. With the match to experimental data and the availability of internal model quantities significant insight can be gained through mapping to reduced order models and their relative importance. We can also demonstrate that starting from simple models has in the past led to the wrong conclusions. The critical new insight presented here is that the QD changes its shape. The quantitative simulation agreement with experiment without any material or geometry parameter adjustment in a general atomistic tool leads us to believe that the era of nano Technology Computer Aided Design (nano-TCAD) is approaching. NEMO 3-D will be released on nanoHUB.org where the community can duplicate and expand on the results presented here through interactive simulations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
