Search

The transformative impact of modern computational paradigms and technologies, such as high-performance computing (HPC), quantum computing, and cloud computing, has opened up profound new opportunities for scientific simulations. Scalable computational chemistry is one beneficiary of this technological progress. The main focus of this paper is on the performance of various quantum chemical formulations, ranging from low-order methods to high-accuracy approaches, implemented in different computational chemistry packages and libraries, such as NWChem, NWChemEx, Scalable Predictive Methods for Excitations and Correlated Phenomena, ExaChem, and Fermi–Löwdin orbital self-interaction correction on Azure Quantum Elements, Microsoft's cloud services platform for scientific discovery. We pay particular attention to the intricate workflows for performing complex chemistry simulations, associated data curation, and mechanisms for accuracy assessment, which is demonstrated with the Arrows automated workflow for high throughput simulations. Finally, we provide a perspective on the role of cloud computing in supporting the mission of leadership computational facilities. [ABSTRACT FROM AUTHOR]

Broadband blue emission in zero-dimensional perovskites has received considerable attention, which is very important for the realization of stable blue-light emitters; however, the underlying formation mechanism remains unclear. Based on first-principles calculations, we have systematically studied the self-trapped excitons (STEs) behavior and luminescence properties in 0D-(DMA)4PbI6 perovskite. Our calculations show that there is a significant difference between the intrinsic STE luminescence mechanism (∼2.51 eV) and experimental observations (∼2.70 eV). In contrast, we found that the iodine vacancy (VI) is energetically accessible and exhibits a shallow charge transition level at ∼2.69 eV (0/+1) above the valence band maximum, which provides the initial local well for the STEs formation. Moreover, the low electronic dimension synergistic Jahn–Teller distortion facilitates the formation of extrinsic excitons self-trapping. Further excited state electronic structure analysis and configuration coordinate diagram calculations confirmed that the broadband blue emission in 0D-(DMA)4PbI6 is the origin of VI-induced extrinsic STEs instead of intrinsic STEs. Therefore, our simulation results rationalize the experimental phenomena and provide important insights into the formation mechanism of STEs in low-dimensional perovskite systems. [ABSTRACT FROM AUTHOR]

Multi-dimensional quantum mechanical simulations are carried out to understand the multi-state and multi-mode vibronic interactions in the first six low-lying viz., X ̃ 2 B 1 , A ̃ 2 A 1 , B ̃ 2 B 2 , C ̃ 2 A 2 , D ̃ 2 A 1 , and E ̃ 2 B 1 electronic states of c-C2H4O·+. Vibronic coupling theory is applied to study interactions among electronic states using symmetry selection rules. A model 6 × 6 diabatic electronic Hamiltonian is constructed. The parameters of the diabatic Hamiltonian are estimated by performing extensive ab initio electronic structure calculations, using the EOM-IP-CCSD method. The nuclear dynamics calculations are performed with both time-independent and time-dependent quantum mechanical methods. The calculated vibronic structures of six electronic states are found to be in excellent agreement with the available experimental findings. Progressions found in the theoretical spectrum are assigned in terms of vibrational modes. It is found that extremely strong vibronic interactions among the X ̃ 2 B 1 – A ̃ 2 A 1 , B ̃ 2 B 2 – C ̃ 2 A 2 , and D ̃ 2 A 1 – E ̃ 2 B 1 electronic states results into highly overlapping vibronic bands due to multiple multi-state conical intersections. The impact of associated nonadiabatic effects on the vibronic structure and dynamics of the mentioned electronic states is examined at length. Interesting comparison is made with the results obtained for the isomeric acetaldehyde radical cation. [ABSTRACT FROM AUTHOR]

Nonadiabatic quantum dynamics are carried out to illustrate the photoionized spectrum of the cyanopropyne (CH3–C≡C–C≡N) as reported in recent experimental measurements [Lamarre et al., J. Mol. Spectrosc. 315, 206 (2015)]. A detailed electronic structure calculation is performed to analyze the topographical details of the first five ionized states, of which three are degenerate states ( X ̃ 2 E , B ̃ 2 E , and C ̃ 2 E) and two are non-degenerate states ( A ̃ 2 A 1 and D ̃ 2 A 1 ). The degenerate E states of the C3V symmetry molecule are prone to Jahn–Teller (JT) instability, and in addition, symmetry allowed A1 − E vibronic coupling, i.e., pseudo-Jahn–Teller (PJT), effects are expected to have a significant impact in the detailed vibronic structure of these electronic states. The JT splittings of X ̃ 2 E and B ̃ 2 E degenerate states are small, whereas it is quite large at three high frequencies in the C ̃ 2 E electronic states. The large energy separation of X ̃ 2 E from the other states and the non-zero PJT coupling of the B ̃ 2 E state with the close-lying A ̃ 2 A 1 state indicate the uncoupled nature of the X ̃ , A ̃ , and B ̃ vibronic bands of C4H3N. The intersection minima of B ̃ and C ̃ states with the D ̃ state nearly coincide with the energetic minimum of D ̃ state. Therefore, the PJT couplings among these states will lead to a strong vibronic interaction to shape the respective band structure. To completely understand the JT and PJT interactions in the photoionized spectrum of C4H3N, the vibronic coupling model Hamiltonian was constructed to perform nuclear dynamics studies for these electronic states. The vibrational progressions in each vibronic band are identified and compared with the available experimental data in the literature. The impacts of JT and PJT effects in the first five ionized states of cyanopropyne are investigated and discussed in detail. [ABSTRACT FROM AUTHOR]

Recent quantum algorithms pertaining to electronic structure theory primarily focus on the threshold-based dynamic construction of ansatz by selectively including important many-body operators. These methods can be made systematically more accurate by tuning the threshold to include a greater number of operators into the ansatz. However, such improvements come at the cost of rapid proliferation of the circuit depth, especially for highly correlated molecular systems. In this work, we address this issue by the development of a novel theoretical framework that relies on the segregation of an ansatz into a dynamically selected core "principal" component, which is, by construction, adiabatically decoupled from the remaining operators. This enables us to perform computations involving the principal component using extremely shallow-depth circuits, whereas the effect of the remaining "auxiliary" component is folded into the energy function via a cost-efficient non-iterative correction, ensuring the requisite accuracy. We propose a formalism that analytically predicts the auxiliary parameters from the principal ones, followed by a suite of non-iterative auxiliary subspace correction techniques with different levels of sophistication. The auxiliary subspace corrections incur no additional quantum resources yet complement an inadequately expressive core of the ansatz to recover a significant amount of electronic correlations. We have numerically validated the resource efficiency and accuracy of our formalism with a number of strongly correlated molecular systems. [ABSTRACT FROM AUTHOR]

In this work, we explore the electronic properties and chemical bonding in the recently discovered mineral zaykovite, the first natural rhodium selenide Rh3Se4. We comprehensively studied the bulk electronic structure, hybridization of rhodium and selenium orbitals, and the influence of spin–orbit interaction on the electronic spectrum, as well as inspected its topological properties. In addition, we investigated the surface electronic structure of zaykovite and revealed the anisotropic Rashba-type spin splitting in the surface states. In addition, using calculations of the phonon spectra and enthalpy of formation, we predicted the family of similar selenides based on other 4d and 5d transition metals such as Ir, Pd, and Pt. The structural and electronic properties of these materials are discussed. [ABSTRACT FROM AUTHOR]

We extend the DeePMD neural network architecture to predict electronic structure properties necessary to perform non-adiabatic dynamics simulations. While learning the excited state energies and forces follows a straightforward extension of the DeePMD approach for ground-state energies and forces, how to learn the map between the non-adiabatic coupling vectors (NACV) and the local chemical environment descriptors of DeePMD is less trivial. Most implementations of machine-learning-based non-adiabatic dynamics inherently approximate the NACVs, with an underlying assumption that the energy-difference-scaled NACVs are conservative fields. We overcome this approximation, implementing the method recently introduced by Richardson [J. Chem. Phys. 158, 011102 (2023)], which learns the symmetric dyad of the energy-difference-scaled NACV. The efficiency and accuracy of our neural network architecture are demonstrated through the example of the methaniminium cation CH2 N H 2 + . [ABSTRACT FROM AUTHOR]

Vanadium dioxide ( VO 2) as a typically strongly correlated material has an attractive near-room-temperature metal–insulator transition (MIT). This transition shifted between metallic rutile and insulating monoclinic phases can be reversed by external thermal stimulation, accompanied by alterations in the lattice structure and electronic energy band structure, which bring about drastic changes in optical, electrical, thermal, and mechanical properties. Based on these unique characteristics, the VO 2 -integrated element can be utilized in the versatile transduction that converts indirect measurable factors into detectable physical quantities, representing a leap forward in sensing technology and facilitating the development of advanced systems toward an unprecedented level of sensitivity and responsiveness. In this perspective, we review the multifunctional advantages of VO 2 , driven by its thermally induced MIT, which transforms its crystal and electronic structures, leading to changes in optical, electrical, thermal, and mechanical properties. Our paper highlights that VO 2 can serve as a high-performance transduction element, leveraging its attractive properties to facilitate the conversion of various variations into measurable signals, thereby enabling diverse advanced applications. To maximize the advantages of VO 2 , we identify the key challenges and opportunities related to this material, offering guidance and recommendations for future research directions. [ABSTRACT FROM AUTHOR]

In this work, we explore the electronic properties and chemical bonding in the recently discovered mineral zaykovite, the first natural rhodium selenide Rh3Se4. We comprehensively studied the bulk electronic structure, hybridization of rhodium and selenium orbitals, and the influence of spin–orbit interaction on the electronic spectrum, as well as inspected its topological properties. In addition, we investigated the surface electronic structure of zaykovite and revealed the anisotropic Rashba-type spin splitting in the surface states. In addition, using calculations of the phonon spectra and enthalpy of formation, we predicted the family of similar selenides based on other 4d and 5d transition metals such as Ir, Pd, and Pt. The structural and electronic properties of these materials are discussed. [ABSTRACT FROM AUTHOR]

Polarization-dependent photoemission electron microscopy (PD-PEEM) exploits spatial variation in the optical selection rules of materials to image domain formation and material organization on the nanoscale. In this Perspective, we discuss the mechanism of PD-PEEM that results in the observed image contrast in experiments and provide examples of a wide range of material domain structures that PD-PEEM has been able to elucidate, including molecular and polymer domains, local electronic structure and defect symmetry, (anti)ferroelectricity, and ferromagnetism. In the end, we discuss challenges and new directions that are possible with this tool for probing domain structure in materials, including investigating the formation of transient ordered states, multiferroics, and the influence of molecular and polymer order and disorder on excited state dynamics and charge transport. [ABSTRACT FROM AUTHOR]

GPU-accelerated on-the-fly nonadiabatic dynamics is enabled by interfacing the linearized semiclassical dynamics approach with the TeraChem electronic structure program. We describe the computational workflow of the "PySCES" code interface, a Python code for semiclassical dynamics with on-the-fly electronic structure, including parallelization over multiple GPU nodes. We showcase the abilities of this code and present timings for two benchmark systems: fulvene solvated in acetonitrile and a charge transfer system in which a photoexcited zinc-phthalocyanine donor transfers charge to a fullerene acceptor through multiple electronic states on an ultrafast timescale. Our implementation paves the way for an efficient semiclassical approach to model the nonadiabatic excited state dynamics of complex molecules, materials, and condensed phase systems. [ABSTRACT FROM AUTHOR]

Electronic structures are fully determined by the exchange–correlation (XC) potential. In this work, we develop a new method to construct reliable XC potentials by properly mixing the exact exchange and the local density approximation potentials in real space. The spatially dependent mixing parameter is derived based on the derivative discontinuity of electron density and is first-principle. We derived the equations for solving the mixing parameter and proposed an approximation to simplify these equations. Based on this approximation, this new method gives reasonable predictions for the ionization energies, fundamental gaps, and singlet–triplet energy differences for various molecular systems. The impact of the approximation on the constructed XC potentials is examined, and it is found that the quality of the XC potentials can be further improved by removing the approximation. This work demonstrates that the derivative discontinuity of electron density is a promising constraint for constructing high-quality XC potentials. [ABSTRACT FROM AUTHOR]

The intrinsic anomalous Hall effect in ferromagnetic materials depends on the spin–orbit-coupling-induced effect in the electronic structure. The presence of anomalous Hall effect can be realized from the existence of the finite Berry curvature. Since the Heusler material Pt2MnGa is ferromagnetic in its collinear state and no report on its transport properties has yet been made, the anomalous Hall effect in this system is investigated using a combined approach of the density functional theory and maximally localized Wannier function. Analysis of the Berry curvature reveals positive and negative peaks along the high symmetry k-points. These peaks together lead to a significant negative anomalous Hall conductivity, which makes the material a promising candidate for applications in spintronics. [ABSTRACT FROM AUTHOR]

CuGBasis is a free and open-source CUDA®/Python library for efficient computation of scalar, vector, and matrix quantities crucial for the post-processing of electronic structure calculations. CuGBasis integrates high-performance Graphical Processing Unit (GPU) computing with the ease and flexibility of Python programming, making it compatible with a vast ecosystem of libraries. We showcase its utility as a Python library and demonstrate its seamless interoperability with existing Python software to gain chemical insight from quantum chemistry calculations. Leveraging GPU-accelerated code, cuGBasis exhibits remarkable performance, making it highly applicable to larger systems or large databases. Our benchmarks reveal a 100-fold performance gain compared to alternative software packages, including serial/multi-threaded Central Processing Unit and GPU implementations. This paper outlines various features and computational strategies that lead to cuGBasis's enhanced performance, guiding developers of GPU-accelerated code. [ABSTRACT FROM AUTHOR]

In this work, we introduce PyCTRAMER, a comprehensive Python package designed for calculating charge transfer (CT) rate constants in disordered condensed-phase systems at finite temperatures, such as organic photovoltaic (OPV) materials. PyCTRAMER is a restructured and enriched version of the CTRAMER (Charge-Transfer RAtes from Molecular dynamics, Electronic structure, and Rate theory) package [Tinnin et al. J. Chem. Phys. 154, 214108 (2021)], enabling the computation of the Marcus CT rate constant and the six levels of the linearized semiclassical approximations of Fermi's golden rule (FGR) rate constant. It supports various types of intramolecular and intermolecular CT transitions from the excitonic states to CT state. Integrating quantum chemistry calculations, all-atom molecular dynamics (MD) simulations, spin-boson model construction, and rate constant calculations, PyCTRAMER offers an automatic workflow for handling photoinduced CT processes in explicit solvent environments and interfacial CT in amorphous donor/acceptor blends. The package also provides versatile tools for individual workflow steps, including electronic state analysis, state-specific force field construction, MD simulations, and spin-boson model construction from energy trajectories. We demonstrate the software's capabilities through two examples, highlighting both intramolecular and intermolecular CT processes in prototypical OPV systems. [ABSTRACT FROM AUTHOR]

Atomic-scale simulations have progressed tremendously over the past decade, largely thanks to the availability of machine-learning interatomic potentials. These potentials combine the accuracy of electronic structure calculations with the ability to reach extensive length and time scales. The i-PI package facilitates integrating the latest developments in this field with advanced modeling techniques thanks to a modular software architecture based on inter-process communication through a socket interface. The choice of Python for implementation facilitates rapid prototyping but can add computational overhead. In this new release, we carefully benchmarked and optimized i-PI for several common simulation scenarios, making such overhead negligible when i-PI is used to model systems up to tens of thousands of atoms using widely adopted machine learning interatomic potentials, such as Behler–Parinello, DeePMD, and MACE neural networks. We also present the implementation of several new features, including an efficient algorithm to model bosonic and fermionic exchange, a framework for uncertainty quantification to be used in conjunction with machine-learning potentials, a communication infrastructure that allows for deeper integration with electronic-driven simulations, and an approach to simulate coupled photon-nuclear dynamics in optical or plasmonic cavities. [ABSTRACT FROM AUTHOR]

Forte is an open-source library specialized in multireference electronic structure theories for molecular systems and the rapid prototyping of new methods. This paper gives an overview of the capabilities of Forte, its software architecture, and examples of applications enabled by the methods it implements. [ABSTRACT FROM AUTHOR]

Colloidal and nanoparticle self-assembly enables the creation of ordered structures with a variety of electronic and photonic functionalities. The outcomes of the self-assembly processes used to synthesize such structures, however, strongly depend on the uniformity of the individual nanoparticles. Here, we explore the simplest form of particle size dispersity—bidispersity—and its impact on the self-assembly process. We investigate the robustness of self-assembling bcc-type crystals via isotropic interaction potentials in binary systems with increasingly disparate particle sizes by determining their terminal size ratio—the most extreme size ratio at which a mixed binary bcc crystal forms. Our findings show that two-well pair potentials produce bcc crystals that are more robust with respect to particle size ratio than one-well pair potentials. This suggests that an improved self-assembly process is accomplished with a second attractive length scale encoded in the particle–particle interaction, which stabilizes the second-nearest neighbor shell. In addition, we document qualitative differences in the process of ordering and disordering: in bidisperse systems of particles interacting via one-well potentials, we observe a breakdown of order prior to demixing, while in systems interacting via two-well potentials, demixing occurs first and bcc continues to form in parts of the droplet down to low size ratios. [ABSTRACT FROM AUTHOR]

We investigate the endofullerene system Ne@C70 by constructing a three-dimensional Potential Energy Surface (PES) describing the translational motion of the Ne atom. This is constructed from electronic structure calculations from a plethora of methods, including MP2, SCS-MP2, SOS-MP2, RPA@PBE, and C(HF)-RPA, which were previously used for He@C60 in Panchagnula et al. [J. Chem. Phys. 160, 104303 (2024)], alongside B86bPBE-25X-XDM and B86bPBE-50X-XDM. The reduction in symmetry moving from C60 to C70 introduces a double well potential along the anisotropic direction, which forms a test of the sensitivity and effectiveness of the electronic structure methods. The nuclear Hamiltonian is diagonalized using a symmetrized double minimum basis set outlined in Panchagnula and Thom [J. Chem. Phys. 159, 164308 (2023)], with translational energies having error bars ±1 and ±2 cm−1. We find no consistency between electronic structure methods as they find a range of barrier heights and minima positions of the double well and different translational eigenspectra, which also differ from the Lennard-Jones (LJ) PES given in Mandziuk and Bačić [J. Chem. Phys. 101, 2126–2140 (1994)]. We find that generating effective LJ parameters for each electronic structure method cannot reproduce the full PES nor recreate the eigenstates, and this suggests that the LJ form of the PES, while simple, may not be best suited to describe these systems. Even though MP2 and RPA@PBE performed best for He@C60, due to the lack of concordance between all electronic structure methods, we require more experimental data in order to properly validate the choice. [ABSTRACT FROM AUTHOR]

We investigate the endofullerene system Ne@C70 by constructing a three-dimensional Potential Energy Surface (PES) describing the translational motion of the Ne atom. This is constructed from electronic structure calculations from a plethora of methods, including MP2, SCS-MP2, SOS-MP2, RPA@PBE, and C(HF)-RPA, which were previously used for He@C60 in Panchagnula et al. [J. Chem. Phys. 160, 104303 (2024)], alongside B86bPBE-25X-XDM and B86bPBE-50X-XDM. The reduction in symmetry moving from C60 to C70 introduces a double well potential along the anisotropic direction, which forms a test of the sensitivity and effectiveness of the electronic structure methods. The nuclear Hamiltonian is diagonalized using a symmetrized double minimum basis set outlined in Panchagnula and Thom [J. Chem. Phys. 159, 164308 (2023)], with translational energies having error bars ±1 and ±2 cm−1. We find no consistency between electronic structure methods as they find a range of barrier heights and minima positions of the double well and different translational eigenspectra, which also differ from the Lennard-Jones (LJ) PES given in Mandziuk and Bačić [J. Chem. Phys. 101, 2126–2140 (1994)]. We find that generating effective LJ parameters for each electronic structure method cannot reproduce the full PES nor recreate the eigenstates, and this suggests that the LJ form of the PES, while simple, may not be best suited to describe these systems. Even though MP2 and RPA@PBE performed best for He@C60, due to the lack of concordance between all electronic structure methods, we require more experimental data in order to properly validate the choice. [ABSTRACT FROM AUTHOR]

Using the density-functional theory approach, including Hubbard U correction, we investigate the defect structures consisting of vanadium (V) atoms embedded in a monolayer silicene. Specifically, we consider V–V atom pairs in antiferromagnetic (AFM), ferromagnetic (FM), and non-magnetic states, which are embedded in substitutional and interstitial sites. We determine the ground-state structures, formation and binding energies, electronic structures, induced magnetization, as well as the spin-exchange coupling between the V–V pair. For the substitutional vanadium atom pair, the stability of the AFM and FM spin configurations depends on the sublattice sites in which the V atoms are sited. When the V pair is located on a similar sublattice site type, the AFM spin alignment is more energetically favored, whereas when the pair is located in a different sublattice site, the FM interactions are more stable. However, the relative stability of the AFM or FM configurations changes rapidly as the separation between the V pair increases. Regarding the interstitial-hole V–V pair configurations, the most stable structure is when the pair is at the nearest-neighbor hole sites and is in an FM alignment. Also, at larger separations, the AFM or FM hole configurations are approximately degenerate in energy. Furthermore, we elucidate on the Ruderman–Kittel–Kasuya–Yosida, direct-exchange, and the superexchange interaction mechanisms in the vanadium-embedded silicene. In addition, we estimate a Curie temperature (Tc) of up to ∼500 K for a silicene structure containing a V pair in the FM spin alignment. Such a high Tc, in addition to the stability of the material, suggests that vanadium-embedded silicene is a potential candidate material for spintronic device applications. [ABSTRACT FROM AUTHOR]

We investigate the effects of solvents on the geometry, absorption spectrum, and first hyperpolarizability of six push–pull molecules, each containing a 4,5-dicyanoimidazole group as an electron acceptor and a N,N-dimethylamino group as an electron donor, with systematically extended π-conjugated systems. Geometry optimizations in dichloromethane, methanol, water, and formamide under normal thermodynamic conditions were performed using the average solvent electrostatic configuration–free energy gradient method, which employs a discrete solvent model. The conformational structure of molecules is moderately affected by the environment, with the π-conjugated system becoming more planar in protic solvents. Solvent effects on the first hyperpolarizability result in marked increases that are in line with the red shifts of the absorption spectrum. The hyperpolarizability of smaller molecules within the set may be significantly influenced by the effects of geometry relaxation in highly polar protic solvents. The results illustrate the role of hydrogen bonds in the structure and electronic properties of push–pull molecules in protic environments. For smaller molecules, hydrogen bonds significantly contribute to enhancing the hyperpolarizability, but the effect of these specific interactions becomes less significant with the length of the π-conjugated system. [ABSTRACT FROM AUTHOR]

MoSi2N4 is a two-dimensional ternary nitride semiconductor that has attracted attention for its excellent mechanical and thermal properties. Theoretical studies predict that zigzag edges of this material can host magnetic edge states and Dirac fermions, but the stability of such edges has not been examined. Here, we present a density functional theory study of the electronic and thermodynamic properties of MoSi2N4 edges. We develop a (partial) ternary phase diagram that identifies a region of chemical potentials within which MoSi2N4 is stable over competing elemental or binary phases. Based on this phase diagram, we determine the thermodynamic stability of several armchair and zigzag edges and elucidate their electronic structures. Bare zigzag edges, predicted to host exotic electronic states, are found to be substantially higher in energy than armchair edges and, thus, unlikely to occur in practice. However, with hydrogen passivation, these zigzag edges can be stabilized relative to their armchair counterparts while retaining metallicity and magnetic order. Our analysis provides a solid thermodynamic basis for further exploration of MoSi2N4 in nanoscale electronics and spintronics. [ABSTRACT FROM AUTHOR]

MoSi2N4 is a two-dimensional ternary nitride semiconductor that has attracted attention for its excellent mechanical and thermal properties. Theoretical studies predict that zigzag edges of this material can host magnetic edge states and Dirac fermions, but the stability of such edges has not been examined. Here, we present a density functional theory study of the electronic and thermodynamic properties of MoSi2N4 edges. We develop a (partial) ternary phase diagram that identifies a region of chemical potentials within which MoSi2N4 is stable over competing elemental or binary phases. Based on this phase diagram, we determine the thermodynamic stability of several armchair and zigzag edges and elucidate their electronic structures. Bare zigzag edges, predicted to host exotic electronic states, are found to be substantially higher in energy than armchair edges and, thus, unlikely to occur in practice. However, with hydrogen passivation, these zigzag edges can be stabilized relative to their armchair counterparts while retaining metallicity and magnetic order. Our analysis provides a solid thermodynamic basis for further exploration of MoSi2N4 in nanoscale electronics and spintronics. [ABSTRACT FROM AUTHOR]

The plane-wave pseudopotential (PW-PP) formalism is widely used for the first-principles electronic structure calculation of extended periodic systems. The PW-PP approach has also been adapted for real-time time-dependent density functional theory (RT-TDDFT) to investigate time-dependent electronic dynamical phenomena. In this work, we detail recent advances in the PW-PP formalism for RT-TDDFT, particularly how maximally localized Wannier functions (MLWFs) are used to accelerate simulations using the exact exchange. We also discuss several related developments, including an anti-Hermitian correction for the time-dependent MLWFs (TD-MLWFs) when a time-dependent electric field is applied, the refinement procedure for TD-MLWFs, comparison of the velocity and length gauge approaches for applying an electric field, and elimination of long-range electrostatic interaction, as well as usage of a complex absorbing potential for modeling isolated systems when using the PW-PP formalism. [ABSTRACT FROM AUTHOR]

Herein, we explore, from a theoretical perspective, the nonradiative photoinduced processes (charge separation and energy transfer) within a family of donor–acceptor supramolecular complexes based on the electron-donor truxene-tetrathiafulvalene (truxTTF) derivative and a series of curved fullerene fragments (buckybowls) of different shapes and sizes (C30H12, C32H12, and C38H14) as electron acceptors that successfully combine with truxTTF via non-covalent interactions. The resulting supramolecular complexes (truxTTF·C30H12, truxTTF·C32H12, and truxTTF·C38H14) undergo charge-separation processes upon photoexcitation through charge-transfer states involving the donor and acceptor units. Despite the not so different size of the buckybowls, they present noticeable differences in the charge-separation efficiency owing to a complex decay post-photoexcitation mechanism involving several low-lying excited states of different natures (local and charge-transfer excitations), all closely spaced in energy. In this intricate scenario, we have adopted a theoretical approach combining electronic structure calculations at (time-dependent) density functional theory, a multistate multifragment diabatization method, the Marcus–Levitch–Jortner semiclassical rate expression, and a kinetic model to estimate the charge separation rate constants of the supramolecular heterodimers. Our outcomes highlight that the efficiency of the photoinduced charge-separation process increases with the extension of the buckybowl backbone. The supramolecular heterodimer with the largest buckybowl (truxTTF·C38H14) displays multiple and efficient electron-transfer pathways, providing a global photoinduced charge separation in the ultrafast time scale in line with the experimental findings. The study reported indicates that modifications in the shape and size of buckybowl systems can give rise to attractive novel acceptors for potential photovoltaic applications. [ABSTRACT FROM AUTHOR]

Modern software engineering of electronic structure codes has seen a paradigm shift from monolithic workflows toward object-based modularity. Software objectivity allows for greater flexibility in the application of electronic structure calculations, with particular benefits when integrated with approaches for data-driven analysis. Here, we discuss different approaches to create deep modular interfaces that connect big-data workflows and electronic structure codes and explore the diversity of use cases that they can enable. We present two such interface approaches for the semi-empirical electronic structure package, DFTB+. In one case, DFTB+ is applied as a library and provides data to an external workflow; in another, DFTB+receives data via external bindings and processes the information subsequently within an internal workflow. We provide a general framework to enable data exchange workflows for embedding new machine-learning-based Hamiltonians within DFTB+ or enabling deep integration of DFTB+ in multiscale embedding workflows. These modular interfaces demonstrate opportunities in emergent software and workflows to accelerate scientific discovery by harnessing existing software capabilities. [ABSTRACT FROM AUTHOR]

As a newly synthesized two-dimensional (2D) carbon material, monolayer quasi-hexagonal phase fullerene (qHP C60) has an excellent electronic structure and low thermal conductivity. qHP C60 attracted significant attention from scientists because it has potential applications in thermoelectric materials. Thermoelectric properties of 2D materials significantly depend on the transport of carriers (such as electrons and phonons), and strain engineering is an essential method for modulating the transport of electrons and phonons in 2D materials. However, the strain engineering method for the modulation of the thermoelectric properties of monolayer qHP C60 has not been reported yet. In the present paper, the first-principles combined with the non-equilibrium Green's function method are used to investigate the ballistic transport properties of electrons and phonons in monolayer qHP C60. The effects of temperature, chemical potential, and biaxial tensile strain on the thermoelectric transport parameters (including conductivity, Seebeck coefficient, power factor, and thermal conductivity) as well as the figure of merit (ZT) of monolayer qHP C60 are presented, compared, discussed, and analyzed. We found that monolayer qHP C60 exhibits anisotropic characteristics in electron and phonon transport properties, showcasing outstanding thermoelectric properties. The distinctive quasi-hexagonal phase fullerene network structure offers a novel platform for exploring innovative 2D thermoelectric materials in research. This study provides crucial theoretical insights to guide the designing and implementation of 2D thermoelectric materials based on fullerenes. [ABSTRACT FROM AUTHOR]

We present an isotropic ab initio (para − H 2 ) 4 four-body interaction potential energy surface (PES). The electronic structure calculations are performed at the correlated coupled-cluster theory level, with single, double, and perturbative triple excitations. They use an atom-centered augmented correlation-consistent double zeta basis set, supplemented by a (3 s 3 p 2 d) midbond function. We use a multilayer perceptron to construct the PES. We apply a rescaling transformation to the output energies during training to improve the prediction of weaker energies in the sample data. At long distances, the interaction energies are adjusted to match the empirically derived four-body dispersion interaction. The four-body interaction energy at short intermolecular separations is net repulsive. The use of this four-body PES, in combination with a first principles pair potential for para − H 2 [J. Chem. Phys. 119, 12551 (2015)] and an isotropic ab initio three-body potential for para − H 2 [J. Chem. Phys. 156, 044301 (2022)], is expected to provide closer agreement with experimental results. [ABSTRACT FROM AUTHOR]

Using resonant Raman spectroscopy with 632.8 nm (1.96 eV) laser, we describe the Raman spectra of single crystals of transition metal dichalcogenides with nominal composition MoSe(2−x)Sx for x = 0.2. Changes in Raman spectra at some regions of the sample indicated non-stoichiometry and, in particular, chalcogenide vacancies. At low temperatures around 77 K, we observed unusual temperature dependent enhancement in the intensity of non-zone center modes as well as overtones and combination modes in Raman spectra. This enhancement in the intensity is correlated to the resonance achieved in the non-stoichiometric regions of the crystal at low temperatures. Observed resonance is attributed to modification in the electronic structure due to defects. Energy dispersive x-ray spectroscopy measurements confirmed chalcogenide vacancies in the crystals. The change in the electronic structure due to defects is also corroborated by photoluminescence spectroscopy measurements. [ABSTRACT FROM AUTHOR]

Nanoscale semiconductors with isolated spin impurities have been touted as promising materials for their potential use at the intersection of quantum, spin, and information technologies. Electron paramagnetic resonance (EPR) studies of spins in semiconducting carbon nanotubes have overwhelmingly focused on spins more strongly localized by sp3-type lattice defects. However, the creation of such impurities is irreversible and requires specific reactions to generate them. Shallow charge impurities, on the other hand, are more readily and widely produced by simple redox chemistry, but have not yet been investigated for their spin properties. Here, we use EPR to study p-doped (6,5) semiconducting single-wall carbon nanotubes (s-SWNTs) and elucidate the role of impurity–impurity interactions in conjunction with exchange and correlation effects for the spin behavior of this material. A quantitative comparison of the EPR signals with phenomenological modeling combined with configuration interaction electronic structure calculations of impurity pairs shows that orbital overlap, combined with exchange and correlation effects, causes the EPR signal to disappear due to spin entanglement for doping levels corresponding to impurity spacings of 14 nm (at 30 K). This transition is predicted to shift to higher doping levels with increasing temperature and to lower levels with increasing screening, providing an opportunity for improved spin control in doped s-SWNTs. [ABSTRACT FROM AUTHOR]

Endohedrally doped clusters form a large category of cage clusters, with unique structures, diverse elemental compositions, and highly tunable electronic structures and physisochemical properties. They have been widely achieved in laboratory and may serve as functional building blocks for assembling new supermolecular structures and devices. In this paper, for the first time, we disclosed the luminescence properties of endohedrally doped group-IV clusters by time-dependent density functional theory calculations. A total of 64 cage clusters have been explored in terms of stability, emission wavelength, and the energy difference between the first excited singlet and triplet states. The key geometric and electronic factors governing the photophysical properties of these cage clusters were unveiled, to provide crucial insights for crafting atomically precise nanoclusters for optical and optoelectronic applications. [ABSTRACT FROM AUTHOR]

The software for chemical interaction networks (SCINE) project aims at pushing the frontier of quantum chemical calculations on molecular structures to a new level. While calculations on individual structures as well as on simple relations between them have become routine in chemistry, new developments have pushed the frontier in the field to high-throughput calculations. Chemical relations may be created by a search for specific molecular properties in a molecular design attempt, or they can be defined by a set of elementary reaction steps that form a chemical reaction network. The software modules of SCINE have been designed to facilitate such studies. The features of the modules are (i) general applicability of the applied methodologies ranging from electronic structure (no restriction to specific elements of the periodic table) to microkinetic modeling (with little restrictions on molecularity), full modularity so that SCINE modules can also be applied as stand-alone programs or be exchanged for external software packages that fulfill a similar purpose (to increase options for computational campaigns and to provide alternatives in case of tasks that are hard or impossible to accomplish with certain programs), (ii) high stability and autonomous operations so that control and steering by an operator are as easy as possible, and (iii) easy embedding into complex heterogeneous environments for molecular structures taken individually or in the context of a reaction network. A graphical user interface unites all modules and ensures interoperability. All components of the software have been made available as open source and free of charge. [ABSTRACT FROM AUTHOR]

The interplay between orientation transition and chiral self-assemblies of para-terphenyl (P3P) molecules on the Cd(0001) surface has been investigated using low temperature scanning tunneling microscopy and density functional theory calculations. Three distinct molecular orientations have been discerned from the self-assembled thin films of P3P. At the low coverage, flat-lying molecules appear in the homochiral domains with the incommensurate registry to the substrate. With the coverage increasing, the incoming molecules are incorporated into the first layer with edge-on orientation and form the self-assembled zigzag chains. The alternative arrangement of zigzag chains with opposite chirality leads to the formation of a c(4 × 2) superstructure, in which the tilted molecules exhibit orientational frustration and fuzzy noises. The analysis of the tunneling spectra reveals that the electronic structure of P3P layers is contingent upon the hybridization between the electronic states of P3P molecules and the Cd(0001) surface. These results provide important insights into the interplay between orientational transition and chiral assembly of P3P molecules on metal substrates. [ABSTRACT FROM AUTHOR]

We optimize the internuclear geometry and electronic structure of a model chiral system to achieve a maximal photoelectron circular dichroism (PECD) in its one-photon ionization by circularly polarized light. The electronic structure calculations are performed by the single center method, while the optimization is done using quantum alchemy employing a Taylor series expansion. Thereby, the effect of bond lengths and uncompensated charge distributions on the chiral response of the model is investigated theoretically in some detail. It is demonstrated that manipulating a chiral asymmetry of the ionic potential may enhance the dichroic parameter (i.e., the PECD) of the randomly oriented model system well beyond β1 = 25%. Furthermore, we demonstrate that quantum alchemy is applicable to PECD despite the unusually strong coupling of spatial and electronic degrees of freedom and discuss the relative impact of the individual degrees of freedom in this model system. We define the necessary conditions for the computational design of PECD for real (non-model) chiral molecules using our approach. [ABSTRACT FROM AUTHOR]

A major difficulty in quantum simulation is the adequate treatment of a large collection of entangled particles, synonymous with electron correlation in electronic structure theory, with coupled cluster (CC) theory being the leading framework for dealing with this problem. Augmenting computationally affordable low-rank approximations in CC theory with a perturbative account of higher-rank excitations is a tractable and effective way of accounting for the missing electron correlation in those approximations. This is perhaps best exemplified by the "gold standard" CCSD(T) method, which bolsters the baseline CCSD with the effects of triple excitations using considerations from many-body perturbation theory (MBPT). Despite this established success, such a synergy between MBPT and the unitary analog of CC theory (UCC) has not been explored. In this work, we propose a similar approach wherein converged UCCSD amplitudes are leveraged to evaluate energy corrections associated with triple excitations, leading to the UCCSD[T] method. In terms of quantum computing, this correction represents an entirely classical post-processing step that improves the energy estimate by accounting for triple excitation effects without necessitating new quantum algorithm developments or increasing demand for quantum resources. The rationale behind this choice is shown to be rigorous by studying the properties of finite-order UCC energy functionals, and our efforts do not support the addition of the fifth-order contributions as in the (T) correction. We assess the performance of these approaches on a collection of small molecules and demonstrate the benefits of harnessing the inherent synergy between MBPT and UCC theories. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The rapid development of deep learning techniques has driven the emergence of a neural network-based variational Monte Carlo (VMC) method (referred to as FermiNet), which has manifested high accuracy and strong predictive power in the electronic structure calculations of atoms, molecules, and some periodic systems. Recently, the implementation of the effective core potential (ECP) scheme has further facilitated more efficient calculations in practice. However, there is still a lack of comprehensive assessments of the ECP's performance under the FermiNet. In this work, we set sail to fill this gap by conducting extensive tests on the first two row elements regarding their atomic, spectral, and molecular properties. Our major finding is that, in general, the qualities of ECPs have been correctly reflected under FermiNet. Two recently built ECP tables, namely, correlation consistent ECP (ccECP) and energy consistent correlated electron pseudopotential (eCEPP), seem to prevail in terms of overall performance. In particular, ccECP performs slightly better on spectral precision and covers more elements, while eCEPP is more systematically built from both shape and energy consistency and better treats the core polarization. On the other hand, the high accuracy of the all-electron calculations is hindered by the absence of relativistic effects as well as the numerical instabilities in some heavier elements. Finally, with further in-depth discussions, we generate possible directions for developing and improving FermiNet in the near future. [ABSTRACT FROM AUTHOR]

This paper reports on the effects of shear rate and interface modeling parameters on the hydrodynamic slip length (LS) for water–graphite interfaces calculated using non-equilibrium molecular dynamics. Five distinct non-bonded solid–liquid interaction parameters were considered to assess their impact on LS. The interfacial force field derivations included sophisticated electronic structure calculation-informed and empirically determined parameters. All interface models exhibited a similar and bimodal LS response when varying the applied shear rate. LS in the low shear rate regime (LSR) is in good agreement with previous calculations obtained through equilibrium molecular dynamics. As the shear rate increases, LS sharply increases and asymptotes to a constant value in the high shear regime (HSR). It is noteworthy that LS in both the LSR and HSR can be characterized by the density depletion length, whereas solid–liquid adhesion metrics failed to do so. For all interface models, LHSR calculations were, on average, ∼28% greater than LLSR, and this slip jump was confirmed using the SPC/E and TIP4P/2005 water models. To address the LS transition from the LSR to the HSR, the viscosity of water and the interfacial friction coefficient were investigated. It was observed that in the LSR, the viscosity and friction coefficient decreased at a similar rate, while in the LSR-to-HSR transition, the friction coefficient decreased at a faster rate than the shear viscosity until they reached a new equilibrium, hence explaining the LS-bimodal behavior. This study provides valuable insights into the interplay between interface modeling parameters, shear rate, and rheological properties in understanding hydrodynamic slip behavior. [ABSTRACT FROM AUTHOR]

Porous, stacked two-dimensional covalent organic frameworks (2D COFs) bearing semiconducting linkers can support directional charge transfer across adjacent layers of the COF. To better inform the current and possible future design rules for enhancing electron and hole transport in such materials, an understanding of how linker selection and functionalization affects interlayer electronic couplings is essential. We report electronic structure simulation and analysis of electronic couplings across adjacent linker units and to encapsulated species in functionalized electroactive 2D COFs. The detailed dependence of these electronic couplings on interlayer interactions is examined through scans along key interlayer degrees of freedom and through configurational sampling from equilibrium molecular dynamics on semiempirical potential energy surfaces. Beyond affirming the sensitivity of the electronic coupling to interlayer distance and orientation, these studies offer guidance toward linker functionalization strategies for enhancing charge carrier transport in electroactive 2D COFs. [ABSTRACT FROM AUTHOR]

In this work, we utilize our recently developed machine learning (ML)-corrected ab initio dispersion (aiD) potential, known as D3-ML, which is based on the comprehensive SAPT10K dataset and relies solely on Cartesian coordinates as input, to address the dispersion deficiencies in second-order Møller−Plesset perturbation theory (MP2) by replacing its problematic dispersion and exchange-dispersion terms with D3-ML. This leads to the development of a new dispersion-corrected MP2 method, MP2+aiD(CCD), which outperforms other spin-component-scaled and dispersion-corrected MP2 methods as well as popular ML models for predicting noncovalent interactions across various datasets, including S66 × 8, NAP6 (containing 6 naphthalene dimers), L7, S12L, DNA−ellipticine, the C60 dimer, and C60[6]CPPA. In addition, MP2+aiD(CCD) exhibits comparable or even superior performance compared to the contemporary ωB97M-V functional. The limited performance of pure ML models for systems outside the training set or larger than those in the training set highlights their instability and unpredictability. Conversely, the outstanding performance and transferability of the hybrid MP2+aiD(CCD) method can be attributed to the fusion of the physical electronic structure method and a data-driven ML model, combining the strengths of both sides. This investigation firmly establishes MP2+aiD(CCD) as one of the most accurate and reliable fifth-order scaling correlated wave function methods currently available for modeling noncovalent interactions, even for large complexes. MP2+aiD(CCD) is expected to be reliably applicable in investigating real-life complexes at the hundred-atom scale. [ABSTRACT FROM AUTHOR]

We present an exchange–correlation approximation in which the Coulomb interaction is split into long- and short-range components and the range separation is determined by a non-empirical density functional. The functional respects important constraints, such as the homogeneous and slowly varying density limits, leads to the correct long-range potential, and eliminates one-electron self-interaction. Our approach is designed for spectroscopic purposes and closely approximates the piecewise linearity of the energy as a function of the particle number. The functional's accuracy for predicting the fundamental gap in generalized Kohn–Sham theory is demonstrated for a large number of systems, including organic semiconductors with a notoriously difficult electronic structure. [ABSTRACT FROM AUTHOR]

The ground state of 3d1 MnO42− was studied by density functional theory (DFT) and complete active space self-consistent field (CASSCF) methods in terms of a variety of molecular point group structures to ascertain the site and Jahn–Teller (JT) distortion effect. Modeling results from UB3LYP/6-31+G(d) calculations with natural bond orbital analysis show the four Mn–O bonds are coordinate covalent. The one-electron matrix elements from CASSCF(AILFT) (ab initio ligand field theory) with a second order perturbation treatment were used to calculate the parameters of the angular overlap model. These allowed the calculation of JT stabilization energies and first and second order JT coupling constants for MnO42− with C2v and D2d symmetry. Absorption spectra and excited state transition energies were calculated assuming the lattice distorted geometry was C2v with time-dependent DFT (TDDFT) using the unrestricted coulomb-attenuating UCAM-B3LYP density functional and with the spin-adapted spin-flip DFT using the UBH&HLYP density functional, both with the AUG-CC-PVTZ basis set. The mean absolute deviation of the calculations from experimental excited states for the ligand field and ligand-to-metal charge transfer (LMCT) bands was better than 0.1 eV. Several new assignments for LMCT excited states were made on the basis of the TDDFT excitation results. [ABSTRACT FROM AUTHOR]

A proof-of-concept study is undertaken to demonstrate the utility of the machine learning combined with the thermodynamic perturbation theory (MLPT) to test the accuracy of electronic structure methods in finite-temperature thermodynamic calculations. As a test example, formic acid dimer is chosen, which is one of the systems included in the popular benchmark set S22 [Jurečka et al., Phys. Chem. Chem. Phys. 8, 1985–1993 (2006)]. Starting from the explicit molecular dynamics and thermodynamic integration performed at the PBE + D2 level, the MLPT is used to obtain fully anharmonic dimerization free and internal energies at the reference quality CCSD(T) level and 19 different density functional approximations, including GGA, meta-GGA, non-local, and hybrid functionals with and without dispersion corrections. Our finite-temperature results are shown to be both qualitatively and quantitatively different from those obtained using the conventional benchmarking strategy based on fixed structures. The hybrid functional HSE06 is identified as the best performing approximate method tested, with the errors in free and internal energies of dimerization being 36 and 41 meV, respectively. [ABSTRACT FROM AUTHOR]

Open-source APOST-3D software features a large number of wavefunction analysis tools developed over the past 20 years, aiming at connecting classical chemical concepts with the electronic structure of molecules. APOST-3D relies on the identification of the atom in the molecule (AIM), and several analysis tools are implemented in the most general way so that they can be used in combination with any chosen AIM. Several Hilbert-space and real-space (fuzzy atom) AIM definitions are implemented. In general, global quantities are decomposed into one- and two-center terms, which can also be further grouped into fragment contributions. Real-space AIM methods involve numerical integrations, which are particularly costly for energy decomposition schemes. The current version of APOST-3D features several strategies to minimize numerical error and improve task parallelization. In addition to conventional population analysis of the density and other scalar fields, APOST-3D implements different schemes for oxidation state assignment (effective oxidation state and oxidation states localized orbitals), molecular energy decomposition schemes, and local spin analysis. The APOST-3D platform offers a user-friendly interface and a comprehensive suite of state-of-the-art tools to bridge the gap between theory and experiment, representing a valuable resource for both seasoned computational chemists and researchers with a focus on experimental work. We provide an overview of the code structure and its capabilities, together with illustrative examples. [ABSTRACT FROM AUTHOR]

We here report the electronic structure calculation of GaInSe3 and GaInTe3 monolayers with the P3m1 (no. 156) space group. The electronic structure and thermoelectric properties of the monolayers are calculated through the Vienna Ab initio Simulation Package and BoltzTraP2 codes. The dynamic and thermodynamic stabilities were verified by calculating their phonon spectra and simulating ab initio molecular dynamics. The monolayers were found to have a direct bandgap, with both PBE + SOC and HSE06 + SOC potentials. The lattice thermal conductivity of GaInTe3 monolayer calculated using Phono3py code shows ultra-low values due to enhanced phonon–phonon scattering. Combining electrical and thermal transport, the values have been evaluated. Importantly, the p-type GaInTe3 has excellent thermoelectric properties at 700 K, with a z T value of 2, indicating that the p-type GaInTe3 has potential application in the field of thermoelectricity. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Recently, some of us reviewed and studied the photoionization dynamics of C60 that are of great interest to the astrochemical community as four of the diffuse interstellar bands (DIBs) have been assigned to electronic transitions in the C60+ cation. Our previous analysis of the threshold photoelectron spectrum (TPES) of C60 [Hrodmarsson et al., Phys. Chem. Chem. Phys. 22, 13880–13892 (2020)] appeared to give indication of D3d ground state symmetry, in contrast to theoretical predictions of D5d symmetry. Here, we revisit our original measurements taking account of a previous theoretical spectrum presented in the work of Manini et al., Phys. Rev. Lett. 91(19), 196402 (2003), obtained within a vibronic model parametrized on density functional theory/local-density approximation electronic structure involving all hg Jahn–Teller active modes, which couple to the 2Hu components of the ground state of the C60+ cation. By reanalyzing our measured TPES of the ground state of the C60 Buckminsterfullerene, we find a striking resemblance to the theoretical spectrum calculated in the work of Manini et al., Phys. Rev. Lett. 91(19), 196402 (2003), and we provide assignments for many of the hg modes. In order to obtain deeper insights into the temperature effects and possible anharmonicity effects, we provide complementary modeling of the photoelectron spectrum via classical molecular dynamics (MD) involving density functional based tight binding (DFTB) computations of the electronic structure for both C60 and C60+. The validity of the DFTB modeling is first checked vs the IR spectra of both species which are well established from IR spectroscopic studies. To aid the interpretation of our measured TPES and the comparisons to the ab initio spectrum we showcase the complementarity of utilizing MD calculations to predict the PES evolution at high temperatures expected in our experiment. The comparison with the theoretical spectrum presented in the work of Manini et al., Phys. Rev. Lett. 91(19), 196402 (2003), furthermore, provides further evidence for a D5d symmetric ground state of the C60+ cation in the gas phase, in complement to IR spectroscopy in frozen noble gas matrices. This not only allows us to assign the first adiabatic ionization transition and thus determine the ionization energy of C60 with greater accuracy than has been achieved at 7.598 ± 0.005 eV, but we also assign the two lowest excited states (2E1u and 2E2u) which are visible in our TPES. Finally, we discuss the energetics of additional DIBs that could be assigned to C60+ in the future. [ABSTRACT FROM AUTHOR]

This study presents a collaborative experimental and theoretical investigation into the structures and electronic properties of niobium-doped germanium clusters. Anion photoelectron spectra for Nb1–2Gen− (n = 3–7) clusters were acquired using 266 nm photon energies, enabling the determination of adiabatic detachment energies and vertical detachment energies. In conjunction with these experimental measurements, density functional theory (DFT) calculations were conducted to validate the experimentally obtained electron detachment energies and elucidate the geometric and electronic structures of each anionic cluster. The agreement between DFT calculations and experimental data establishes a solid foundation for assessing the structural evolution and electronic properties of niobium-doped germanium clusters. It is noted that both neutral and anionic clusters exhibit predominantly similar overall structural characteristics, with the exception of Nb2Ge6− and Nb2Ge6. Furthermore, this investigation emphasizes the exceptional chemical stability of the D3d symmetric chair-shaped structure in Nb2Ge6−, providing insights into its bonding characteristics. [ABSTRACT FROM AUTHOR]

Quantum-chemical subsystem and embedding methods require complex workflows that may involve multiple quantum-chemical program packages. Moreover, such workflows require the exchange of voluminous data that go beyond simple quantities, such as molecular structures and energies. Here, we describe our approach for addressing this interoperability challenge by exchanging electron densities and embedding potentials as grid-based data. We describe the approach that we have implemented to this end in a dedicated code, PyEmbed, currently part of a Python scripting framework. We discuss how it has facilitated the development of quantum-chemical subsystem and embedding methods and highlight several applications that have been enabled by PyEmbed, including wave-function theory (WFT) in density-functional theory (DFT) embedding schemes mixing non-relativistic and relativistic electronic structure methods, real-time time-dependent DFT-in-DFT approaches, the density-based many-body expansion, and workflows including real-space data analysis and visualization. Our approach demonstrates, in particular, the merits of exchanging (complex) grid-based data and, in general, the potential of modular software development in quantum chemistry, which hinges upon libraries that facilitate interoperability. [ABSTRACT FROM AUTHOR]