Antonio Monari

Theoretical modeling of large molecular systems. Advances in the local self consistent field method for mixed quantum mechanics/molecular mechanics calculations.

Molecular mechanics methods can efficiently compute the macroscopic properties of a large molecular system but cannot represent the electronic changes that occur during a chemical reaction or an electronic transition. Quantum mechanical methods can accurately simulate these processes, but they require considerably greater computational resources. Because electronic changes typically occur in a limited part of the system, such as the solute in a molecular solution or the substrate within the active site of enzymatic reactions, researchers can limit the quantum computation to this part of the system. Researchers take into account the influence of the surroundings by embedding this quantum computation into a calculation of the whole system described at the molecular mechanical level, a strategy known as the mixed quantum mechanics/molecular mechanics (QM/MM) approach. The accuracy of this embedding varies according to the types of interactions included, whether they are purely mechanical or classically electrostatic. This embedding can also introduce the induced polarization of the surroundings. The difficulty in QM/MM calculations comes from the splitting of the system into two parts, which requires severing the chemical bonds that link the quantum mechanical subsystem to the classical subsystem. Typically, researchers replace the quantoclassical atoms, those at the boundary between the subsystems, with a monovalent link atom. For example, researchers might add a hydrogen atom when a C-C bond is cut. This Account describes another approach, the Local Self Consistent Field (LSCF), which was developed in our laboratory. LSCF links the quantum mechanical portion of the molecule to the classical portion using a strictly localized bond orbital extracted from a small model molecule for each bond. In this scenario, the quantoclassical atom has an apparent nuclear charge of +1. To achieve correct bond lengths and force constants, we must take into account the inner shell of the atom: for an sp(3) carbon atom, we consider the two core 1s electrons and treat that carbon as an atom with three electrons. This results in an LSCF+3 model. Similarly, a nitrogen atom with a lone pair of electrons available for conjugation is treated as an atom with five electrons (LSCF+5). This approach is particularly well suited to splitting peptide bonds and other bonds that include carbon or nitrogen atoms. To embed the induced polarization within the calculation, researchers must use a polarizable force field. However, because the parameters of the usual force fields include an average of the induction effects, researchers typically can obtain satisfactory results without explicitly introducing the polarization. When considering electronic transitions, researchers must take into account the changes in the electronic polarization. One approach is to simulate the electronic cloud of the surroundings by a continuum whose dielectric constant is equal to the square of the refractive index. This Electronic Response of the Surroundings (ERS) methodology allows researchers to model the changes in induced polarization easily. We illustrate this approach by modeling the electronic absorption of tryptophan in human serum albumin (HSA).

Toward a Quantitative Assessment of Electronic Transitions’ Charge-Transfer Character

We hereby report studies devoted to a topological descriptor of photoinduced electronic charge density variation. Our novel index, symbolized as ϕS, consists in the detachment and attachment densities overlap, where the detachment density physically depicts the electron density removed from the ground state of a molecule during the transition while the attachment density consists in the rearranged density in the excited state. Our method provides a simple and efficient way to quantitatively evaluate how easy the charge-separation is made upon the chromophores light absorption. Furthermore, this model can be applied for instance to address a comment on new push-pull dyes charge-transfer ability in order to assess their potentiality as candidates for light absorption-based devices. Moreover, the ϕS assessment allows us to perform some methodological diagnostic tests concerning the use of long-range corrected exchange-correlation functional in a time-dependent density functional theory (TDDFT) framework. This paper relates the ϕS descriptors mathematical foundations from various perspectives (detachment/attachment densities or natural transition orbitals), together with its application to several types of chromophores. Connections and divergences with a formerly proposed index are finally evidenced.

3,5-Bis(ethynyl)pyridine and 2,6-bis(ethynyl)pyridine spanning two Fe(Cp*)(dppe) units: role of the nitrogen atom on the electronic and magnetic couplings.

The role of the nitrogen atom on the electronic and magnetic couplings of the mono-oxidized and bi-oxidized pyridine-containing complex models [2,6-{Cp(dpe)Fe-C≡C-}(2)(NC(5)H(3))](n+) and [3,5-{Cp(dpe)Fe-C≡C-}(2)(NC(5)H(3))](n+) is theoretically tackled with the aid of density-functional theory (DFT) and multireference configuration interaction (MR-CI) calculations. Results are analyzed and compared to those obtained for the reference complex [1,3-{Cp*(dppe)Fe-C≡C-)}(2)(C(6)H(4))](n+). The mono-oxidized species show an interesting behavior at the borderline between spin localization and delocalization and one through-bond communication path among the two involving the central ring, is favored. Investigation of the spin state of the dicationic complexes indicates ferromagnetic coupling, which can differ in magnitude from one complex to the other. Very importantly, electronic and magnetic properties of these species strongly depend not only upon the location of the nitrogen atom in the ring versus that of the organometallic end-groups but also upon the architectural arrangement of one terminus, with respect to the other and/or vis-à-vis the central ring. To help validate the theoretical results, the related families of compounds [1,3-{Cp*(dppe)Fe-C≡C-)}(2)(C(6)H(4))](n+), [2,6-{Cp*(dppe)Fe-C≡C-}(2)(NC(5)H(3))](n+), [3,5-{Cp*(dppe)Fe-C≡C-}(2)(NC(5)H(3))](n+) (n = 0-2) were experimentally synthesized and characterized. Electrochemical, spectroscopic (infrared (IR), Mössbauer), electronic (near-infrared (NIR)), and magnetic properties (electron paramagnetic resonance (EPR), superconducting quantum interference device (SQUID)) are discussed and interpreted in the light of the theoretical data. The set of data obtained allows for many strong conclusions to be drawn. A N atom in the long branch increases the ferromagnetic interaction between the two Fe(III) spin carriers (J > 500 cm(-1)), whereas, when placed in the short branch, it dramatically reduces the magnetic exchange in the di-oxidized species (J = 2.14(5) cm(-1)). In the mixed-valence compounds, when the N atom is positioned on the long branch, the intermediate excited state is higher in energy than the different ground-state conformers and the relaxation process provides exclusively the Fe(II)/Fe(III) localized system (H(ab) ≠ 0). Positioning the N atom on the short branch modifies the energy profile and the diabatic mediating state lies just above the reactant and product diabatic states. Consequently, the LMCT transition becomes less energetic than the MMCT transition. Here, the direct coupling does not occur (H(ab) = 0) and only the coupling through the bridge (c) and the reactant (a) and product (b) diabatic states is operating (H(ac) = H(bc) ≠ 0).

Spectral properties of polypyridyl ruthenium complexes intercalated in DNA: theoretical insights into the surrounding effects of [Ru(dppz)(bpy)2]2+.

The UV/Visible absorption properties of a polypyridyl ruthenium complex upon intercalation on DNA are studied at the mixed quantum mechanics molecular mechanics level of theory. Vertical excitation transitions are computed by time dependent density functional theory. Particular emphasis is put on the different levels at which the macromolecular environment is treated, and in particular on the analysis of the effect of mechanical, electrostatic and polarizable embedding. We show that with the highest level of theory the experimental absorption wavelengths are reproduced with a difference of only 2 or 3 nm for the low energy bands. The systematic analysis of the individual vertical transitions allows us to get much more insights into the role played by the environment, in particular, in metal to ligand and intra ligand charge transfer transitions that can lead to the production of DNA oxidative lesions exploitable in phototherapy.

Theoretical study of new ruthenium-based dyes for dye-sensitized solar cells.

Two relevant, recently reported, ruthenium-based complexes to be used as sensitizers in Grätzel photovoltaic cells are theoretically studied. The UV/vis absorption spectra have been computed within the time-dependent density functional theory formalism. The obtained excitation energies are compared with the experimental results, and the nature of the transition is analyzed in terms of the electronic density. A preliminary study on the performance of different functionals against the equation of motion coupled cluster is performed on a smaller model system.

Improved Treatment of Surrounding Effects: UV/vis Absorption Properties of a Solvated Ru(II) Complex.

The UV/vis and circular-dichroism spectra of a bis-bipyridinyl ruthenium complex are computed at the density functional theory level and the time dependent density functional level of theory. The effects of the solvent, here water, have been taken into account, by polarizable continuum methods and by a hybrid quantum-mechanics/molecular-mechanics approach combined with molecular dynamics. The effects of the solvent have been decomposed in geometric, electrostatic, and polarization of the environment. The principal transitions have been analyzed by means of natural transition orbitals.

A QM/MM Study of the Absorption Spectrum of Harmane in Water Solution and Interacting with DNA: The Crucial Role of Dynamic Effects

We present a time-dependent density functional theory computation of the absorption spectra of one β-carboline system: the harmane molecule in its neutral and cationic forms. The spectra are computed in aqueous solution. The interaction of cationic harmane with DNA is also studied. In particular, the use of hybrid quantum mechanics/molecular mechanics methods is discussed, together with its coupling to a molecular dynamics strategy to take into account dynamic effects of the environment and the vibrational degrees of freedom of the chromophore. Different levels of treatment of the environment are addressed starting from purely mechanical embedding to electrostatic and polarizable embedding. We show that a static description of the spectrum based on equilibrium geometry only is unable to give a correct agreement with experimental results, and dynamic effects need to be taken into account. The presence of two stable noncovalent interaction modes between harmane and DNA is also presented, as well as the associated absorption spectrum of harmane cation.

Interaction of Palmatine with DNA: An Environmentally Controlled Phototherapy Drug

Palmatine is one of the four main protoberberine alkaloids and is largely employed in pharmacy and medicine as a versatile drug with considerable biological activities. More recently, palmatine has been proposed as a promising DNA phototherapy drug, notably due to its ability to produce in situ singlet oxygen only when interacting with DNA. The fine mechanisms of palmatine-DNA interactions as well as its complicated photophysics are not yet fully understood. In this paper, we identify via molecular dynamic techniques two stable interaction modes between palmatine and B-DNA, namely insertion and minor groove binding, whose structural and electronic bases are analyzed and rationalized. These two competitive modes share the same UV-vis signature and estimated binding free energies, and thus they may indeed coexist. By using hybrid quantum mechanics/molecular mechanics protocols coupled to molecular dynamics, we analyze palmatine excited state properties in water solution and in interaction with DNA. The environmentally controlled production of singlet oxygen is thus rationalized in terms of the competition between local and charge-transfer excited states.

A theoretical study of linear beryllium chains: Full configuration interaction

We present a full configuration interaction study of Be(N) (N=2,3,4,5) linear chains. A comparative study of the basis-set effect on the reproduction of the energy profile has been reported. In particular, the 3s1p, 4s2p, 4s2p1d, 5s3p2d, and 5s3p2d1f bases were selected. For the smallest chains (i.e., Be(2) and Be(3)), smaller basis sets give dissociative energy profiles, so large basis set is demanded for the reproduction of equilibrium minima in the structures. For Be(4) and Be(5) linear chains, the energy profiles show a minimum also by using the smallest basis sets, but the largest ones give a much stronger stabilization energy. For all the structures, two spin states have been studied: the singlet and the triplet. It is shown that the energy separation of the two states, in the equilibrium region, is small and decays exponentially with respect to the number of atoms in the chain. Finally an interpolative technique allowing for the estimation of the long-chain parameters from shorter ones is presented.

Resolving the Benzophenone DNA-Photosensitization Mechanism at QM/MM Level.

Benzophenone, the parent of the diarylketone family, is a versatile compound commonly used as a UV blocker. It may also trigger triplet-based DNA photosensitization. Therefore, benzophenone is involved in DNA photodamage induction. In the absence of experimentally resolved structure, the mechanism of DNA damage production remains elusive. Employing a hybrid quantum mechanics/molecular mechanics approach, here we address the spin transfer mechanism between this drug and proximal thymine, that is, the DNA nucleobase most prone to suffer triplet damages.