Born-Oppenheimer Molecular Dynamics Trajectories Research Articles

Pd8 Cluster: Too Small to Melt? A BOMD Study.

The question of whether a solid-liquid phase transition occurs in small clusters poses a fundamental challenge. In this study, we attempt to elucidate this phenomenon through a thorough examination of the thermal behavior and structural stability of Pd8 clusters employing ab initio simulations. Initially, a systematic global search is carried out to identify the various isomers of the Pd8 cluster. This is accomplished by employing an ab initio basin-hopping algorithm and using the PBE/SDD scheme integrated in the Gaussian code. The resulting isomers are further refined through reoptimization using the deMon2k package. To ensure the structural firmness of the lowest-energy isomer, we calculated normal modes. The structural stability as a function of temperature is analyzed through the Born-Oppenheimer molecular dynamics (BOMD) approach. Multiple BOMD trajectories at distinct simulated temperatures are examined with data clustering analysis to determine cluster isomers. This analysis establishes a connection between the potential energy landscape and the simulated temperature. To address the question of cluster melting, canonical parallel-tempering BOMD runs are performed and analyzed with the multiple-histogram method. A broad maximum in the heat capacity curve indicates a melting transition between 500 and 600 K. To further examine this transition, the mean-squared displacement and the pair-distance distribution function are calculated. The results of these calculations confirm the existence of a solid-liquid phase transition, as indicated by the heat capacity curve.

Carbene Cross-Linking in Gas-Phase Peptide Ion Scaffolds.

Scaffolds consisting of a peptide, a phthalate linker, and a 4,4-azipentyl group were synthesized and used to study intramolecular peptide-carbene cross-linking in gas-phase cations. Carbene intermediates were generated by UV-laser photodissociation at 355 nm of the diazirine ring in mass-selected ions, and the cross-linked products were detected and quantified by collision-induced dissociation tandem mass spectrometry (CID-MSn, n = 3-5). Peptide scaffolds containing Ala and Leu residues with a C-terminal Gly gave 21-26% yields of cross-linked products, while the presence of the Pro and His residues decreased the yields. Experiments using hydrogen-deuterium-hydrogen exchange, carboxyl group blocking, and analysis of CID-MSn spectra of reference synthetic products revealed that a significant fraction of cross-links involved the Gly amide and carboxyl groups. Interpretation of the cross-linking results was aided by Born-Oppenheimer molecular dynamics (BOMD) and density functional theory calculations that allowed us to establish the protonation sites and conformations of the precursor ions. Analysis of long (100 ps) BOMD trajectories was used to count close contacts between the incipient carbene and peptide atoms, and the counting statistics was correlated with the results of gas-phase cross-linking.

For every retinylidene a force field: Using machine learning to improve molecular simulations of rhodopsins.

G-protein-coupled receptors (GPCRs) are a superfamily of transmembrane proteins that perform diverse signaling roles in organisms, including humans. This makes them an attractive target for pharmaceuticals, with 134 different GPCRs being currently targeted by roughly 700 drugs. Here we study the retinylidene cofactor of GPCR rhodopsin that is responsible for initiating phototransduction. There are many distinct retinal forms—isomers, rotamers, and protonation states—that all have unique electronic landscapes. However, current molecular dynamics force field libraries are incomplete in that they ignore these differences and treat each retinylidene form identically. It is unknown what effect this has on current molecular dynamics simulations of retinal proteins. Our work addresses this problem by using machine learning methods to generate unique force fields for each of forty retinylidene forms. This approach improved on traditional parameterization schemes by using Born-Oppenheimer molecular dynamics trajectories to sample the retinylidene energetic landscapes. To validate our methods, the new force fields were then used in all-atom molecular dynamics simulations of the entire rhodopsin-membrane-water system. From these simulations we were able to recover hydrogen relaxation rates that closely matched experimental solid-state deuterium (2H) NMR measurements. Comparison of force field parameters showed a small but—given the highly allosteric nature of rhodopsin—potentially significant difference between retinylidene forms. We expect this work will improve simulations of many different types of rhodopsins, and that it will lead to an increased appreciation of the complexity of the retinylidene energetic landscape. Work supported by US NIH (R01EY026041) and US NSF (CHE 1904125, MCB 1817862) and Ministry of Education and Science of the Russian Federation (Project FSRM-2020-0006).

Noncanonical Isomers of Nucleoside Cation Radicals: An Ab Initio Study of the Dark Matter of DNA Ionization.

Cation radicals of DNA nucleosides, 2'-deoxyadenosine, 2'-deoxyguanosine, 2'-deoxycytidine, and 2'-deoxythymidine, can exist in standard canonical forms or as noncanonical isomers in which the charge is introduced by protonation of the nucleobase, whereas the radical predominantly resides in the deoxyribose moiety. Density functional theory as well as correlated ab initio calculations with coupled clusters (CCSD(T)) that were extrapolated to the complete basis set limit showed that noncanonical nucleoside ion isomers were thermodynamically more stable than their canonical forms in both the gas phase and as water-solvated ions. This indicated the possibility of exothermic conversion of canonical to noncanonical forms. The noncanonical isomers were calculated to have very low adiabatic ion-electron recombination energies (REad) for the lowest-energy isomers 2'-deoxy-(N-3H)adenos-1'-yl (4.74 eV), 2'-deoxy-(N-7H)guanos-1'-yl (4.66 eV), 2'-deoxy-(N-3H)cytid-1'-yl (5.12 eV), and 2'-deoxy-5-methylene-(O-2H)uridine (5.24 eV). These were substantially lower than the REad value calculated for the canonical 2'-deoxyadenosine, 2'-deoxy guanosine, 2'-deoxy cytidine, and 2'-deoxy thymidine cation radicals, which were 7.82, 7.46, 8.14, and 8.20 eV, respectively, for the lowest-energy ion conformers of each type. Charge and spin distributions in noncovalent cation-radical dA⊂dT and dG⊂dC nucleoside pairs and dAT, dCT, dTC, and dGC dinucleotides were analyzed to elucidate the electronic structure of the cation radicals. Born-Oppenheimer molecular dynamics trajectory calculations of the dinucleotides and nucleoside pairs indicated rapid exothermic proton transfer from noncanonical T+· to A in both dAT+· and dA⊂dT+·, leading to charge and radical separation. Noncanonical T+· in dCT+· and dTC+· initiated rapid proton transfer to cytosine, whereas the canonical dCT+· dinucleotide ion retained the cation radical structure without isomerization. No spontaneous proton transfer was found in dGC+· and dG⊂dC+· containing canonical neutral and noncanonical ionized deoxycytidine.

Automatic Approach to Explore the Multireaction Mechanism for Medium-Sized Bimolecular Reactions via Collision Dynamics Simulations and Transition State Searches

We develop a broadly applicable computational method for the automatic exploration of the bimolecular multireaction mechanism. The current methodology mainly involves the high-energy Born-Oppenheimer molecular dynamics (BOMD) simulation and the successive reaction pathway construction. Several computational tricks are introduced, which include the selection of the reactive regions based on the electronic structure calculations and the employment of the virtual collision dynamics simulations with the monitoring of the atomic distance before the BOMD simulation. These prescreening steps largely reduce the number of trajectories in the BOMD simulations and significantly save the computational cost. The hidden Markov model combined with the modified atomic connectivity matrix is used for the detection of reaction events in each BOMD trajectory. Starting from several geometries close to the reaction events, the further intermediate optimization and transition state searches are conducted. The proposed method allows us to build the complicated multireaction mechanism of medium-sized bimolecular systems automatically. Here, we examine the feasibility and efficiency of the current method by its performance in searching the mechanisms of two prototype reactions in environmental science, which are the penicillin G anion + H2O and penicillin G anion + OH radical reactions. The result indicates that the proposed theoretical method is a powerful protocol for the automatic search of the bimolecular reaction mechanisms for medium-sized compounds.

First-principle study of the structures, growth pattern, and properties of (Pt3Cu)n, n = 1-9, clusters.

In this work, a first-principles systematic study of (Pt3Cu)n, n = 1-9, clusters was performed employing the linear combination of Gaussian-type orbital auxiliary density functional theory approach. The growth of the clusters has been achieved by increasing the previous cluster by one Pt3Cu unit at a time. To explore in detail the potential energy surface of these clusters, initial structures were obtained from Born-Oppenheimer molecular dynamics trajectories generated at different temperatures and spin multiplicities. For each cluster size, several dozens of structures were optimized without any constraints. The most stable structures were characterized by frequency analysis calculations. This study demonstrates that the obtained most stable structures prefer low spin multiplicities. To gain insight into the growing pattern of these systems, average bond lengths were calculated for the lowest stable structures. This work reveals that the Cu atoms prefer to be together and to localize inside the cluster structures. Moreover, these systems tend to form octahedra moieties in the size range of n going from 4 to 9 Pt3Cu units. Magnetic moment per atom and spin density plots were obtained for the neutral, cationic, and anionic ground state structures. Dissociation energies, ionization potential, and electron affinity were calculated, too. The dissociation energy and the electron affinity increase as the number of Pt3Cu units grows, whereas the ionization potential decreases.

Ab Initio QM/MM Simulation of Ferrocene Homogeneous Electron-Transfer Reaction.

Here, we demonstrate the feasibility of hybrid computational methods to predict the homogeneous electron exchange between the ferrocene and its oxidized (ferrocenium) state. The free energy for ferrocene oxidation was determined from thermodynamic cycles and implicit solvation strategies within density functional theory (DFT) methods leading to no more than 15% of deviation (in the range of 0.1-0.2 eV) when compared to absolute redox free energies obtained experimentally. Reorganization energy, as defined according to the Marcus theory of electron-transfer rate, was obtained by sampling the vertical ionization/electron affinity energies using hybrid quantum/classical (QM/MM) Born-Oppenheimer molecular dynamics trajectories. Calculated reorganization energies show a subtle but noteworthy dependence with the nature and the localization of the compensating countercharge. We concluded that the adopted hybrid computational strategy, to simulate homogeneous redox reactions, was successfully demonstrated and it further permits applications in more complex systems (required in daily life applications), where the electron transfer occurs heterogeneously.

Resolving the ultrafast intersystem crossing in a bimetallic platinum complex.

Bimetallic platinum complexes have interesting luminescent properties and feature long-lasting vibrational coherence and ultrafast intersystem crossing (ISC) after photoexcitation. Ultrafast triplet formation is driven by very strong spin-orbit coupling in these platinum (II) systems, where relativistic theoretical approaches beyond first-order perturbation theory are desirable. Using a fully variational relativistic theoretical method recently developed by the authors, we investigate the origins of ultrafast ISC in the [Pt(ppy) (μ-tBu2pz)]2 complex (ppy = phenylpyridine, pz = pyrazolate). Spin-orbit coupling values, evaluated along a Born-Oppenheimer molecular dynamics trajectory, are used to propagate electronic populations in time. Using this technique, we estimate ultrafast ISC rates of 15-134 fs in this species for the possible ISC pathways into the three low-lying triplet states.

Photodissociative Cross-Linking of Diazirine-Tagged Peptides with DNA Dinucleotides in the Gas Phase.

Non-covalent complexes of DNA dinucleotides dAA, dAT, dGG, dGC, and dCG with diazirine-tagged Cys-Ala-Gln-Lys peptides were generated as singly charged ions in the gas phase. Laser photodissociation at 355nm of the diazirine ring in the gas-phase complexes created carbene intermediates that underwent covalent cross-linking to the dinucleotides. The dinucleotides differed in the cross-linking yields, ranging from 27 to 36% for dAA and dAT up to 90-98% for dGG, dGC, and dCG. Collision-induced dissociation tandem mass spectrometry (CID-MS3) of the cross-linked conjugates revealed that fragmentation occurred chiefly in the dinucleotide moieties, resulting in a loss of a nucleobase and backbone cleavages. The CID-MS3 spectra further revealed that cross-links were primarily formed in the 3'-nucleotides for the dAT, dGC, and dCG combinations. Gas-phase and solution structures of dGG complexes with S-tagged CAQK were investigated by Born-Oppenheimer molecular dynamics (BOMD) and density functional theory calculations. The low free-energy complexes had zwitterionic structures in which the peptide was protonated at the N-terminus and in the Lys residue whereas the carboxyl or dGG phosphate were deprotonated, corresponding to the respective (Cys+, Lys+, COO-)+ and (Cys+, Lys+, phosphate-)+ protomeric types. Both types preferred structures in which the peptide N-terminal cysteine carrying the S-photo-tag was aligned with the 3'-guanine moiety. BOMD trajectories at 310K were analyzed for close contacts of the incipient peptide carbene with the positions in dGG that pointed to frequent contacts with the N-1, NH2, and N-7 atoms of 3'-guanine, in agreement with the cross-linking results. Carbene insertion to the guanine N-1-H and NH2 bonds was calculated by density functional and Møller-Plesset perturbational theory to be 350-380kJmol-1 exothermic. Based on calculations, we proposed a mechanism for the carbene reaction with guanine starting with an exothermic attack at N-7 to form a dipolar intermediate that can close an aziridine ring in another exothermic reaction, forming a stable covalent cross link. Graphical Abstract.

Photodissociative Cross-Linking of Non-covalent Peptide-Peptide Ion Complexes in the Gas Phase.

We report a gas-phase UV photodissociation study investigating non-covalent interactions between neutral hydrophobic pentapeptides and peptide ions incorporating a diazirine-tagged photoleucine residue. Phenylalanine (Phe) and proline (Pro) were chosen as the conformation-affecting residues that were incorporated into a small library of neutral pentapeptides. Gas-phase ion-molecule complexes of these peptides with photo-labeled pentapeptides were subjected to photodissociation. Selective photocleavage of the diazirine ring at 355nm formed short-lived carbene intermediates that underwent cross-linking by insertion into H-X bonds of the target peptide. The cross-link positions were established from collision-induced dissociation tandem mass spectra (CID-MS3) providing sequence information on the covalent adducts. Effects of the amino acid residue (Pro or Phe) and its position in the target peptide sequence were evaluated. For proline-containing peptides, interactions resulting in covalent cross-links in these complexes became more prominent as proline was moved towards the C-terminus of the target peptide sequence. The photocross-linking yields of phenylalanine-containing peptides depended on the position of both phenylalanine and photoleucine. Density functional theory calculations were used to assign structures of low-energy conformers of the (GLPMG + GLL*LK + H)+ complex. Born-Oppenheimer molecular dynamics trajectory calculations were used to capture the thermal motion in the complexes within 100ps and determine close contacts between the incipient carbene and the H-X bonds in the target peptide. This provided atomic-level resolution of potential cross-links that aided spectra interpretation and was in agreement with experimental data. Graphical Abstract ᅟ.

Ab initio studies on the photodissociation dynamics of the 1,1-difluoroethyl radical.

Born-Oppenheimer molecular dynamics trajectory calculations at the HCTH147/6-31G** level of theory simulate the dissociation dynamics of photolytically excited 1,1-difluoroethyl radicals. EOMCCSD/AUG-cc-pVDZ calculations show that an excitation energy of 94.82 kcal/mol is necessary to initiate photodissociation reactions. In contrast to photodissociation dynamics of ethyl radicals where a large discrepancy between actual dissociation rates and rates that are predicted by statistical rate theories, we find reaction rates of 5.1 × 1011 s-1 for the dissociation of an H atom, which is in perfect accord with what is predicted by Rice-Ramsperger-Kassel-Marcus (RRKM) calculations and there is no indication of any nonstatistical effects. However, our trajectory calculations show a much larger fraction of C-C bond breakage reaction of 56% occurring than that expected by RRKM (only 16%).

Aqueous Solvation of SmI2: A Born-Oppenheimer Molecular Dynamics Density Functional Theory Cluster Approach.

We report the results of Born-Oppenheimer molecular dynamics (BOMD) simulations on the aqueous solvation of the SmI2 molecule at room temperature using the cluster microsolvation approach including 32 water molecules. The electronic structure calculations were done using the M062X hybrid exchange-correlation functional in conjunction with the 6-31G** basis sets for oxygen and hydrogen. For the iodine and samarium atoms the Stuttgart-Köln relativistic effective-core potentials were utilized with their associated valence basis sets. Starting from the optimized geometry of SmI2 embeded in the microsolvation environment, we find a swift substitution of the iodine ions by eight tightly bound water molecules around Sm(II). Through the Sm-O radial distribution function and the evolution of the Sm-O distances, the present study predicts a first rigid Sm(II) solvation shell from 2.6 to 3.4 Å, whose integration leads to a coordination number of 8.4 water molecules, and a second softer solvation sphere from 3.5 to ca. 6 Å. The Sm(II)-O radial distribution function is in excellent agreement with that reported for Sr2+ from EXAFS studies, a fact that can be explained because Sr2+ and Sm2+ have almost identical ionic radii (ca. 1.26 Å) and coordination numbers: 8 for Sr2+ and 8.4 for Sm2+. The theoretical EXAFS spectrum was obtained from the BOMD trajectory and is discussed in the light of the experimental spectra for Sm(III). Once microsolvation is achieved, no water exchange events were found to occur around Sm2+, in agreement with the experimental data for Eu2+ (which has a nearly identical charge-to-ionic radius relation as Sm2+), where the mean residence time of a water molecule in [Eu(H2O)8]2+ is known to be ca. 230 ps.

Stability and Structures of Silver Subnanometer Clusters in EMT Zeolite with Maximum Aluminum Content

The stability of small silver clusters, Agn with n = 3, 4, 6, and 8, supported in a sodalite of EMT zeolite with a maximum possible aluminum content (Si/Al = 1) was studied with the density-functional-theory-based cluster approach. The results showed that various silver clusters, partially or completely reduced, can be stabilized in this zeolite with the maximum negative defects. Computed formation energies for two reduction reactions reveal that completely reduced Ag clusters can be stabilized through hydrogenation, giving some preference to the formation of hydrogenated Ag6 species. The hydrogenated neutral AgnHm/EMT clusters remain stable at T = 300 K, as obtained from the analysis of the Born–Oppenheimer molecular dynamics trajectories. Ion exchange of Na+ by Ag+ cations is obtained to be exothermic with energy gain of 0.5 eV per Ag cation. The formation of Ag–O ionic bonds, known to stabilize metal clusters in zeolitic frameworks, is found to be preponderant in the structures with fewer than six silv...

Efficient calculation of nuclear spin-rotation constants from auxiliary density functional theory.

The computation of the spin-rotation tensor within the framework of auxiliary density functional theory (ADFT) in combination with the gauge including atomic orbital (GIAO) scheme, to treat the gauge origin problem, is presented. For the spin-rotation tensor, the calculation of the magnetic shielding tensor represents the most demanding computational task. Employing the ADFT-GIAO methodology, the central processing unit time for the magnetic shielding tensor calculation can be dramatically reduced. In this work, the quality of spin-rotation constants obtained with the ADFT-GIAO methodology is compared with available experimental data as well as with other theoretical results at the Hartree-Fock and coupled-cluster level of theory. It is found that the agreement between the ADFT-GIAO results and the experiment is good and very similar to the ones obtained by the coupled-cluster single-doubles-perturbative triples-GIAO methodology. With the improved computational performance achieved, the computation of the spin-rotation tensors of large systems or along Born-Oppenheimer molecular dynamics trajectories becomes feasible in reasonable times. Three models of carbon fullerenes containing hundreds of atoms and thousands of basis functions are used for benchmarking the performance. Furthermore, a theoretical study of temperature effects on the structure and spin-rotation tensor of the H(12)C-(12)CH-DF complex is presented. Here, the temperature dependency of the spin-rotation tensor of the fluorine nucleus can be used to identify experimentally the so far unknown bent isomer of this complex. To the best of our knowledge this is the first time that temperature effects on the spin-rotation tensor are investigated.

Dynamic Molecular Graphs: “Hopping” Structures

This work aims to contribute to the discussion about the suitability of bond paths and bond-critical points as indicators of chemical bonding defined within the theoretical framework of the quantum theory of atoms in molecules. For this purpose, we consider the temporal evolution of the molecular structure of [Fe{C(CH2 )3 }(CO)3 ] throughout Born-Oppenheimer molecular dynamics (BOMD), which illustrates the changing behaviour of the molecular graph (MG) of an electronic system. Several MGs with significant lifespans are observed across the BOMD simulations. The bond paths between the trimethylenemethane and the metallic core are uninterruptedly formed and broken. This situation is reminiscent of a "hopping" ligand over the iron atom. The molecular graph wherein the bonding between trimethylenemethane and the iron atom takes place only by means of the tertiary carbon atom has the longest lifespan of all the considered structures, which is consistent with the MG found by X-ray diffraction experiments and quantum chemical calculations. In contrast, the η(4) complex predicted by molecular-orbital theory has an extremely brief lifetime. The lifespan of different molecular structures is related to bond descriptors on the basis of the topology of the electron density such as the ellipticities at the FeCH2 bond-critical points and electron delocalisation indices. This work also proposes the concept of a dynamic molecular graph composed of the different structures found throughout the BOMD trajectories in analogy to a resonance hybrid of Lewis structures. It is our hope that the notion of dynamic molecular graphs will prove useful in the discussion of electronic systems, in particular for those in which analysis on the basis of static structures leads to controversial conclusions.

Effect of molecular vibrations on the MD/QC-simulated absorption spectra

Effect of molecular vibrations on the absorption spectra simulated via a sequential approach combining molecular dynamics (MD) with quantum-chemical calculations has been investigated. Simulated spectra have been obtained from the time-dependent density functional theory results averaged over series of molecular geometries retrieved from Born–Oppenheimer MD trajectories. Distributions of bond lengths have been analyzed and related to the features of calculated spectra. For NVE simulations of small systems, absorption spectra exhibit bimodal bandshape as a result of classical treatment of vibrations. For NVE trajectories of larger systems or simulations in the NVT ensemble calculated absorption bands are symmetric, however, they may not agree with the results of Franck–Condon analysis. These results are practical manifestations of effects predicted theoretically from general principles. Consequences for the modeling of absorption spectra have been discussed. © 2013 Wiley Periodicals, Inc.

A density functional study of Rh13

A density functional study was performed for the Rh13 cluster using the linear combination of Gaussian-type orbitals density functional theory (LCGTO-DFT) approach. The calculations employed both the local density approximation (LDA) as well as the generalized gradient approximation (GGA) in combination with a quasi-relativistic effective core potential (QECP). Initial structures for the geometry optimization were taken along Born–Oppenheimer molecular dynamics (BOMD) trajectories. The BOMD trajectories were performed at different temperatures and considered different potential energy surfaces (PES). As a result, several hundred isomers of the Rh13 cluster in different spin multiplicities were optimized with the aim to determine the lowest energy structures. All geometry optimizations were performed without any symmetry restriction. A vibrational analysis was performed to characterize these isomers. Structural parameters, relative stability energy, harmonic frequencies, binding energy, and most relevant Kohn–Sham (KS) molecular orbitals are reported. The obtained results are compared with available data from the literature. This study predicts a low symmetry biplanarlike structure as the ground-state structure of Rh13 with 11 unpaired electrons. This isomer was first noticed by inspection of first-principle Born–Oppenheimer molecular dynamics (BOMD) simulations between 300 and 600 K. This represents the most extensive theoretical study on the ground-state structure of the Rh13 cluster and underlines the importance of BOMD simulations to fully explore the PES landscapes of complicated systems.

An extrapolation method for the efficient calculation of molecular response properties within Born–Oppenheimer molecular dynamics

The calculation of molecular response properties in dynamic molecular systems is a major challenge that requires sampling over many steps of, e.g., Born-Oppenheimer molecular dynamics (BO-MD) simulations. We present an extrapolation scheme to accelerate such calculations for multiple steps within BO-MD trajectories or equivalently within other sampling methods of conformational space. The extrapolation scheme is related to the one introduced by Pulay and Fogarasi [Chem. Phys. Lett., 2004, 386, 272] for self-consistent field (SCF) energy calculations. We extend the extrapolation to the quantities within our density matrix-based Laplace-transformed coupled perturbed SCF (DL-CPSCF) method that allows for linear-scaling calculations of response properties for large molecular systems. Here, we focus on the example of calculating NMR chemical shifts for which the number of required DL-CPSCF iterations reduces by roughly 40-70%.

Coupling quantum interpretative techniques: another look at chemical mechanisms in organic reactions.

A cross ELF-NCI analysis is tested over prototypical organic reactions. The synergetic use of ELF and NCI enables the understanding of reaction mechanisms since each method can respectively identify regions of strong and weak electron pairing. Chemically intuitive results are recovered and enriched by the identification of new features. Non covalent interactions are found to foresee the evolution of the reaction from the initial steps. Within NCI, no topological catastrophe is observed as changes are continuous to such an extent that future reaction steps can be predicted from the evolution of the initial NCI critical points. Indeed, strong convergences through the reaction paths between ELF and NCI critical points enable to identify key interactions at the origin of the bond formation. VMD scripts enabling the automatic generation of movies depicting the cross NCI/ELF analysis along a reaction path (or following a Born-Oppenheimer molecular dynamics trajectory) are provided as S.I.

Charge controlled changes in the cluster and spin dynamics of Sc3N@C80(CF3)2: the flexible spin density distribution and its impact on ESR spectra

Cluster and spin dynamics of a Sc(3)N@C(80)(CF(3))(2) derivative are studied by DFT in different charge states, from -3 to +1. For the neutral Sc(3)N@C(80)(CF(3))(2), static DFT computations of many cluster conformers as well as Born-Oppenheimer molecular dynamics (BOMD) show that addition of two CF(3) groups to Sc(3)N@C(80) significantly changes dynamics of the Sc(3)N cluster: instead of free rotation as in Sc(3)N@C(80), the cluster in Sc(3)N@C(80)(CF(3))(2) exhibits only hindered motions. Similar cluster dynamics is found in the mono- and trianions of Sc(3)N@C(80)(CF(3))(2), while free rotation of the cluster is found in the cation. In the radical species, motions of the cluster dramatically change spin-density distribution. Spin populations of the metal atoms and the carbon cage are followed along the BOMD trajectories to reveal the details of the spin-flow. (45)Sc ESR hyperfine coupling constants integrated over BOMD trajectories are found to be substantially different from the results of static DFT computations, which emphasizes that cluster dynamics should be taken into account for reliable predictions of spectroscopic properties.