Experimentally, determining the absolute valence and conductance band positions can be difficult. While DFT often underestimates band gaps, experimentally the band gap is a far more accessible measurement that could be used to find the band position of the valence or conductance band if one is k...
I want to understand the electronic band structure diagram of the following image, corresponding to $\text{MoS}_2$ (TMD): I read about DFT (density functional theory). DFT is based on solving the Schrodinger equation for a set of atoms. Through the Born Oppenheimer approach, it is possible to de...
2D materials (such as graphene) are mostly exfoliated from their 3D bulk counterpart. I am wondering how we can calculate work function for a 2D material? Normally electrostatic potential in the out-of-plane direction is required to obtain E$_{\rm{vac}}$, which is later used in the definition of ...
The Raman spectrum of a material shows certain 'active' modes for which we get experimental signatures. Phonon calculations are mostly used to calculate vibrational modes. How can we determine the nature of these modes from a phonon calculation (e.g., if the mode is an in-plane or out-of-plane or...
From Statistical Mechanics obtain an expression for the heat capacity of a solid is given by1: \begin{equation} \tag{1} {C_V}\left( T \right) = k\int\limits_0^\infty {\frac{{{{\left( {uv} \right)}^2}{e^{uv}}}}{{{{\left( {{e^{uv}} - 1} \right)}^2}}}g\left( v \right)dv} \end{equation} with $u=h/kT...
The famous Hohenberg-Kohn theorems say that there is a one-to-one mapping between the many-body Hamiltonian, $\mathcal{H}$, of a solid and its ground-state electron density $\rho(\mathbf{r})$. As far as I understand, this also means that all the properties of the ground-state wavefunction are enc...
A fractal is, accordingly Oxford English Dictionary: A curve or geometric figure, each part of which has the same statistical character as the whole. Fractals are useful in modeling structures (such as eroded coastlines or snowflakes) in which similar patterns recur at progressively smaller scal...
I have two main modeling research lines: one related to structural biology (including rational drug design, de novo design, polymorphism, etc.) and the other one related to condensed matter (nanostructures, crystals, etc.). Any time I submit a manuscript in the biological area, either the editors...
I am trying to understand the symmetry elements of space group number 194 (P$6_3$/mmc), which is hexagonal and has 24 symmetry operations. In the table of symmetry operations it says that it has three glide planes, which I marked in the figure (see blue rectangles). However, more than three glide...
I am trying to synthesize MoSe$_2$ crystals in a furnace with two temperature zones. The problem is that the oven is not very wide, so the temperature of the hottest zone influences the temperature of the cold zone. This is visible in the Excel chart I obtained from the practice, in which I place...
Studying magnetic systems, two frequently used approximations are the Heisenberg and Ising models (a discussion about these approximations can be read here): \begin{equation} \tag{Heisenberg} \hat{H}_H=-\sum_{\langle i j\rangle}J\hat{S}_i\hat{S}_j \end{equation} \begin{equation} \tag{Ising} \hat{...
This question is related to other one here in the MatterModelingSE: Is it possible to calculate/estimate the value of the J parameter to be used in the Heisengerg/Ising hamiltonians? Here J is the exchange interaction parameters between two nearest-neighbor spins using in Heisenberg and Ising h...
Does anyone know if magnetism and topological insulating behavior coexist in a material? If yes, can someone refer to a recent work?
I have recently been working with non-collinear antiferromagnetic systems. I understand that, in order to obtain spin density, one must obtain the electron density in one configuration of spin and subtract the electron density in the other configuration. I am unsure as to how to go about this in ...
I work with magnetic materials and I haven't found a way to automatically plot structures with their corresponding magnetic moments. For this specific project I am using VASP. I am trying to plot something like this:
For many applications with heavy metals, pseudopotentials can be used to include some amount of relativistic effects. But for what sort of systems does it become necessary to actually use a relativistic method (e.g Dirac Hartree Fock)? Is this there active research into material modeling with a r...
Since the effect of spin-orbit coupling plays an important role in many transition metal complexes, what are the common methods to incorporate the effect of spin-orbit coupling?
In Jensen's Introduction to Computational Chemistry it says that the total non-relativistic Hamiltonian operator, transformed to the center of mass system, can be written, in atomic units, as $$ \hat{H}=\hat{T}_N+\hat{H}_e+\hat{H}_{mp} $$ where $\hat{T}_N$ is the kinetic energy of the nuclei, $\...
I'm new to DFT. Does the DFT take into account the spin-orbit interaction? On the one hand, this is a relativistic effect, perhaps DFT does not take it into account. However, on the other hand, I read that DFT is an extremely accurate method and the error of this method lies only in the inaccura...
I have used Vanderbilt USPP for generating a pseudopotential that does not include the spin-orbit term. But I want to have the spin-orbit coupling term included in the pseudopotential, to do research on a van der Waal (VDW) hetero-structure of 2D materials. How do I use Vanderbilt USPP to gen...
I am working on DFT code having plane wave along xy axis and bspline in z direction. For calculating the properties of TMD materials, spin orbit coupling must be included but i did not find fully relativistic uspp pseudo-potentials of atoms that can be generated from vanderbilt code. Then I get i...
I have been assigned a project on the generation of a force-field of Al (aluminum) for a LAMMPS potential. I have no idea how to proceed and what to do next. Can anyone please let me know how should I proceed further?
I would like to know what are the different types of pseudopotentials, the pro and cons, and what properties can/cannot be calculated with them?
If you do a bibliography search about the codes used in matter modeling, you will find a large list. I played with a few of them (both free and commercial) and one of my findings is that the properties they can calculate are different, even being important properties. I am aware that they used di...
I've read in some texts that GW-BSE quasiparticle calculations have only been implemented on Norm-conserving Pseudopotentials (and not Ultra-soft or PAW PPs). Is there a conceptual reason for this? I can't seem to locate a reliable source but here it is mentioned that the YAMBO code (used for GW-...
I am currently trying to figure out how to compute band structures for my system, using the hybrid functional HSE06. I'm doing this on Quantum Espresso. As I understand, there are a handful of ways to do this: 1) Generate maximally localized wannier functions (MLWFs) with the HSE functional turn...
Regarding the calculation of binding energy of molecular absorption on graphene, for which this has been discussed before, the process of calculating $E_{AB}$ is fairly straight forward, i.e., it is the single point energy of the relaxed complex AB. However, the calculation of $E_A$ and $E_B$ can...
I want to calculate the density of states of GaAs conduction band (preferably each valley like Gamma, L and X). Is there any software available where I can easily generate DOS? Has the DOS of GaAs been studied previously?
I am interested in studying material heat transport using Density Functional Theory. Searching, I found a work1 about using ab initio equilibrium molecular dynamics (AIMD) as implemented in QuantumEspresso (QE) software. I already used/tested QE for other simpler calculations and I found it very ...
I am facing problems in calculating the phonon band structure for bilayer 2D systems. While the band structure mostly shows positive frequencies for the 1 layer system, going to 2 layers makes some frequencies imaginary. Moreover, the choice of van der Waals interaction also comes into account fo...
Does it make any sense to calculate edge states for topological semimetals while they don't have any global bands?
Can someone explain how to verify the presence of Rashba effect from the band structure calculations?
The image below refers to a phenomenon that occurs in TMDs (transition metal dichalcogenides) that allowed the development of valleytronics. Why are there separate bands of different colors in this image (some with red on top and blue on bottom and others with blue on top and red on bottom)? Fig...
I want to draw the energy ($E$) diagrams for a simple cubic cell of parameter $a$, where each atom provides two electrons for the almost free electron levels for planes [100], [110] and [111]. I calculate $k$ for the first Brillouin zone: $k_{\text{ first Brillouin zone}}=\frac{\pi}{a}=\frac{3,14...
Many problems in Computational Physics need the use of random number generators. When studying magnetic materials using the Heisenberg/Ising hamiltonians (related questions about them can be seen here, here and here), one of the must used tool is the Metropolis method1. The Metropolis method can ...
Previous questions (here, here and here) were about modeling magnetic homogeneous systems using Heisenberg and/or Ising hamiltonians: \begin{equation} \tag{Heisenberg} \hat{H}_H=-\sum_{\langle i j\rangle}J\hat{S}_i\hat{S}_j \end{equation} \begin{equation} \tag{Ising} \hat{H}_I=-\sum_{\langle ij\r...
I find it surprisingly difficult to find researches/papers on systematic "many-body interaction" extensions of the Ising model. Can somebody tell me a good review/article etc on this matter that goes through possible extensions of the Ising model in a somewhat systematic way? The usual $s_is_j$ i...
In the paper that introduced "Self-learning MC" (an ML-inspired MC technique, as I understand) PhysRevB.95.041101 (https://journals.aps.org/prb/abstract/10.1103/PhysRevB.95.041101) the authors consider a many-body Ising model as an example to show the efficiency of their algorithm. The model look...
I would like to demonstrate that the two-band Hamiltonian is given by the following formula: \begin{equation} \hat{H}_0=at(\tau k_x\hat{\sigma}_x+k_y\hat{\sigma}_y)+\frac{\Delta}{2}\hat{\sigma}_z-\lambda\tau\frac{\hat{\sigma}_z-1}{2}\hat{s}_z \tag{1} \end{equation} where $\tau=\pm1$ is the valley...
I have been unable to figure out what the Kohn anomaly is. What exactly is it, and why is it important?
I came across the concept of using phonons to establish a material's dynamical stability, based on whether or not imaginary frequencies are present in its phonon band structure. What I am struggling with is how to determine the phase transitions based on a phonon band structure with such imaginar...
I read in David Vanderbilt's book named "Berry Phases in Electronic Structure Theory - Electric Polarization, Orbital Magnetization and Topological Insulators" the definition of Berry curvature: "Berry curvature $\Omega(\mathbf{\lambda})$ is simply defined as the Berry phase per unit area in ($\l...
I am trying to understand the Berry phase through the evolution of a system that evolves adiabatically. Schrodinger's equation is: \begin{equation} H(\lambda)|n(\lambda)\rangle=E_n|n(\lambda)\rangle \tag{1} \end{equation} where $n$ labels the eigenstates. If $\lambda$ doesn't change with time, th...
The wave function for the adiabatic approach is as follows: \begin{equation} |\psi(t)\rangle=e^{i\phi(\lambda(t))}e^{-i\gamma(t)}|n(t)\rangle \tag{1} \end{equation} where $e^{i\phi(\lambda(t))}$ is the geometric phase or Berry phase and $e^{-i\gamma(t)}$ is dynamic phase. Berry phase has the foll...
What does the quantum anomalous Hall phenomenon consist of? I am studying the transition metal dichalcogenides (TMDs) and I have seen webinars and articles that said that these materials exhibited the anomalous quantum Hall effect related to the curvature of Berry, which changes the speed of elec...
I want to represent parallel transport as in the following figure, but next to this figure I would like to place a diagram with the two arrows (of the initial state and of the final state that are present at the top of the spherical surface) and the angle between them, to see better and put more ...
With the advent of more computational power than ever in the recent years, interest in in silico design of interesting compounds has grown as well. I am wondering about the state of the art for the case of designing novel materials: What methods are commonly used (DFT, semi-empirical, more / le...
How can someone properly decide on an appropriate level of theory MM, QM... What are the checklists throughout the experiment (in-silico)? How to decide a proper control case/group? For example from personal experience, it happened to me to work without a control (ground truth). I have realized t...
High-throughput density functional theory (DFT) calculations are used to screen for new materials and conduct fundamental research in materials science and materials innovation. It involves computations on tens of thousands of compounds, and such a scale demands unique calculation and data manage...
The Materials Genome Initiative held a 2017 workshop, that led to this published report: de Pablo et al., "New frontiers for the materials initiative," npj Computational Materials 5, 41 (2019). In part, this report lists various successes of the initiative, due to complementary theoretical, compu...
There are several pages where you can find scripts/simulations to generate the first Brillouin zone for square and hexagonal 2D lattices. I wonder if there is a tool to generate the Brillouin for other 2D lattices like the tiles presented here and here. PS: I am aware that not all the tiles can...
When we calculate the band structure of certain material, we only have to calculate the value along the high symmetry point which enclose the Irreducible Brillouin Zone. Why the information lie in the IBZ is enough for us to know the information of entire BZ ? How to find the IBZ given a reciproc...
For the choice of k-mesh we have few options like: Monkhorst-Pack, Chadi-Cohen, Gamma-centered mesh, etc. What are the different types of k-mesh and what are their pros and cons?
A Brillouin zone is defined as a Wigner~Secitz primitive cell in the reciprocal lattice [1,2]. The construction of the first Brillouin zone for two different 2D lattices are shown below: The construction procedure for 3D lattices is basically the same as for 2D lattices. In the image below, ...
Since the beginning of the century, Wien2K has proven to be a very powerful player in computational condensed matter and materials physics. Wien2K is an all-electron periodic DFT code based on the Augmented Planewave + Local Orbital [APW+lo] method, and also widely known as Full-Potential - Linea...
Is there any connection between charge density wave and band inversion? Or is there any system in which band inversion follows the same mechanism as the CDW? I am just trying to find a connection between them and I will really appreciate it if someone can help me.
I used VASP to do the band structure calculation on Graphene. Originally I would expect that the Dirac cone can be observed in the K point. If we watched qualitatively, we can see this behavior. But when I zoom in the region where the conduction band touches with the valences band, a bandgap can...
I'm trying to simulate a graphene flake with its edge C atoms saturated by H atoms, in a temperature ramp from 300 K to 1600 K, using LAMMPS and the AIREBO potential proposed by Stuart et al. in 2000. The problem is when the temperature is close to 500 K the C-H bonds start breaking and at 1600 K...
I saw a lot of references say that the Dirac cone in graphene is protected by inversion and time reversal symmetries. How can one understand this statement? How can one show explicitly that the gapless state will be destroyed if the Hamiltonian violates either one of them? I know that inversion s...
I am working with graphene growth on metal substrates. But I need some sample codes to know if I am doing it in the right way because running multiple simulations in LAMMPS takes a lot of time. Can anyone help me in that? It will be really helpful. Any sample file regarding graphene growth will d...
Prior to the discovery of graphene, the Mermin-Wagner theorem was used to argue that purely two-dimensional materials would not be stable as two-dimensional order would show logarithmic divergences at long range and fluctuations would spontaneously melt the crystal. Yet, graphene exists and it ha...
Most periodic density functional theory (DFT) codes use plane-wave basis sets in conjunction with three-dimensional periodic boundary conditions. In contrast, for molecular systems of finite size, Gaussian basis sets are often used. The former are quite efficient for periodic systems, but more ex...
I'm interested in looking at current flow across a nanoscale junction, specifically a pair of electrodes linked by a molecular bridge. How is this sort of problem typically approached? I'm vaguely familiar with the idea of using Green's functions for this purpose, but don't know much in detail b...
Taking into account the fact that the theory of quantum gravity does not exist and the QED calculations are not possible for most realistic chemical systems, what levels of accuracy can we expect from a theoretical calculation on simple (small) materials? Examples that come to mind are: Simplest...
A few years back, I performed a few calculations of metallic slabs using VASP and the Atomic Simulation Environment (ASE). I'm by no means an expert, but I noticed that the calculations were very CPU intensive e.g. ~24 cpu-hours (in parallel) to get a single energy point and gradient calculation....
In order to do DFT calculations of metallic alloys, the start point is a supercell, whose atoms are changed to match the desired stoichiometry. Gold and Silver, for example, both have FCC structure. In Quantum-ESPRESSO, we can both tell to the code the Bravais lattice (ibrav parameter) and specif...
I tried to recover a .CIF file (Crystallographic Information File) with the crystal structure for the compound discussed in a paper by Reuter et. al.[1], strontium hydroxide octahydrate - $\ce{Sr(OH)2·8H2O}$ - from Cambridge Crystallographic Data Centre (CCDC), but couldn't find it. The paper its...
I am studying a material which was first discovered in the early 1960s. It was reported then that it crystallized in the hexagonal $P6_3/mmc$ space group. Recently the material was synthesized again and found that it posses $P6_3mc$ symmetry. My own DFT calculations and another published work fin...
I am using VASP package to model a material that has fractional site occupancies. An appropriate supercell (cif) was generated using supercell package and converted to POSCAR using VESTA. Then the structural relaxation was performed. Now I have the relaxed structure (CONTCAR) and I want to analyz...
Transition metal dichalcogenides (TMDs) are composed of three atomic planes and often two atomic species: a metal and two chalcogens. The honeycomb, hexagonal lattice has three fold symmetry and can permit mirror plane symmetry and / or inversion symmetry. [16] In the macroscopic bulk crystal, or...
I am aware that there exists a flowchart about how to determine the point group of a molecule with a given geometry, such as given in the diagram here: https://chem.libretexts.org/Bookshelves/Physical_and_Theoretical_Chemistry_Textbook_Maps/Book%3A_Symmetry_(Vallance)/03._Symmetry_Classification_...
Taking a bulk WS2 (bilayer) structure and relaxing it with ISIF=3 and IVDW=10 (DFT-D2 method) using the VASP code. Then opening this structure with VESTA software. Are the shown a and c lattice constants are the optimized ones? or there is something else that must be done to extract them?
The method that I am following in creating a hexagonal lattice POSCAR file is by downloading it from https://materialsproject.org/ website. Now, I want to learn how to create it manually starting from given information in an article. Could someone of you give me a detailed method to do that? Plea...
I have seen that in many articles of 2D materials calculations that VASP is the most used one among DFT codes. Could please explain to me why?
I am seeking some models or references on how to couple or what to take into consideration while coupling first-order phase transition (for magnetic systems) and second-order phase transition regarding order and disorder in these systems. I am attaching a link for one such system https://www.natu...
The Jahn-Teller effect arises from the symmetry breaking (distortion) of the local environment for the crystal field of 3d transitional metal ions.The distortion result in energy saving of electronic energy from the splitting of the degeneracy energy levels of octahedron / tetrahedral environment...
next day → last day (1384 days later) »