Energy Landscapes of Carbon Clusters from Tight-Binding Quantum Potentials

Infrared bands of neutral gas-phase carbon clusters in a broad spectral range

The identification of species in the interstellar medium requires precise and molecule-specific spectroscopic information in the laboratory framework, in broad spectral ranges and under conditions relevant to interstellar environments. In this work, we measure the gas-phase infrared spectra of neutral carbon clusters, CN (N = 6-11), in a molecular beam. The CN distribution is formed by photofragmentation of C60 molecules, concurrently showing a top-down formation mechanism. A broad spectral range in the infrared between 500-3200 cm-1 (20-3.125 μm) is investigated. We observe strong bands between 5 and 6 μm, in conjunction with novel features in the 3 μm region. Density functional theory calculations reveal that these short wavelength modes correspond to combination bands with significant infrared intensity. Moreover, we identify the N ≤ 10 clusters as linear, while C11 adopts a ring configuration, placing the linear-to-ring transition at N = 11 under our molecular beam conditions. The linearity of C10 is discussed based on the formation pathway from larger clusters in energetic conditions. Given the vast and very precise infrared information already been released from the James Webb Space Telescope mission, this infrared spectroscopic data set in conjunction with information on formation mechanisms is of major relevance for identifying neutral carbon clusters in astronomical environments.

Multilevel Framework for Analysis of Protein Folding Involving Disulfide Bond Formation

In this study, a three-layered multicenter ONIOM approach is implemented to characterize the naive folding pathway of bovine pancreatic trypsin inhibitor (BPTI). Each layer represents a distinct level of theory, where the initial layer, encompassing the entire protein, is modeled by a general all-atom force-field GFN-FF. An intermediate electronic structure layer consisting of three multicenter fragments is introduced with the state-of-the-art semiempirical tight-binding method GFN2-xTB. Higher accuracy, specifically addressing the breaking and formation of the three disulfide bonds, is achieved at the innermost layer using the composite DFT method r2SCAN-3c. Our analysis sheds light on the structural stability of BPTI, particularly the significance of interlinking disulfide bonds. The accuracy and efficiency of the multicenter QM/SQM/MM approach are benchmarked using the oxidative formation of cystine. For the folding pathway of BPTI, relative stabilities are investigated through the calculation of free energy contributions for selected intermediates, focusing on the impact of the disulfide bond. Our results highlight the intricate trade-off between accuracy and computational cost, demonstrating that the multicenter ONIOM approach provides a well-balanced and comprehensive solution to describe electronic structure effects in biomolecular systems. We conclude that multiscale energy landscape exploration provides a robust methodology for the study of intriguing biological targets.

Fast Screening of Minimum Energy Crossing Points with Semiempirical Tight-Binding Methods

The investigation of photochemical processes is a highly active field in computational chemistry. One research direction is the automated exploration and identification of minimum energy conical intersection (MECI) geometries. However, due to the immense technical effort required to calculate nonadiabatic potential energy landscapes, the routine application of such computational protocols is severely limited. In this study, we will discuss the prospect of combining adiabatic potential energy surfaces from semiempirical quantum mechanical calculations with specialized confinement potential and metadynamics simulations to identify S₀/T₁ minimum energy crossing point (MECP) geometries. It is shown that MECPs calculated at the GFN2-xTB level can provide suitable approximations to high-level S₀/S₁ab initio conical intersection geometries at a fraction of the computational cost. Reference MECIs of benzene are studied to illustrate the basic concept. An example application of the presented protocol is demonstrated for a set of photoswitch molecules.

Effects of orbital angles on the modeling of conjugated systems with curvature

Models with angle corrections give well predictions of both neutral and charged fullerenes. The integrals of nonparallel orbitals explain why angle features of designed and deep-learning models are necessary to describe conjugated systems.