Density Functional Theory Formalism Research Articles

Advancing molecular modeling and reverse vaccinology in broad-spectrum yellow fever virus vaccine development

Yellow fever outbreaks are prevalent, particularly in endemic regions. Given the lack of an established treatment for this disease, significant attention has been directed toward managing this arbovirus. In response, we developed a multiepitope vaccine designed to elicit an immune response, utilizing advanced immunoinformatic and molecular modeling techniques. To achieve this, we predicted B- and T-cell epitopes using the sequences from all structural (E, prM, and C) and nonstructural proteins of 196 YFV strains. Through comprehensive analysis, we identified 10 cytotoxic T-lymphocyte (CTL) and 5T-helper (Th) epitopes that exhibited overlap with B-lymphocyte epitopes. These epitopes were further evaluated for their affinity to a wide range of human leukocyte antigen system alleles and were rigorously tested for antigenicity, immunogenicity, allergenicity, toxicity, and conservation. These epitopes were linked to an adjuvant (β\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\beta$$\\end{document}-defensin) and to each other using ligands, resulting in a vaccine sequence with appropriate physicochemical properties. The 3D structure of this sequence was created, improved, and quality checked; then it was anchored to the Toll-like receptor. Molecular Dynamics and Quantum Mechanics/Molecular Mechanics simulations were employed to enhance the accuracy of docking calculations, with the QM portion of the simulations carried out utilizing the density functional theory formalism. Moreover, the inoculation model was able to provide an optimal codon sequence that was inserted into the pET-28a( +) vector for in silico cloning and could even stimulate highly relevant humoral and cellular immunological responses. Overall, these results suggest that the designed multi-epitope vaccine can serve as prophylaxis against the yellow fever virus.

Current density functional framework for spin-orbit coupling: Extension to periodic systems.

Spin-orbit coupling induces a current density in the ground state, which consequently requires a generalization for meta-generalized gradient approximations. That is, the exchange-correlation energy has to be constructed as an explicit functional of the current density, and a generalized kinetic energy density has to be formed to satisfy theoretical constraints. Herein, we generalize our previously presented formalism of spin-orbit current density functional theory [Holzer et al., J. Chem. Phys. 157, 204102 (2022)] to non-magnetic and magnetic periodic systems of arbitrary dimension. In addition to the ground-state exchange-correlation potential, analytical derivatives such as geometry gradients and stress tensors are implemented. The importance of the current density is assessed for band gaps, lattice constants, magnetic transitions, and Rashba splittings. In the latter, the impact of the current density may be larger than the deviation between different density functional approximations.

Predicting Solvation Free Energies from the Minnesota Solvation Database Using Classical Density Functional Theory Based on the PC-SAFT Equation of State.

We critically assess the capabilities of classical density functional theory (DFT) based on the perturbed-chain statistical associating fluid theory (PC-SAFT) equation of state to predict the solvation free energies of small molecules in various hydrocarbon solvents. We compare DFT results with experimental data from the Minnesota solvation database and utilize statistical methods to analyze the accuracy of our approach, as well as its weaknesses. The mean absolute error of the solvation free energies is 3.7 kJ mol-1 for n-alkane solvents, ranging from pentane to hexadecane, with 473 solute-solvent systems. For solvents consisting of cyclic hydrocarbons (cyclohexane, benzene, toluene, and ethylbenzene) with 245 solute-solvent systems, we report a slightly larger mean absolute error of 4.2 kJ mol-1. We identify three possible sources of errors: (i) the neglect of solute-solvent and solvent-solvent Coulomb interactions, which limits the applicability of PC-SAFT DFT to nonpolar and weakly polar molecules; (ii) the solute's Lennard-Jones parameters supplied by the general AMBER force field, which are not parametrized toward solvation free energies; and (iii) the application of the Lorentz-Berthelot combining rules to the dispersive interactions between a segment of the PC-SAFT solvent and a Lennard-Jones interaction site of the solute. The approach is more accurate than standard implementations of phenomenological models in common chemistry software packages, which exhibit mean absolute errors larger than 9.12 kJ mol-1, even though newer phenomenological models achieve a mean absolute error of about 2 kJ mol-1. PC-SAFT DFT is more computationally efficient than state of the art explicit molecular simulations in combination with free energy perturbation methods. It is predictive with respect to solvation free energies, i.e., the input for the model is the (element-specific) molecular force field, the solute configuration from molecular dynamics simulations, and the (substance-specific) PC-SAFT parameters. The PC-SAFT parametrization uses pure-component data and does not require experimental solvation free energies. The PC-SAFT equation of state, without applying a DFT formalism, can also be used to calculate solvation free energies, provided that the PC-SAFT parameters for the solute are available. A large number of substances was recently parametrized by members of our group (Esper, T.; Bauer, G.; Rehner, P.; Gross, J. Ind. Eng. Chem. Res. 2023, 62), which enables a comparison to the DFT approach for 103 substances. Accurate results are obtained from the PC-SAFT equation of state with an MAE below 2.51 kJ mol-1. The DFT approach does not require PC-SAFT parameters for the solute and can be applied to all solutes that can be represented by the molecular force field.

Realization of IR Photodetector based on AGNRs operating in the Atmospheric Windows

Nowadays, the importance of Atmospheric windows in optoelectronic industries is not hidden from anyone. On the other hand, using graphene nanoribbons widely has been applied in optical devices. So, armchair graphene nanoribbons (AGNRs) from three different families as IR photodetectors are proposed and modeled, in this work. Our results show that the energy gap of AGNRs has been in the 1–3 μm and 8–14 μm atmospheric windows in the pure state and with nitrogen (N)/boron (B) substitutional doping. Also, responsivity and specific detectivity (D*) have been calculated, by introducing a simple structure as an IR photodetector. The energy gap and dark current (Jdark) calculations have been done by density functional theory (DFT) and non-equilibrium Green's function (NEGF) formalism, respectively. The specific detectivity peaks on the order of 109 cmHz1/2/W for 9-AGNR and 10-AGNR at 1–3 μm atmospheric window have been concluded. Furthermore, our results show that 11-AGNR with B dopants can be used to design dual-band detectors in the 1–3 μm and 8–14 μm atmospheric windows.

Boosting the Performance of Diketopyrrolopyrrole-Triphenylamine-Based Organic Solar Cells via π-Linker Engineering.

The design and development of novel and efficient donor-π-acceptor (D-π-A) type conjugated systems has attracted substantial interest in the field of organic electronics owing to their intriguing properties. In this paper, we have designed seven new and efficient D-π-A type conjugated systems (M1-M7) by a variety of π-linkers with triphenylamine (TPA) as the electron donor and diphenyldiketopyrrolopyrrole (DPP) as the electron acceptor using density functional theory (DFT) formalism for organic solar cells (OSCs). The π-linker has been substituted between the donor and acceptor for efficient electron transfer. Here, our primary focus is on narrowing the highest occupied molecular orbital-lowest unoccupied molecular orbital gaps, electronic transition, charge transfer rate, reorganization energies, and the theoretical power conversion efficiencies (PCEs). Our study reveals that the designed compounds exhibit excellent charge transfer rates. The absorption properties of the compounds have been examined using the time-dependent density functional theory (TD-DFT) method. The TD-DFT study shows that compound M2 possesses the highest absorption maxima with a maximum bathochromic shift. For a better understanding of the electron transport process of our designed compounds, we have designed donor/acceptor (D/A) blends, and each of the developed blends (FREA/M1-M7) can encourage charge carrier separation. According to the photovoltaic performance of the D/A blends, compound FREA-M2, which has a theoretical PCE of 16.53%, is the most appealing choice for use in OSCs. We expect that by thoroughly examining the relationship between structure, characteristics, and performance, this work will serve as a roadmap for future research and development of TPA-DPP-based photovoltaic materials.

Theoretical investigation of fluorine-decorated gallium nitride monolayers with vacancy and doping defects

In this manuscript, through the formalism of density functional theory (DFT), we have investigated fluorine-decorated gallium nitride monolayers’ structural stability and electronic properties. We also have introduced vacancy and chemical substitution defects in the decorated structures to expand our investigation, since it is well-known that introducing these defects induces the appearance of notable electronic behaviors. Our results show that adding fluorine to the gallium nitride monolayer increases stability compared to the pristine GaN. In addition, it is possible to verify that the inserting defects in the GaN monolayer fluoride induce remarkable modifications in the energy gap values. Finally, we have shown that an essential property of GaN monolayer is that after fluorine adsorption, such nanostructures become more resistant to redox.

Modeling the near-field effect on molecular excited states using the discrete interaction model/quantum mechanical method.

Strong light-matter interactions significantly modify the optical properties of molecules in the vicinity of plasmonic metal nanoparticles. Since the dimension of the plasmonic cavity approaches that of the molecules, it is critical to explicitly describe the nanoparticle junctions. In this work, we use the discrete interaction model/quantum mechanical (DIM/QM) method to model the coupling between the plasmonic near-field and molecular excited states. DIM/QM is a combined electrodynamics/quantum mechanical model that uses an atomistic description of the nanoparticle. We extend the DIM/QM method to include the local field effects in the sum-over-state formalism of time-dependent density functional theory. As a test of the method, we study the interactions between small organic chromophores and metal nanoparticles. In particular, we examine how the inclusion of multiple electronic transitions and intermolecular interactions modify the coupling between molecules and nanoparticles. Using the sum-over-state formalism of DIM/QM, we show that two-state models break down when the plasmon excitation is detuned from the molecular excitations. To gain further insight, we compare the simple coupled-dipole model (CDM) with the DIM/QM model. We find that CDM works well for simple systems but fails when going beyond the single molecule or single nanoparticle cases. We also find that the coupling depends strongly on the site of the nanoparticle in which the chromophore couples to. Our work suggests the importance of explicitly describing the cavity to capture the atomistic level local field environment in which the molecule strongly couples to.

A New Era for Gold Nitrate Nanoparticle Applications

Gold nanoparticles are included in the group of metallic nanoparticles, and their application has become possible recently because of new optical equipment, electronic devices, and probes capable of manufacturing, detecting, and identifying biomolecules of medical interest. Among the range of nanoparticles, one stands out: the gold nanoparticle. It doesn't possess toxic characteristics and can be administered in vivo. In this study, the structural, optical and electronic properties of gold nitrate were investigated using the density functional theory formalism, considering the gradient generalized approximation. Finally, we observe its application in nanobiotechnology, specifically photothermal applications.

Pressure-induced tuning of physical properties in high-throughput metal halide MSn2Br5 (M = K, Cs) perovskites for optoelectronic applications.

The physical properties of the ferromagnetic oxide perovskites MSn2Br5 (M = K, Cs) were thoroughly examined using the GGA + PBE formalism of density functional theory. The investigation includes a comprehensive characterization of these materials under hydrostatic pressures ranging up to 25 GPa. Our work represents the first theoretical framework for exploring the behavior of MSn2Br5 (M = K, Cs) under pressure, providing valuable insights into their properties. To ensure the thermodynamic and mechanical stability of the studied compounds, we justified their stability through the analysis of formation energy and Born stability criteria. Furthermore, we conducted a thorough examination of the mechanical features of MSn2Br5 (M = K, Cs) based on various parameters, such as elastic constants, elastic moduli, the Kleinman parameter, the machinability index, and the Vickers hardness. Pugh's ratio and Poisson's ratio data show a ductile behavior for both compounds under stress. Moreover, our analysis of the refractive index suggests that both materials hold significant potential as candidates for ultrahigh-density optical data storage devices, particularly when subjected to appropriate laser irradiation. This finding opens up exciting possibilities for utilizing MSn2Br5 (M = K, Cs) in advanced optical technologies.

Ga and In-based hybrid halide perovskites as an alternative to Pb: a first principles study.

Lead-free hybrid halide perovskites have gained much attention in the field of photovoltaics due to their non-toxicity, stability and unique photo-physical properties. Sn and Ge-based ABX3 perovskites have been widely studied due to their similar electronic properties to Pb-based materials. However, the unstable oxidation state of Sn is a major challenge for the commercialization of this class of materials. To overcome this problem, here, we have designed a series of novel Ga and In-based A3B2X9-type perovskite materials incorporating the methylammonium (MA) organic cation in the A site and I- as the halide ion in the X site. In this regard, we have investigated different structural, electronic, optical and photovoltaic properties by employing the density functional theory formalism. The formation of a stable three dimensional perovskite structure is determined by the observed values of tolerance factor (TF) and octahedral factor (μ). The observed negative values of formation enthalpy manifest that our studied materials are also thermodynamically stable. The obtained band gap values reveal that our designed perovskite materials can act as semiconducting materials for application in photovoltaics. We have also investigated the optical properties of our studied materials and the observed values of dielectric function and absorption coefficient in the visible range of the electromagnetic spectrum indicate their excellent photo absorption. The observed theoretical power conversion efficiency (PCE) values reveal that (MA)3In2I9 (13.82%) and (MA)3 (Ga.50In.50)2I9 (12.8%) can be chosen as potential candidates for application in perovskite-based photovoltaics. This research provides a pathway for the development of less toxic and efficient semiconducting materials, offering exciting prospects for their utilization in optoelectronics and contributing to the ongoing efforts to advance sustainable energy technologies.

Computational study of lattice and thermoelectric properties of Rb2LiTlF6 double perovskite

We applied the density functional theory formalism to compute the lattice thermal conductivity, transport and thermoelectric properties of Rb2LiTlF6 double perovskite in the cubic phase. To gain some insights on lattice related dynamic behaviour of this material, we determined its phonon lifetimes and Grüneisen parameters within the considered frequency range. This quaternary compound has a very low value of lattice thermal conductivity, with the predicted value of 0.687 W/mK at 300 K. The highest value of dimensionless figure of merit was determined to be 5.73 × 10−3 at 800 K when hole concentration is 1020 cm−3, this low ZT value could be attributed to the wide band gap nature of this compound. Our results will serve as a benchmark for future experimental and theoretical research on lattice and thermoelectric properties of this material.

1T’-[formula omitted] hybrid monolayer as a novel magnetic material: A first principles study

The advent of two-dimensional (2D) materials, such as transition metal dichalcogenides (2D-TMD), has led to an extensive amount of interest amongst scientists and engineers alike in search of major breakthroughs in the electronic, magnetic, and thermodynamic properties of 2D materials, generating interest in novel device applications. This paper discusses a theoretical study of the structural, electronic, magnetic, stability, phonon dispersion, and thermodynamic properties of the monoclinic RuWTe2 monolayer, a type of TMD in its 1T’ phase. The study uses computational modeling, employing the Density Functional Theory (DFT) and Density Functional Perturbation Theory (DFPT) formalisms, using GGA approximation and pseudopotentials to replace nucleus electrons in each atomic species. Just for bandgap estimation, it was used a single point energy of HSE06 hybrid functional which resulted in a small gap semiconductor (≈ 0.4eV). Even after adding the DFT+U correction to PBE functional, it presented magnetic properties with a trend toward ferrimagnetism. It presents structural and energetic stability through calculations of cohesive, formation energy, and phonon dispersion, since it does not present virtual phonon frequency. Through thermodynamics, the free energy (F) indicates that the synthesis process for 1T’-RuWTe2 would be spontaneous even at low temperatures. These results demonstrate that the monoclinic 1T’-RuWTe2 hybrid monolayer is a promising candidate suitable for applications in magnetic and thermal devices.

A self-adaptive first-principles approach for magnetic excited states

The profound impact of excited magnetic states on the intricate interplay between electron and lattice behaviors in magnetic materials is a topic of great interest. Unfortunately, despite the significant strides that have been made in first-principles methods, accurately tracking these phenomena remains a challenging and elusive task. The crux of the challenge that lies before us is centered on the intricate task of characterizing the magnetic configuration of an excited state, utilizing a first-principle approach that is firmly rooted in the ground state of the system. We propose a versatile self-adaptive spin-constrained density functional theory formalism. By iteratively optimizing the constraining field alongside the electron wave function during energy minimization, we are able to obtain an accurate potential energy surface that captures the longitudinal and transverse variations of magnetization in itinerant ferromagnetic Fe. Moreover, this technique allows us to identify the subtle coupling between magnetic moments and other degrees of freedom by tracking energy variation, providing new insights into the intricate interplay between magnetic interactions, electronic band structure, and phonon dispersion curves in single-layered CrI3\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$\\mathrm{CrI} _{3}$\\end{document}. This new methodology represents a significant breakthrough in our ability to probe the complex and multifaceted properties of magnetic systems.

The convexity condition of density-functional theory.

It has long been postulated that within density-functional theory (DFT), the total energy of a finite electronic system is convex with respect to electron count so that 2Ev[N0] ≤ Ev[N0 - 1] + Ev[N0 + 1]. Using the infinite-separation-limit technique, this Communication proves the convexity condition for any formulation of DFT that is (1) exact for all v-representable densities, (2) size-consistent, and (3) translationally invariant. An analogous result is also proven for one-body reduced density matrix functional theory. While there are known DFT formulations in which the ground state is not always accessible, indicating that convexity does not hold in such cases, this proof, nonetheless, confirms a stringent constraint on the exact exchange-correlation functional. We also provide sufficient conditions for convexity in approximate DFT, which could aid in the development of density-functional approximations. This result lifts a standing assumption in the proof of the piecewise linearity condition with respect to electron count, which has proven central to understanding the Kohn-Sham bandgap and the exchange-correlation derivative discontinuity of DFT.

The effect of oxygen atom substitution in the sulphur sites of the bulk ZnS: A DFT study

Systematic first-principles calculations were performed, using the density functional theory (DFT) formalism to investigate the structural electronic and optical properties of ZnS1-xOx alloys in wurtzite phases. The electronic band structure shows that the fundamental band gap is direct at Γ point. The calculated band gap bowing disorder shows a largest contribution of the charge exchange. We demonstrate the formation of an isolated intermediate electronic band structure which may be attributed to coupling between impurity O-2p states in the conduction band. The optical parameters like dielectric function, optical reflectivity, absorption coefficient, transmittance spectra, and refractive index are performed. Moreover, the intermediate electronic level leads to active absorption of light and low transmittance below the band gap which can generate a significant photocurrent in ZnS1-xOx semiconductors.

Electric field effect on [Formula: see text] clusters for applications in MOSFETs and DSSCs: a DFT study.

The structural, electronic, non-linear optical (NLO) and spectral properties of [Formula: see text] clusters with [Formula: see text] have been studied under the influence of an external electric field using density functional theory (DFT). The effect of variation in the Hf:Ti ratio on different properties of clusters is investigated. The motivation to study [Formula: see text] clusters lies in the fact that HfTiO thin films have wide applications in various optoelectronic and photovoltaic devices. So, it will be interesting to study the effect of electric field on [Formula: see text] clusters with the variation in the number of Hf, Ti and O atoms. It is observed that out of all the clusters, [Formula: see text] and [Formula: see text] are the most stable clusters with high values of binding energy and HOMO-LUMO gap. The application of an external electric field on these most stable clusters distorts their geometry and their HOMO-LUMO gap decreases, dipole moment and polarizability increases as the electric field is increased from 0 a.u. to 340 x[Formula: see text] a.u. The applied electric field increases the polar character of clusters due to electron cloud deformation and hence, increases the reactivity of the clusters, thus making these clusters suitable for electrocatalytic reactions. The electric field controlled high values of dielectric constant makes these clusters suitable to be used in the oxide layer of metal oxide semiconductor field effect transistors (MOSFETs) with better capacitance. Under the effect of an electric field, the absorption peaks of UV-VIS spectra gets red-shifted. Due to the tuning of absorption spectra from ultraviolet to visible region, [Formula: see text] clusters can be thought of as a good replacement for [Formula: see text] in dye-sensitized solar cells (DSSCs). The computational study of [Formula: see text] clusters has been performed using DFT. For the ground state of [Formula: see text] clusters, the optimization and frequency calculations have been performed using hybrid B3LYP (Becke three-parameter exchange functional combined with Lee, Yang and Parr correlation functional) functional with LANL2DZ (Los Alamos National Laboratory effective core potentials with Double Zeta atomic set) basis set under hybrid-GGA (generalized gradient approximation). Optimization and frequency calculations have been performed in each case. The excited state calculations have been carried out within time-dependent DFT formalism for a total of 50 excited states. The computational chemistry software package Gaussian 16 along with its graphical interface Gaussview have been employed for the current study of [Formula: see text] clusters.

Deciphering Interactions between Potential Inhibitors and the Plasmodium falciparum DHODH Enzyme: A Computational Perspective.

Malaria is a parasitic disease that, in its most severe form, can even lead to death. Insect-resistant vectors, insufficiently effective vaccines, and drugs that cannot stop parasitic infestations are making the fight against the disease increasingly difficult. It is known that the enzyme dihydroorotate dehydrogenase (DHODH) is of paramount importance for the synthesis of pyrimidine from the Plasmodium precursor, that is, for its growth and reproduction. Therefore, its blockade can lead to disruption of the parasite's life cycle in the vertebrate host. In this scenario, PfDHODH inhibitors have been considered candidates for a new therapy to stop the parasitic energy source. Given what is known, in this work, we applied molecular fractionation with conjugated caps (MFCC) in the framework of the quantum formalism of density functional theory (DFT) to evaluate the energies of the interactions between the enzyme and the different triazolopyrimidines (DSM483, DMS557, and DSM1), including a complex carrying the mutation C276F. From these results, it was possible to identify the main features of each system, focusing on the wild-type and mutant PfDHODH and examining the major amino acid residues that are part of the four complexes. Our analysis provides new information that can be used to develop new drugs that could prove to be more effective alternatives to present antimalarial drugs.

Double nanorings: Structural and optoelectronic properties based on first principles calculations

In this article, we report the structural and electronic properties of the double nanoring using the Density Functional Theory (DFT) formalism. The double nanorings were built with atoms of Boron (B), Aluminium (Al), Gallium (Ga), Indium (In), Thallium (Tl), Nitrogen (N), Phosphorus (P), Arsenic (As) or Bismuth (Bi). We calculated 19 new structures, being that these molecules showed stability for good convergence criteria. All calculated vibrational frequencies are positive for any molecule studied, indicating that all these structures are in a stable state. The calculated HOMO-LUMO gaps (HLgap) converged to values between 0.244 eV and 2.064 eV, showing promising electronic properties, except for the B8N8 nanoring that reached a HLgap of 4.739 eV, qualifying it as an insulating. The optical absorptions were also obtained for the double nanorings using the Time-Dependent Density Functional Theory (TD-DFT), and their dependence on the wavelength were analyzed, where four of these molecules absorb on the visible region (Al8N8, Al8P8, B8P8, Ga8N8 and Ga8P8).

What can lattice DFT teach us about real-space DFT?

In this paper we establish a connection between density functional theory (DFT) for lattice models and common real-space DFT. We consider the lattice DFT description of a two-level model subject to generic interactions in Mermin's DFT formulation in the grand canonical ensemble at finite temperature. The case of only density-density and Hund's rule interaction studied in earlier work is shown to be equivalent to an exact-exchange description of DFT in the real-space picture. In addition, we also include the so-called pair-hopping interaction which can be treated analytically and, crucially, leads to non-integer occupations of the Kohn-Sham (KS) levels even in the limit of zero temperature. Treating the hydrogen molecule in a minimal basis is shown to be equivalent to our two-level lattice DFT model. By means of the fractional occupations of the KS orbitals (which, in this case, are identical to the many-body ones) we reproduce the results of full configuration interaction, even in the dissociation limit and without breaking the spin symmetry. Beyond the minimal basis, we embed our HOMO-LUMO model into a standard DFT calculation and, again, obtain results in overall good agreement with exact ones without the need of breaking the spin symmetry.

A multi-GPU implementation of Real-Time Time-Dependent Auxiliary Density Functional Theory for the investigation of nanosystems irradiations

We report a new Multi-GPU (Graphical Processor Unit) implementation of real-time time-dependent Auxiliary Density Functional Theory (DFT) for simulations of attosecond electronic dynamics in molecular systems subjected to strong perturbations. Our code relies on the Kohn-Sham formalism of DFT and has been implemented in the deMon2k Fortran code. We expand single-particle wave functions (i.e molecular orbitals) as linear combinations of Gaussian-type-orbitals centered on atoms. The density matrix propagation is carried out on GPU while the Kohn-Sham potential is operated on CPUs (Central Processor Unit) with the help of variationally fitted densities. We propose a parallelization strategy using the MAGMA/CUDA libraries to calculate the exponential of dense Hermitian matrices entering the mathematical definition of the propagator, here using Taylor expansions. We report performance benchmarks on water droplets and on fullerenes (C50 to C540). They show a clear advantage of GPU over CPU (using the Scalapack library). The benchmarks also show the benefit of using more than one GPU for systems comprised of up to more than 10,000 basis functions. There, a speed-up of almost 40 between pure 40 CPU and four 4 GPU is obtained. Attosecond electron dynamics simulation in molecular systems comprised of several thousands of electrons becomes amenable to routine simulations in our code. We assess the accuracy of the GPU implementation considering various applications, namely, the calculation of extreme UV absorption spectra with non-Hermitian dynamics, the response of C180 to an electric perturbation, and finally the irradiation of a DNA/protein complex by a 0.4 MeV proton. The results demonstrate the robustness of the implementation. This work also paves the way for future even more efficient implementations. Program SummaryProgram title: deMon2k-RT-TDDFTCPC Library link to program files: https://doi.org/10.17632/whwc5zkst5.1Developer's repository link: DOI: 10.5281/zenodo.8301468Licensing provisions: CeCILL Free Software License Agreement v2.1Programming language: FortranExternal routines/libraries: MPI, BLAS/LAPACK, CUDA, MAGMASupplementary material: User manual, installation instructions, test cases. Nature of problem: Electron dynamics within molecular systems is amenable to simulation with Real-Time Time-Dependent Density Functional Theory (RT-TD-DFT). The methodology consists in propagating on the attosecond time scale the single particle electronic wave functions on the underlying Kohn-Sham potential [1,2]. In deMon2k-RT-TDDFT, the propagation is coupled to the Auxiliary DFT framework for fast evaluation of the potential [3]. A remaining computational bottleneck is the matrix exponential calculation of the Hermitian Hamiltonian matrix. This task is achieved in deMon2k via diagolanization or expansions formulas (Taylor, Chebychev or Baker-Campbell-Hausdorff). The operation is cumbersome on central-processor-unit machines for large matrices (e.g. > 5000 × 5000), which corresponds to molecular systems comprising around 2,000 electrons. Alternatives are needed to overcome this bottleneck.Solution method: We propose to carry out the matrix exponentiation on graphical-processor-units thanks to the MAGMA[4]/CUDA libraries. The code permits a drastic reduction of the computational cost, opening the attosecond dynamic simulations of large molecular systems comprised thousands of electrons with controlled accuracy. The method is tested here for Taylor expansions. Additional comments including restrictions and unusual features: The code described in this article opens the realization of RT-TD-ADFT simulations. The code is linked to the Master version 6.1.6 of deMon2k (http://www.demon-software.com), which must be pre-installed by the user.