18 March 2020 1 2K Report

I'm implementing the Occupation Matrix Control scheme (OMC) in VASP DFT+U calculations via direcct input of the Occupation Matrix. This requires an external input file (OCCMATRIX) containing the occupation matrix and the tag OCCEXT = 1 in the INCAR file in order for the system to read the file. I'm studying the occupation of f - orbitals (L=3) in uranium mononitride (UN) therefore a 7x7 diagonal matrix (positions 10 to 16) is required. No off-diagonal elements are required as UN is collinear. There are 21 possible settings for placing 2 electrons in 7f orbitals - in the example below I have the (0000011) configuration so the positions 15 and 16 are set to 1 for electron occupancy. The OCCMATRIX file is as follows:

4 No of atoms to be specified

1 3 1 Atom No (in POSCAR), L, spin(2) or not (1)

spin component 1 Not really important what is here !

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000

spin component 2

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

Blank line between atoms

2 3 1 Atom No (in POSCAR), L, spin(2) or not (1)

spin component 1 Not really important what is here !

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000

spin component 2

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

Blank line between atoms

3 3 2 Atom No (in POSCAR), L, spin(2) or not (1)

spin component 2 Spins switched to give the AFM arrangement

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000

spin component 1

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

Blank line between atoms

4 3 2 Atom No (in POSCAR), L, spin(2) or not (1)

spin component 2 Spins switched to give the AFM arrangement

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000

spin component 1

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

From the job output I always get the error message:

"Error in OCCEXT format. Format is in the occupation matrix documentation"

Despite I followed step-by-step the example in the attached occupation matrix documentation (w/o the right hand notes) I always get the same error message.

I've also tried, for the AFM order on atoms 3 and 4, not to switch the matrix position and leaving matrix 1 before matrix 2 but with no results.

Any help would be much appreciated!!!

Many Thanks.

More Nicola Zagni's questions See All
How to compare means of continuous data with zero-inflation and skew?

I have two groups (unpaired) of plant biomass measurements and want to know whether there is a test that can be used to determine statistical significance where assumptions are not violated. The...

18 February 2021 1,258 4 View

Fluo 3-AM in pancreatic beta-cells?

Dear all, I'm trying to understand if a molecule of my interest (capable of inducing glucose-stimulated insulin secretion) determines an increase in intracellular calcium levels in the presence of...

16 February 2021 453 2 View

Are there commercially available fluorescence-tagged Influenza A virus?

30 November 2020 9,445 1 View

What tests for: 1. differences between continuous variables when grouped with different categorical ones; 2. significance of potential relationships?

Hello. I was wondering if anyone had suggestions in terms of what exploratory analysis to carry out for my dissertation with the following data. I measured N=26 participants pre and post an...

09 August 2020 5,256 5 View

Wavy looking bands In western blots?

Hello everybody, from some time, I keep having waving bands in the lower part of my gel, like the proteins are not running on a flat front. I cast my own gels and run them at 90V at the beginning...

26 May 2020 937 3 View

SURVEY - how people changed attitude toward green spaces & urban forests during the COVID19 quarantine?

https://vub.fra1.qualtrics.com/jfe/form/SV_5bSsVfI1g1VgI7j While shaking our society and turning our way of living upside down, the COVID19 pandemic is also putting the relation between nature...

15 May 2020 2,339 8 View

Glycerol topology/structure for GROMOS87/Berger forcefield?

Hello all, I need the topology /structure file of the glycerol for the GROMOS87 forcefield. Is some repository available online or do someone know how can I derive them by myself (maybe some tool...

12 March 2020 8,858 2 View

How to interpret moderating effect in PLS?

I am running a Partial Least Squadre Path Modelling in R and I have a Moderating effect (explained by one scalar variable) between my exogenous x and my endogenous y. The only variable that...

25 February 2020 8,022 4 View

Recombinant protein for antibody neutralization?

I need to understand if an antibody I'm using is specific or not. The antibody is not validated. In immunoblotting, the antibody detects a prominent band at the right molecular weight and two...

29 April 2019 4,752 4 View

How do I isolate LC nucleus for cell culture?

I need to isolate LC for primary cell culture? How should I do? What is better: use frozen tissue or fresh tissue?

21 January 2019 9,383 3 View

Similar questions and discussions
How to analyze a particular region of my simulation box in GROMACS?

Hi, I am using GROMACS 5.1.4 to simulate water in contact with a surface. I am only interested to calculate some properties of water near the surface. Let's say at a distance of x nm from the...

03 March 2021 9,989 2 View

Does any one know how i can prepare liver tissue with transmission electron microscope with good fixation?

i have problem in preparation of liver tissue with bad cell organelles and how i can adjust osmolarity

03 March 2021 8,890 1 View

How to increase a model accuracy ?

dear community, my model is based feature extraction from non stationary signals using discrete Wavelet Transform and then using statistical features then machine learning classifiers in order to...

03 March 2021 6,994 5 View

I need to image mesoporous silica nanoparticles that are in ~100 nm size and do high resolution TEM to look at mesopores. what grids should I use?

Need to image mesoporous silica nanoparticles using the TEM. Also, need high resolution TEM images to see the mesoporous structure. Kindly suggest what kind of grids to use. Thanks, Shatadru

02 March 2021 1,787 2 View

Arena simulation - how do i model the following model in ARENA (see description) (new to ARENA)?

Hello, I'm fairly new to ARENA and I have the following problem: I have two queues in which, in each queue, a different order type arrives (both having an arrival rate of exponential 10). Order...

02 March 2021 2,431 3 View

Why almost every element would be metallized in sufficiently high pressure and low temperature?

(This statement is from wikipedia, BTW). What is so special about metallic bond that is limiting high-P low-T state of a matter (as long as individual atoms exist, not white dwarf of neutron star...

02 March 2021 3,309 2 View

Shall I skip the measurement model in SmartPLS?

Hello everyone, I had run a mediation model in SmartPLS using questionnaires that were already validated among my statistical population. A professor who edited my paper -though she is not...

01 March 2021 3,052 3 View

How to classify strictly upper triangular matrices by similarity?

Let F be a field. Consider U, the set of n times n strictly upper triangular matrices in F. For X, Y in U, we call them similar if there exists some S, which is non-degenerate and upper...

01 March 2021 2,957 8 View

Is CD44 degraded in lysosome together with hyaluronan after CD44 mediated endocytosis of hyaluronic acid?

Also when RHAMM binds hyaluronic acid, they get internalized, will RHAMM also be degraded? Or both CD44 and RHAMM will be transferred back to the cell membrane? Asking for breast cancer cell line...

01 March 2021 8,169 2 View

Why is the sum of the elements in a row in a Q matrix equal to zero?

When explaining substitution models, the substitutions are expressed as Q matrix. Why is the sum of the elements in a row zero?

28 February 2021 3,864 3 View

Our partners will collect data and use cookies for ad personalization and measurement. Learn how we and our ad partner Google, collect and use data. Agree & close