Show Posts

This section allows you to view all posts made by this member. Note that you can only see posts made in areas you currently have access to.


Messages - Julian Schneider

Pages: 1 ... 8 9 [10] 11
136
General Questions and Answers / Re: Building 6H-SiC Question
« on: September 11, 2015, 14:50 »
BTW. We'll soon have a brand new tutorial on vapor deposition simulations on the SiC(0001) surface..

137
General Questions and Answers / Re: Building 6H-SiC Question
« on: September 11, 2015, 14:47 »
The bond stub is still there, so it looks as if the structure is still intact, only that the one silicon atom is outside the box. Try repeating the cell in A- and B-direction, or using Bulk Tools > Wrap to wrap the atoms inside the box, and check if the bonds still is broken.
I don't see any of these problems when I use the cif-file structure and do the same cleaving operations.


I've been reading,


Ogasawara N, Norimatsu W, Irle S, Kusunoki M. Growth mechanisms and selectivity for graphene or carbon nanotube formation on SiC (0001¯): A density-functional tight-binding molecular dynamics study. Chemical Physics Letters [serial online]. March 18, 2014;595-596:266-271. Available from: ScienceDirect, Ipswich, MA. Accessed September 11, 2015.


They use 5x5x.5 unit cells, I wonder if that is a coincidence?


I guess that's just the choice the authors made to obtain a supercell of a suitable size.

138
We have quite a lot of predefined potentials for SiC. You can simply use one of those. Send your SiC configuration to the ScriptGenerator, add a New Calculator, double-click it to open its settings, and select ATK-Classical, a list of possible potentials will be available. I'd recommend you use one of the Tersoff potentials, e.g. Tersoff_SiC_2005 or Tersoff_SiC_1989.

139
General Questions and Answers / Re: Building 6H-SiC Question
« on: September 11, 2015, 08:56 »

In the diagram the structure looks similar to " 6H SiC Example.png " in my original post. I'm guessing the structure in Jakse was not really meant to represent the physical appearance of the material from one face, but perhaps just a 2-d representation of the stacking sequence?


Yes, you are right. From the picture in the paper it could be both 6H-SiC or 2H-Sic.

Guessing from the side view in your picture " 6H SiC Example.png", it might rather be 2H-SiC (cf. attached picture), but it is not 6H-SiC.

140
General Questions and Answers / Re: Building 6H-SiC Question
« on: September 10, 2015, 13:12 »

I'm confused, doesn't 6H-SiC have 12 atoms in a unit cell? There's only 6 in the coordinates you gave for 6H.
If you import the .cif file, attached in my previous post, into the Builder (via Add > From Files), you will see that it has 12 atoms, although the cif file has only 6 coordinates. The remaining coordinates are deduced from the symmetry.


Also, different sources seem to get different structures for the unit cell, is there more than one?


The two structures in the pictures might be equivalent within periodic boundary conditions, but it is really difficult to tell just from the flat picture.

My knowledge in solid state physics is limited. I'm trying to build 6H-SiC to run tests on the surfaces. I'm not sure if I'm doing this correctly. I start with the Hexagonal SiC in the data base, it consists of 4 atoms. I repeat it and then move atoms around to get the structure of the 6H-SiC unit cell (12 atoms) by going off schematics I've found from various sources. Then I repeat the 6H-SiC structure 5X5X1. I attached images of two of my stash items with the unit cell and the 5X5X1, I also attached two images I looked off of to help me.

Basically my question is...is this right?



If you take the unit cell from the .cif file and cleave the (0001) surface (using the Surface (Cleaver) plugin, and selecting a carbon atom as the topmost terminating atom) then you get the surface as in the picture (see attached figure). Then you can start introducing the modifications.


141
General Questions and Answers / Re: Tremolox Package
« on: September 2, 2015, 15:04 »
You are right, the parameters for EAM and ReaxFF are in stored in separate files. These files are shipped with VNL, so you don't need to generate anything, and you can just use the predefined potentials as they are.

In case you want to modify the predefined EAM or ReaxFF-potentials, you find the files,  if you go into the directory where VNL is installed, and then search for the filename from here. They will be in a subdirectory. You should never directly modify any of these files, as the corresponding predefined potential will be changed or lost without any notice!!!
Instead you should make a copy of this file in the same directory and modify the copied file.
If you then want to employ this new potential file, you simply have to specify the new file name instead of the old one in the EAMsetflPotential constructor (or ReaxFF constructor).
Still, you should be extremely careful when you make any changes in the VNL-installation directory!!

142
General Questions and Answers / Re: Tremolox Package
« on: September 2, 2015, 09:20 »
For most of the potentials, the parameters for the predefined potentials are in the source code itself.
You can, however, easily access and modify the parameters, by setting up a calculator with the desired potential in the Script Generator, and then setting "Script detail" in the "Global IO" settings of the Script Generator to "Show Defaults". If you then send the script to the Editor, a detailed python script with all potential functions and parameters is generated, where the parameters can easily be modified.
See also the tutorial "Adding, combining, and modifying classical potentials" (http://www.quantumwise.com/publications/tutorials/item/848-adding-combining-and-modifying-classical-potentials).

143
General Questions and Answers / Re: Tremolox Package
« on: August 31, 2015, 12:53 »
Generating new potential parameters is a generally a very difficult procedure and it is far from being a black-box procedure. Although, we are working on a tool that will facilitate fitting of new potentials, the current version of VNL does not include such functionality. If there is no parameter set for the desired combination of elements available in the literature, you need to do the fitting yourself.
The common approach is the following: First you generate reference training configurations of your system and calculate representative properties using a higher-level method (e.g. DFT), or use experimental values. These can be lattice constants, elastic constants, energy-displacement curves, atomic forces upon displacement, NEB-energy-curves etc. The reference date should resemble the situation you are finally interested in, as close as possible, so when you are interested in elastic properties, elastic constants should be reproduced by the potential, if you are interested in surface reactions, then bond breaking barriers (e.g. via NEB) should be included in the training set etc. Then you have to choose a potential form which might suit the general material class, e.g. embedded atom for metals and alloys, Tersoff for covalent crystals, or ionic potentials for oxides. Ideally, you can start from combing existing potentials, then you fit some or all potential parameters by hand or using an optimization algorithm (e.g. from numpy/scipy) until the reference properties calculated with the new parameter set match the reference properties sufficiently.
However, you should be aware that fitting potentials requires quite some experience in using classical potentials at the first place.

144
That is nothing to worry about.
ATK-Classical actually comprises 3 calculators: TremoloX, Brenner, and EMT, the last two essentially for historic reasons. When you open the New Calculator widget and select ATK-Classical, VNL checks all of these 3 calculators if they have potentials for the current configuration. If one of these 3 doesn't have parameters for the required elements, then this info message is shown, but hat doesn't mean that there are no parameters available in general.
In your case you can still choose the e.g. the suggested Tersoff potential (from the TremoloXCalculator) without any problem, and simply ignore that message.

145
General Questions and Answers / Re: Building 6H-SiC Question
« on: August 10, 2015, 09:08 »
It looks right, but it's difficult to tell exactly just from the picture.
I have attached a .cif file of the SiC-6H (http://rruff.geo.arizona.edu/AMS/minerals/Moissanite) structure to which you can compare your structure.

146
1) An instead of using finite temperature difference between the phonon occupancies of the reservoirs, the differential ( nL (TL) - nR(TR) ) / (TL - T R) is evaluated as dn(T)/dT at the specified temperature.

2) As it says in the manual M is a diagonal matrix with matrix elements corresponding to the masses of the atoms.

147
(3) The PartialCharges are zero if the classical potential you selected does not use coulomb interactions to model the electrostatics of the system. In this case no partial charges on the atoms are defined and that is what the PartialCharges analysis object tells you.

148
For the Stillinger-Weber and Tersoff-potentials the default values usually work well. You don't need to specify the maximum interaction range in the DynamicalMatrix parameters.

In your case, I'd rather suspect that there is something odd with the structure. Negative frequencies mean that there are soft modes present in the system, either by atoms sitting in a local saddle point.
You could try to calculate the phonon bandstructure (without checking Accoustic sum rule and Symmetrize) and then the location of the negative frequencies might give you some hints. If the negative frequencies occur at gamma, then you might need to increase the maximum interaction range, if the negative frequencies occur elsewhere, then it is likely that your structure is not in a local minimum.

149
Just as an information, if you want to use constraints together with the Image-dependent pair potential or the Halgren-Lipscomb method in ATK2014, you can tag the group of constrained atoms using the tag name 'constrained'. Then the NEB builder will automatically recognize these tagged atoms and keep them fixed.
In ATK2015, the NEB builder will have a proper constraints widget similar to OptimizeGeometry or MolecularDynamics.

150
The maximum interaction range is an optional parameter that can be set in the DynamicalMatrix analysis object. When calculating the dynamical matrix, this parameter specifies the distance around a displaced central atom (i) in which the derivative of the forces on atoms (j)  are taken into account to calculate d^2 U / dr_i dr_j .
If this parameter is not specified, the default values are used (which are the covalent distances multiplied by a fuzz factor of 4.0).  In most cases this works well, but when e.g. a classical potential is used, which has a large cutoff radius (e.g. CoulombDSF (http://www.quantumwise.com/documents/manuals/latest/ReferenceManual/index.html/ref.coulombdsf.html) with a default of 9 Angstrom), then the maximum interaction range parameter has to be increased to be larger than this cutoff radius, i.e. larger than 9 Angstrom, to avoid discontinuities.
You can find out if large cutoff values are used by your potential, if you send the script to the Editor with the "Script details" parameter set to "Show defaults". Then all the details of the classical potential are reflected in the script.
Which exact potential are you using?

Apart from that, negative frequencies might also occur if you are sitting exactly at a local maximum (e.g. due to symmetries in your structure). In this case, small random displacements to the atoms followed by another accurate geometry optimization should lead you to a minimum.

Pages: 1 ... 8 9 [10] 11