Cleavage of an Aromatic C-C Bond in Ferrocene by Insertion of an Iridium Nitrido Nitrogen Atom.

The intermolecular cleavage of C-C bonds is a rare event. Herein, we report on a late transition-metal terminal nitrido complex, which upon oxidation undergoes insertion of the nitrido nitrogen atom into the aromatic C-C bond of ferrocene. This reaction path was confirmed through 15N and deuterium isotope labeling experiments of the nitrido complex and ferrocenium, respectively. Cyclic voltammetry and UV/vis spectroscopy monitoring of the reaction revealed that oxidation is the initial step, yielding the tentative radical cationic nitrido complex, which is experimentally supported by extended X and Q-band electron paramagnetic resonance (EPR) and ENDOR, UV/vis, vT 1H NMR, and vibrational spectroscopic data. Density functional theory (DFT) and multireference calculations of this highly reactive intermediate revealed an S = 1/2 ground state. The high reactivity can be traced to the increased electrophilicity in the oxidized complex. Based on high-level PNO-UCCSD(T) calculations and UV/vis kinetic measurements, it is proposed that the reaction proceeds by initial electrophilic exo attack of the nitrido nitrogen atom at the cyclopentadienyl ring and consecutive ring expansion to a pyridine ring.

Three port logic gate using forward volume spin wave interference in a thin yttrium iron garnet film.

We demonstrate a logic gate based on interference of forward volume spin waves (FVSWs) propagating in a 54 nm thick, 100 µm wide yttrium iron garnet waveguide grown epitaxially on a garnet substrate. Two FVSWs injected by coplanar waveguides were made to interfere constructively and destructively by varying their phase difference, showing an XNOR logic function. The reflected and resonant waves generated at the edges of the waveguide were suppressed using spin wave absorbers. The observed isolation ratio was 19 dB for a magnetic field of ~2.80 kOe ( = 223 kA m-1) applied perpendicular to the film. The wavelength and device length were ~8.9 µm and ~53 µm, respectively. Further, the interference state of the SWs was analyzed using three-dimensional radio frequency simulations.

Multistate Multiconfiguration Quantum Chemical Computation of the Two-Photon Absorption Spectra of Bovine Rhodopsin.

Recently, progress in IR sources has led to the discovery that humans can detect infrared (IR) light. This is hypothesized to be due to the two-photon absorption (TPA) events promoting the retina dim-light rod photoreceptor rhodopsin to the same excited state populated via one-photon absorption (OPA). Here, we combine quantum mechanics/molecular mechanics and extended multiconfiguration quasi-degenerate perturbation theory calculations to simulate the TPA spectrum of bovine rhodopsin (Rh) as a model for the human photoreceptor. The results show that the TPA spectrum of Rh has an intense S0 → S1 band but shows also S0 → S2 and S0 → S3 transitions whose intensities, relative to the S0 → S1 band, are significantly increased when compared to the corresponding bands of the OPA spectrum. In conclusion, we show that IR light in the 950 nm region can be perceived by rod photoreceptors, thus supporting the two-photon origin of the IR perception. We also found that the same photoreceptor can perceive red (i.e., close to 680 nm) light provided that TPA induces population of S2.

The role of Snell's law for a magnonic majority gate.

In the fifty years since the postulation of Moore's Law, the increasing energy consumption in silicon electronics has motivated research into emerging devices. An attractive research direction is processing information via the phase of spin waves within magnonic-logic circuits, which function without charge transport and the accompanying heat generation. The functional completeness of magnonic logic circuits based on the majority function was recently proved. However, the performance of such logic circuits was rather poor due to the difficulty of controlling spin waves in the input junction of the waveguides. Here, we show how Snell's law describes the propagation of spin waves in the junction of a Ψ-shaped magnonic majority gate composed of yttrium iron garnet with a partially metallized surface. Based on the analysis, we propose a magnonic counterpart of a core-cladding waveguide to control the wave propagation in the junction. This study has therefore experimentally demonstrated a fundamental building block of a magnonic logic circuit.

Demonstration of a robust magnonic spin wave interferometer.

Magnonics is an emerging field dealing with ultralow power consumption logic circuits, in which the flow of spin waves, rather than electric charges, transmits and processes information. Waves, including spin waves, excel at encoding information via their phase using interference. This enables a number of inputs to be processed in one device, which offers the promise of multi-input multi-output logic gates. To realize such an integrated device, it is essential to demonstrate spin wave interferometers using spatially isotropic spin waves with high operational stability. However, spin wave reflection at the waveguide edge has previously limited the stability of interfering waves, precluding the use of isotropic spin waves, i.e., forward volume waves. Here, a spin wave absorber is demonstrated comprising a yttrium iron garnet waveguide partially covered by gold. This device is shown experimentally to be a robust spin wave interferometer using the forward volume mode, with a large ON/OFF isolation value of 13.7 dB even in magnetic fields over 30 Oe.

Communication: An efficient approach to compute state-specific nuclear gradients for a generic state-averaged multi-configuration self consistent field wavefunction.

We present a new, very efficient semi-numerical approach for the computation of state-specific nuclear gradients of a generic state-averaged multi-configuration self consistent field wavefunction. Our approach eliminates the costly coupled-perturbed multi-configuration Hartree-Fock step as well as the associated integral transformation stage. The details of the implementation within the Firefly quantum chemistry package are discussed and several sample applications are given. The new approach is routinely applicable to geometry optimization of molecular systems with 1000+ basis functions using a standalone multi-core workstation.

Shape of Multireference, Equation-of-Motion Coupled-Cluster, and Density Functional Theory Potential Energy Surfaces at a Conical Intersection.

We report and characterize ground-state and excited-state potential energy profiles using a variety of electronic structure methods along a loop lying on the branching plane associated with a conical intersection (CI) of a reduced retinal model, the penta-2,4-dieniminium cation (PSB3). Whereas the performance of the equation-of-motion coupled-cluster, density functional theory, and multireference methods had been tested along the excited- and ground-state paths of PSB3 in our earlier work, the ability of these methods to correctly describe the potential energy surface shape along a CI branching plane has not yet been investigated. This is the focus of the present contribution. We find, in agreement with earlier studies by others, that standard time-dependent DFT (TDDFT) does not yield the correct two-dimensional (i.e., conical) crossing along the branching plane but rather a one-dimensional (i.e., linear) crossing along the same plane. The same type of behavior is found for SS-CASPT2(IPEA=0), SS-CASPT2(IPEA=0.25), spin-projected SF-TDDFT, EOM-SF-CCSD, and, finally, for the reference MRCISD+Q method. In contrast, we found that MRCISD, CASSCF, MS-CASPT2(IPEA=0), MS-CASPT2(IPEA=0.25), XMCQDPT2, QD-NEVPT2, non-spin-projected SF-TDDFT, and SI-SA-REKS yield the expected conical crossing. To assess the effect of the different crossing topologies (i.e., linear or conical) on the PSB3 photoisomerization efficiency, we discuss the results of 100 semiclassical trajectories computed by CASSCF and SS-CASPT2(IPEA=0.25) for a PSB3 derivative. We show that for the same initial conditions, the two methods yield similar dynamics leading to isomerization quantum yields that differ by only a few percent.

C20H4(C4F8)3: a fluorine-containing annulated corannulene that is a better electron acceptor than C60.

At sixes and sevens: The reaction of corannulene with 35 equivalents of 1,4-C4F8I2 is an efficient and a relatively selective process that yields two main products in which six H atoms are substituted with three C4F8 moieties that form six- and seven-membered rings. Low-temperature photoelectron spectroscopy showed the electron affinity of the major isomer (shown) exceeds that of C60 (2.74±0.02 and 2.689±0.008 eV, respectively).

Mapping the Excited State Potential Energy Surface of a Retinal Chromophore Model with Multireference and Equation-of-Motion Coupled-Cluster Methods.

The photoisomerization of the retinal chromophore of visual pigments proceeds along a complex reaction coordinate on a multidimensional surface that comprises a hydrogen-out-of-plane (HOOP) coordinate, a bond length alternation (BLA) coordinate, a single bond torsion and, finally, the reactive double bond torsion. These degrees of freedom are coupled with changes in the electronic structure of the chromophore and, therefore, the computational investigation of the photochemistry of such systems requires the use of a methodology capable of describing electronic structure changes along all those coordinates. Here, we employ the penta-2,4-dieniminium (PSB3) cation as a minimal model of the retinal chromophore of visual pigments and compare its excited state isomerization paths at the CASSCF and CASPT2 levels of theory. These paths connect the cis isomer and the trans isomer of PSB3 with two structurally and energetically distinct conical intersections (CIs) that belong to the same intersection space. MRCISD+Q energy profiles along these paths provide benchmark values against which other ab initio methods are validated. Accordingly, we compare the energy profiles of MRPT2 methods (CASPT2, QD-NEVPT2, and XMCQDPT2) and EOM-SF-CC methods (EOM-SF-CCSD and EOM-SF-CCSD(dT)) to the MRCISD+Q reference profiles. We find that the paths produced with CASSCF and CASPT2 are topologically and energetically different, partially due to the existence of a "locally excited" region on the CASPT2 excited state near the Franck-Condon point that is absent in CASSCF and that involves a single bond, rather than double bond, torsion. We also find that MRPT2 methods as well as EOM-SF-CCSD(dT) are capable of quantitatively describing the processes involved in the photoisomerization of systems like PSB3.

Phase transitions, magnetotransport and magnetocaloric effects in a new family of quaternary Ni-Mn-In-Z Heusler alloys.

The magnetic, magnetotransport, and magnetocaloric properties near compound phase transitions in Ni50Mn35In14Z (Z = In, Ge, Al), and Ni48Co2Mn35In15 Heusler alloys have been studied using VSM and SQUID magnetometers (at magnetic fields (H) up to 5 T), four-probe method (at H = 0.005-1.5 T), and an adiabatic magnetocalorimeter (for H changes up to deltaH = 1.8 T), respectively. The martensitic transformation (MT) is accompanied by large magnetoresistance (up to 70%), a significant change in resistivity (up to 200%), and a sign reversal of the ordinary Hall effect coefficient, all related to a strong change in the electronic spectrum at the MT. The field dependences of the Hall resistance are complex in the vicinity of the MT, indicating a change in the relative concentrations of the austenite and martensite phases at strong fields. Negative and positive changes in adiabatic temperatures of about -2 K and +2 K have been observed in the vicinity of MT and Curie temperatures, respectively, for deltaH = 1.8 T.

Dynamic Electron Correlation Effects on the Ground State Potential Energy Surface of a Retinal Chromophore Model.

The ground state potential energy surface of the retinal chromophore of visual pigments (e.g., bovine rhodopsin) features a low-lying conical intersection surrounded by regions with variable charge-transfer and diradical electronic structures. This implies that dynamic electron correlation may have a large effect on the shape of the force fields driving its reactivity. To investigate this effect, we focus on mapping the potential energy for three paths located along the ground state CASSCF potential energy surface of the penta-2,4-dieniminium cation taken as a minimal model of the retinal chromophore. The first path spans the bond length alternation coordinate and intercepts a conical intersection point. The other two are minimum energy paths along two distinct but kinetically competitive thermal isomerization coordinates. We show that the effect of introducing the missing dynamic electron correlation variationally (with MRCISD) and perturbatively (with the CASPT2, NEVPT2, and XMCQDPT2 methods) leads, invariably, to a stabilization of the regions with charge transfer character and to a significant reshaping of the reference CASSCF potential energy surface and suggesting a change in the dominating isomerization mechanism. The possible impact of such a correction on the photoisomerization of the retinal chromophore is discussed.

Extended multi-configuration quasi-degenerate perturbation theory: the new approach to multi-state multi-reference perturbation theory.

The distinctive desirable features, both mathematically and physically meaningful, for all partially contracted multi-state multi-reference perturbation theories (MS-MR-PT) are explicitly formulated. The original approach to MS-MR-PT theory, called extended multi-configuration quasi-degenerate perturbation theory (XMCQDPT), having most, if not all, of the desirable properties is introduced. The new method is applied at the second order of perturbation theory (XMCQDPT2) to the 1(1)A(')-2(1)A(') conical intersection in allene molecule, the avoided crossing in LiF molecule, and the 1(1)A(1) to 2(1)A(1) electronic transition in cis-1,3-butadiene. The new theory has several advantages compared to those of well-established approaches, such as second order multi-configuration quasi-degenerate perturbation theory and multi-state-second order complete active space perturbation theory. The analysis of the prevalent approaches to the MS-MR-PT theory performed within the framework of the XMCQDPT theory unveils the origin of their common inherent problems. We describe the efficient implementation strategy that makes XMCQDPT2 an especially useful general-purpose tool in the high-level modeling of small to large molecular systems.

Photoisomerization around a fulvene double bond: coherent population transfer to the electronic ground state?

Photoisomerization around a central fulvene-type double bond is known to proceed through a conical intersection at the perpendicular geometry. The process is studied with an indenylidene-dihydropyridine model compound, allowing the use of visible excitation pulses. Transient absorption shows that 1) stimulated emission shifts to the red and loses oscillator strength on a 50 fs timescale, and 2) bleach recovery is highly nonexponential and not affected by solvent viscosity or methyl substitution at the dihydropyridine ring. Quantum-chemical calculations are used to explain point 1 as a result of initial elongation of the central C=C bond with mixing of S(2) and S(1) states. From point 2 it is concluded that internal conversion of S(1)→S(0) does not require torsional motion to the fully perpendicular state. The S(1) population appears to encounter a sink on the torsional coordinate before the conical intersection is reached. Rate equations cannot model the observed ground-state recovery adequately. Instead the dynamics are best described with a strongly damped oscillatory contribution, which could indicate coherent S(1)-S(0) population transfer.

Gas phase absorption studies of photoactive yellow protein chromophore derivatives.

Photoabsorption spectra of deprotonated trans p-coumaric acid and two of its methyl substituted derivatives have been studied in gas phase both experimentally and theoretically. We have focused on the spectroscopic effect of the location of the two possible deprotonation sites on the trans p-coumaric acid which originate to either a phenoxide or a carboxylate. Surprisingly, the three chromophores were found to have the same absorption maximum at 430 nm, in spite of having different deprotonation positions. However, the absorption of the chromophore in polar solution is substantially different for the distinct deprotonation locations. We also report on the time scales and pathways of relaxation after photoexcitation for the three photoactive yellow protein chromophore derivatives. As a result of these experiments, we could detect the phenoxide isomer within the deprotonated trans p-coumaric acid in gas phase; however, the occurrence of the carboxylate is uncertain. Several computational methods were used simultaneously to provide insights and assistance in the interpretation of our experimental results. The calculated excitation energies S(0)-S(1) are in good agreement with experiment for those systems having a negative charge on a phenoxide moiety. Although our augmented multiconfigurational quasidegenerate perturbation theory calculations agree with experiment in the description of the absorption spectrum of anions with a carboxylate functional group, there are some puzzling disagreements between experiment and some calculational methods in the description of these systems.

Modeling photoabsorption of the asFP595 chromophore.

The fluorescent protein asFP595 is a promising photoswitchable biomarker for studying processes in living cells. We present the results of a high level theoretical study of photoabsorption properties of the model asFP595 chromophore molecule in biologically relevant protonation states: anionic, zwitterionic, and neutral. Ground state equilibrium geometry parameters are optimized in the PBE0/(aug)-cc-pVDZ density functional theory approximation. An augmented version of multiconfigurational quasidegenerate perturbation theory (aug-MCQDPT2) following the state-averaged CASSCF/(aug)-cc-pVDZ calculations is used to estimate the vertical S0-S1 excitation energies for all chromophore species. An accuracy of this approach is validated by comparing the computed estimates of the S0-S1 absorption maximum of the closely related chromophore from the DsRed protein to the known experimental value in the gas phase. An influence of the CASSCF active space on the aug-MCQDPT2 excitation energies is analyzed. The zwitterionic form of the asFP595 chromophore is found to be the most sensitive to a particular choice and amount of active orbitals. This observation is explained by the charge-transfer type of the S0-S1 transition involving the entire conjugated pi-electron system for the zwitterionic protonation state. According to the calculation results, the anionic form in the trans conformation is found to possess the most red-shifted absorption band with the maximum located at 543 nm. The bands of the zwitterionic and neutral forms are considerably blue-shifted compared to those of the anionic form. These conclusions are at variance with the results obtained in the TDDFT approximation for the asFP595 chromophore. The absorption wavelengths computed in the aug-MCQDPT2/CASSCF theory are as follows: 543 (535), 470 (476), and 415 (417) nm for the anionic, zwitterionic, and neutral forms of the trans and cis (in parentheses) isomers of the asFP595 chromophore. These data can be used as a reference for further theoretical studies of the asFP595 chromophore in different media and for modeling photoabsorption properties of the asFP595 fluorescent protein.

An opsin shift in rhodopsin: retinal S0-S1 excitation in protein, in solution, and in the gas phase.

We considered a series of model systems for treating the photoabsorption of the 11-cis retinal chromophore in the protonated Schiff-base form in vacuum, solutions, and the protein environment. A high computational level, including the quantum mechanical-molecular mechanical (QM/MM) approach for solution and protein was utilized in simulations. The S0-S1 excitation energies in quantum subsystems were evaluated by means of an augmented version of the multiconfigurational quasidegenerate perturbation theory (aug-MCQDPT2) with the ground-state geometry parameters optimized in the density functional theory PBE0/cc-pVDZ approximation. The computed positions of absorption bands lambdamax, 599(g), 448(s), and 515(p) nm for the gas phase, solution, and protein, respectively, are in excellent agreement with the corresponding experimental data, 610(g), 445(s), and 500(p) nm. Such consistency provides a support for the formulated qualitative conclusions on the role of the chromophore geometry, environmental electrostatic field, and the counterion in different media. An essentially nonplanar geometry conformation of the chromophore group in the region of the C14-C15 bond was obtained for the protein, in particular, owing to the presence of the neighboring charged amino acid residue Glu181. Nonplanarity of the C14-C15 bond region along with the influence of the negatively charged counterions Glu181 and Glu113 are found to be important to reproduce the spectroscopic features of retinal chromophore inside the Rh cavity. Furthermore, the protein field is responsible for the largest bond-order decrease at the C11-C12 double bond upon excitation, which may be the reason for the 11-cis photoisomerization specificity.

Anomalous Faraday effect of a system with extraordinary optical transmittance.

It is shown theoretically that the Faraday rotation becomes anomalously large and exhibits extraordinary behavior near the frequencies of the extraordinary optical transmittance through optically thick perforated metal film with holes filled with a magneto-optically active material. This phenomenon is explained as result of strong confinement of the evanescent electromagnetic field within magnetic material, which occurs due to excitation of the coupled plasmon-polaritons on the opposite surfaces of the film.