Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “excited state methods”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Manipulating Two-Photon Absorption of Molecules through Efficient Optimization of Entangled Light

We report how the unique temporal and spectral features of pulsed entangled photons from a parametric downconversion source can be utilized for manipulating electronic excitations through the optimization of their spectral phase. A new comprehensive optimization protocol based on Bayesian optimization has been developed in this work to selectively excite electronic states accessible by two-photon absorption. Using our optimization method, the entangled two-photon absorption probability for a thiophene dendrimer can be enhanced by up to a factor of 20, while classical light turns out to be nonoptimizable. Moreover, the optimization involving photon entanglement enables selective excitation that would not be possible otherwise. In addition to optimization, we have explored entangled two-photon absorption in the small entanglement time limit showing that entangled light can excite molecular electronic states that are vanishingly small for classical light. As a result, we demonstrate these opportunities with an application to a thiophene dendrimer.

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

Real time evolution for ultracompact Hamiltonian eigenstates on quantum hardware

In this work we present a detailed analysis of variational quantum phase estimation (VQPE), a method based on real-time evolution for ground and excited state estimation on near-term hardware. We derive the theoretical ground on which the approach stands, and demonstrate that it provides one of the most compact variational expansions to date for solving strongly correlated Hamiltonians. At the center of VQPE lies a set of equations, with a simple geometrical interpretation, which provides conditions for the time evolution grid in order to decouple eigenstates out of the set of time evolved expansion states, and connects the method to the classical filter diagonalization algorithm. Further, we introduce what we call the unitary formulation of VQPE, in which the number of matrix elements that need to be measured scales linearly with the number of expansion states, and we provide an analysis of the effects of noise which substantially improves previous considerations. The unitary formulation allows for a direct comparison to iterative phase estimation. Our results mark VQPE as both a natural and highly efficient quantum algorithm for ground and excited state calculations of general many-body systems. We demonstrate a hardware implementation of VQPE for the transverse field Ising model. Further, we illustrate its power on a paradigmatic example of strong correlation (Cr2 in the SVP basis set), and show that it is possible to reach chemical accuracy with as few as ~50 timesteps.

Klymko, Katherine↗
DOE OSTI2026

High-Performance Semiempirical Excited-State Molecular Dynamics Powered by Graphics Processing Units

Here, this Letter introduces excited-state molecular dynamics in PYSEQM, a GPU-accelerated semiempirical quantum chemistry engine implemented in PyTorch. The new module enables Born–Oppenheimer molecular dynamics (BOMD) using configuration-interaction singles and random phase approximation for excited states, allowing long trajectories and large statistical ensembles to be simulated efficiently on a single GPU. We also implement an extended Lagrangian excited-state BOMD (XL-ESMD) scheme that propagates auxiliary electronic variables, enabling relaxed ground and excited-state convergence thresholds without compromising energy conservation. The excited-state BOMD implementation scales smoothly from small chromophores to a nearly 900-atom dendrimer (taking 6.5 s per MD step). PYSEQM also supports batched execution, allowing many geometries or trajectories to be evaluated in a single GPU launch, substantially increasing throughput and making ensemble-based protocols routine. As a demonstration, we compute absorption, emission, and infrared spectra from trajectories propagated on the ground and first excited states. The XL-ESMD scheme yields identical spectra at significantly lower computational cost, establishing the role of extended Lagrangian based dynamics for efficient excited-state BOMD simulations. Beyond raw performance, PYSEQM’s PyTorch foundation provides automatic differentiation for forces, efficient GPU batching, and seamless interfacing with machine learning models. These capabilities position PYSEQM as a practical platform for machine learning-augmented excited-state dynamics and lay the foundation for future data-driven nonadiabatic excited-state dynamics modeling of ultrafast spectroscopic probes.

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

Alternant Hydrocarbon Diradicals as Optically Addressable Molecular Qubits

High-spin molecules allow for bottom-up qubit design and are promising platforms for magnetic sensing and quantum information science. Optical addressability of molecular electron spins has also been proposed in first-row transition-metal complexes via optically detected magnetic resonance (ODMR) mechanisms analogous to the diamond-nitrogen-vacancy color center. However, significantly less progress has been made on the front of metal-free molecules, which can deliver lower costs and milder environmental impacts. At present, most luminescent open-shell organic molecules are π-diradicals, but such systems often suffer from poor ground-state open-shell characters necessary to realize a stable ground-state molecular qubit. In this work, we use alternancy symmetry to selectively minimize radical–radical interactions in the ground state, generating π-systems with high diradical characters. We call them m-dimers, referencing the need to covalently link two benzylic radicals at their meta carbon atoms for the desired symmetry. Through a detailed electronic structure analysis, we find that the excited states of alternant hydrocarbon m-diradicals contain important symmetries that can be used to construct ODMR mechanisms leading to ground-state spin polarization. Furthermore, the molecular parameters are set in the context of a tris(2,4,6-trichlorophenyl)methyl (TTM) radical dimer covalently tethered at the meta position, demonstrating the feasibility of alternant m-diradicals as molecular color centers.

Chemistry↗
DOE OSTI2026

Decay spectroscopy in the neutron-rich Mo–Ru–Pd region

A leading challenge of nuclear-structure research is to understand the properties of nuclides of extreme isospin. Experiments at radioactive-ion-beam facilities, such as the Facility for Rare Isotope Beams in the US, may answer key questions that address diverse topics including fundamental nuclear physics, stellar nucleosynthesis and nuclear applications. The neutron-rich Mo–Ru–Pd (Z = 42 – 46) nuclides are hypothesised to exhibit triaxial-oblate deformation. We performed an experiment with the Facility for Rare Isotope Beams Decay Station initiator (FDSi) to study the structure and decay properties of nuclides in this region. Over 100 different nuclides have been identified in a preliminary analysis of the data. This work presents a first look at several examples between Rb (Z = 37) and Ag (Z = 47). Performance of the FDSi and methods developed to measure ground-state and excited-state lifetimes are presented, and plans for future work are also discussed.

Allmond, James [ORNL] (ORCID:0000000165338721)↗
DOE OSTI2021

Investigation of secondary γ -ray angular distributions using the N 15 ( p , α 1 γ ) C * 12 reaction

The observation of secondary γ-rays provides an alternative method of measuring cross sections that populate excited final states in nuclear reactions. The angular distributions of these γ-rays also provide information on the underlying reaction mechanism. Despite a large amount of data of this type in the literature, publicly available R-matrix codes do not have the ability to calculate these types of angular distributions. In this report, the mathematical formalism derived in Brune and deBoer is implemented in the R-matrix code AZURE2 and calculations are compared with previous data from the literature for the 15 N( p, α1γ ) 12 C* reaction. In addition, new measurements, made at the University of Notre Dame Nuclear Science Laboratory using the Hybrid Array of Gamma Ray Detectors (HAGRiD), are reported that span an energy range from E p = 0.88 to 4.0 MeV. Excellent agreement between data and the phenomenological fit is obtained up to the limit of the previous fit at E p = 2.0 MeV and the R-matrix fit is extended from E x ≈ 13.5 up to E x ≈ 15.3 MeV, where 15 N+p and 12 C+α reactions are fit simultaneously for the first time. An excellent reproduction of the 15 N( p, α1γ ) 12 C* and 12 C(α, α) 12 C data is achieved, but inconsistencies and difficulty in fitting other data is encountered and discussed.

6 ≤ A ≤ 19↗
DOE OSTI2025

Characterizing the Excited States and Electronic Absorption Spectra of Small Alkylperoxy (RO 2 •) and Hydroperoxy (•QOOH) Radicals

Organic peroxy (ROO•) and hydroperoxy (•QOOH) radicals are key reactive intermediates that are formed via the oxidation of volatile organic compounds during combustion or in the Earth’s atmosphere. Their primary fate is continued unimolecular decay or bimolecular chemistry, the relative branching for which is heavily structure- and temperature-dependent. This article outlines a combined single- and multi-reference quantum chemical study to characterize the near-UV accessible electronically excited states of the prototypical ROO• and •QOOH intermediates, tert-butyl peroxy and hydroperoxy-tert-butyl radicals—the ground-state chemistries of which have been well studied both experimentally and computationally. Additionally, we simulate the electronic absorption profiles of these ROO• and •QOOH intermediates with a variety of multi- and single-reference methods. The results show an interesting conformer dependence on the electronically excited-state character and electronic absorption maxima of •QOOH. The results show promise for electronic absorption spectroscopy to be used as a selected probe for determining •QOOH conformers. Additionally, electronic absorption may contribute to the daytime removal of long-lived •QOOH intermediates formed in the troposphere. We expect that our studies will motivate experiments on the electronic absorption spectra of experimentally achievable ROO• and •QOOH.

excited states↗
DOE OSTI2021

CAS without SCF—Why to use CASCI and where to get the orbitals

The complete active space self-consistent field (CASSCF) method has seen broad adoption due to its ability to describe the electronic structure of both the ground and excited states of molecules over a broader swath of the potential energy surface than is possible with the simpler Hartree–Fock approximation. However, it also has a reputation for being unwieldy, computationally costly, and un-black-box. Here, we discuss a class of alternatives, complete active space configuration interaction (CASCI) methods, paying particular attention to their application to electronic excited states. The goal of this Perspective is fourfold. First, we argue that CASCI is not merely an approximation to CASSCF, in that it can be designed to have important qualitative advantages over CASSCF. Second, we present several insights drawn from our experience experimenting with different schemes for computing orbitals to be employed in CASCI. Third, we argue that CASCI is well suited for application to nanomaterials. Lastly, we reason that, with the rise in new low-scaling approaches for describing multireference systems, there is a greater need than ever to develop new methods for defining orbitals that provide an efficient and accurate description of both static correlation and electronic excitations in a limited active space.

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

Adaptive variational quantum eigensolvers for highly excited states

Highly excited states of quantum many-body systems are central objects in the study of quantum dynamics and thermalization that challenge classical computational methods due to their volume-law entanglement content. In this work, we explore the potential of variational quantum algorithms to approximate such states. We propose an adaptive variational quantum eigensolver (VQE) for excited states (X) that self-generates a variational ansatz for arbitrary eigenstates of a many-body Hamiltonian H by attempting to minimize the energy variance with respect to H. We benchmark the method by applying it to an Ising spin chain with integrable and nonintegrable regimes, where we calculate various quantities of interest, including the total energy, magnetization density, and entanglement entropy. We also compare the performance of adaptive VQE-X to an adaptive variant of the folded-spectrum method. For both methods, we find a strong dependence of the algorithm's performance on the choice of operator pool used for the adaptive construction of the ansatz. In particular, an operator pool including long-range two-body gates accelerates the convergence of both algorithms in the nonintegrable regime. Here, we also study the scaling of the number of variational parameters with system size, finding that an exponentially large number of parameters may be necessary to approximate individual highly excited states. Nevertheless, we argue that these methods lay a foundation for the use of quantum algorithms to study finite-energy-density properties of many-body systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Lifetimes of two- and four-neutron-hole states in 204 Pb and 206 Pb

The centroid-shift method was used to determine half-lives of excited levels in 204 Pb and 206 Pb. For the states with 𝐼 𝜋 = 12 + , 17 − , and 19 − in 204 Pb, values of 𝑇 1/2 = 4.9⁢(3), 3.9(3), and 0.9(3) ns, respectively, were determined. The half-life of the 𝐼 𝜋 = 9 − level in 206 Pb was measured to be 1.1(3) ns. In 204 Pb, the 12 + state arises from two neutron holes in the 𝑖 13/2 subshell, while the 17 − and 19 − levels are built on 𝜈⁡(𝑖$^{−3}_{13/2}$, 𝑝$^{−1}_{1/2}$) and 𝜈⁡(𝑖$^{−3}_{13/2}$, 𝑓$^{−1}_{5/2}$) four-neutron-hole configurations, respectively. The 9 − state in 206 Pb results from a 𝜈⁡(𝑖$^{−1}_{13/2}$, 𝑓${−1}_{5/2}$) excitation. As a result, the transition probabilities for the decay of these isomeric levels are reasonably consistent with the expected single-particle values, and these states are described well by shell-model calculations using the kshell code with the KHHE interaction.

190 ≤ A ≤ 219↗
DOE OSTI2020

Binding and excitations in Si x H y molecular systems using quantum Monte Carlo

In this work, we present high-accuracy correlated calculations of small SixHy molecular systems in both the ground and excited states. We employ quantum Monte Carlo (QMC) together with a variety of many-body wave function approaches based on basis set expansions. The calculations are carried out in a valence-only framework using recently derived correlation consistent effective core potentials. Our primary goal is to understand the fixed-node diffusion QMC errors in both the ground and excited states with single-reference trial wave functions. Using a combination of methods, we demonstrate the very high accuracy of the QMC atomization energies being within ≈0.07 eV or better when compared with essentially exact results. By employing proper choices for trial wave functions, we have found that the fixed-node QMC biases for total energies are remarkably uniform ranging between 1% and 3.5% with absolute values at most ≈0.2 eV across the systems and several types of excitations such as singlets and triplets as well as low-lying and Rydberg-like states. Our results further corroborate that Si systems, and presumably also related main group IV and V elements of the periodic table (Ge, Sn, etc), exhibit some of the lowest fixed-node biases found in valence-only electronic structure QMC calculations.

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

Lifetime measurements of 0 + states in Er 168 with the Doppler-shift attenuation method

Here, the lowest-lying shape oscillations of deformed nuclei have been described as quadrupole in nature (λ = 2), resulting in two types of vibrations or oscillations: β vibrations with oscillations along the symmetry axis (K π = 0 + ) and γ vibrations breaking axial symmetry with a projection of K π = 2 + on the symmetry axis. The γ vibration seems to be well characterized as the first K π = $2^+_1$ (or $2^+_γ$) band in deformed nuclei and exhibits a systematic behavior across the region. The nature of the K π = 0 + excitations, however, has remained poorly understood and has been open to debate for some decades. The goal of this work is to understand the nature of 0 + states observed in 168 Er through measurements of the lifetimes of these states and to determine if they are consistent with oscillations built on a deformed ground state, the minima of other coexisting shapes, single-particle states, or a mixture of effects. Lifetimes of excited states in the 168 Er nucleus were measured with the Doppler Shift Attenuation Method (DSAM) and the inelastic neutron scattering reaction, (n, n'γ), at the University of Kentucky Accelerator Laboratory. Numerous 0 + states had been observed by the (p, t) reaction. We confirm the 0 + states at 1217.2, 1421.5, 1833.6, 2364.9, 2392.1, and 2643.0 keV in 168 Er. We could not, however, support the previous assignments of 0 + levels at 2114.1, 2200.6, 2572.5, and 2617.4 keV. We report measured lifetimes for six confirmed 0 + excitations and additional members of 0 + bands. The results for 168 Er show that it is the third excited K π = 0 + ($0^+_4$) excitation that carries the collective strength and, therefore, the potential to be an oscillation on the ground state. This result is similar to the case in 166 Er, where it was also the $0^+_4$ state that exhibited greater collectivity than the first excited K π = 0 + band. The Delaroche et al. prediction for a collective K π = 0 + band is at E T =1.818 MeV which corresponds the third excited K π = 0 + band.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Kondo Lattice Model of Magic-Angle Twisted-Bilayer Graphene: Hund’s Rule, Local-Moment Fluctuations, and Low-Energy Effective Theory

For this work, we apply a generalized Schrieffer-Wolff transformation to the extended Anderson-like topological heavy fermion (THF) model for the magic-angle (θ=1.05°) twisted bilayer graphene (MATBLG) [Phys. Rev. Lett. 129, 047601 (2022)], to obtain its Kondo lattice limit. In this limit localized $\mathcal{f}$ electrons on a triangular lattice interact with topological conduction $\mathcal{c}$ electrons. By solving the exact limit of the THF model, we show that the integer fillings ν=0, ±1, ±2 are controlled by the heavy $\mathcal{f}$ electrons, while ν=±3 is at the border of a phase transition between two $\mathcal{f}$-electron fillings. For ν=0, ±1, ±2, we then calculate the Ruderman-Kittel-Kasuya-Yosida (RKKY) interactions between the $\mathcal{f}$ moments in the full model and analytically prove the SU(4) Hund’s rule for the ground state which maintains that two $\mathcal{f}$ electrons fill the same valley-spin flavor. Our (ferromagnetic interactions in the) spin model dramatically differ from the usual Heisenberg antiferromagnetic interactions expected at strong coupling. We show the ground state in some limits can be found exactly by employing a positive semidefinite “bond-operators” method. We then compute the excitation spectrum of the $\mathcal{f}$ moments in the ordered ground state, prove the stability of the ground state favored by RKKY interactions, and discuss the properties of the Goldstone modes, the (reason for the accidental) degeneracy of (some of) the excitation modes, and the physics of their phase stiffness. We develop a low-energy effective theory for the $\mathcal{f}$ moments and obtain analytic expressions for the dispersion of the collective modes. We discuss the relevance of our results to the spin-entropy experiments in TBG.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Quantum-enhanced photoprotection in neuroprotein architectures emerges from collective light-matter interactions

Background Superradiance is the phenomenon of many identical quantum systems absorbing and/or emitting photons collectively at a higher rate than any one system can individually. This phenomenon has been studied analytically in idealized distributions of electronic two-level systems (TLSs), each with a ground and excited state, as well as numerically in realistic photosynthetic nanotubes and cytoskeletal architectures. Methods Superradiant effects are studied here in idealized toy model systems and realistic biological mega-networks of tryptophan (Trp) molecules, which are strongly fluorescent amino acids found in many proteins. Each Trp molecule acts as a chromophore absorbing in the ultraviolet spectrum and can be treated approximately as a TLS, with its 1 L a excited singlet state; thus, organized Trp networks can exhibit superradiance. Such networks are found, for example, in microtubules, actin filaments, and amyloid fibrils. Microtubules and actin filaments are spiral-cylindrical protein polymers that play significant biological roles as primary constituents of the eukaryotic cytoskeleton, while amyloid fibrils have been targeted in a variety of neurodegenerative diseases. We treat these proteinaceous Trp networks as open quantum systems, using a non-Hermitian Hamiltonian to describe interactions of the chromophore network with the electromagnetic field. We numerically diagonalize the Hamiltonian to obtain its complex eigenvalues, where the real part is the energy and the imaginary part is its associated enhancement rate. We also consider multiple realizations of increasing static disorder in either the site energies or the decay rates. Results We obtained the energies and enhancement rates for realistic microtubules, actin filament bundles, and amyloid fibrils of differing lengths, and we use these values to calculate the quantum yield, which is the ratio of the number of photons emitted to the number of photons absorbed. We find that all three of these structures exhibit highly superradiant states near the low-energy portion of the spectrum, which enhances the magnitude and robustness of the quantum yield to static disorder and thermal noise. Conclusion The high quantum yield and stable superradiant states in these biological architectures may play a photoprotective rolein vivo, downconverting energetic ultraviolet photons—absorbed from those emitted by reactive free radical species—to longer, safer wavelengths and thereby mitigating biochemical stress and photophysical damage. Contrary to conventional assumptions that quantum effects cannot survive in large biosystems at high temperatures, our results suggest that macropolymeric collectives of TLSs in microtubules, actin filaments, and amyloid fibrils exhibit increasingly observable and robust effects with increasing length, up to the micron scale, due to quantum coherent interactions in the single-photon limit. Superradiant enhancement and high quantum yield exhibited in neuroprotein polymers could thus play a crucial role in information processing in the brain, the development of neurodegenerative diseases such as Alzheimer’s and related dementias, and a wide array of other pathologies characterized by anomalous protein aggregates.

Physics↗
DOE OSTI2021

Using projection operators with maximum overlap methods to simplify challenging self-consistent field optimization

Maximum overlap methods are effective tools for optimizing challenging ground- and excited-state wave functions using self-consistent field models such as Hartree-Fock and Kohn-Sham density functional theory. Nevertheless, such models have shown significant sensitivity to the user-defined initial guess of the target wave function. Here in this work, a projection operator framework is defined and used to provide a metric for non-aufbau orbital selection in maximum-overlap-methods. The resulting algorithms, termed the Projection-based Maximum Overlap Method (PMOM) and Projection-based Initial Maximum Overlap Method (PIMOM), are shown to perform exceptionally well when using simple user-defined target solutions based on occupied/virtual molecular orbital permutations. This work also presents a new metric that provides a simple and conceptually convenient measure of agreement between the desired target and the current or final SCF results during a calculation employing a maximum-overlap method.

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

Quantum Solver of Contracted Eigenvalue Equations for Scalable Molecular Simulations on Quantum Computing Devices

The accurate computation of ground and excited states of many-fermion quantum systems is one of the most consequential, contemporary challenges in the physical and computational sciences whose solution stands to benefit significantly from the advent of quantum computing devices. Existing methodologies using phase estimation or variational algorithms have potential drawbacks such as deep circuits requiring substantial error correction or non-trivial high-dimensional classical optimization. In this work, we introduce a quantum solver of contracted eigenvalue equations, the quantum analogue of classical methods for the energies and reduced density matrices of ground and excited states. The solver does not require deep circuits or difficult classical optimization and achieves an exponential speed-up over its classical counterpart. We demonstrate the algorithm though computations on both a quantum simulator and two IBM quantum processing units.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Strategies to Obtain Reliable Energy Landscapes from Embedded Multireference Correlated Wavefunction Methods for Surface Reactions

Embedded correlated wavefunction (ECW) theory is a powerful tool for studying ground- and excited-state reaction mechanisms and associated energetics in heterogeneous catalysis. Several factors are important to obtaining reliable ECW energies, critically the construction of consistent active spaces (ASs) along reaction pathways when using a multireference correlated wavefunction (CW) method that relies on a subset of orbital spaces in the configuration interaction expansion to account for static electron correlation, e.g., complete AS self-consistent field theory, in addition to the adequate partitioning of the system into a cluster and environment, as well as the choice of a suitable basis set and number of states included in excited-state simulations. Here, in this work, we conducted a series of systematic studies to develop best-practice guidelines for ground- and excited-state ECW theory simulations, utilizing the decomposition of NH 3 on Pd(111) as an example. We determine that ECW theory results are relatively insensitive to cluster size, the aug-cc-pVDZ basis set provides an adequate compromise between computational complexity and accuracy, and that a fixed-clean-surface approximation holds well for the derivation of the embedding potential. Additionally, we demonstrate that a merging approach, which involves generating ASs from the molecular fragments at each configuration, is preferable to a creeping approach, which utilizes ASs from adjacent structures as an initial guess, for the generation of consistent potential energy curves involving open-d-shell metal surfaces, and, finally, we show that it is essential to include bands of excited states in their entirety when simulating excited-state reaction pathways.

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

COLUMBUS─An Efficient and General Program Package for Ground and Excited State Computations Including Spin–Orbit Couplings and Dynamics

The COLUMBUS program system provides the tools for performing high-level multireference (MR) computations, including the multireference configuration interaction (MRCI) method and its multireference averaged quadratic coupled cluster (MR-AQCC) extension, allowing computations on a wide range of fascinating atomic and molecular systems, including the treatment of open-shells and complicated excited state phenomena. The inclusion of spin−orbit coupling (SOC) directly within the MRCI step enables the description of systems containing heavy elements, such as lanthanides and actinides, whose properties are strongly influenced by SOC. Analytic energy gradients and nonadiabatic couplings at the correlated MRCI level provide the foundation for a variety of dynamics studies, giving insight into ultrafast photochemistry. New and ongoing method developments in COLUMBUS include the computation of spin densities, improved descriptions of ionic states, enhancements to the AQCC method, and the porting of COLUMBUS to graphical processing units (GPUs). New external interfaces enable an enhanced description of electronic resonances and molecules in strong laser fields. This work highlights these new developments while providing a detailed account of the diverse applications of COLUMBUS in recent years.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗