Ab Initio Configuration Interaction Calculations Research Articles

Exploring the photochemistry of OAlOH: Photodissociation pathways and electronic spectra.

This study was focused on the photochemistry of OAlOH and three possible pathways, which were studied with high-level multireference configuration interaction abinitio calculations. We computed cuts of the six-dimensional potential energy surfaces for the ground, the lowest singlet and triplet excited states, and probed the photodissociation mechanisms and the stabilities. The OAlOH electronic spectrum, with an energy reaching 7.15eV, contained four prominent peaks. Photodissociation to AlO, OH, and AlOH constituted a plausible mechanism within the deep-UV range (λ = 250.4nm). Our data indicated the photostability of OAlOH in the near-UV‒Vis region, so detection with laser-induced fluorescence is possible. Fluorescence and phosphorescence may occur upon excitation at 363.5nm. The roles of OAlOH in the photochemical reactions of Al-bearing molecules in the upper atmosphere and VY Canis Majoris are discussed.

Quantum study of the CH3+ photodissociation in full-dimensional neural network potential energy surfaces.

C H 3 + , a cornerstone intermediate in interstellar chemistry, has recently been detected for the first time by using the James Webb Space Telescope. The photodissociation of this ion is studied here. Accurate explicitly correlated multi-reference configuration interaction abinitio calculations are done, and full-dimensional potential energy surfaces are developed for the three lower electronic states, with a fundamental invariant neural network method. The photodissociation cross section is calculated using a full-dimensional quantum wave packet method in heliocentric Radau coordinates. The wave packet is represented in angular and radial grids, allowing us to reduce the number of points physically accessible, requiring to push up the spurious states appearing when evaluating the angular kinetic terms, through projection technique. The photodissociation spectra, when employed in astrochemical models to simulate the conditions of the Orion bar, result in a lesser destruction of CH3+ compared to that obtained when utilizing the recommended values in the kinetic database for astrochemistry.

A High-level Ab Initio Study of the Destruction of Methanimine under UV Radiation

The photodecomposition of methanimine (CH2NH) in the interstellar medium through several possible pathways is investigated by means of high-level multireference configuration interaction ab initio calculations. Among these pathways are photodissociation pathways involving hydrogen-atom elimination from both the CH2 and NH groups, and fragmentation into CH2 and NH. Potential-energy curves for the ground and several excited electronic states, as well as nonadiabatic couplings between them, are calculated. Possible dissociation mechanisms are discussed for the different pathways. It is found that the minimum excitation energy required for methanimine dissociation is above 7 eV. By using a two-dimensional representation of methanimine, CH2NH → CHNH2 isomerization is explored as an additional methanimine decomposition pathway. Hydrogen-atom elimination from the CH2 group is also investigated along the isomerization pathway. The results show that the isomerization proceeds by overcoming a transition state that in the first two excited states would require excitation energies similar to or somewhat lower than the typical minimum energies needed for breaking the molecule through the fragmentation pathways. Therefore, CH2NH → CHNH2 isomerization can effectively contribute to methanimine decomposition, competing efficiently with the photodissociation pathways. The radiation content present in the interstellar medium makes possible the occurrence of all the pathways studied.

Probing the quadrupole transition strength of C15 via deuteron inelastic scattering

Deuteron elastic scattering from $^{15}\mathrm{C}$ and inelastic scattering reactions to the first excited state of $^{15}\mathrm{C}$ were studied using a radioactive beam of $^{15}\mathrm{C}$ in inverse kinematics. The scattered deuterons were measured using HELIOS. The elastic scattering differential cross sections were analyzed using the optical model. A matter deformation length ${\ensuremath{\delta}}_{d}=1.04(11)\phantom{\rule{0.16em}{0ex}}\mathrm{fm}$ has been extracted from the differential cross sections of inelastic scattering to the first excited state. The ratio of neutron and proton matrix elements ${M}_{n}/{M}_{p}=3.6(4)$ has been determined from this quadrupole transition. Neutron effective charges and core-polarization parameters of $^{15}\mathrm{C}$ were determined and discussed. Results from ab initio no-core configuration interaction calculations were also compared with the experimental observations. This result supports a moderate core decoupling effect of the valence neutron in $^{15}\mathrm{C}$ similarly to its isotone $^{17}\mathrm{O}$, in line with the interpretation of other neutron-rich carbon isotopes.

Precision half-life determination for the β+ emitter N13

A new precision half-life measurement of $^{13}\mathrm{N}$ has been conducted using the TwinSol $\ensuremath{\beta}$-counting station at the University of Notre Dame. The measured value of ${t}_{1/2}^{\mathrm{new}}=597.05(19)\phantom{\rule{4pt}{0ex}}\mathrm{s}$ differs from the previous world value by about $2.8\ensuremath{\sigma}$. An evaluation of the $^{13}\mathrm{N}$ half-life results in a ${t}_{1/2}^{\mathrm{world}}=597.19(22)\phantom{\rule{4pt}{0ex}}\mathrm{s}$. Updated standard model predictions for the Fermi to Gamow-Teller mixing ratio $\ensuremath{\rho}$ and its associated correlation parameters have been calculated using the new $^{13}\mathrm{N}$ world half-life in preparation for a future measurement of the mixing ratio. Finally, an ab initio no-core configuration interaction (NCCI) calculation for the $B(\mathrm{GT})$ of this decay, carried out using the Daejeon16 interaction, has been performed, revealing the need for higher-order chiral corrections.

Ab initio estimation of E2 strengths in Li8 and its neighbors by normalization to the measured quadrupole moment

For electric quadrupole ($E2$) observables, which depend on the large-distance tails of the nuclear wave function, ab initio no-core configuration interaction (NCCI) calculations converge slowly, making meaningful predictions challenging to obtain. Nonetheless, the calculated values for different $E2$ matrix elements, particularly those involving levels with closely-related structure (e.g., within the same rotational band) are found to be robustly proportional. This observation suggests that a known value for one observable may be used to determine the overall scale of $E2$ strengths, and thereby provide predictions for others. In particular, we demonstrate that meaningful predictions for $E2$ transitions may be obtained by calibration to the ground-state quadrupole moment. We test this approach for well-measured low-lying $E2$ transitions in $^7$Li and $^9$Be, then provide predictions for transitions in $^8$Li and $^9$Li. In particular, we address the $2^+\rightarrow1^+$ transition in $^8$Li, for which the reported measured strength exceeds ab initio Green's function Monte Carlo (GFMC) predictions by over an order of magnitude.

Natural orbitals for theab initiono-core configuration interaction approach

Ab initio no-core configuration interaction (NCCI) calculations for the nuclear many-body problem have traditionally relied upon an antisymmetrized product (Slater determinant) basis built from harmonic oscillator orbitals. The accuracy of such calculations is limited by the finite dimensions which are computationally feasible for the truncated many-body space. We therefore seek to improve the accuracy obtained for a given basis size by optimizing the choice of single-particle orbitals. Natural orbitals, which diagonalize the one-body density matrix, provide a basis which maximizes the occupation of low-lying orbitals, thus accelerating convergence in a configuration-interaction basis, while also possibly providing physical insight into the single-particle structure of the many-body wave function. We describe the implementation of natural orbitals in the NCCI framework and examine the nature of the natural orbitals thus obtained, the properties of the resulting many-body wave functions, and the convergence of observables. After taking ${}^{3}\mathrm{He}$ as an illustrative testbed, we explore aspects of NCCI calculations with natural orbitals for the ground state of the $p$-shell neutron halo nucleus ${}^{6}\mathrm{He}$.

Thermodynamics stability, electronic structures and spectroscopic properties of defects and Ce3+ ions in Y2O3

Thermodynamic properties, electronic structures and spectroscopic properties of defects and Ce3+ in Y2O3 are studied by using the hybrid density functional theory associated with multi-reference configuration interaction ab-initio calculations. Thermodynamic transition energy levels of the easily generated oxygen vacancies in the host are analyzed according to HSE06-calculated formation energies, which may be conducive to interpretations of the persistent luminescence (PersL) of Y2O3-based phosphors. Besides, the locations of impurity states (caused by VO and Ce3+) in energy bands are obtained from derived density of states. Moreover, energies and oscillator strengths of 4f1 → 5d1−5 transitions of Ce3+ ions (at Y1 and Y2 sites) calculated from the CASSCF/CASPT2/RASSI−SO method agree reasonably well with experimental excitation spectra of Y2O3: Ce3+ phosphors, achieving the assignment of excitation spectra. The presented calculations can be applied to identify luminescent centers in Ce3+-doped phosphors and reveals possible native defects and their roles in the PersL of phosphors.

Singlet Fission Rate: Optimized Packing of a Molecular Pair. Ethylene as a Model.

A procedure is described for unbiased identification of all π-electron chromophore pair geometry choices that locally maximize the rate of conversion of a singlet exciton into a singlet biexciton (triplet pair), using a simplified version of the diabatic frontier orbital model of singlet fission (SF). The resulting approximate optimal geometries provide insight and are expected to represent useful starting points for searches by more advanced methods. The general procedure is illustrated on a pair of ethylenes as the simplest model of a π-electron system, but it is applicable to pairs of much larger molecules, with dozens of non-hydrogen atoms, and not necessarily planar. We first examine the value of |TA|2, the square of the electronic matrix element for SF with initial excitation fully localized on partner A, on a grid of several billion geometries within the six-dimensional space of physically realizable possibilities. Several of the optimized pair geometries are somewhat unexpected, but all are found to follow the qualitative guidance proposed earlier. In the neighborhood of each local maximum of |TA|2, consideration of mixing with charge-transfer configurations and of excitonic interaction between partners A and B determines the SF energy balance and yields squared matrix elements |T*|2 and |T**|2 for the lower and upper excitonic states S* and S**, respectively. Assuming Boltzmann populations of these states, the geometry is further optimized to maximize k, the sum of the SF rates obtained from Marcus theory, and this reorders the suitable geometries substantially. At 87 pair geometries, the |T*|2 and |T**|2 values are compared with those obtained from high-level ab initio nonorthogonal configuration interaction calculations and found to follow the same trend. Finally, the biexciton binding energy at the optimized geometries is calculated. Altogether, 13 significant local maxima of SF rate for a pair of ethylenes are identified in the physically relevant part of space that avoids molecular interpenetration in the hard-sphere approximation. The three best geometries are twist-stacked, slip-stacked, and L-shaped. The maxima occur at the (five-dimensional) surfaces of seven six-dimensional "parent" regions of space centered at physically inaccessible geometries at which the calculated SF rate is very large but the two ethylenes interpenetrate. The results are displayed in interactive graphics. The computer code ("Simple") written for these calculations is flexible in that it permits a choice of performing the search for local maxima in six dimensions on |TA|2, |T*|2, or k. It is available as freeware at https://cloud.uochb.cas.cz/simple .

SCCS−radical: Renner-Teller effect and spin-orbit coupling in the X 2Πu electronic state

SCCS- was detected by laser-induced fluorescence spectroscopy in 2003 (M. Nakajima, Y. Yoneda, Y. Sumiyoshi, T. Nagata, Y. Endo, J. Chem. Phys. 119 (2003) 7805) and its spectrum was analyzed and the results presented together with ab initio calculations data. Symmetrical stretching vibrations were assigned in both the ground X 2?u and the first excited A 2?g electronic states, but no data about spin?orbit splitting and bending vibrational modes were given. In the present work, the vibronic levels in the ground electronic state of SCCS- were calculated by means of a model developed by Peric and coworkers for the treatment of the Renner?Teller effect in any-atomic linear species in its variational form (M. Mitic, R. Rankovic, M. Milovanovic, S. Jerosimic, M. Peric, Chem. Phys. 464 (2016) 55) using the ab initio multireference configuration interaction calculations for obtaining potential energy curves in the Born?Oppenheimer approximation. Additionally, the spin?orbit splitting in the ground state was investigated taking into account the interaction with the first excited state, and the energies obtained from combined treatment of vibronic and spin?orbit interactions in the ground state are reported. Finally, based on the present results, several assignments of unidentified bands reported by Nakajima et al. are proposed.

Full-Dimensional Ab Initio Potential Energy Surface and Vibrational Energy Levels of Li₂H.

We built a full-dimensional analytical potential energy surface of the ground electronic state of Li2H from ca. 20,000 ab initio multi-reference configuration interaction calculations, including core–valence correlation effects. The surface is flexible enough to accurately describe the three dissociation channels: Li (2s 2S) + LiH (1Σ+), Li2 (1Σg+) + H (1s 2S) and 2Li (2s 2S) + H (1s 2S). Using a local fit of this surface, we calculated pure (J = 0) vibrational states of Li2H up to the barrier to linearity (ca. 3400 cm−1 above the global minimum) using a vibrational self-consistent field/virtual state configuration interaction method. We found 18 vibrational states below this barrier, with a maximum of 6 quanta in the bending mode, which indicates that Li2H could be spectroscopically observable. Moreover, we show that some of these vibrational states are highly correlated already ca. 1000 cm−1 below the height of the barrier. We hope these calculations can help the assignment of experimental spectra. In addition, the first low-lying excited states of each B1, B2 and A2 symmetry of Li2H were characterized.

Theoretical study of electronic properties and isotope effects in the UV absorption spectrum of disulfur

The electronic structures of triplet S2 ground and excited states are studied by ab initio molecular orbital and configuration interaction calculation. Potential energy curves correlated with S(3P)+S(3P) and S(3P)+S(1D) at the dissociation limit are evaluated, and electronic terms for a total of 11 states are assigned. Transition dipole moments, as a function of internuclear distance, are determined for two allowed transitions to B″3Πu and B3Σu- excited states. The total absorption cross-sections are computed to estimate isotope-fractionation constants, ε, for four most common isotopologues: 32S32S, 32S33S, 32S34S, and 32S36S by quantum close-coupling (R-matrix) expansion approach and they are found to lie in a mostly opaque to competing absorbers spectral window. We suggest that the photochemistry and isotopic effects of S2 are of significant importance and provide data showing high sensitivity of mass-independent fractionation to excitation wavelength. Zero-point energy based constants εZPE are estimated as well to compare with the obtained isotope effects and two modes for MIF are present in three-isotope plots; large isotopic effects were observed for both 36S and 33S with an excitation wavelength-dependent fluctuation.

An investigation of the sites occupied by atomic barium in solid xenon-A 2D-EE luminescence spectroscopy and molecular dynamics study.

A detailed characterisation of the luminescence recorded for the 6p 1P1-6s 1S0 transition of atomic barium isolated in annealed solid xenon has been undertaken using two-dimensional excitation-emission (2D-EE) spectroscopy. In the excitation spectra extracted from the 2D-EE scans, two dominant thermally stable sites were identified, consisting of a classic, three-fold split Jahn-Teller band, labeled the blue site, and an unusual asymmetric 2 + 1 split band, the violet site. A much weaker band has also been identified, whose emission is strongly overlapped by the violet site. The temperature dependence of the luminescence for these sites was monitored revealing that the blue site has a non-radiative channel competing effectively with the fluorescence even at 9.8 K. By contrast, the fluorescence decay time of the violet site was recorded to be 4.3 ns and independent of temperature up to 24 K. The nature of the dominant thermally stable trapping sites was investigated theoretically with Diatomics-in-Molecule (DIM) molecular dynamics simulations. The DIM model was parameterized with ab initio multi-reference configuration interaction calculations for the lowest energy excited states of the Ba⋅Xe pair. The simulated absorption spectra are compared with the experimental results obtained from site-resolved excitation spectroscopy. The simulations allow us to assign the experimental blue feature spectrum to a tetra-vacancy trapping site in the bulk xenon fcc crystal-a site often observed when trapping other metal atoms in rare gas matrices. By contrast, the violet site is assigned to a specific 5-atom vacancy trapping site located at a grain boundary.

Natural orbital description of the halo nucleus 6He

Ab initio calculations of nuclei face the challenge of simultaneously describing strong short-range internucleon correlations and the long-range properties of weakly-bound halo nucleons. Natural orbitals, which diagonalize the one-body density matrix, provide a basis which is better matched to the physical structure of the many-body wave function. We demonstrate that the use of natural orbitals significantly improves convergence for ab initio no-core configuration interaction calculations of the neutron halo nucleus 6He, relative to the traditional oscillator basis.

RETRACTED: Theoretical study of electronic properties and isotope effects in the UV absorption spectrum of disulfur

Abstract The electronic structures of triplet S2 ground and excited states are studied by ab initio molecular orbital and configuration interaction calculation. Potential energy curves correlated with S ( 3 P ) + S ( 3 P ) and S ( 3 P ) + S ( 1 D ) at the dissociation limit are evaluated, and electronic terms for a total of 11 states are assigned. Transition dipole moments, as a function of internuclear distance, are determined for two allowed transitions to B ″ 3 Π u and B 3 Σ u - excited states. The total absorption cross-sections are computed to estimate isotope-fractionation constants, e , for four most common isotopologues: 32 S 32 S , 32 S 33 S , 32 S 34 S , and 32 S 36 S by quantum close-coupling (R-matrix) expansion approach and they are found to lie in a mostly opaque to competing absorbers spectral window. We suggest that the photochemistry and isotopic effects of S2 are of significant importance and provide data showing high sensitivity of mass-independent fractionation to excitation wavelength. Zero-point energy based constants e ZPE are estimated as well to compare with the obtained isotope effects and two modes for MIF are present in three-isotope plots; large isotopic effects were observed for both 36 S and 33 S with an excitation wavelength-dependent fluctuation.

Femtosecond real-time probing of reactions MMXVII: The predissociation of sodium iodide in the A 0+ state

We revisit the femtosecond transition-state spectroscopy of sodium iodide taking advantage of modern lasers and pulse-shaping to better map the low-lying electronic states, some forming predissociative wells through curve crossings. We also carry out high-level ab initio multi-reference configuration interaction calculations including spin-orbit coupling terms and using large correlation-consistent basis sets to arrive at accurate ground- and excited-state potential energy curves of NaI. Density matrix calculations employing vibrational wave functions determined from the ab initio X 0+ and A 0+ potentials are used to simulate time dependent wave packet dynamics of NaI pumped to the A 0+ state.

Magnetic Excitations and Electronic Interactions in Sr_{2}CuTeO_{6}: A Spin-1/2 Square Lattice Heisenberg Antiferromagnet.

Sr_{2}CuTeO_{6} presents an opportunity for exploring low-dimensional magnetism on a square lattice of S=1/2 Cu^{2+} ions. We employ abinitio multireference configuration interaction calculations to unravel the Cu^{2+} electronic structure and to evaluate exchange interactions in Sr_{2}CuTeO_{6}. The latter results are validated by inelastic neutron scattering using linear spin-wave theory and series-expansion corrections for quantum effects to extract true coupling parameters. Using this methodology, which is quite general, we demonstrate that Sr_{2}CuTeO_{6} is an almost ideal realization of a nearest-neighbor Heisenberg antiferromagnet but with relatively weak coupling of 7.18(5)meV.

Study of the electronic and rovibronic structure of the X ²Σ⁺, A ²Π, and B ²Σ⁺ states of AlO.

The electronic structure of the X (2)Σ(+), A (2)Π, and B (2)Σ(+) states of aluminum monoxide (AlO) are studied via ab initio multi-reference configuration interaction calculations. Core correlation corrections, several basis sets, and active space choices are considered. Angular momentum and spin-orbit coupling terms are obtained at different levels of theory. The resulting ab initio curves are used to solve the associated rovibronic problem for the total angular momentum J up to 112.5 and then also refined by fitting to the experimental wavenumbers available in the literature, reproducing them with the root-mean-square error of 0.07 cm(-1). Theoretical rovibronic energy levels of AlO in its X (2)Σ(+), A (2)Π, and B (2)Σ(+) electronic states are presented including those from the X - B blue-green system.

Halo nucleiHe6andHe8with the Coulomb-Sturmian basis

The rapid Gaussian falloff of the oscillator functions at large radius makes them poorly suited for the description of the asymptotic properties of the nuclear wave function, a problem which becomes particularly acute for halo nuclei. We consider an alternative basis for ab initio no-core configuration interaction (NCCI) calculations, built from Coulomb-Sturmian radial functions, allowing for realistic exponential falloff. NCCI calculations are carried out for the neutron halo nuclei 6,8He, as well as the baseline case 4He, with the JISP16 nucleon-nucleon interaction. Estimates are made for the root-mean-square radii of the proton and matter distributions.

Potential energy surfaces and quasiclassical trajectory study of the O + H2(+)→ OH(+) + H, OH + H(+) proton and hydrogen atom transfer reactions and isotopic variants (D2(+), HD(+)).

The rate constants (k; T: 200-900 K) and cross-sections (σ; Ecol: 0.010-0.50 eV) of the O + H2(+)→ OH(+) + H (1), OH + H(+) (2) reactions, which occur on the ground (1(2)A'') and first excited (1(2)A') potential energy surfaces (PESs), respectively, were investigated for the first time, considering also the rate constants for D2(+) and HD(+). Ab initio multireference configuration interaction calculations were performed on both barrierless PESs (where the minimum energy path involves the insertion of the O atom into the middle of the H2(+) bond), and suitable analytical expressions were developed for the first time and used in quasiclassical trajectory (QCT) calculations. k(1) ≈ 3k(2) independent of the isotopic variant, k(H2(+)) > k(HD(+)) > k(D2(+)) for , and the intermolecular and intramolecular isotopic effects are essentially independent of T. Comparison with the Langevin-Gioumousis-Stevenson (LGS) simple capture model shows that these results are similar to the QCT ones, especially for ; and the isotopic effects are coincident with the QCT ones for both reactions. For O + H2(+), σ(1) ≈ 3σ(2) at Ecol≤ 0.10 eV, and σ(1) = 1.5σ(2) at 0.40 and 0.50 eV. The larger value of σ(1(2)A'') with respect to σ(1(2)A') arises from the larger value of bmax(1(2)A'') with respect to bmax(1(2)A'), and this results from the more attractive character of the former PES. Besides, the reaction probabilities are quite large [0.78-0.98 (1(2)A'') and 0.78-0.93 (1(2)A')], and the decreasing trend of both cross-sections as Ecol increases arises from the barrierless character of both PESs. We expect that these results (in particular, the competition between proton transfer and hydrogen atom transfer) will encourage experimentalists to carry out investigations on this interesting reaction.