SCALEE - Revision history

Karsai at 14:05, 16 October 2024

2024-10-16T14:05:03Z

← Older revision		Revision as of 14:05, 16 October 2024
Line 13:		Line 13:

	----		----
	[[Category:INCAR tag~~]][[Category:Molecular dynamics~~]][[Category:Advanced molecular-dynamics sampling~~]][[Category:Thermodynamic integration~~]]		[[Category:INCAR tag]][[Category:Advanced molecular-dynamics sampling]]

Karsai at 14:04, 16 October 2024

2024-10-16T14:04:16Z

@@ Line 4: / Line 4: @@
 ----
 The tag {{TAG|SCALEE}} sets the coupling parameter <math>\lambda</math> and hence controls the Hamiltonian of the calculation.
-By default {{TAG|SCALEE}}=1 and the scaling of the energies and forces via the coupling constant is internally skipped in the code. To enable the scaling {{TAG|SCALEE}}<math>\ne</math>1 has to be specified. A VASP calculation outputs the integrand for a given coupling constant at every molecular dynamics step. How to choose the ensemble size and carry out the integration is described in the main text and especially in the supplementary information of reference {{cite|dorner:PRL:2018}}.
+By default {{TAG|SCALEE}}=1 and the scaling of the energies and forces via the coupling constant is internally skipped in the code. To enable the scaling {{TAG|SCALEE}}<math>\ne</math>1 has to be specified.
-Two possible options are available for the reference system:
+More information using this tag is given [[Thermodynamic integration calculations|here]].
 == Related tags and articles ==
 {{TAG|VCAIMAGES}}, {{TAG|IMAGES}}, {{TAG|NCORE IN IMAGE1}}, {{TAG|PHON_NSTRUCT}}, {{TAG|IBRION}}
 ----
 [[Category:INCAR tag]][[Category:Molecular dynamics]][[Category:Advanced molecular-dynamics sampling]][[Category:Thermodynamic integration]]

Karsai at 09:36, 15 October 2024

2024-10-15T09:36:49Z

← Older revision		Revision as of 09:36, 15 October 2024
Line 42:		Line 42:
	----		----

	[[Category:INCAR tag]][[Category:Molecular dynamics]][[Category:Thermodynamic integration]]		[[Category:INCAR tag]][[Category:Molecular dynamics]][[Category:Advanced molecular-dynamics sampling]][[Category:Thermodynamic integration]]

Ftran at 14:26, 7 October 2024

2024-10-07T14:26:26Z

← Older revision		Revision as of 14:26, 7 October 2024
Line 28:		Line 28:
	::<math> F = -\frac{1}{\beta} \mathrm{ln} \left[ \frac{V^{N}}{\Alpha^{3N} N!} \right] </math>		::<math> F = -\frac{1}{\beta} \mathrm{ln} \left[ \frac{V^{N}}{\Alpha^{3N} N!} \right] </math>

	~~whre~~ <math>V</math> is the volume of the system, <math>N</math> is the number of particles in the system and <math>\Alpha</math> is the de Broglie wavelength. The Stirling approximation applies in principle only in the limes of infinitely many particles. In reference {{cite\|dorner:PRL:2018}} the exact ideal gas equation was used since it helped to speed up the convergence of the final free energy of liquid Si with respect to the system size.		where <math>V</math> is the volume of the system, <math>N</math> is the number of particles in the system and <math>\Alpha</math> is the de Broglie wavelength. The Stirling approximation applies in principle only in the limes of infinitely many particles. In reference {{cite\|dorner:PRL:2018}} the exact ideal gas equation was used since it helped to speed up the convergence of the final free energy of liquid Si with respect to the system size.

	*Harmonic solid:		*Harmonic solid:

Huebsch at 14:25, 7 April 2022

2022-04-07T14:25:55Z

Karsai at 08:19, 3 April 2020

2020-04-03T08:19:37Z

← Older revision		Revision as of 08:19, 3 April 2020
Line 8:		Line 8:


	The free energy of a fully interacting system can be written as the sum of a non-interacting reference system and the difference of the fully interacting system and the non-interacting system		The free energy of a fully interacting system can be written as the sum of the free energy a non-interacting reference system and the difference in the free energy of the fully interacting system and the non-interacting system

	<math> F_{1} = F_{0} + \Delta F_{0\rightarrow 1} </math>.		<math> F_{1} = F_{0} + \Delta F_{0\rightarrow 1} </math>.

Karsai at 08:18, 3 April 2020

2020-04-03T08:18:07Z

← Older revision		Revision as of 08:18, 3 April 2020
Line 32:		Line 32:
	<math> F = -\frac{1}{\beta} \mathrm{ln} \left[ \frac{V^{N}}{\Alpha^{3N} N!} \right] </math>		<math> F = -\frac{1}{\beta} \mathrm{ln} \left[ \frac{V^{N}}{\Alpha^{3N} N!} \right] </math>

	whre <math>V</math> is the volume of the system, <math>N</math> is the number of particles in the system and <math>\Alpha</math> is the de Broglie wavelength.		whre <math>V</math> is the volume of the system, <math>N</math> is the number of particles in the system and <math>\Alpha</math> is the de Broglie wavelength. The Stirling approximation applies in principle only in the limes of infinitely many particles. In reference {{cite\|dorner:PRL:2018}} the exact ideal gas equation was used since it helped to speed up the convergence of the final free energy of liquid Si with respect to the system size.

	*Harmonic solid:		*Harmonic solid:

Karsai at 08:14, 3 April 2020

2020-04-03T08:14:11Z

← Older revision		Revision as of 08:14, 3 April 2020
Line 16:		Line 16:
	<math> \Delta F_{0\rightarrow 1} = \int\limits_{0}^{1} d\lambda \langle U_{1}(\lambda) - U_{0}(\lambda) \rangle_{\lambda} </math>.		<math> \Delta F_{0\rightarrow 1} = \int\limits_{0}^{1} d\lambda \langle U_{1}(\lambda) - U_{0}(\lambda) \rangle_{\lambda} </math>.

	Here <math>U_{1}(\lambda)</math> and <math>U_{0}(\lambda)</math> describe the potential energies of a fully-interacting and a non-interacting reference system, respectively. The coupling strength of the systems is controlled via the coupling parameter <math>\lambda</math>. The notation <math>\langle \ldots \rangle_{\lambda}</math> denotes an ensemble average of a system driven by the following classical Hamiltonian		Here <math>U_{1}(\lambda)</math> and <math>U_{0}(\lambda)</math> describe the potential energies of a fully-interacting and a non-interacting reference system, respectively. The coupling strength of the systems is controlled via the coupling parameter <math>\lambda</math>. It is neccessary that the connection of the two systems via the coupling constant is reversible. The notation <math>\langle \ldots \rangle_{\lambda}</math> denotes an ensemble average of a system driven by the following classical Hamiltonian

	<math> H_{\lambda}= \lambda H_{1} + (1-\lambda) H_{0} </math>.		<math> H_{\lambda}= \lambda H_{1} + (1-\lambda) H_{0} </math>.

Karsai at 08:10, 3 April 2020

2020-04-03T08:10:33Z

← Older revision		Revision as of 08:10, 3 April 2020
Line 10:		Line 10:
	The free energy of a fully interacting system can be written as the sum of a non-interacting reference system and the difference of the fully interacting system and the non-interacting system		The free energy of a fully interacting system can be written as the sum of a non-interacting reference system and the difference of the fully interacting system and the non-interacting system

	<math> F_{1} = F_{0} + \Delta F_{1\rightarrow 0} </math>.		<math> F_{1} = F_{0} + \Delta F_{0\rightarrow 1} </math>.

	Using thermodynamic integration the free energy difference between the two systems is written as		Using thermodynamic integration the free energy difference between the two systems is written as

	<math> \Delta F_{1\rightarrow 0} = \int\limits_{0}^{1} d\lambda \langle U_{1}(\lambda) - U_{0}(\lambda) \rangle_{\lambda} </math>.		<math> \Delta F_{0\rightarrow 1} = \int\limits_{0}^{1} d\lambda \langle U_{1}(\lambda) - U_{0}(\lambda) \rangle_{\lambda} </math>.

	Here <math>U_{1}(\lambda)</math> and <math>U_{0}(\lambda)</math> describe the potential energies of a fully-interacting and a non-interacting reference system, respectively. The coupling strength of the systems is controlled via the coupling parameter <math>\lambda</math>. The notation <math>\langle \ldots \rangle_{\lambda}</math> denotes an ensemble average of a system driven by the following classical Hamiltonian		Here <math>U_{1}(\lambda)</math> and <math>U_{0}(\lambda)</math> describe the potential energies of a fully-interacting and a non-interacting reference system, respectively. The coupling strength of the systems is controlled via the coupling parameter <math>\lambda</math>. The notation <math>\langle \ldots \rangle_{\lambda}</math> denotes an ensemble average of a system driven by the following classical Hamiltonian

Karsai at 08:09, 3 April 2020

2020-04-03T08:09:55Z

← Older revision		Revision as of 08:09, 3 April 2020
Line 7:		Line 7:


	~~Using thermodynamic integration the free energy difference between two systems is written as~~

	<math> \Delta F = \int\limits_{0}^{1} d\lambda \langle U_{1}(\lambda) - U_{0}(\lambda) \rangle_{\lambda} </math>.		The free energy of a fully interacting system can be written as the sum of a non-interacting reference system and the difference of the fully interacting system and the non-interacting system

			<math> F_{1} = F_{0} + \Delta F_{1\rightarrow 0} </math>.

			Using thermodynamic integration the free energy difference between the two systems is written as

			<math> \Delta F_{1\rightarrow 0} = \int\limits_{0}^{1} d\lambda \langle U_{1}(\lambda) - U_{0}(\lambda) \rangle_{\lambda} </math>.

	Here <math>U_{1}(\lambda)</math> and <math>U_{0}(\lambda)</math> describe the potential energies of a fully-interacting and a non-interacting reference system, respectively. The coupling strength of the systems is controlled via the coupling parameter <math>\lambda</math>. The notation <math>\langle \ldots \rangle_{\lambda}</math> denotes an ensemble average of a system driven by the following classical Hamiltonian		Here <math>U_{1}(\lambda)</math> and <math>U_{0}(\lambda)</math> describe the potential energies of a fully-interacting and a non-interacting reference system, respectively. The coupling strength of the systems is controlled via the coupling parameter <math>\lambda</math>. The notation <math>\langle \ldots \rangle_{\lambda}</math> denotes an ensemble average of a system driven by the following classical Hamiltonian

← Older revision		Revision as of 14:25, 7 April 2022
Line 5:		Line 5:
	----		----
	A detailed description of calculations using thermodynamic integration within VASP is given in the supplemental information of reference {{cite\|dorner:PRL:2018}} ('''caution''': the tag ''ISPECIAL''=0 used in that reference is not valid anymore, instead the tag {{TAG\|PHON_NSTRUCT}}=-1 is used).		A detailed description of calculations using thermodynamic integration within VASP is given in the supplemental information of reference {{cite\|dorner:PRL:2018}} ('''caution''': the tag ''ISPECIAL''=0 used in that reference is not valid anymore, instead the tag {{TAG\|PHON_NSTRUCT}}=-1 is used).



	The free energy of a fully interacting system can be written as the sum of the free energy a non-interacting reference system and the difference in the free energy of the fully interacting system and the non-interacting system		The free energy of a fully interacting system can be written as the sum of the free energy a non-interacting reference system and the difference in the free energy of the fully interacting system and the non-interacting system

	<math> F_{1} = F_{0} + \Delta F_{0\rightarrow 1} </math>.		::<math> F_{1} = F_{0} + \Delta F_{0\rightarrow 1} </math>.

	Using thermodynamic integration the free energy difference between the two systems is written as		Using thermodynamic integration the free energy difference between the two systems is written as
Line 18:		Line 16:
	Here <math>U_{1}(\lambda)</math> and <math>U_{0}(\lambda)</math> describe the potential energies of a fully-interacting and a non-interacting reference system, respectively. The coupling strength of the systems is controlled via the coupling parameter <math>\lambda</math>. It is neccessary that the connection of the two systems via the coupling constant is reversible. The notation <math>\langle \ldots \rangle_{\lambda}</math> denotes an ensemble average of a system driven by the following classical Hamiltonian		Here <math>U_{1}(\lambda)</math> and <math>U_{0}(\lambda)</math> describe the potential energies of a fully-interacting and a non-interacting reference system, respectively. The coupling strength of the systems is controlled via the coupling parameter <math>\lambda</math>. It is neccessary that the connection of the two systems via the coupling constant is reversible. The notation <math>\langle \ldots \rangle_{\lambda}</math> denotes an ensemble average of a system driven by the following classical Hamiltonian

	<math> H_{\lambda}= \lambda H_{1} + (1-\lambda) H_{0} </math>.		::<math> H_{\lambda}= \lambda H_{1} + (1-\lambda) H_{0} </math>.


	The tag {{TAG\|SCALEE}} sets the coupling parameter <math>\lambda</math> and hence controls the Hamiltonian of the calculation.		The tag {{TAG\|SCALEE}} sets the coupling parameter <math>\lambda</math> and hence controls the Hamiltonian of the calculation.
	By default {{TAG\|SCALEE}}=1 and the scaling of the energies and forces via the coupling constant is internally skipped in the code. To enable the scaling {{TAG\|SCALEE}}<math>\ne</math>1 has to be specified. A VASP calculation outputs the integrand for a given coupling constant at every molecular dynamics step. How to choose the ensemble size and carry out the integration is described in the main text and especially in the supplementary information of reference {{cite\|dorner:PRL:2018}}.		By default {{TAG\|SCALEE}}=1 and the scaling of the energies and forces via the coupling constant is internally skipped in the code. To enable the scaling {{TAG\|SCALEE}}<math>\ne</math>1 has to be specified. A VASP calculation outputs the integrand for a given coupling constant at every molecular dynamics step. How to choose the ensemble size and carry out the integration is described in the main text and especially in the supplementary information of reference {{cite\|dorner:PRL:2018}}.


	Two possible options are available for the reference system:		Two possible options are available for the reference system:
Line 30:		Line 26:
	By default the thermodynamic integration is carried out from the ideal gas to the fully interacting case (in the case when no {{TAG\|DYNMATFULL}} is present in the calculation folder). Usually the Stirling approximation is used for the free energy of the ideal gas written as		By default the thermodynamic integration is carried out from the ideal gas to the fully interacting case (in the case when no {{TAG\|DYNMATFULL}} is present in the calculation folder). Usually the Stirling approximation is used for the free energy of the ideal gas written as

	<math> F = -\frac{1}{\beta} \mathrm{ln} \left[ \frac{V^{N}}{\Alpha^{3N} N!} \right] </math>		::<math> F = -\frac{1}{\beta} \mathrm{ln} \left[ \frac{V^{N}}{\Alpha^{3N} N!} \right] </math>

	whre <math>V</math> is the volume of the system, <math>N</math> is the number of particles in the system and <math>\Alpha</math> is the de Broglie wavelength. The Stirling approximation applies in principle only in the limes of infinitely many particles. In reference {{cite\|dorner:PRL:2018}} the exact ideal gas equation was used since it helped to speed up the convergence of the final free energy of liquid Si with respect to the system size.		whre <math>V</math> is the volume of the system, <math>N</math> is the number of particles in the system and <math>\Alpha</math> is the de Broglie wavelength. The Stirling approximation applies in principle only in the limes of infinitely many particles. In reference {{cite\|dorner:PRL:2018}} the exact ideal gas equation was used since it helped to speed up the convergence of the final free energy of liquid Si with respect to the system size.
Line 37:		Line 33:
	If the file {{TAG\|DYNMATFULL}} exists in the calculation directory and {{TAG\|SCALEE}}<math>\ne</math>1, the second order Hessian matrix is added to the force and thermodynamic integration from a harmonic model to a fully interacting system is carried out. The {{TAG\|DYNMATFULL}} file stores the eigenmodes and eigenvalues from diagonalizing the dynamic matrix. This file is written by a previous calculation using the {{TAG\|INCAR}} tags {{TAG\|IBRION}}=6 and {{TAG\|PHON_NSTRUCT}}=-1.		If the file {{TAG\|DYNMATFULL}} exists in the calculation directory and {{TAG\|SCALEE}}<math>\ne</math>1, the second order Hessian matrix is added to the force and thermodynamic integration from a harmonic model to a fully interacting system is carried out. The {{TAG\|DYNMATFULL}} file stores the eigenmodes and eigenvalues from diagonalizing the dynamic matrix. This file is written by a previous calculation using the {{TAG\|INCAR}} tags {{TAG\|IBRION}}=6 and {{TAG\|PHON_NSTRUCT}}=-1.

	== Related ~~Tags~~ and ~~Sections~~ ==		== Related tags and articles ==
	{{TAG\|VCAIMAGES}}, {{TAG\|IMAGES}}, {{TAG\|NCORE IN IMAGE1}}, {{TAG\|PHON_NSTRUCT}}, {{TAG\|IBRION}}		{{TAG\|VCAIMAGES}}, {{TAG\|IMAGES}}, {{TAG\|NCORE IN IMAGE1}}, {{TAG\|PHON_NSTRUCT}}, {{TAG\|IBRION}}

Line 46:		Line 42:
	----		----

	[[Category:INCAR]][[Category:Molecular ~~Dynamics~~]][[Category:Thermodynamic integration]]		[[Category:INCAR tag]][[Category:Molecular dynamics]][[Category:Thermodynamic integration]]