Search NASA⌕ Search

SEARCH · Search NASA

Results for “CASSCF”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

Internal Conversion between Bright (1 1 B $^+_u$) and Dark (21 A $^-_g$) States in s- trans -Butadiene and s- trans -Hexatriene

Internal conversion (IC) between the two lowest singlet excited states, 1 1 B$_u^+$ and 2 1 A$_g^–$, of s-trans-butadiene and s-trans-hexatriene is investigated using a series of single- and multi- reference wave function and density functional theory (DFT) methodologies. Three independent types of the equation-of-motion coupled-cluster (EOMCC) theory capable of providing an accurate and balanced description of one- as well as two-electron transitions, abbreviated as δ-CR-EOMCC(2,3), DIP-EOMCC(4h2p){N o }, and DEA-EOMCC(4p2h){N u } or DEA-EOMCC(3p1h,4p2h){N u }, consistently predict that the 1 1 B$_u^+$/2 1 A$_g^–$ crossing in both molecules occurs along the bond length alternation coordinate. However, the analogous 1 1 B$_u^+$ and 2 1 A$_g^–$ potentials obtained with some multireference approaches, such as CASSCF and MRCIS(D), as well as with the linear-response formulation of time-dependent DFT (TDDFT), do not cross. Hence, caution needs to be exercised when studying the low-lying singlet excited states of polyenes with conventional multiconfigurational methods and TDDFT. The multistate many-body perturbation theory methods, such as XMCQDPT 2 , do correctly reproduce the curve crossing. Among the simplest and least expensive computational methodologies, the DFT approaches that incorporate the contributions of doubly excited configurations, abbreviated as MRSF (mixed reference spin-flip) TDDFT and SSR( 4,4 ), accurately reproduce our best EOMCC results. This is highly promising for nonadiabatic molecular dynamics simulations in larger systems.

36 MATERIALS SCIENCE↗

Excited States of Crystalline Point Defects with Multireference Density Matrix Embedding Theory

Accurate and affordable methods to characterize the electronic structure of solids are important for targeted materials design. Embedding-based methods provide an appealing balance in the trade-off between cost and accuracy-particularly when studying localized phenomena. Here, we use the density matrix embedding theory (DMET) algorithm to study the electronic excitations in solid-state defects with a restricted open-shell Hartree-Fock (ROHF) bath and multireference impurity solvers, specifically, complete active space self-consistent field (CASSCF) and n-electron valence state second-order perturbation theory (NEVPT2). In this work, we apply the method to investigate the electronic excitations in an oxygen vacancy (OV) on a MgO(100) surface and find absolute deviations within 0.05 eV between DMET using the CASSCF/NEVPT2 solver, denoted as CAS-DMET/NEVPT2-DMET, and the nonembedded CASSCF/NEVPT2 approach. Next, we establish the practicality of DMET by extending it to larger supercells for the OV defect and a neutral silicon vacancy in diamond where the use of nonembedded CASSCF/NEVPT2 is extremely expensive.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Local Excitations of a Charged Nitrogen Vacancy in Diamond with Multireference Density Matrix Embedding Theory

Here, we investigate the negatively charged nitrogen-vacancy center in diamond using periodic density matrix embedding theory (pDMET). To describe the strongly correlated excited states of this system, the complete active space self-consistent field (CASSCF) followed by n-electron valence state second-order perturbation theory (NEVPT2) was used as the impurity solver. Since the NEVPT2-DMET energies show a linear dependence on the inverse of the size of the embedding subspace, we performed an extrapolation of the excitation energies to the nonembedding limit using a linear regression. The extrapolated NEVPT2-DMET first triplet–triplet excitation energy is 2.31 eV and that for the optically inactive singlet–singlet transition is 1.02 eV, both in agreement with the experimentally observed vertical excitation energies of ∼2.18 eV and ∼1.26 eV, respectively. This is the first application of pDMET to a charged periodic system and the first investigation of the NV – defect using NEVPT2 for periodic supercell models.

36 MATERIALS SCIENCE↗

Exploring Ligand-to-Metal Charge-Transfer States in the Photo-Ferrioxalate System Using Excited-State Specific Optimization

The photo-ferrioxalate system (PFS), [Fe(III)(C 2 O 4 )] 3– , more than an exact chemical actinometer, has been extensively applied in wastewater and environment treatment. Despite many experimental efforts to improve clarity, important aspects of the mechanism of ferrioxalate photolysis are still under debate. In this paper, we employ the recently developed WΓ-CASSCF to investigate the ligand-to-metal charge-transfer states that are key to ferrioxalate photolysis. This investigation provides a qualitative picture of these states and key potential energy surface features related to the photolysis. Furthermore, our theoretical results are consistent with the prompt charge-transfer picture seen in recent experiments and clarify some features that are not visible in experiments. Two ligand-to-metal charge-transfer states contribute to the photolysis of ferrioxalate, and the avoided crossing barrier between them is low compared with the initial photoexcitation energy. Our data also clarify that one Fe–O bond cleaves first, followed by the C–C bond and the other Fe–O bond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spectroscopic, Magnetic, and Computational Investigations on a Series of Rhenium(III) Cyclopentadienide β-diketiminate Halide and Pseudohalide Complexes

Here the low-valent rhenium(I) salt Na[Re(η 5 -Cp)(BDI)] (BDI = N,N'-bis(2,6-diisopropylphenyl)-3,5-dimethyl-β-diketiminate) was shown to react with various halide and chalcogenolate reagents, leading to the isolation of a series of rhenium(III)-halide and -pseudohalide complexes: Re(X)(η 5 -Cp)(BDI) (1-X, X = F, Cl, Br, I; 2, X = OTf) and Re(ER)(η 5 -Cp)(BDI) (3-ER, ER = SBn, SeBn, TePh). The 1 H NMR spectra of these complexes displayed sharp resonances shifted several ppm from typical diamagnetic regions, as well as distinct chemical shift trends down both the halide (with the exception of F) and chalcogenolate series, with both Cp and BDI backbone peaks moving downfield with increasing atomic number. Subsequent magnetic susceptibility measurements of the rhenium(III) halides 1-Cl, 1-Br, and 1-I indicated that these complexes display temperature-independent paramagnetism (TIP), with χ TIP values of 7.41(44) × 10 -4 to 1.50(51) × 10 -3 cm 3 mol -1 . Multireference complete active space self-consistent field (CASSCF) computations incorporating spin-orbit state mixing revealed small energy separations (1.7-3.0 kcal/mol) between thermally isolated ground states and the first excited spin-orbit states for the rhenium halides, confirming that low-lying excited states are responsible for the observed TIP behavior.

36 MATERIALS SCIENCE↗

Elucidating the Role of Halides and Iron during Radiolysis-Driven Oxidative Etching of Gold Nanocrystals Using Liquid Cell Transmission Electron Microscopy and Pulse Radiolysis

Graphene liquid cell transmission electron microscopy (TEM) has enabled the observation of a variety of nanoscale transformations. Yet understanding the chemistry of the liquid cell solution and its impact on the observed transformations remains an important step toward translating insights from liquid cell TEM to benchtop chemistry. Gold nanocrystal etching can be used as a model system to probe the reactivity of the solution. FeCl 3 has been widely used to promote gold oxidation in bulk and liquid cell TEM studies, but the roles of the halide and iron species have not been fully elucidated. In this work, we observed the etching trajectories of gold nanocrystals in different iron halide solutions. We observed an increase in gold nanocrystal etch rate going from Cl - - to Br - - to I - -containing solutions. Furthermore, this is consistent with a mechanism in which the dominant role of halides is as complexation agents for oxidized gold species. Additionally, the mechanism through which FeCl 3 induces etching in liquid cell TEM remains unclear. Ground-state bleaching of the Fe(III) absorption band observed through pulse radiolysis indicates that iron may react with Cl 2 · - radicals to form an oxidized transient species under irradiation. Complete active space self-consistent field (CASSCF) calculations indicate that the FeCl 3 complex is oxidized to an Fe species with an OH radical ligand. Together our data indicate that an oxidized Fe species may be the active oxidant, while halides modulate the etch rate by tuning the reduction potential of gold nanocrystals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Uranium(II) Arene Complex That Acts as a Uranium(I) Synthon

We report two-electron reduction of the amidate-supported U(III) mono(arene) complex U(TDA) 3 (2) with KC 8 yields the anionic bis(arene) complex [K[2.2.2]cryptand][U(TDA) 2 ] (3) (TDA = N-(2,6-di-isopropylphenyl)pivalamido). EPR spectroscopy, magnetic susceptibility measurements, and calculations using DFT as well as multireference CASSCF methods all provide strong evidence that the electronic structure of 3 is best represented as a 5f 4 U(II) metal center bound to a monoreduced arene ligand. Reactivity studies show 3 reacts as a U(I) synthon by behaving as a two-electron reductant toward I 2 to form the dinuclear U(III)–U(III) triiodide species [K[2.2.2]cryptand][(UI(TDA) 2 ) 2 (μ-I)] (6) and as a three-electron reductant toward cycloheptatriene (CHT) to form the U(IV) complex [K[2.2.2]cryptand][U(η 7 -C 7 H 7 )(TDA) 2 (THF)] (7). The reaction of 3 with cyclooctatetraene (COT) generates a mixture of the U(III) anion [K[2.2.2]cryptand][U(TDA) 4 ] (1-crypt) and U(COT) 2 , while the addition of COT to complex 2 instead yields the dinuclear U(IV)–U(IV) inverse sandwich complex [U(TDA) 3 ] 2 (μ-η 8 :η 3 -C 8 H 8 ) (8). Two-electron reduction of the homoleptic Th(IV) amidate complex Th(TDA) 4 (4) with KC 8 gives the mono(arene) complex [K[2.2.2]cryptand][Th(TDA) 3 (THF)] (5). The C–C bond lengths and torsion angles in the bound arene of 5 suggest a direduced arene bound to a Th(IV) metal center; this conclusion is supported by DFT calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanism of Nitrogen-Carbon Bond Formation from Iron(IV) Disilylhydrazido Intermediates during N 2 Reduction

We recently reported a reaction sequence that activates C–H bonds in simple arenes as well as the N–N triple bond in N 2 , delivering the aryl group to N 2 to form a new N–C bond ( Nature 2020 , 584, 221). This enables the transformation of abundant feedstocks (arenes and N 2 ) into N-containing organic compounds. The key N–C bond forming step occurs upon partial silylation of N 2 . However, the pathway through which reduction, silylation, and migration occurred was unknown. Here, in this work, we describe synthetic, structural, magnetic, spectroscopic, kinetic, and computational studies that elucidate the steps of this transformation. N 2 must be silylated twice at the distal N atom before aryl migration can occur, and sequential silyl radical and silyl cation addition is a kinetically competent pathway to a formally iron(IV)–NN(SiMe 3 ) 2 intermediate that can be isolated at low temperature. Kinetic studies show its first-order conversion to the migrated product, and DFT calculations indicate a concerted transition state for migration. The electronic structure of the formally iron(IV) intermediate is examined using DFT and CASSCF calculations, which reveal contributions from iron(II) and iron(III) resonance forms with oxidized NNSi 2 ligands. The depletion of electron density from the Fe-coordinated N atom makes it electrophilic enough to accept the incoming aryl group. This new pathway for the N–C bond formation offers a method for functionalizing N 2 using organometallic chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Conical Intersection Accessibility Dictates Brightness in Red Fluorescent Proteins

Red fluorescent protein (RFP) variants are highly sought after for in-vivo imaging since longer wavelengths improve depth and contrast in fluorescence imaging. However, the lower energy emission wavelength usually correlates with a lower fluorescent quantum yield compared to their green emitting counterparts. To guide the rational design of bright variants, we have theoretically assessed two variants (mScarlet and mRouge) which are reported to have very different brightness. Using an α-CASSCF QM/MM framework (chromophore and all protein residues within 6 Å of it in the QM region, for a total of more than 450 QM atoms), we identify key points on the ground and first excited state potential energy surfaces. The brighter variant mScarlet has a rigid scaffold, and the chromophore stays largely planar on the ground state. The dimmer variant mRouge shows more flexibility and can accommodate a pre-twisted chromophore conformation which provides easier access to conical intersections. Notably, the main difference between the variants lies in the intersection seam regions, which appear largely inaccessible in mScarlet but partially accessible in mRouge. This observation is mainly related with changes in the cavity charge distribution, the hydrogen-bonding network involving the chromophore and a key ARG/THR mutation (which changes both charge and steric hindrance).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electrically tuned hyperfine spectrum in neutral Tb( II )(Cp iPr5 ) 2 single-molecule magnet

Molecular spin qubits with long spin coherence time as well as non-invasive operation methods on such qubits are in high demand. It was shown that both molecular electronic and nuclear spin levels can be used as qubits. In solid state systems with dopants, an electric field was shown to effectively change the spacing between the nuclear spin qubit levels when the electron spin density is high at the nucleus of the dopant. Inspired by such solid-state systems, we propose that divalent lanthanide (Ln) complexes with an unusual electronic configuration of Ln 2+ have a strong interaction between the Ln nuclear spin and the electronic degrees of freedom, which renders electrical tuning of the interaction. Furthermore, as an example, we study electronic structure and hyperfine interaction of the 159 Tb nucleus in a neutral Tb(II)(Cp iPr5 ) 2 single-molecule magnet (SMM), which exhibits unusually long magnetization relaxation time, using the complete active space self-consistent field (CASSCF) method with spin–orbit interaction included within the restricted active space state interaction (RASSI). Our calculations show that the low-energy states arise from 4f 8 (6s,5d z 2 ) 1 , 4f 8 (5d x 2 -y 2 ) 1 , and 4f 8 (5d xy ) 1 configurations. We compute the hyperfine interaction parameters and the electronic–nuclear spectrum within our multiconfigurational approach. We find that the hyperfine interaction is about one order of magnitude greater than that for Tb(III)Pc 2 SMMs. This stems from the strong Fermi contact interaction between the Tb nuclear spin and the electron spin density at the nucleus that originates from the occupation of the (6s,5d) orbitals. We also uncover that the response of the Fermi contact term to electric field results in electrical tuning of the electronic–nuclear level separations. This hyperfine Stark effect may be useful for applications of molecular nuclear spins for quantum computing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Actinide arene-metalates: ion pairing effects on the electronic structure of unsupported uranium–arenide sandwich complexes

Addition of [UI 2 (THF) 3 (μ-OMe)] 2 ·THF (2·THF) to THF solutions containing 6 equiv. of K[C 14 H 10 ] generates the heteroleptic dimeric complexes [K(18-crown-6)(THF) 2 ] 2 [U(η 6 -C 14 H 10 )(η 4 -C 14 H 10 )(μ-OMe)] 2 ·4THF (1 18C6 ·4THF) and {[K(THF) 3 ][U(η 6 -C 14 H 10 )(η 4 -C 14 H 10 )(μ-OMe)]} 2 (1 THF ) upon crystallization of the products in THF in the presence or absence of 18-crown-6, respectively. Both 1 18C6 ·4THF and 1 THF are thermally stable in the solid-state at room temperature; however, after crystallization, they become insoluble in THF or DME solutions and instead gradually decompose upon standing. X-ray diffraction analysis reveals 1 18C6 ·4THF and 1 THF to be structurally similar, possessing uranium centres sandwiched between bent anthracenide ligands of mixed tetrahapto and hexahapto ligation modes. Yet, the two complexes are distinguished by the close contact potassium-arenide ion pairing that is seen in 1 THF but absent in 1 18C6 ·4THF, which is observed to have a significant effect on the electronic characteristics of the two complexes. Structural analysis, SQUID magnetometry data, XANES spectral characterization, and computational analyses are generally consistent with U(IV) formal assignments for the metal centres in both 1 18C6 ·4THF and 1 THF , though noticeable differences are detected between the two species. For instance, the effective magnetic moment of 1 THF (3.74 μ B ) is significantly lower than that of 1 18C6 ·4THF (4.40 μ B ) at 300 K. Furthermore, the XANES data shows the U L III -edge absorption energy for 1 THF to be 0.9 eV higher than that of 1 18C6 ·4THF, suggestive of more oxidized metal centres in the former. Of note, CASSCF calculations on the model complex {[U(η 6 -C 14 H 10 )(η 4 -C 14 H 10 )(μ-OMe)] 2 } 2– (1*) shows highly polarized uranium–arenide interactions defined by π-type bonds where the metal contributions are primarily comprised by the 6d-orbitals (7.3 ± 0.6%) with minor participation from the 5f-orbitals (1.5 ± 0.5%). These unique complexes provide new insights into actinide–arenide bonding interactions and show the sensitivity of the electronic structures of the uranium atoms to coordination sphere effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Using a multistate mapping approach to surface hopping to predict the ultrafast electron diffraction signal of gas-phase cyclobutanone

Using the recently developed multistate mapping approach to surface hopping (multistate MASH) method combined with SA(3)-CASSCF(12,12)/aug-cc-pVDZ electronic structure calculations, the gas-phase isotropic ultrafast electron diffraction (UED) of cyclobutanone is predicted and analyzed. After excitation into the n-3s Rydberg state (S2), cyclobutanone can relax through two S2/S1 conical intersections, one characterized by compression of the CO bond and the other by dissociation of the α–CC bond. Subsequent transfer into the ground state (S0) is then achieved via two additional S1/S0 conical intersections that lead to three reaction pathways: α ring-opening, ethene/ketene production, and CO liberation. The isotropic gas-phase UED signal is predicted from the multistate MASH simulations, allowing for a direct comparison to the experimental data. This work, which is a contribution to the cyclobutanone prediction challenge, facilitates the identification of the main photoproducts in the UED signal and thereby emphasizes the importance of dynamics simulations for the interpretation of ultrafast experiments.

Chemistry↗

Ultrafast electron diffraction of photoexcited gas-phase cyclobutanone predicted by ab initio multiple cloning simulations

We present the result of our calculations of ultrafast electron diffraction (UED) for cyclobutanone excited into the S2 electronic state, which is based on the non-adiabatic dynamics simulations with the Ab Initio Multiple Cloning (AIMC) method with the electronic structure calculated at the SA(3)-CASSCF(12,12)/aug-cc-pVDZ level of theory. The key features in the UED pattern were identified, which can be used to distinguish between the reaction pathways observed in the AIMC dynamics, although there is a significant overlap between representative signals due to the structural similarity of the products. The calculated UED pattern can be compared with the experiment.

Chemistry↗

Excited electronic states of Na 2 and K 2 : The potential for long-lived “reservoir” states leading to collision induced population inversions

Potential energy curves (PECs) for the spin-free (ΛS) and spin–orbit (Ω) states associated with the four lowest-lying dissociation channels of Na 2 and K 2 were calculated at the SA-CASSCF/SO-CASPT2/aug-cc-pwCVQZ-DK level. The PECs of Na 2 were consistent with the experimental data and with the FS-CCSD (2,0) calculations, reproducing the double-well and the “shelf” character for some of the potentials of the excited states. For K 2 , the PECs behaved in a similar way and the spectroscopic parameters for the ground and the excited states are in good agreement with the available experimental values. The dissociation energy of K 2 was predicted to be D e = 4454 cm −1 , within an agreement of 5 cm −1 with the experiments. For Na 2 , D e = 5789 cm −1 compared to the experimental value of 6022 cm −1 . The inclusion of spin–orbit coupling effects resulted in avoided crossings, which affect the PECs. Spin–orbit changes the predicted curves for some excited Ω states arising from ΛS states that overlap each other, affecting their associated vibrational frequencies and bond distances. Here, the current studies of the low-lying states in K 2 reveal a similar structure to those of Na 2 , which suggests the accessibility of long-lived energy storing reservoir states and possible population inversions in K 2 following prior experimental work on the reaction of halogen atoms with Na 3 to produce excited states of Na 2 .

Ab-initio methods↗

Fast methods for multisite charge transfer. Processes II. Analytic nuclear gradients and nonadiabatic dynamics for cCASSCF(1,n) and cCASSCF(2n-1,n) wavefunctions

In this work we derive and implement analytic nuclear gradients and derivative couplings for a constrained complete active space self-consistent field with a small active space designed to model electron or hole transfer. Using a Lagrangian formalism, we are able to differentiate both the CASSCF energy and the constraint (which is required for smooth surfaces over a wide range of parameter space), and the resulting efficient algorithm can be immediately applied to nonadiabatic dynamics simulations of charge transfer processes. Here, we run initial surface-hopping simulations of a proton coupled electron transfer event for a phenoxyl–phenol system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Many roads to the seam: How conformational flexibility drives nonadiabatic relaxation in a prototypical tetrapyrrolic chromophore

Large and structurally flexible chromophores pose challenges for in silico modeling of photodeactivation due to the many vibrational modes that can funnel the system toward energy degeneracy. In this work, we examine how the multiple degrees of freedom in biliverdin, a prototypical tetrapyrrolic chromophore, cooperate to drive access to the S 1 /S 0 intersection seam in vacuo. We begin by mapping the ground-state potential energy surface to identify representative biliverdin conformers relevant to photoexcitation. We then use a CASSCF-based framework to map the excited-state landscape and characterize the intersection seam, identifying distinct conical-intersection types. Finally, we employ ab initio multiple spawning to resolve the dynamical pathways by which the system accesses these regions. DFT potential-energy and free-energy mappings indicate that, although several conformers are relevant, the “locked-helix” ZsZsZs conformer predominates in the ground state. The intersection seam comprises numerous geometrically distinct regions characterized by varying degrees and combinations of dihedral torsion, pyramidalization, and bond-length alternation. Yet only select regions lie within energetic reach, and moderate barriers separate them from the S 1 minimum. Nonadiabatic dynamics combined with multivariate analyses show that, despite extensive mode coupling during deactivation that guides the system toward multiple regions of the seam, a single dihedral torsion, together with bond-length alternation, predominantly drives energy degeneracy. This work offers new insight into biliverdin’s intrinsic photochemical response and underscores a general feature of flexible chromophores: many modes may participate during photorelaxation, but only a limited subset ultimately dictates seam accessibility.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Activation of methane by U + studied by guided ion beam tandem mass spectrometry and quantum chemistry

Reaction pathways of all products formed in the U + + CH 4 (CD 4 ) reaction were explored as a function of kinetic energy using guided ion beam tandem mass spectrometry and quantum chemical calculations. UH + , UC + , UCH + , UCH 2 + , and UCH 3 + (and their perdeuterated analogues) are formed in endothermic reactions. In both systems, the UCH 2 + (UCD 2 + ) dehydrogenated product was the dominant product in the low-energy region, whereas the UH + (UD + ) hydride product became predominant at high energies. The kinetic energy behavior of the various products is consistent with a common intermediate of H–U + –CH 3 (D–U + –CD 3 ). Here, the kinetic energy dependence of all product cross sections was modeled to obtain experimental bond dissociation energies at 0 K (in eV): D 0 (U + –H) = 2.42 ± 0.10, D 0 (U + –C) = 3.95 ± 0.12, D 0 (U + –CH) = 4.91 ± 0.09, D 0 (U + –CH 2 ) = 4.11 ± 0.04, and D 0 (U + –CH 3 ) = 2.41 ± 0.09. Quantum chemical calculations using the UCCSD(T) and UB3LYP approaches with the cc-pwCVXZ-PP basis set with MDF-60 pseudopotential for U + and the aug-cc-pCVXZ and aug-cc-pVXZ (X = T, Q) basis set for carbon and hydrogen, respectively, validate the experimental bond dissociation energies and outline the potential energy surface for all reactions observed. In addition, spin–orbit corrections of the bond energies for all products were calculated at a CASSCF-CASPT2-RASSI level.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗