Search NASA⌕ Search

SEARCH · Search NASA

Results for “Electronic wave function”

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 1,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Controlling Molecular Structure and Spin with Multiconfigurational Quantum Chemistry (Final Technical Report)

For many first-row transition metal complexes, structure-property relationships can be obtained from high-level molecular geometry optimizations and subsequent electronic structure studies. The project funded under this award utilized newly implemented fully internally contracted (FIC) nuclear gradients for extended multi-state (XMS) complete active space second-order multireference perturbation theory (CASPT2) to explore geometry changes in first-row transition metal coordination complexes. At the start of the project, only one full geometry optimization using FIC-CASPT2 analytical gradients had been reported (J. Chem. Theory Comput., 2016, 12 (8), 3781), and many open-questions regarding the performance and achievable accuracy in applying such computations to larger complexes persisted. Key deliverables in the report include the implementation of the numerical Hessian and subsequent vibrational analysis in the BAGEL program package. This includes both full Hessian vibrational analysis (FHVA) and a partial Hessian vibrational analysis (PHVA). The latter of which allows us to reduce the significant cost of computing the vibrational frequencies in a large molecule by focusing on the modes of highest interest. The gradient code also proved useful in the development of an approach to improve the Hubbard-U correction by removing bias in predicting spin-splitting in Fe(II) complexes via plane-wave density functional theory (DFT). Finally, we showed for three families of complexes (spin-crossover (SCO) complexes, metallocorroles complexes, and those with metal-metal bonds) that CASPT2 can result in good molecular geometries and established best practices in undertaking this work. We are now using this approach in collaborative projects with experimental groups, which would not have been possible at the start of this project. A common theme also arose through this work that static and dynamic correlation must be recovered in a balanced way to yield quantitative results. By systematically varying how both effects were included, we were able to make key chemical insights. This work supported the training of two postdoctoral scholars, one graduate student, and one undergraduate student in applying high-level wavefunction based methods to challenging systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1999

Modeling of THz Lasers Based on Intersubband Transitions in Semiconductor Quantum Wells

In semiconductor quantum well structures, the intersubband energy separation can be adjusted to the terahertz (THz) frequency range by changing the well width and material combinations. The electronic and optical properties of these nanostructures can also be controlled by an applied dc electric field. These unique features lead to a large frequency tunability of the quantum well devices. In the on-going project of modeling of the THz lasers, we investigate the possibility of using optical pumping to generate THz radiation based on intersubband transitions in semiconductor quantum wells. We choose the optical pumping because in the electric current injection it is difficult to realize population inversion in the THz frequency range due to the small intersubband separation (4-40 meV). We considered both small conduction band offset (GaAs/AlGaAs) and large band offset (InGaAs/AlAsSb) quantum well structures. For GaAs/AlGaAs quantum wells, mid-infrared C02 lasers are used as pumping sources. For InGaAs/AlAsSb quantum wells, the resonant intersubband transitions can be excited by the near-infrared diode lasers. For three- and four-subband quantum wells, we solve the pumpfield-induced nonequilibrium distribution function for each subband of the quantum well system from a set of rate equations that include both intrasubband and intersubband relaxation processes. Taking into account the coherent interactions between pump and THz (signal) waves, we calculate the optical gain for the THz field. The gain arising from population inversion and stimulated Raman processes is calculated in a unified manner. A graph shows the calculated THz gain spectra for three-subband GaAs/AlGaAs quantum wells. We see that the coherent pump and signal wave interactions contribute significantly to the gain. The pump intensity dependence of the THz gain is also studied. The calculated results are shown. Because of the optical Stark effect and pump-induced population redistribution, the maximum THz gain saturates at larger pump intensities.

Liu, Ansheng↗
DOE OSTI2021

Computing Edge States without Hard Truncation

We present a numerical method which accurately computes the discrete spectrum and associated bound states of semi-infinite Hamiltonians which model electronic “edge” states localized at boundaries of one- and two-dimensional crystalline materials. The problem is nontrivial since arbitrarily large finite “hard” (Dirichlet) truncations of the Hamiltonian in the infinite bulk direction tend to produce spurious bound states partially supported at the truncation. Our method, which overcomes this difficulty, is to compute the Green's function of the semi-infinite Hamiltonian by imposing an appropriate boundary condition in the bulk direction; then, the spectral data is recovered via Riesz projection. We demonstrate our method's effectiveness by studies of edge states at a graphene zig-zag edge in the presence of defects modeled both by a discrete tight-binding model and a continuum PDE model under finite difference discretization. Our method may also be used to study states localized at domain wall-type edges in one- and two-dimensional materials where the edge Hamiltonian is infinite in both directions; we demonstrate this for the case of a tight-binding model of distinct honeycomb structures joined along a zig-zag edge. Here, we expect our method to be useful for designing novel devices based on precise wave-guiding by edge states.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Accuracy of XAS theory for unraveling structural changes of adsorbates: CO on Ni(100)

Studying surface reactions using ultrafast optical pump and X-ray probe experiments relies on accurate calculations of X-ray spectra of adsorbates for correct identification of the spectral signatures and their dynamical evolution. We show that experimental XAS can be well reproduced for the different binding sites in a static prototype system CO/Ni(100) at a standard DFT (GGA) level of theory, using a plane-wave basis and pseudopotentials. This validates its utility in analyzing ultrafast X-ray probe experiments. The accuracy of computed relative core level binding energies is about 0.2 eV, representing a lower limit for which spectral features can be resolved with this method. We also show that the commonly used Z +1 approximation gives very good core binding energy shifts overall. However, we find a discrepancy for CO adsorbed in the hollow site, which we assign to the significantly stronger hybridization in hollow bonding than in on-top.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Constraints for stellar electron-capture rates on 86 Kr via the 86 Kr(t, 3 He +γ) 86 Br reaction and the implications for core-collapse supernovae

Background: In the late stages of stellar core collapse just prior to core bounce, electron captures on medium-heavy nuclei drive deleptonization. Therefore, simulations require the use of accurate reaction rates. Nuclei with neutron number near N=50 above atomic number Z=28 play an important role. Rates presently used in astrophysical simulations rely primarily on a relatively simple single-state approximation. In order to improve the accuracy of the astrophysical simulations, experimental data are needed to test the electron-capture rates and to guide the development of better theoretical models and astrophysical simulations. Purpose: The purpose of the present work was to measure the Gamow-Teller transition strength from 86 Kr to 86 Br, to derive the stellar electron-capture rates based on the extracted strengths, and to compare the derived rates with rates based on shell-model and quasiparticle random-phase approximation (QRPA) Gamow-Teller strengths calculations, as well as the single-state approximation. An additional purpose was to test the impact of using improved electron-capture rates on the late evolution of core-collapse supernovae. Method: The Gamow-Teller strengths from 86 Kr were extracted from the 86 Kr(t, 3 He +γ) charge-exchange reaction at 115MeV/u. The electron-capture rates were calculated as a function of stellar density and temperature. Besides the case of 86 Kr, the electron-capture rates based on the QRPA calculations were calculated for 78 additional isotopes near N=50 above Z=28. The impact of using these rates instead of those based on the single-state approximation is studied in a spherically symmetrical simulation of core collapse just prior to bounce. Results: The derived electron-capture rates on 86 Kr from the experimental Gamow-Teller strength distribution are much smaller than the rates estimated based on the single-state approximation. Rates based on Gamow-Teller strengths estimated in shell-model and QRPA calculations are more accurate. The core-collapse supernova simulation with electron-capture rates based on the QRPA calculations indicate a significant reduction in the deleptonization during the collapse phase. Conclusions: It is important to utilize microscopic theoretical models that are tested by experimental data to constrain and estimate Gamow-Teller strengths and derived electron-capture rates for nuclei near N=50 that are inputs for astrophysical simulations of core-collapse supernovae and their multimessenger signals, such as the emission of neutrinos and gravitational waves.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2011

Solar Wind Strahl Observations and Their Implication to the Core-Halo Formation due to Scattering

A study of the kinetic properties of the strahl electron velocity distribution functions (VDF?s) in the solar wind is presented. This study focuses on the mechanisms that control and regulate the electron VDF?s and the stability of the strahl electrons in the solar wind; mechanisms that are not yet well understood. Various parameters are investigated such as the strahl-electron density, temperature anisotropy, and electron heat-flux. These parameters are used to investigate the stability of the strahl population. The analysis check for whether the strahl electrons are constrained by some instability (e.g., the whistler or KAW instabilities), or are maintained by other types of processes. The electron heat-flux and temperature anisotropy are determined by modeling of the 3D-VDF?s from which the moments properties of the various populations are obtained. The results of this study have profound implication on the current hypothesis about the probable formation of the solar wind halo electrons produced from the scattering of the strahl population. This hypothesis is strengthened by direct observations of the strahl electrons being scattered into the core-halo in an isolated event. The observation implies that the scattering of the strahl is not a continuous process but occurs in bursts in regions where conditions for wave growth providing the scattering are optimum. Sometimes, observations indicate that the strahl component is anisotropic (Tper/Tpal approx. 2). This provides a possible free energy source for the excitation of whistler waves as a possible scattering mechanism, however this condition is not always observed. The study is based on high time resolution data from the Cluster/PEACE electron spectrometer.

Vinas, Adolfo F.↗
DOE OSTI2022

Sign-free determinant quantum Monte Carlo study of excitonic density orders in a two-orbital Hubbard-Kanamori model

While excitonic instabilities in multiorbital systems recently have come under scrutiny in a variety of transition-metal compounds, understanding emergence of these instabilities from strong electronic interactions has remained a challenge. In this work, we present a sign-problem-free determinant quantum Monte Carlo study of excitonic density orders in a half-filled two-orbital Hubbard-Kanamori model with broken orbital degeneracy, which accounts for the role of Hund's coupling in transition-metal compounds. For strong inverted (negative) Hund's exchange, we find numerical evidence for the emergence of excitonic density order, with competition between anti-ferro-orbital order and Q = (π, π) excitonic density order as a function of orbital splitting and Hund's coupling. While inverted Hund's coupling stabilizes a spin-singlet excitonic density phase for weak orbital splitting, positive Hund's coupling favors a spin-triplet excitonic density phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

f-Orbital Mixing in the Octahedral f 2 Compounds UX 6 2– [X = F, Br, Cl, I] and PrCl 6 3–

Understanding how interactions between the f orbitals and ligand orbitals in lanthanide and actinide systems affect their physical properties is the central issue in f-element chemistry. A wide variety of approaches including both theoretical and experimental tools have been used to study these relationships. Among the most widely used tools has been crystal field theory (CFT), which bridges theory and experiment in that it is a model based largely on atomic theory that is parametrized using experimental data. Crystal field theory is quite accurate for the lanthanides, due in part to the highly contracted nature of the 4f orbitals. For actinides, crystal field theory is less accurate, potentially due to the treatment of orbital mixing. In CFT, orbital mixing is handled implicitly by allowing the electron repulsion parameters (Slater F k parameters) and the spin–orbit coupling constant to vary. As a result, orbital mixing in CFT is isotropic in that the F k parameters and the spin–orbit coupling constant affect all f orbitals equally. This approximation works well for the lanthanides due to the limited degree of orbital mixing in these complexes. In actinide complexes, the 5f orbitals have greater overlap with the ligand orbitals, and this approximation is less accurate than in the lanthanides. In this work, we report a modification of CFT that includes the effect of orbital mixing on electron repulsion and spin–orbit coupling for each f orbital. The model is applied to the tetravalent uranium hexahalide dianions and PrCl 6 3– for which the energies of many low-lying excited states are known. The new model generally fits the data as well the traditional CFT although with fewer parameters. However, the new model does not fit the data better than the more complex CFT models of Faucher and co-workers. The results of the model show in detail how changes in overlap and orbital energies influence the energies of the bonding and antibonding orbitals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Nematicity and superconductivity: Competition versus cooperation

Electronic nematic behavior has been identified and studied in iron-based superconductors for some time, particularly in the well-known BaFe 2 As 2 system, where it is well known to compete with superconductivity. On the other hand, it has been shown recently that FeSe displays a negligible effect of nematicity on superconductivity near the superconducting transition, and actual cooperation between the two orders when the system is doped with S. Recently it has also been proposed that LiFeAs undergoes a nematic transition in the superconducting state itself. Generally, we expect superconductivity to be anisotropic when it coexists with nematic order, but it is not clear under what circumstances the two orders compete or cooperate, nor how the anisotropy of the superconducting state correlates with that in the nematic state. To address this, we study a simple mean-field model of a d-wave Pomeranchuk instability together with a mixed s,d pairing interaction, and identify when nematicity is enhanced or suppressed by superconductivity. Here, we show that the competition or cooperation depends significantly on the distortion of the Fermi surface due to nematicity relative to the anisotropy of the superconducting gap function. Further, we discuss the implications of our results for FeSe and LiFeAs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Optimization of transformer ratio and beam loading in a plasma wakefield accelerator with a structure-exploiting algorithm

Plasma-based acceleration has emerged as a promising candidate as an accelerator technology for a future linear collider or a next-generation light source. We consider the plasma wakefield accelerator (PWFA) concept where a plasma wave wake is excited by a particle beam and a trailing beam surfs on the wake. For a linear collider, the energy transfer efficiency from the drive beam to the wake and from the wake to the trailing beam must be large, while the emittance and energy spread of the trailing bunch must be preserved. One way to simultaneously achieve this when accelerating electrons is to use longitudinally shaped bunches and nonlinear wakes. In the linear regime, there is an analytical formalism to obtain the optimal shapes. In the nonlinear regime, however, the optimal shape of the driver to maximize the energy transfer efficiency cannot be precisely obtained because currently no theory describes the wake structure and excitation process for all degrees of nonlinearity. In addition, the ion channel radius is not well defined at the front of the wake where the plasma electrons are not fully blown out by the drive beam. We present results using a novel optimization method to effectively determine a current profile for the drive and trailing beam in PWFA that provides low energy spread, low emittance, and high acceleration efficiency. We parameterize the longitudinal beam current profile as a piecewise-linear function and define optimization objectives. For the trailing beam, the algorithm converges quickly to a nearly inverse trapezoidal trailing beam current profile similar to that predicted by the ultrarelativistic limit of the nonlinear wakefield theory. For the drive beam, the beam profile found by the optimization in the nonlinear regime that maximizes the transformer ratio also resembles that predicted by linear theory. Furthermore, the current profiles found from the optimization method provide higher transformer ratios compared with the linear ramp predicted by the relativistic limit of the nonlinear theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Imaging Nanoscale Energy Transport and Conversion with Ultrafast Electron Microscopy (Final Technical Report)

Light-matter interactions are ubiquitous in nature and reside at the heart of innumerable technologies. The cascade of processes that occur when a material absorbs a photon of light are exceedingly complex and are interwoven in both space and time, rendering precise determination of the atomic-scale and ultrafast mechanisms immensely challenging. The advent of methods for generating short pulses of light several decades ago led to major advances in understanding the initial moments of light absorption and the resultant effects, though directly interrogating the response of the atoms within the material continued to prove challenging. More recently, methods for generating ultrashort pulses of X-rays and fast electrons have opened the way to probing photoinduced structural dynamics of a wide range of matter in multiple phases. An especially promising laboratory-scale method is ultrafast electron microscopy (UEM), wherein the modalities of conventional transmission electron microscopes (imaging, diffraction, spectroscopy) are extended into the femtosecond temporal regime. Here, the ultrafast imaging and diffraction modalities of UEM were used to study the transient structural aspects of photoexcitation of semiconducting materials (e.g., spatially-resolved electron-phonon coupling, excitation and emission of acoustic phonons, and discrete nanoscale scattering processes). The project had three main objectives: (1) determination of the excitation mechanisms of dense, hypersonic charge-carrier waves and the spatially mediated means by which they couple to the lattice via coherent phonon emission, especially with UEM imaging, wherein effects of nanoscale structural and morphological features on the coupling and relaxation dynamics are expected to dictate nucleation sites and preferred wave vectors, (2) elucidation of photoinduced acoustic-phonon seeding, emergence, propagation, and decay over nanoscale crystal regions and especially with respect to local strain fields and atomic-scale disruptions in lattice order, and (3) interwoven with the first two were efforts aimed at realizing combined angstrom-femtosecond spatiotemporal imaging with UEM. The outcomes and impacts of this project were the generation of new knowledge with respect to fundamental light-matter interactions and, especially, the spatially-mediated excitation and evolution of the structural response of materials following coherent photoexcitation. Importantly, the spatial and temporal resolutions of the UEM imaging modalities used are well-suited for such studies and enable spatially-resolved mechanisms to be determined as a function of atomic order, structural features, and morphology. In addition, ultrafast crystallographic measurements were used to correlate real- and reciprocal-space dynamics in order to determine atomic-scale preferential wave vectors and ultrafast scattering mechanisms, especially as dictated by specimen boundary conditions. The obtained results, as detailed in peer-reviewed publications and presentations, illustrate the importance of ultrafast, angstrom-scale real-space imaging for developing a comprehensive understanding of energy transport and conversion in materials.

47 OTHER INSTRUMENTATION↗
NASA NTRS1994

Neutron and electromagnetic emissions during the 1990 May 24 solar flare

In this paper, we are primarily concerened with the solar neutron emission during the 1990 May 24 flare, utilizing the counting rate of the Climax neutron monitor and the time profiles of hard X-rays and gamma-rays obtained with the GRANAT satellite (Pelaez et al., 1992; Talon et al., 1993; Terekhov et al., 1993). We compare the derived neutron injection function with macroscopic parameters of the flare region as obtained from the H-alpha and microwave observations made at the Big Bear Solar Observatory (BBSO) and the Owens Valley Radio Observatory (OVRO) respectively. Our results are summarized as folows: (1) to explain the neutron monitor counting rate and 57.5-110 MeV and 2.2 MeV gamma-ray time profiles, we consider a two-component neutron injection function, Q(E, t). (2) From the H-alpha observations, we find a relatively small loop of length approximattely equal to 2 x 10(exp 4) km, which may be regarded as the source for the fast-decaying component of gamma-rays (57.5-110 MeV) and for the first component of neutron emission. From microwave visibility and the microwave total power spectrum we postulate the presence of a rather big loop (approximately equal to 2 x 10(exp 5) km), which we regard as being responsible for the slow-decaying component of the high-energy emission. We show how the neutron and gamma-ray emission data can be explained in terms of the macroscopic parameters derived from the H-alpha and microwave observations. (3) The H-alpha observations also reveal the presence of a fast mode MHD shock (the Moreton wave) which precedes the microwave peak by 20-30 s and the peak of gamma-ray intensity by 40-50 s. From this relative timing and the single-pulsed time profiles of both radiations, we can attribute the whole event as due to a prompt acceleration of both electrons and protons by the shock and subsequent deceleration of the trapped particles while they propagate inside the magnetic loops.

Kocharov, L. G.↗
DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ∼1/Np with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly non-equilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2014

Energetic Electrons in Dipolarization Events: Spatial Properties and Anisotropy

Using the electromagnetic fields of an MHD simulation of magnetotail reconnection, flow bursts, and dipolarization, we further investigate the acceleration of electrons to suprathermal energies. Particular emphasis is on spatial properties and anisotropies as functions of energy and time. The simulation results are compared with Time History of Events and Macroscale Interactions during Substorms observations. The test particle approach successfully reproduces several observed injection features and puts them into a context of spatial maps of the injection region(s): a dominance of perpendicular anisotropies farther down the tail and closer to the equatorial plane, an increasing importance of parallel anisotropy closer to Earth and at higher latitudes, a drop in energy fluxes at energies below approximately 10 keV, coinciding with the plasma density drop, together with increases at higher energy, a triple peak structure of flux increases near 0 deg, 90 deg, and 180 deg, and a tendency of flux increases to extend to higher energy closer to Earth and at lower latitudes. We identified the plasma sheet boundary layers and adjacent lobes as a main source region for both increased and decreased energetic electron fluxes, related to the different effects of adiabatic acceleration at high and low energies. The simulated anisotropies tend to exceed the observed ones, particularly for perpendicular fluxes at high energies. The most plausible reason is that the MHD simulation lacks the effects of anisotropy-driven microinstabilities and waves, which would reduce anisotropies.

anisotropy↗
NASA NTRS2016

Impacts of Spontaneous Hot Flow Anomalies on the Magnetosheath and Magnetopause

Spacecraft observations and global hybrid (kinetic ions and fluid electrons) simulations have demonstrated that ion dissipation processes at the quasi-parallel bow shock are associated with the formation of structures called spontaneous hot flow anomalies (SHFAs). Previous simulations and recent spacecraft observations have also established that SHFAs result in the formation of magnetosheath filamentary structures(MFS). In this paper we demonstrate that in addition to MFS, SHFAs also result in the formation of magnetos heath cavities that are associated with decreases in density, velocity, and magnetic field and enhancements in temperature. We use the results of a global MHD run to determine the change in the magnetosheath properties associated with cavities due to ion kinetic effects. The results also show the formation of regions of high flow speed called magnetosheath jets whose properties as a function of solar wind Mach number are described in this study. Comparing the properties of the simulated magnetosheath cavities and jets to past spacecraft observations provides good agreement in both cases. We also demonstrate that pressure variations associated with cavities and SHFAs in the sheath result in a continuous sunward and anti sunward magnetopause motion. This result is consistent with previous suggestions that SHFAs may be responsible for the generation of ion cyclotron waves and precipitation of ring current protons in the outer magnetosphere.

Omidi, N.↗
DOE OSTI2022

Accurate localization of Kosterlitz-Thouless-type quantum phase transitions for one-dimensional spinless fermions

We investigate the charge-density wave (CDW) transition for one-dimensional spinless fermions at half band filling with nearest-neighbor electron transfer amplitude t and interaction V. The model is equivalent to the anisotropic XXZ Heisenberg model for which the Bethe Ansatz provides an exact solution. For V>V c =2t, the CDW order parameter and the single-particle gap are finite but exponentially small, as is characteristic for a Kosterlitz-Thouless transition. It is notoriously difficult to locate such infinite-order phase transitions in the phase diagram using approximate analytical and numerical approaches. Second-order Hartree-Fock theory is qualitatively applicable for all interaction strengths, and predicts the CDW transition to occur at V$^{(2)}_{c,2}$≈1.5t. Second-order Hartree Fock theory is almost variational because the density of quasiparticle excitations is small. We apply the density-matrix renormalization group (DMRG) for periodic boundary conditions for system sizes up to 514 sites, which permits a reliable extrapolation of all physical quantities to the thermodynamic limit, apart from the critical region. We investigate the ground-state energy, the gap, the order parameter, the momentum distribution, the quasiparticle density, and the density-density correlation function to locate V c from the DMRG data. In conclusion, tracing the breakdown of the Luttinger liquid and the peak in the quasiparticle density at the band edge permits us to reproduce V c with an accuracy of one percent.

1-dimensional spin chains↗
DOE OSTI2018

An MS-CASPT2 Calculation of the Excited Electronic States of an Axial Difluoroborondipyrromethene (BODIPY) Dimer

The previously reported (Duman et al., J. Org. Chem. 2012, 77, 4516) calculated state energies of monomeric difluoroborondipyrromethene (BODIPY) and its axial dimer would suggest that these dyes are promising candidates for singlet fission, and the dimer was computed to have an unusual low-lying doubly excited state. We find that these results were affected by the use of an imbalanced active space in multireference calculations and are not correct. Multistate complete-active-space second-order perturbation theory (MS-CASPT2/cc-pVDZ) calculations using an [8,8] (8 electrons in 8 orbitals) active space for the monomer and a [16,16] active space for the dimer reproduce quite well the observed excitation energies of the S 1 states of both, and yield T 1 excitation energies well in excess of half of the S 1 excitation energies. We conclude that neither BODIPY monomer nor its axial dimer would permit exothermic singlet fission and are not worthy of investigation as potentially useful candidates, and that the unusual low-energy doubly excited states of the dimer were artifacts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Nature of Oxygen Adsorption on Defective Carbonaceous Materials

Here, plane-wave density functional theory has been used to study oxygen adsorption on graphene, graphite, and (12,0) zigzag single-walled carbon nanotubes with and without Stone–Wales (SW) and single-vacancy (SV) defects to understand the role of defects on carbonaceous material reactivity. Atomic oxygen adsorption leads to the formation of an epoxide on defect-free graphene and graphite and an ether on the exterior wall of carbon nanotubes and SW-defected materials. O 2 chemisorption is endothermic on defect-free graphene and graphite and slightly exothermic on defect-free nanotubes. O 2 chemisorption energies are predicted to be –1.1 to –1.4 eV on an SW defect and –6.0 to –8.0 eV on an SV defect. An SW defect lowers the energy barriers by 0.90 and 0.50 eV for O 2 chemisorption on graphene and nanotubes, respectively. The formation of a C–O–O–C group is important for O 2 dissociation on defect-free and SW-defected materials. The energy barrier is less than 0.30 eV on an SV defect. The more reactive SW defect toward O adsorption on graphene is mostly due to the strained defective carbon atoms being able to donate more electrons to an O to form an ether. The larger 2s character in the hybrid orbitals in an ether than in an epoxide makes the ether C–O bond stronger. Stronger C–O binding on an SW-defective carbon nanotube than on a defect-free nanotube is in part due to more flexibility of the defect to release the epoxide ring strain to form an ether.

25 ENERGY STORAGE↗