Controlling the dynamics of many-electron systems from first principles: a combination of optimal control and time-dependent density-functional theory

Modeling of High-Harmonic Generation in the C60 Fullerene Using Ab Initio, DFT-Based, and Semiempirical Methods

We report calculations of the high-harmonic generation spectra of the C60 fullerene molecule carried out by employing a diverse set of real-time time-dependent quantum chemical methods. All methodologies involve expanding the propagated electronic wave function in bases consisting of the ground and singly excited time-independent eigenstates obtained through the solution of the corresponding linear-response equations. We identify the correlation and exchange effect in the spectra by comparing the results from methods relying on the Hartree-Fock reference determinant with those obtained using approaches based on the density functional theory with different exchange-correlation functionals. The effect of the full random-phase approximation treatment of the excited electronic states is also analyzed and compared with the configuration interaction singles and the Tamm-Dancoff approximation. We also showcase the fact that the real-time extension of the semiempirical method INDO/S can be effectively applied for an approximate description of laser-driven dynamics in large systems.

Bayesian-Based Hybrid Method for Rapid Optimization of NV Center Sensors

NV centers are among the most promising platforms in the field of quantum sensing. Magnetometry based on NV centers, especially, has achieved concrete development in areas of biomedicine and medical diagnostics. Improving the sensitivity of NV center sensors under wide inhomogeneous broadening and fieldamplitude drift is a crucial issue of continuous concern that relies on the coherent control of NV centers with high average fidelity. Quantum optimal control (QOC) methods provide access to this target; nevertheless, the high time consumption of current methods due to the large number of needful sample points as well as the complexity of the parameter space has hindered their usability. In this paper, we propose the Bayesian estimation phase-modulated (B-PM) method to tackle this problem. In the case of the state transforming of an NV center ensemble, the B-PM method reduced the time consumption by more than 90% compared with the conventional standard Fourier basis (SFB) method while increasing the average fidelity from 0.894 to 0.905. In the AC magnetometry scenario, the optimized control pulse obtained with the B-PM method achieved an eight-fold extension of coherence time T2 compared with the rectangular π pulse. Similar application can be made in other sensing situations. As a general algorithm, the B-PM method can be further extended to the open- and closed-loop optimization of complex systems based on a variety of quantum platforms.

Many-electron dynamics in laser-driven molecules: wavefunction theory vs. density functional theory

With recent experimental advances in laser-driven electron dynamics in polyatomic molecules, the need arises for their reliable theoretical modelling. Among efficient, yet fairly accurate methods for many-electron dynamics are Time-Dependent Configuration Interaction Singles (TD-CIS) (a Wave Function Theory (WFT) method), and Real-Time Time-Dependent Density Functional Theory (RT-TD-DFT), respectively. Here we compare TD-CIS combined with extended Atomic Orbital (AO) bases, TD-CIS/AO, with RT-TD-DFT in a grid representation of the Kohn-Sham orbitals, RT-TD-DFT/Grid. Possible ionization losses are treated by complex absorbing potentials in energy space (for TD-CIS/AO) or real space (for RT-TD-DFT), respectively. The comparison is made for two test cases: (i) state-to-state transitions using resonant lasers (π-pulses), i.e., bound electron motion, and (ii) large-amplitude electron motion leading to High Harmonic Generation (HHG). Test systems are a H2 molecule and cis- and trans-1,2-dichlorethene, C2H2Cl2, (DCE). From time-dependent electronic energies, dipole moments and from HHG spectra, the following observations are made: first, for bound state-to-state transitions enforced by π-pulses, TD-CIS nicely accounts for the expected population inversion in contrast to RT-TD-DFT, in agreement with earlier findings. Secondly, when using laser pulses under non-resonant conditions, dipole moments and lower harmonics in HHG spectra are obtained by TD-CIS/AO which are in good agreement with those obtained with RT-TD-DFT/Grid. Deviations become larger for higher harmonics and at low laser intensities, i.e., for low-intensity HHG signals. We also carefully test effects of basis sets for TD-CIS/AO and grid size for RT-TD-DFT/Grid, different exchange-correlation functionals in RT-TD-DFT, and absorbing boundaries. Finally, for the present examples, TD-CIS/AO is observed to be at least an order of magnitude more computationally efficient than RT-TD-DFT/Grid.

Extreme nonlinear strong-field photoemission from carbon nanotubes

Strong-field photoemission produces attosecond (10-18 s) electron pulses that are synchronized to the waveform of the incident light. This nonlinear photoemission lies at the heart of current attosecond technologies. Here we report a new nonlinear photoemission behaviour-the nonlinearity in strong-field regime sharply increases (approaching 40th power-law scaling), making use of sub-nanometric carbon nanotubes and 800 nm pulses. As a result, the carrier-envelope phase sensitive photoemission current shows a greatly improved modulation depth of up to 100% (with a total modulation current up to 2 nA). The calculations reveal that the behaviour is an interplay of valence band optical-field emission with charge interaction, and the nonlinear dynamics can be tunable by changing the bandgap of carbon nanotubes. The extreme nonlinear photoemission offers a new means of producing extreme temporal-spatial resolved electron pulses, and provides a new design philosophy for attosecond electronics and photonics.

Role of initial magnetic disorder: A time-dependent ab initio study of ultrafast demagnetization mechanisms

Despite more than 20 years of development, the underlying physics of the laser-induced demagnetization process is still debated. We present a fast, real-time time-dependent density functional theory (rt-TDDFT) algorithm together with the phenomenological atomic Landau-Lifshitz-Gilbert model to investigate this problem. Our Hamiltonian considers noncollinear magnetic moment, spin-orbit coupling (SOC), electron-electron, electron-phonon, and electron-light interactions. The algorithm for time evolution achieves hundreds of times of speedup enabling calculation of large systems. Our simulations yield a demagnetization rate similar to experiments. We found that (i) the angular momentum flow from light to the system is not essential and the spin Zeeman effect is negligible. (ii) The phonon can play a role but is not essential. (iii) The initial spin disorder and the self-consistent update of the electron-electron interaction play dominant roles and enhance the demagnetization to the experimentally observed rate. The spin disorder connects the electronic structure theory with the phenomenological three-temperature model.

Exact Time-Dependent Exchange-Correlation Potential in Electron Scattering Processes

We identify peak and valley structures in the exact exchange-correlation potential of time-dependent density functional theory that are crucial for time-resolved electron scattering in a model one-dimensional system. These structures are completely missed by adiabatic approximations that, consequently, significantly underestimate the scattering probability. A recently proposed nonadiabatic approximation is shown to correctly capture the approach of the electron to the target when the initial Kohn-Sham state is chosen judiciously, and it is more accurate than standard adiabatic functionals but ultimately fails to accurately capture reflection. These results may explain the underestimation of scattering probabilities in some recent studies on molecules and surfaces.

Multielectron dissociative ionization of methane and formaldehyde molecules with optimally tailored intense femtosecond laser pulses

The multielectron dissociative ionization of CH₄ and CH₂O molecules has been investigated using optimum convolution of different dual tailored short laser pulses. Based on three dimensional molecular dynamics simulations and TDDFT approach, the dissociation probability is enhanced by designing the dual chirped-chirped laser pulses and chirped-ordinary laser pulses for formaldehyde molecule. However, it is interesting to notice that the sensitivity of enhanced dissociation probability into different tailored laser pulses is not significant for methane molecule. In this presented modifications, time variation of bond length, velocity, time dependent electron localization function and evolution of the efficient occupation states are presented to analyze the time evolution of molecular dynamics. This work is proved to be a potential way to reduce the controlling costs with a currently available pulse shaping technology.

Time-Local Equation for the Exact Optimized Effective Potential in Time-Dependent Density Functional Theory

A long-standing challenge in the time-dependent density functional theory is to efficiently solve the exact time-dependent optimized effective potential (TDOEP) integral equation derived from orbital-dependent functionals, especially for the study of nonadiabatic dynamics in time-dependent external fields. In this Letter, we formulate a completely equivalent time-local TDOEP equation that admits a unique real-time solution in terms of time-dependent Kohn-Sham and effective memory orbitals. The time-local formulation is numerically implemented, with the incorporation of exponential memory loss to address the unaccounted for correlation component in the exact-exchange-only functional, to enable the study of the many-electron dynamics of a one-dimensional hydrogen chain. It is shown that the long time behavior of the electric dipole converges correctly and the zero-force theorem is fulfilled in the current implementation.

Selective photo-dissociative ionization of methane molecule with TDDFT study

Three dimensional calculation of control dynamics for finding the optimized laser filed is formulated using an iterative method and time-dependent density functional approach. An appropriate laser pulse is designed to control the desired products in the dissociation of methane molecule. The tailored laser pulse profile, eigenstate distributions and evolution of the efficient occupation numbers are predicted and exact energy levels of this five-atomic molecule is obtained. Dissociation rates of up to 78%, 80%, 90%, and 82% for CH₂⁺, CH⁺, C⁺ and C⁺⁺ are achieved. Based on the present approach one can reduce the controlling costs.

Theoretical Shaping of Femtosecond Laser Pulses for Molecular Photodissociation with Control Techniques Based on Ehrenfest's Dynamics and Time-Dependent Density Functional Theory

The combination of nonadiabatic Ehrenfest-path molecular dynamics (EMD) based on time-dependent density functional theory (TDDFT) and quantum optimal control formalism (QOCT) was used to optimize the shape of ultra-short laser pulses to achieve photodissociation of a hydrogen molecule and the trihydrogen cation H3 (+) . This work completes a previous one [A. Castro, ChemPhysChem, 2013, 14, 1488-1495], in which the same objective was achieved by demonstrating the combination of QOCT and TDDFT for many-electron systems on static nuclear potentials. The optimization model, therefore, did not include the nuclear movement and the obtained dissociation mechanism could only be sequential: fast laser-assisted electronic excitation to nonbonding states (during which the nuclei are considered to be static), followed by field-free dissociation. Here, in contrast, the optimization was performed with the QOCT constructed on top of the full dynamic model comprised of both electrons and nuclei, as described within EMD based on TDDFT. This is the first numerical demonstration of an optimal control formalism for a hybrid quantum-classical model, that is, a molecular dynamics method.

Real-space grids and the Octopus code as tools for the development of new simulation approaches for electronic systems

Real-space grids are a powerful alternative for the simulation of electronic systems. One of the main advantages of the approach is the flexibility and simplicity of working directly in real space where the different fields are discretized on a grid, combined with competitive numerical performance and great potential for parallelization. These properties constitute a great advantage at the time of implementing and testing new physical models. Based on our experience with the Octopus code, in this article we discuss how the real-space approach has allowed for the recent development of new ideas for the simulation of electronic systems. Among these applications are approaches to calculate response properties, modeling of photoemission, optimal control of quantum systems, simulation of plasmonic systems, and the exact solution of the Schrödinger equation for low-dimensionality systems.

Numerical density-to-potential inversions in time-dependent density functional theory

We present an unconditionally stable TDDFT inversion algorithm based on a constrained partial differential equation optimization framework and capable of recovering highly singular potentials.

Time-Dependent Electronic Populations in Fragment-Based Time-Dependent Density Functional Theory

Conceiving a molecule as being composed of smaller molecular fragments, or subunits, is one of the pillars of the chemical and physical sciences and leads to productive methods in quantum chemistry. Using a fragmentation scheme, efficient algorithms can be proposed to address problems in the description of chemical bond formation and breaking. We present a formally exact time-dependent density functional theory for the electronic dynamics of molecular fragments with a variable number of electrons. This new formalism is an extension of previous work [Phys. Rev. Lett. 111, 023001 (2013)]. We also introduce a stable density-inversion method that is applicable to time-dependent and ground-state density functional theories and their extensions, including those discussed in this work.

Adaptive hybrid optimal quantum control for imprecisely characterized systems

Optimal quantum control theory carries a huge promise for quantum technology. Its experimental application, however, is often hindered by imprecise knowledge of the input variables, the quantum system's parameters. We show how to overcome this by adaptive hybrid optimal control, using a protocol named Ad-HOC. This protocol combines open- and closed-loop optimal control by first performing a gradient search towards a near-optimal control pulse and then an experimental fidelity estimation with a gradient-free method. For typical settings in solid-state quantum information processing, adaptive hybrid optimal control enhances gate fidelities by an order of magnitude, making optimal control theory applicable and useful.