Search

The atom-centered density matrix propagation (ADMP) method is an extended Lagrangian approach to ab initio molecular dynamics, which includes the density matrix in an orthonormalized atom-centered Gaussian basis as additional, fictitious, electronic degrees of freedom, classically propagated along with the nuclear ones. A high adiabaticity between the nuclear and electronic subsystems is mandatory in order to keep the trajectory close to the Born–Oppenheimer (BO) surface. In this regard, the fictitious electronic mass μ, being a symmetric, nondiagonal matrix in its most general form, represents a free parameter, exploitable to optimize the propagation of the electronic density. Although mass-weighting schemes in ADMP exist, a systematic procedure to define an optimal value of the fictitious masses is not available yet. In this work, in order to rationally evaluate the electronic mass, fictitious electronic normal modes are defined through the diagonalization of the Hessian of the electronic density matrix. If the same frequency is imposed on all such modes (compatible with the chosen integration time step), then the corresponding μmatrix can be calculated and then employed for the following propagation. Analysis of several ADMP test simulations reveals that such Hessian-based mass-weighting approach is able to ensure, together with a 0.1/0.2 fs time steps, a high separation between the (real) nuclear and the (fictitious) electronic frequencies, which determines a high adiabaticity. This high, unprecedented, accuracy in the propagation leads, in turn, to low errors in the estimated nuclear vibrational frequencies, making the ADMP method totally comparable to a fully converged BO molecular dynamics simulation but more computationally efficient. This work, therefore, contributes to a further development of the ADMP ab initio molecular dynamics method, aimed at improving its accuracy through a more rational evaluation of the fictitious electronic mass parameter.

In this work, we computed and analyzed, by means of density-based descriptors, the real-time evolution of both the locally excited (LE) and charge-transfer (CT) excited states for the planar and twisted conformations of the DMABN (4-(N,N-dimethylamino)benzonitrile) molecule using real-time time-dependent density functional theory (DFT) and three different exchange–correlation energy functionals (EXC) belonging to the same family (the PBE one). Our results based on the analysis of density-based descriptors show that the underlying EXC modifies the evolution in time of the density. In particular, comparing the frequency of density reorganization computed with the three functionals (PBE, PBE0, and LC-PBE), we found that the frequency of electronic interconversion of the individual determinants involved during the dynamics increases from PBE to PBE0 and to LC-PBE. This allows us to show that there is a correlation between the delocalization of the electronic density and the frequency of reorganization. In particular, the greater the mean hole–electron distance during the dynamics, the lower is the frequency of density reorganization.

Modern, nonlinear, time-resolved spectroscopic techniques have opened new doors for investigating the intriguing but complex world of photoinduced ultrafast out-of-equilibrium phenomena and charge dynamics. The interaction between light and matter introduces an additional dimension, where the complex interplay between electronic and vibrational dynamics needs the most advanced theoretical–computational protocols to be fully understood on the molecular scale. In this study, we showcase the capabilities of ab initio molecular dynamics simulation integrated with a multiresolution wavelet protocol to carefully investigate the excited-state relaxation dynamics in a noncovalent complex involving tetramethylbenzene (TMB) and tetracyanoquinodimethane (TCNQ) undergoing charge transfer (CT) upon photoexcitation. Our protocol provides an accurate description that facilitates a direct comparison between transient vibrational analysis and time-resolved spectroscopic signals. This molecular level perspective enhances our understanding of photorelaxation processes confined in the adiabatic regime and offers an improved interpretation of vibrational spectra. Furthermore, it enables the quantification of anharmonic vibrational couplings between high- and low-frequency modes, specifically the TCNQ “rocking” and “bending” modes. Additionally, it identifies the primary vibrational mode that governs the adiabaticity between the ground state and the CT state. This comprehensive understanding of photorelaxation processes holds significant importance in the rational design and precise control of more efficient photovoltaic and sensor devices.

Here is presented the ultrafast hole–electron dynamics of photoinduced metal to ligand charge-transfer (MLCT) states in a Ru(II) complex, [Ru(dcbpy)2(NCS)2]4–(dcbpy = 4,4′-dicarboxy-2,2′-bipyridine), a photoactive molecule employed in dye sensitized solar cells. Via cutting-edge computational techniques, a tailored computational protocol is here presented and developed to provide a detailed analysis of the electronic manifold coupled with nuclear vibrations to better understand the nonradiative pathways and the resulting overall dye performances in light-harvesting processes (electron injection). Thus, the effects of different vibrational modes were investigated on both the electronic levels and charge transfer dynamics through a theoretical-computational approach. First, the linear response time-dependent density functional (LR-TDDFT) formalism was employed to characterize excitation energies and spacing among electronic levels (the electronic layouts). Then, to understand the ultrafast (femtosecond) charge dynamics on the molecular scale, we relied on the nonperturbative mean-field quantum electronic dynamics via real-time (RT-) TDDFT. Three vibrational modes were selected, representative for collective nuclear movements that can have a significant influence on the electronic structure: two involving NCS–ligands and one involving dcbpy ligands. As main results, we observed that such MLCT states, under vibrational distortions, are strongly affected and a faster interligand electron transfer mechanism is observed along with an increasing MLCT character of the adiabatic electronic states approaching closer in energy due to the vibrations. Such findings can help both in providing a molecular picture of multidimensional vibro-electronic spectroscopic techniques, used to characterize ultrafast coherent and noncoherent dynamics of complex systems, and to improve dye performances with particular attention to the study of energy or charge transport processes and vibronic couplings.

The interplay between light absorption and the molecular environment has a central role in the observed photophysics of a wide range of photoinduced chemical and biological phenomena. The understanding of the interplay between vibrational and electronic transitions is the focus of this work, since it can provide a rationale to tune the optical properties of charge transfer (CT) materials used for technological applications. A clear description of these processes poses a nontrivial challenge from both the theoretical and experimental points of view, where the main issue is how to accurately describe and probe drastic changes in the electronic structure and the ultrafast molecular relaxation and dynamics. In this work we focused on the intermolecular CT reaction that occurs upon photon absorption in a π-stacked model system in dichloromethane solution, in which the 1-chloronaphthalene (1ClN) acts as the electron donor and tetracyanoethylene (TCNE) is the electron acceptor. Density functional theory calculations have been carried out to characterize both the ground-state properties and more importantly the low-lying CT electronic transition, and excellent agreement with recently available experimental results [Mathies, R. A.; et al. J. Phys. Chem. A2018, 122(14), 3594] was obtained. The minima of the ground state and first singlet excited state have been accurately characterized in terms of spatial arrangements and vibrational Raman frequencies, and the CT natures of the first two low-lying electronic transitions in the absorption spectra have been addressed and clarified too. Finally, by modeling the possible coordination sites of the TCNE electron acceptor with respect to monovalent ions (Na+, K+) in an implicit solution of acetonitrile, we find that TCNE can accommodate a counterion in two different arrangements, parallel and orthogonal to the C═C axis, leading to the formation of a contact ion pair. The nature of the counterion and its relative position entail structural modifications of the TCNE radical anion, mainly the central C═C and C≡N bonds, compared to the isolated case. An important red shift of the C═C stretching frequency was observed when the counterion is orthogonal to the double bond, to a greater extent for Na+. On the contrary, in the second case, where the counterion ion lies along the internuclear C═C axis, we find that K+polarizes the electron density of the double bond more, resulting in a greater red shift than with Na+.

Defects within nanodiamonds have enabled a variety of quantum sensors based on optical properties that are sensitive to pressure, optical/electric fields, and the configuration of both electronic and nuclear of spins. The presence of dopants can introduce midgap states influencing their optical properties. Complex defects based on nitrogen aggregates have been shown to form in nitrogen-containing diamonds; as such, it is important to differentiate these systems from pure diamond and single nitrogen vacancy (NV) systems. Here we report the effect of N4VxB- and H4-aggregate defects on the infrared vibrational, optical, and X-ray absorption spectroscopies of nanodiamonds. It is found that the presence of these polyatomic nitrogen-aggregate defects introduces unique vibrational responses, as well as electronic levels giving rise to unique optical and X-ray absorption features.

Diamonds approaching the nanoscale have the potential for use as probe materials as their optical properties can be sensitive to optical/electric fields, mechanical stress/pressure, and the configuration of nuclear spins. The surface of nanodiamonds impacts their optical properties and sensing capabilities, and examining the nanodiamond surface through X-ray absorption can give insights into molecular surface structures. Here, quantum dot models with varying amounts of surface carbon passivation are prepared, optimized, and compared. The loss of the diamond sp3lattice is examined by investigating the bond length and tetrahedral character of the carbons comprising nanodiamonds for the appearance of aromatic sp2surface domains. Electronic transitions in the carbon K-edge region, using the energy-specific time-dependent density functional theory method, as well as vibrational spectra are computed from the optimized models. The surface reorganization is shown to affect the electronic characteristics of the nanodiamond. As a result, there is a distinct absorption peak in the carbon K-edge region, along with stretching modes in the vibrational spectra, that can be correlated to the nature of the surface hybridization of the nanodiamond.

Projected Hartree–Fock (PHF) theory can restore important symmetries to broken symmetry wave functions. Variation after projection (VAP) implementations make it possible to deliberately break and then restore a given symmetry by directly minimizing the projected energy expression. This technique can be applied to any symmetry that can be broken from relaxing constraints on single Slater determinant wave functions. For instance, generalized Hartree–Fock (GHF) wave functions are eigenfunctions of neither Ŝznor S2. By relaxing these constraints, the wave function can explore a larger variational space and can reach lower energies than more constrained HF solutions. We have implemented spin-projected GHF (SGHF), which retains many of the advantages of breaking symmetry while also being a spin eigenfunction, with some notable improvements over previous implementations. Our new algorithm involves the formation of new intermediate matrices not previously discussed in the literature. Discretization of the necessary integration over the rotation group SO(3) is also accomplished much more efficiently using Lebedev grids. A novel scheme to incrementally build rotated Fock matrices is also introduced and compared with more standard approaches.

Photoinduced charge transfer (CT) excited states and their relaxation mechanisms can be highly interdependent on the environment effects and the consequent changes in the electronic density. Providing a molecular interpretation of the ultrafast (subpicosecond) interplay between initial photoexcited states in such dense electronic manifolds in condensed phase is crucial for improving and understanding such phenomena. Real-time time-dependent density functional theory is here the method of choice to observe the charge density, explicitly propagated in an ultrafast time domain, along with all time-dependent properties that can be easily extracted from it. A designed protocol of analysis for real-time electronic dynamics to be applied to time evolving electronic density related properties to characterize both in time and in space CT dynamics in complex systems is here introduced and validated, proposing easy to be read cross-correlation maps. As case studies to test such tools, we present the photoinduced charge-transfer electronic dynamics of 5-benzyluracil, a mimic of nucleic acid/protein interactions, and the metal-to-ligand charge-transfer electronic dynamics in water solution of [Ru(dcbpy)2(NCS)2]4–, dcbpy = (4,4′-dicarboxy-2,2′-bipyridine), or “N34–”, a dye sensitizer for solar cells. Electrostatic and explicit ab initio treatment of solvent molecules have been compared in the latter case, revealing the importance of the accurate modeling of mutual solute–solvent polarization on CT kinetics. We observed that explicit quantum mechanical treatment of solvent slowed down the charge carriers mobilities with respect to the gas-phase. When all water molecules were modeled instead as simpler embedded point charges, the electronic dynamics appeared enhanced, with a reduced hole–electron distance and higher mean velocities due to the close fixed charges and an artificially increased polarization effect. Such analysis tools and the presented case studies can help to unveil the influence of the electronic manifold, as well as of the finite temperature-induced structural distortions and the environment on the ultrafast charge motions.

We simulated the intrinsic reaction path of the Green Fluorescent Protein (GFP) proton shuttle in both the ground state (S0) and first singlet excited state (S1), accounting for the main energetic and steric effects of the protein in a convenient model including the chromophore, the crystallographic water, and the residues directly involved in the proton transfer event. We adopted density functional theory (DFT) and time-dependent density functional theory (TD-DFT) levels to define the potential energy surfaces of the two electronic states, and we compared results obtained by the Damped Velocity Verlet and the Hessian-based Predictor–Corrector integrators of the intrinsic reaction coordinate, which gave a comparable and consistent picture of the mechanism. We show that, at S1, the GFP proton transfer becomes favored, with respect to S0, as suggested by the experimental evidence. As an important finding, this change is strictly related to the rearrangement of the hydrogen bond network composing the reaction path, which, in S1, relaxes to a tighter and planar configuration, as a consequence of the photoinduced relaxation in the GFP chromophore structure, thus prompting more effectively for the proton shuttle. Therefore, we give an unprecedented direct proof of the key role played by the photoinduced structural relaxation of the GFP on the chromophore photoacidity, validating, in particular, the hypothesis of Fang and co-workers [Nature2009, 462, 200].

Ion diffusion in semiconductor nanocrystals, or quantum dots (QDs), has gained recognition in recent years as a crucial process for advancing both energy storage and, more generally, the postsynthetic p-type doping chemistry of these materials. In this report, we present first an energetic analysis of group I cations (H+, Li+, and Na+) diffusion in (MX)84–QDs, with M = Zn, Cd and X = S, Se. The bound solutions to the corresponding one-dimensional nuclear Schrödinger equation were solved for these systems, relying on the discrete variable representation method. From this vantage, the quantum nature of the intercalating ion can be revealed. Evidence for the importance of including quantum effects in the treatment of these diffusion processes is presented, both with the density of energy eigenstates of the intercalating ion and from a comparison of the standard deviation in the population distribution of the intercalating ion to the lattice spacings of its host material. Results suggest that the use of classical mechanics for simulations of the ion diffusion processes in these and other related materials can be a questionable practice for the smallest group I cations. Trends devised herein can be useful to help guide the development of new experimental approaches to postsynthetic doping of semiconductor nanocrystals, and in designing electrode materials for next generation electrochemical energy storage devices.

Tractable methods for studying the molecular dynamics of chemical processes driven by electronic nonadiabaticity are highly sought after to provide insight into, for example, photochemical reaction mechanisms, molecular collisions, and thermalized electronic band structures. Starting from the time-dependent Schrödinger equation for a many-body system, a direct ab initiotrajectory surface-hopping (TSH) method relying on an analytical treatment of nonadiabatic couplings between electronic states is developed in this work. An approach that combines time-dependent perturbation theory and explicit time evolution via TSH to expedite calculation of nonadiabatic transition rates, namely, meta-surface-hopping dynamics, is presented, and an extrapolatory approach using time-dependent perturbation theory for recovering unbiased transition rates is assessed. The meta-surface-hopping method is applied to the problem of estimating nonradiative relaxation rates of a photoexcited iminium ion, CH2NH2+, and evidence for internal consistency of the combined dynamics/perturbation theory approach is presented.

Charge-transfer dynamics and interligand electron transfer (ILET) phenomena play a pivotal role in dye-sensitizers, mostly represented by the Ru-based polypyridyl complexes, for TiO2and ZnO-based solar cells. Starting from metal-to-ligand charge-transfer (MLCT) excited states, charge dynamics and ILET can influence the overall device efficiency. In this letter, we focus on N34–dye ( [Ru(dcbpy)2(NCS)2]4–, dcbpy = 4,4′-dicarboxy-2,2′-bipyridine) to provide a first direct observation with high time resolution (<20 fs) of the ultrafast electron exchange between bpy-like ligands. ILET is observed in water solution after photoexcitation in the ∼400 nm MLCT band, and assessment of its ultrafast time-scale is here given through a real-time electronic dynamics simulation on the basis of state-of-the-art electronic structure methods. Indirect effects of water at finite temperature are also disentangled by investigating the system in a symmetric gas-phase structure. As main result, remarkably, the ILET mechanism appears to be based upon a purely electronic evolution among the dense, experimentally accessible, MLCT excited states manifold at ∼400 nm, which rules out nuclear–electronic couplings and proves further the importance of the dense electronic manifold in improving the efficiency of dye sensitizers in solar cell devices.