Terahertz pulse induced transitions between ionic and neutral phases and electronic polarization reversal in TTF-CA

We have investigated the terahertz (THz)-pulse-induced dynamics of tetrathiafulvalene-$p$-chloranil near the boundary between the ionic and neutral phases with the use of exact numerical calculations of an extended Hubbard model coupled with lattice motion. For the ionic phase, when the applied electric field of the THz pulse opposes the electronic contribution to the electric polarization (electronic polarization) ${\overline{P}}_{\mathrm{el}}$ of the ground state and the maximum amplitude of electric field is greater than a threshold value, the THz-pulse excited state changes as ${\mathrm{I}}_{\mathrm{A}}\ensuremath{\rightarrow}\mathrm{N}\ensuremath{\rightarrow}{\mathrm{I}}_{\mathrm{B}}\ensuremath{\rightarrow}\mathrm{N}\ensuremath{\rightarrow}{\mathrm{I}}_{\mathrm{A}}\ensuremath{\rightarrow}\mathrm{N}\ensuremath{\rightarrow}\ensuremath{\cdots}\phantom{\rule{4pt}{0ex}}({\mathrm{I}}_{\mathrm{A}}$ ground-state case) or ${\mathrm{I}}_{\mathrm{B}}\ensuremath{\rightarrow}\mathrm{N}\ensuremath{\rightarrow}{\mathrm{I}}_{\mathrm{A}}\ensuremath{\rightarrow}\mathrm{N}\ensuremath{\rightarrow}{\mathrm{I}}_{\mathrm{B}}\ensuremath{\rightarrow}\mathrm{N}\ensuremath{\rightarrow}\ensuremath{\cdots}\phantom{\rule{4pt}{0ex}}({\mathrm{I}}_{\mathrm{B}}$ ground-state case), where N shows the neutral state, ${\mathrm{I}}_{\mathrm{A}}\phantom{\rule{4pt}{0ex}}({\mathrm{I}}_{\mathrm{B}})$ shows the ionic state with ${\overline{P}}_{\mathrm{el}}l0\phantom{\rule{4pt}{0ex}}({\overline{P}}_{\mathrm{el}}g0)$, and the phase of the bond-length alternation of ${\mathrm{I}}_{\mathrm{A}}$ is opposite to that of ${\mathrm{I}}_{\mathrm{B}}$. For the neutral phase, when the maximum amplitude of the electric field is greater than a threshold value, the THz-pulse excited state changes as $\mathrm{N}\ensuremath{\rightarrow}{\mathrm{I}}_{\mathrm{B}}\ensuremath{\rightarrow}\mathrm{N}\ensuremath{\rightarrow}{\mathrm{I}}_{\mathrm{A}}\ensuremath{\rightarrow}\mathrm{N}\ensuremath{\rightarrow}{\mathrm{I}}_{\mathrm{B}}\ensuremath{\rightarrow}\ensuremath{\cdots}$ (positive electric field case) or $\mathrm{N}\ensuremath{\rightarrow}{\mathrm{I}}_{\mathrm{A}}\ensuremath{\rightarrow}\mathrm{N}\ensuremath{\rightarrow}{\mathrm{I}}_{\mathrm{B}}\ensuremath{\rightarrow}\mathrm{N}\ensuremath{\rightarrow}{\mathrm{I}}_{\mathrm{A}}\ensuremath{\rightarrow}\ensuremath{\cdots}$ (negative electric-field case). The phase transitions and electronic polarization reversal are driven by time variation of the lattice order parameter, which indicates the magnitude and phase of the bond-length alternation, and the lattice motion is induced by THz-pulse excitation through the electron--lattice coupling. Transitions between the ionic and neutral phases occur and electronic polarization reverses on a picosecond timescale together with the realizable magnitude of the THz pulse both in the ionic and neutral phases near the phase boundary.