An efficient method for calculating dynamical hyperpolarizabilities using real-time time-dependent density functional...

In this paper we present a time-domain time-dependent density functional theory (TDDFT) approach to calculate frequency-dependent polarizability and hyperpolarizabilities. In this approach, the electronic degrees of freedom are propagated within the density matrix based TDDFT framework using the efficient modified midpoint and unitary transformation algorithm. We use monochromatic waves as external perturbations and apply the finite field method to extract various orders of the time-dependent dipole moment. By fitting each order of time-dependent dipole to sinusoidal waves with harmonic frequencies, one can obtain the corresponding (hyper)polarizability tensors. This approach avoids explicit Fourier transform and therefore does not require long simulation time. The method is illustrated with application to the optically active organic molecule para-nitroaniline, of which the frequency-dependent polarizability α(-ω; ω), second-harmonic generation β(-2ω; ω, ω), optical rectification β(0; -ω, ω), third-harmonic generation γ(-3ω; ω, ω, ω), and degenerate four-wave mixing γ(-ω; ω, ω, -ω) are calculated.