Search

We present a real-time time-dependent four-component Dirac-Kohn-Sham (RT-TDDKS) implementation based on the BERTHA code. This new implementation takes advantage of modern software engineering, including the prototyping techniques. The software design follows a three step approach: (i) the prototype implementation of a time-propagation algorithm in nonrelativistic real-time TDDFT within the Psi4NumPy framework, which provides a suitable environment for the creation of a clear, readable, and easy to test reference code in Python, (ii) the design of an original Python application programming interface for the relativistic four-component code BERTHA (PyBERTHA), which has an efficient computational kernel for relativistic integrals written in FORTRAN, and (iii) the porting of the time-propagation scheme enveloped within the Psi4NumPy framework to PyBERTHA. The propagation scheme consequently resides in a single readable Python computer code that is easy to maintain and in which the key quantities, such as the Dirac-Kohn-Sham and dipole matrices, can be accessed directly from the PyBERTHA module. For linear algebra operations (matrix-matrix multiplications and diagonalization) we use the highly optimized procedures implemented in the popular NumPy library. The overhead introduced by the Python interface to BERTHA is almost negligible (less than 1% evaluated on the SCF procedure), and the interoperability between different programming languages (FORTRAN, C, and Python) does not affect the numerical stability of the time-propagation scheme. Our new RT-TDDKS implementation has been employed to investigate the stability of the time-propagation procedure in combination with a density-fitting algorithm (both for the Coulomb and for the exchange-correlation matrix construction), which are employed in BERTHA to speed up the Dirac-Kohn-Sham matrix evaluation. On the basis of systematic calculations, employing several density-fitting basis sets of increasing accuracy, we showed that quantitative agreement can be achieved in combination with extended-fitting basis sets, with an error in the Coulomb energy below 1 μ-hartree. Convergence of the transition energies increasing of quality of the fitting basis sets has been also observed. Our data suggest that the error in the Coulomb energy may also represent a good estimate of the fitting basis set quality for real-time electron dynamics simulations. Further, we study the applicability of the RT-TDDKS method in combination with both weak- and extreme strong-field regime. Numerical results of excited-state transitions for the Group 12 atoms are reported and compared with a previous real-time Dirac-Kohn-Sham implementation (Repisky et al. J. Chem. Theory Comput. 2015, 11, 980-991). Finally, calculations of high harmonic generation in the hydrogen molecule and Au dimer have been also carried out. We were able to generate high harmonics with relatively well-defined peaks up to the 21st and 13th order in the case of H 2 and Au 2 , respectively. Our findings show that the four-component structure of the Dirac-Kohn-Sham Hamiltonian provides a suitable theoretical framework, with no intrinsic unfavorable features, to study molecules in the strong-field regime.

We have recently introduced a simple yet powerful tool for analyzing quantitatively the coordination bond in terms of the donation and back-donation constituents of the Dewar-Chatt-Duncanson model. Our approach is based on the decomposition, via natural orbitals for chemical valence (NOCV), of the so-called charge-displacement (CD) function into additive chemically meaningful components (Bistoni et al. J. Chem. Phys. 2015, 142, 084112 ). The method, referred to as NOCV/CD, provides clear-cut measures of donation and back-donation charge flows following bond formation, and its robustness has been demonstrated by a tight correlation of the related charge-transfer estimates with experimental observables. In this paper we extend the NOCV/CD analysis scheme to the four-component relativistic framework, which includes spin-orbit coupling variationally. This formalism is incorporated into a recently developed, highly efficient parallel version of the relativistic Dirac-Kohn-Sham (DKS) program BERTHA (Rampino et al. J. Chem. Theory Comput. 2014, 10, 3766-3776 ). We test the accuracy and numerical stability of this new implementation through the analysis of the convergence properties of the basis sets employed to expand the DKS spinor solution and those used to linearize the electronic density in the density-fitting algorithm speeding up the evaluation of the DKS matrix. An illustration of NOCV/CD analysis in the relativistic framework is also given through a study of the metal-carbonyl coordination bond in a series of [M-CO] + (M = Cu, Ag, Au) complexes of group 11 metals, where relativistic effects, including spin-orbit coupling, are found to play an important role.

We propose a new complete memory-distributed algorithm, which significantly improves the parallel implementation of the all-electron four-component Dirac-Kohn-Sham (DKS) module of BERTHA (J. Chem. Theory Comput. 2010, 6, 384). We devised an original procedure for mapping the DKS matrix between an efficient integral-driven distribution, guided by the structure of specific G-spinor basis sets and by density fitting algorithms, and the two-dimensional block-cyclic distribution scheme required by the ScaLAPACK library employed for the linear algebra operations. This implementation, because of the efficiency in the memory distribution, represents a leap forward in the applicability of the DKS procedure to arbitrarily large molecular systems and its porting on last-generation massively parallel systems. The performance of the code is illustrated by some test calculations on several gold clusters of increasing size. The DKS self-consistent procedure has been explicitly converged for two representative clusters, namely Au20 and Au34, for which the density of electronic states is reported and discussed. The largest gold cluster uses more than 39k basis functions and DKS matrices of the order of 23 GB.

Observations of long-lived coherences in photosynthetic light-harvesting complexes utilize short pulses with broad spectral bandwidths to coherently excite multiple transitions and coherent superpositions. In order to identify the role that such quantum effects might play in efficient energy transfer, however, an alternative approach is required. We have developed a technique for two-color photon echo spectroscopy to selectively excite the pathway of interest and measure its evolution in the absence of any other excitation. We use this technique to excite a coherence pathway in phycocyanin-645 from cryptophyte algae and measure the dynamics of this coherence. A decoherence time of 500 fs was measured, and clear signatures for strong coupling between the electronic states and phonon modes were observed, allowing coherent coupling between otherwise nonresonant transitions. This provides detailed experimental evidence of the long-lived coherences and the nature of the quantum mechanical interactions between electronic states and phonon modes in phycocyanin-645 from cryptophyte marine algae.

We show that all-electron relativistic four-component Dirac-Kohn-Sham (DKS) computations, using G-spinor basis sets and state-of-the-art density fitting algorithms, can be efficiently parallelized and applied to large molecular systems, including large clusters of heavy atoms. The performance of the parallel implementation of the DKS module of the program BERTHA is illustrated and analyzed by some test calculations on several gold clusters up to Au32, showing that calculations with more than 25 000 basis functions (i.e., DKS matrices on the order of 10 GB) are now feasible. As a first application of this novel implementation, we investigate the interaction of the atom Hg with the Au20 cluster.

Spectroscopic constants, including dissociation energies, harmonic and anharmonic vibrational frequencies, and dipole moments, are calculated for the complete alkali auride series (LiAu, NaAu, KAu, RbAu, CsAu). The four-component formulation of relativistic density functional theory has been employed in this study, using the G-spinor basis sets implemented recently in the program BERTHA. The performance of four standard nonrelativistic density functionals employed is investigated by comparing the results with the best available theoretical and experimental data. The present work provides the first theoretical predictions on the molecular properties of RbAu. The intermetallic bond that occurs in the alkali auride series is highly polar and is characterized by a large charge transfer from the alkali metals to gold. The extent of this electron transfer has been investigated using several different charge analysis methods, enabling us to reach some general conclusions on their relative performance. We further report a detailed analysis of the topological properties of relativistic electron density in the bonding region, discussing the features of this approach which characterize the nature of the chemical bond. We have also computed the fully relativistic density for the alkali halides MBr and MI (M = Li, Na, K, Rb, and Cs). The comparative study shows that, on the basis of several topological properties and the variation in bond lengths, the gold atom behaves similarly to a halogen intermediate between Br and I.