this is for holding javascript data
Alfredo A. Correa edited section_Computational_and_Theoretical_Details__.tex
over 8 years ago
Commit id: 3f7330833e4775af8867df908653b49b3a96a636
deletions | additions
diff --git a/section_Computational_and_Theoretical_Details__.tex b/section_Computational_and_Theoretical_Details__.tex
index ca12573..4683033 100644
--- a/section_Computational_and_Theoretical_Details__.tex
+++ b/section_Computational_and_Theoretical_Details__.tex
...
In this work we employed the formalism of TDDFT coupled with Ehrenfest molecular dynamics (EMD) \cite{Gross_1996,Calvayrac_2000,Mason_2007,Alonso_2008,Andrade_2009} to simulate the collision processes between the target electrons and the ion (proton). In TDDFT-EMD method, the dynamics of the electrons are treated quantum mechanically described by TDDFT and the nuclei are point particles treated classically using EMD. The strength of this method is used to calculate the $S_\text{e}$ for metals. %We compared our results with those contained in \textsc{SRIM} database for the case of proton in $\mathrm{Cu}$.
The energy transfered to the electrons of the target ($\mathrm{Cu}$) atom due to a constant velocity moving proton is
carefully monitored.
The energy loss of the proton is negligible hence total energy of the system is not conserved. This is because at At the time scales of the simulations, the large mass of the proton guarantees a negligible decline in its velocity.
For simplicity the proton is forced to move at constant speed hence total energy of the system is not conserved, but the total energy excess of the system correlates with the stopping power.
As the proton moves, the time-dependent Kohn-Sham (TDKS) equation \cite{Runge_1984} describes electronic density and energy of the system due to the dynamics of effective single particle states under the external potential generated by the proton and the crystal of Cu nuclei. The TDKS equation can be written
as (Hartree atomic units are used here): as:
\begin{equation}
\mathrm i\hbar\frac\partial{\partial t}\psi_i(\textbf r, t) = \left\{-\frac{\hbar^2\nabla^2}{2m} + V_\text{KS}[n(t)](\{\mathbf R_i(t)\}_i, \mathbf r, t)\right\}\psi_i(\textbf r, t)