Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision		 Previous revision
Last revisionBoth sides next revision
nuclear-electronic_orbital_method [2024/01/29 13:06] – [Dumping] rmatalhasecknuclear-electronic_orbital_method [2024/01/29 13:29] – [Bibliography] rmatalhaseck
Line 96: Line 96:
   * **''NEORD'', //number//** sets the start for the fast rotational update of the orbitals in the local version   * **''NEORD'', //number//** sets the start for the fast rotational update of the orbitals in the local version
   * **''NOBLOCKDIAG''** disables the block diagonalization of the nuclear starting guess (this is generally not recommended!!)   * **''NOBLOCKDIAG''** disables the block diagonalization of the nuclear starting guess (this is generally not recommended!!)
 +  * **''NEOMIXBAS''** enables the use of user-defined mixed basis sets (see example for use)
 ===== Adaptive NEO ===== ===== Adaptive NEO =====
  
Line 162: Line 162:
 </code> </code>
  
-The second example shows the input of a **''LDF-NEO-RHF''** computation of the same molecule starting from a prior RHF calculation. +The second example shows the input of a **''LDF-NEO-RHF''** computation of the same molecule starting from a prior RHF calculation. In this example a [[dump_density_or_orbital_values_cube|cube]] file is requested. This will output the quantum nuclei density.
  
 <code> <code>
Line 201: Line 201:
 {cube,nuclear.cube;density,2102.2} {cube,nuclear.cube;density,2102.2}
 </code> </code>
-In the last example a [[dump_density_or_orbital_values_cube|cube]] file is requested. This will output the quantum nuclei density.+ 
 +The following example shows a NEO calculation, where a user-defined mixed basis set is used. Thereby, the electronic basis set at the quantum nuclei is larger than for regular hydrogen atoms. The use of the **''NEOMIXBAS''** requires the additional definition of the **''elebas''** and **''elefit''** basis sets as shown below. 
 + 
 +<code> 
 +memory,50,m 
 +gdirect 
 +nosym 
 + 
 +geometry={ 
 +3 
 + 
 +H1  -3.5008791    1.2736107    0.7596000 
 +H2  -4.9109791    1.2967107    0.1521000 
 +O   -3.9840791    1.3301107   -0.0574000 
 +} 
 + 
 +charge=0 
 + 
 +basis={ 
 +default=cc-pvtz 
 +H1=cc-pv5z 
 + 
 +set,nucbas 
 +default=neo-basis 
 +H1=pb4-f2 
 + 
 +set,nucfit 
 +default=neo-basis 
 +H1=10s10p10d10f 
 + 
 +set,elebas 
 +default=cc-pvtz 
 +H1=cc-pv5z 
 + 
 +set,elefit,context=jkfit 
 +default=cc-pvtz 
 +H1=cc-pv5z 
 +} 
 + 
 +qnuc,H1 
 + 
 +{df-neo-rhf,maxdis=10,maxit=1000,df_basis=elefit 
 +neoit,100 
 +neothre,1.d-6 
 +neothrie,1.d-8 
 +neothrin,1.d-8 
 +neothrd,1.d-8 
 +neothrg,1.d-8 
 +neoatden 
 +neomixbas 
 +} 
 +</code> 
 + 
 +The example below shows the input for an adaptive NEO calculation, where the nuclear basis function centers convergence is set below 1E-5 bohr and a damping factor of 0.5 is applied. 
 + 
 +<code> 
 +memory,50,m 
 +gdirect 
 +nosym 
 + 
 +geometry={ 
 +3 
 + 
 +H1  -3.5008791    1.2736107    0.7596000 
 +H2  -4.9109791    1.2967107    0.1521000 
 +O   -3.9840791    1.3301107   -0.0574000 
 +} 
 + 
 +charge=0 
 + 
 +basis={ 
 +default=cc-pvdz 
 + 
 +set,nucbas 
 +default=neo-basis 
 +H1=pb4-f2 
 + 
 +set,nucfit 
 +default=neo-basis 
 +H1=10s10p10d10f 
 +} 
 + 
 +qnuc,H1 
 + 
 +{df-neo-rhf,maxdis=10,maxit=500,df_basis=cc-pvdz 
 +neoit,100 
 +neothre,1.d-6 
 +neothrie,1.d-8 
 +neothrin,1.d-8 
 +neothrd,1.d-8 
 +neothrg,1.d-8 
 +adaptive 
 +adthres,1.d-5 
 +addump,0.5 
 +} 
 +</code>
 ===== Bibliography ===== ===== Bibliography =====
  
Line 216: Line 311:
 ===(L)DF-NEO-RHF=== ===(L)DF-NEO-RHF===
  
-Lukas Hasecke, and Ricardo A. Mata [[https://doi.org/10.21203/rs.3.rs-3231458/v1|Nuclear quantum effects made accessible: local-density fitting in multicomponent methods]] //Research Square// **2023** preprint. +Lukas Hasecke, and Ricardo A. Mata [[https://doi.org/10.1021/acs.jctc.3c01055|Nuclear Quantum Effects Made Accessible: Local Density Fitting in Multicomponent Methods]] //J. Chem. Theory Comput.// **2023** //19// (22), 8223–8233.