Basis Set Research Articles

DFT and Monte Carlo simulation studies of potential corrosion inhibition properties of some basic heterocyclic compounds

ABSTRACT This study utilises a combination of density functional theory (DFT) approach, 6-31G (d, p) basis set, and Monte Carlo (MC) simulations. The goal is to understand the relationship between the structure and properties of heterocyclic inhibitors and how they interact with Fe (110) surfaces. Furan, pyrrole, thiophene, pyridine, pyridazine, pyrimidine, indole, benzofuran, carbazole, quinolone, isoquinoline, and imidazole are the chemicals that were looked at in this study. The DFT calculations provide important insights into various properties of the inhibitors, such as molecular structure, dipole moments, electrostatic potential maps, Reduced Density Gradient (RDG), UV spectroscopy, and thermal characteristics. Carbazole has the smallest energy difference (4.789 eV) and the highest softness (0.209 eV−1), It also has the most chemical activity, Furthermore, the Monte Carlo models demonstrate that these derivatives have favourable adsorption energies on the Fe (110) surface. Among them, carbazole shows significant inhibitory potential. The study confirms that the carbazole compound with the biggest negative adsorption energy (−95.495) actually has the smallest energy difference and fits with the MC simulation for adsorption on Fe (110) when it is not charged. Carbazole has much stronger adsorption than oxazole, showing how the structure of a molecule affects how well an inhibitor works.

Electronic, thermodynamic, NLO, and spectroscopic properties of OAFLC compounds via DFT: A focus on 1F6B, 1F6Bi, 1F6BiN, and 1F6BN

This study presents a comprehensive quantum mechanical investigation of four novel liquid crystalline compounds: 1F6B, 1F6Bi, 1F6BiN, and 1F6BN. Using Density Functional Theory (DFT) with the B3LYP functional and 6-311G(d,p) basis set, we optimized the molecular geometries and analyzed their electronic and vibrational properties. Atomic Polar Tensor (APT) charge analysis showed significant charge distributions. Electrostatic potential (ESP) mapping identified potential reactive sites, while Electron Localization Function (ELF) and Localized Orbital Locator (LOL) analyses depicted the electron density distributions. Thermodynamic parameters indicated the compounds’ vibrational stability and responsiveness to thermal changes. Non-linear optical (NLO) properties were characterized by significant dipole moments and polarizabilities, suggesting potential in optoelectronic applications. Frontier molecular orbital analysis indicated stability with HOMO-LUMO gaps ranging from 4.15 to 4.64 eV. UV–Vis spectra revealed key electronic transitions, while Raman spectroscopy detailed various vibrational modes. The study offers valuable insights into the molecular behavior of these compounds, laying the groundwork for their future applications in material science and advanced technologies.

Structural analysis and theoretical studies of 2-(2-chlorobenzylidene) malononitrile (CS) “a riot controlling agent”

In solvent methanol, malononitrile and 2-chlorobenzaldehyde were stirred in the presence of base piperidine to synthesise the target chemical 2-(2-chlorobenzylidene) malononitrile (CS). The crystal structure of the target compound C10H5ClN2 has been elucidated through X-ray crystallographic analysis. It crystallises in the monoclinic crystal system within the P21 space group, containing four molecules per asymmetric unit. The unit cell dimensions are a = 3.8825(3) Å, b = 21.0202(18) Å, c = 21.3481(18) Å, and β = 95.096(4)°. The phenyl and malononitrile groups, which make up the molecular components, are each planar but slanted towards one another at a dihedral angle of 12.7° (4). Besides this, target molecule is also stabilized by hydrogen bonding, halogen bonding and π stacking interactions. Solid-state structure has been analysed through Hirshfeld surface analysis, including the evaluation of the different energy frameworks, indicating that the molecular sheets are primarily formed by hydrogen bonds and halogen bonds and the stabilization is dominated via the electrostatic energy contribution. Fingerprint plots were carried out to study the relative contribution of noncovalent interactions in the target compound CS. Further, the density functional theory calculations were employed using B3LYP hybrid functional with 6-311++G level basic set to optimize the structural coordinates in gas phase and in different solvents (chloroform, water, and methanol). The chemically active regions of the target compound were identified by the analysis of molecular electrostatic potential surface. Further, from PASS analysis and literature reports, molecular docking of title compounds was carried out with human transient receptor potential ankyrin 1 ion channel (hTRPA1) (PDB-ID 6PQO). ΔGb (docking energy) and inhibition constant was calculated to be −7.53 kcal/mol and 3.33 μM respectively, with the target ion channel.

Spectroscopic, quantum chemical and molecular docking studies of N-((anthracen-9-yl) methylene-2-methyl benzenamine): a potent probe to lipid bilayer systems

ABSTRACT A new synthesized Schiff base, N-((anthracen-9-yl) methylene-2-methyl benzenamine) (AMMB), was prepared from o-toluidine and 9-anthracenecarboxaldehyde. The structure of the compound was confirmed by single crystal X-ray diffraction and characterized through various analytical tools such as high resolution mass, 1H and 13C NMR, infrared, UV-vis and elemental composition analysis. The molecular geometries of AMMB were optimized by means of DFT/B3LYP/6311-G++ (d,p) basis set. The optimized geometric parameters of the compound were concordant with the single crystal structure. The optimized structure was utilized in comparison of experimental FT-IR with calculated frequency calculation. TD-DFT analysis was performed to compute the electronic absorption of the compound. SCF-GIAO was employed to compute 1H and 13C NMR, which was consistent with the experimental results. Additionally, NLO, NBO, MEP and global reactivity parameters were also computed using the same level of theory. Hirshfield surface analysis, fingerprint plot and energy frameworks were presented to explore the intermolecular interactions within the crystal. Molecular docking study was carried out to predict the probable binding site of AMMB into lipid bilayers such as DMPC, DPPC, POPC, DLPE, DOPC and POPE. The molecular docking analysis of AMMB showed complex formation with binding affinity energy ranging from −9.3 to −3.4 kcal/mol.

Design, synthesis, characterization, and theoretical calculations, along with in silico and in vitro antimicrobial proprieties of new isoxazole-amide conjugates

Abstract Functionalized isoxazoles provide valuable structural motifs, opening up a wide range of uses in the medicinal, pharmacological, and pharmaceutical fields. Within this scope, an efficient approach has been adopted to synthesize a novel series of functionalized isoxazole derivatives, starting from aza-aurone, providing reproducible access to the desired isoxazoles in excellent yields. All synthesized compounds were structurally elucidated through the use of various spectroscopic techniques and mass spectrometry. The derivatives generated were screened for their antimicrobial potential against the fungus Candida albicans as well as three bacterial strains. The results show that almost all of the tested isoxazole derivatives were found to be significantly potent against the fungus C. albicans. The functionalized isoxazoles were also computed using the Gaussian software package with the 6-31++G(d,p) basis set at B3LYP, HF, and M062X levels, and their chemical activities were compared. Moreover, the molecular docking studies of tested isoxazole compounds were performed against the C. albicans receptor. The results suggest that the newly synthesized compounds exhibit docking scores ranging from −10.29 to −15.08 kcal/mol, revealing a high affinity for the target enzyme (5V5Z). Lastly, drug similarity studies and ADMET (absorption, distribution, metabolism, excretion, and toxicity) properties assessments indicate that isoxazole derivatives have favorable absorption, distribution, and metabolism properties associated with a proven lack of toxicity.

Molecular docking, DFT and antiproliferative properties of 4-(3,4-dimethoxyphenyl)-3-(4-methoxyphenyl)-1-phenyl-1H-pyrazolo[3,4-b]pyridine as potent anticancer agent with CDK2 and PIM1 inhibition potency.

Due to the limited effeteness and safety concerns associated with current cancer treatments, there is a pressing need to develop novel therapeutic agents. 4-(3,4-Dimethoxyphenyl)-3-(4-methoxyphenyl)-1-phenyl-1H-pyrazolo[3,4-b]pyridine (3) was synthesized and Initially screened on 59 cancer cell lines showed promising anticancer activity, so, it was chosen for a 5-dose experiment by the NCI/USA. The GI50 values ranged from 1.04 to 8.02 μM on the entire nine panels (57 cell lines), with a GI50 of 2.70 μM for (MG-MID) panel, indicating an encouraging action. To further explore the molecular attributes of compound 3, we optimized its structure using DFT with the B3LYP/6-31 + + G(d,p) basis set. We have considered vibrational analysis, bond lengths and angles, FMOs, and MEP for the structure. Additionally, pharmacokinetic assessments were conducted using various in-silico platforms to evaluate the compound safety. A molecular modeling study created a kinase profile on 44 different kinases. This allowed us to study our compound's binding affinity to these kinases and compare it to the co-crystallized one. Our findings revealed compound 3 exhibited better binding for half of the tested kinases, suggesting its potential as a multi-kinase inhibitor. To further validate our computational results, we tested compound 3 for its inhibitory effects on CDK2 and PIM1. Compound 3 exhibited an IC50 of 0.30 µM for CDK2 inhibition, making it five times less active than Roscovitine, which has an IC50 of 0.06 µM. However, compound 3 demonstrated slightly better inhibition of PIM1 compared to Staurosporine. These findings suggest that compound 3 is a promising anticancer agent with the potential for further development into a highly active compound.

Atomistic first-principles modeling of single donor spin-qubit

Using an impurity atom in crystal silicon as a spin-1/2 qubit has been made experimentally possible recently where the impurity atom acts as a quantum dot (QD). Quantum transport in and out of such a donor QD occurs in the sequential tunneling regime where a physical quantity of importance is the charging (addition) energy, which measures the energy necessary for adding an electron into the donor QD. In this work, we present a first-principles method to quantitatively predict the addition energy of the donor QD. Using density functional theory (DFT), we determine the impurity states that serve as the basis set for subsequent exact diagonalization calculation of the many-body states and energies of the donor QD. Due to the large effective Bohr radius of the conduction electrons in Si, very large supercells containing more than 10 000 atoms must be used to obtain accurate results. For the donor QD of a phosphorus impurity in bulk Si, the combined DFT and exact diagonalization predicts the first addition energy to be 53 meV, in good agreement with the corresponding experimental value. For the donor QD of an arsenic impurity in Si, the first addition energy is predicted to be 44.2 meV. The calculated many-body wave functions provide a vivid electronic picture of the donor QD.

A static quantum embedding scheme based on coupled cluster theory.

We develop a static quantum embedding scheme that utilizes different levels of approximations to coupled cluster (CC) theory for an active fragment region and its environment. To reduce the computational cost, we solve the local fragment problem using a high-level CC method and address the environment problem with a lower-level Møller-Plesset (MP) perturbative method. This embedding approach inherits many conceptual developments from the hybrid second-order Møller-Plesset (MP2) and CC works by Nooijen [J. Chem. Phys. 111, 10815 (1999)] and Bochevarov and Sherrill [J. Chem. Phys. 122, 234110 (2005)]. We go beyond those works here by primarily targeting a specific localized fragment of a molecule and also introducing an alternative mechanism to relax the environment within this framework. We will call this approach MP-CC. We demonstrate the effectiveness of MP-CC on several potential energy curves and a set of thermochemical reaction energies, using CC with singles and doubles as the fragment solver, and MP2-like treatments of the environment. The results are substantially improved by the inclusion of orbital relaxation in the environment. Using localized bonds as the active fragment, we also report results for N=N bond breaking in azomethane and for the central C-C bond torsion in butadiene. We find that when the fragment Hilbert space size remains fixed (e.g., when determined by an intrinsic atomic orbital approach), the method achieves comparable accuracy with both a small and a large basis set. Additionally, our results indicate that increasing the fragment Hilbert space size systematically enhances the accuracy of observables, approaching the precision of the full CC solver.

A qubit-efficient variational selected configuration-interaction method

Abstract Finding the ground-state energy of molecules is an important and challenging computational problem for which quantum computing can potentially find efficient solutions. The variational quantum eigensolver (VQE) is a quantum algorithm that tackles the molecular groundstate problem and is regarded as one of the flagships of quantum computing. Yet, to date, only very small molecules were computed via VQE, due to high noise levels in current quantum devices. Here we present an alternative variational quantum scheme that requires significantly less qubits than VQE. The reduction in the qubit number allows for shallower circuits to be sufficient, rendering the method more resistant to noise. The proposed algorithm, termed variational quantum selected-configuration-interaction (VQ-SCI), is based on: (a) representing the target groundstate as a superposition of Slater determinant configurations, encoded directly upon the quantum computational basis states; and (b) selecting a-priory only the most dominant configurations. This is demonstrated through a set of groundstate calculations of the H2, LiH, BeH2, H2O, NH3 and C2H4 molecules in the sto-3g basis set, performed on IBM quantum devices. We show that the VQ-SCI reaches the full configuration interaction energy within chemical accuracy using the lowest number of qubits reported to date. Moreover, when the SCI matrix is generated ‘on the fly’, the VQ-SCI requires exponentially less memory than classical SCI methods. This offers a potential remedy to a severe memory bottleneck problem in classical SCI calculations. Finally, the proposed scheme is general and can be straightforwardly applied for finding the groundstate of any Hermitian matrix, outside the chemical context.

Post-CCSD(T) Thermochemistry of Chlorine Fluorides as a Challenging Test Case for Evaluating Density Functional Theory and Composite Ab Initio Methods.

Quantum chemistry plays a key role in exploring the chemical properties of highly reactive chlorine polyfluorides compounds (ClFn). Here, we investigate the thermochemical properties of ClFn species (n = 2-6) by means of high-level thermochemical procedures approximating the CCSDT(Q) and CCSDTQ5 energies at the complete basis set limit. We consider total atomization energies (TAEs), Cl-F bond dissociation energies (BDEs), F2 elimination energies (F2 elim.), ionization potentials (IPs), and electron affinities (EAs). The TAEs have significant contributions from post-CCSD(T) correlation effects. The higher-order triple excitations, CCSDT-CCSD(T), are negative and amount to -0.338 (ClF2), -0.727 (ClF3), -0.903 (ClF4), -1.335 (ClF5), and -1.946 (ClF6) kcal/mol. However, the contributions from quadruple (and, where available, also quintuple) excitations are much larger and positive and amount to +1.335 (ClF2), +1.387 (ClF3), +2.367 (ClF4), +2.399 (ClF5), and +3.432 (ClF6) kcal/mol. Thus, the contributions from post-CCSD(T) excitations exceed the threshold of chemical accuracy in nearly all cases. Due to their increasing hyper-valency and multireference character, the ClFn series provides an interesting and challenging test case for both density functional theory and low-level composite ab initio procedures. Here, we highlight the limitations in achieving overall chemical accuracy across all DFT and most composite ab initio procedures.

Effect of Fluorine Substitution on Magnetizabilities: Insights from Density Functional Theory and Benchmark Coupled-Cluster Calculations.

The quantitative calculation of the magnetizability tensor of fluorine-containing molecules has been a longstanding challenge for quantum chemistry. We report a benchmark study on the effect of fluorine substitution on the magnetizability (both isotropic and anisotropic) of 24 small closed-shell molecules using coupled-cluster (CCSD(T)) theory, large basis sets, and London atomic orbitals. By extrapolation to a complete basis set (CBS) result, we establish the CCSD(T)/CBS limit and take the opportunity to assess the performance of various density functional theory approximations. Correcting for zero-point vibration further allows us to directly compare theory with (gas-phase) experimental data. We revisit Flygare's hypothesis on molecular magnetizabilities, which relates the successive replacement of hydrogen by fluorine atoms to changes in the magnetizability anisotropy. Some exceptions to this hypothesis are documented, and we rationalize these observations on the basis of hybridization on carbon.

Computational Exploration of the Impact of Low‐Spin and High‐Spin Ground State on the Chelating Ability of Dimethylglyoxime Ligand on Dihalo Transition Metal: A QTAIM, EDA, and CDA Analysis

ABSTRACTWe have explored the chelation of dimethylglyoxime ligand to divalent (ndx: x = 6, 7, 8) transition metal (TM) cations in two media (gas phase and water) at the B3LYP//LANL2DZ/6–311+G (d,p) and B3LYP/def2‐TZVP level at lower multiplicity and higher multiplicity states. Majority of the 18 optimized halide (chloride and bromide) complexes prefer square planar configuration. The correlations discerned between the experimental structural data and their estimated counterparts demonstrate a good credibility for complexes at lower multiplicity state. The basis set superposition errors (BSSEs) estimated is very small which reflects the fact that the choice of different basis sets (B3LYP//LANL2DZ/6–311+G (d,p)) introduces a slight bias in the calculation of energies. The ADMP (atom‐centered density matrix propagation) simulations in water on chloride complexes indicate the irreversible nature of these M—N dissociation in trajectory simulation process. This fact explains our exclusive focus on the examination of the [glyoxime ligand]…[MX2] interactions. In addition, the solvation of (3d and 4d) transition metal chloride complexes causes a sensitive augmentation of the metal ion affinity (MIA) with an average of 0.29 and 0.24 kcal/mol. In both multiplicity states, the topological parameters have illustrated that the M—N and M—X bonds are typical metal–ligand in both media. The average ΔEorb/ΔEsteric ratio equal to 0.45 and 0.11 in gas phase and water, respectively, reveals the predominance of the contributions from non‐covalent bonding interactions (NCI) compared to those of covalent bonding. But, the maximal value equal to 6.760 is obtained for bromide rhodium complex in water. NBO analysis in both media highlights the fact that a more pronounced ionic character is observed for the majority of the chloride complexes at both spin multiplicity states because of their higher retained charges on the metal atom. For [dimethylglyoxime]…[MX2] interaction (X = Cl and Br), the charge decomposition analysis demonstrates that the lowest value of the d/b ratio is found for the chloride platinum complex at lower multiplicity state in water. This is a proof of its strong relativistic effects.

Molecular docking and spectroscopic exploration (FT‐IR, FT‐Raman) with HOMO–LUMO, ADMET, and 5,7‐dihydroxy‐2‐(4‐hydroxyphenyl) chroman‐4‐one against lung cancer: A potential uptake B‐RAF inhibitor

AbstractNaringenin has been proven to inhibit cell proliferation, which has anticancer properties. The role of Naringenin molecule in lung cancer and its processes are yet unknown. Naringenin is chemically called as 5,7‐dihydroxy‐2‐(4‐hydroxyphenyl) chroman‐4‐one, which was optimized geometrical parameters analysis such as bond length, bond angle and torsion angle were analyzed from Naringenin by utilizing Gaussian 09 W program. Characterizations of Naringenin analyzed by B3LYP density functional theory with basis set 6–311G (d,p). The energy value of Naringenin molecules are analyzed with (ground state) HOMO value −4.616, LUMO value for −0.169 (first excited state) energies and also predicted by energy cap value is 4.446. ADMET and drug likeness of the title compound was predicted, it's qualified with the Lipinski's rule of five. Naringenin molecular structural changes, distribution and also its reactive site investigated with MEP (Molecular Electrostatic Potential) analyses spans from −0.118eO to 0.118eO. The density of state (DOS) used to know molecular orbital contribution for selected compound. It was determined that Naringenin molecule in the active site of B‐RAF inhibitor (PDB code: 4MNF) has a binding energy of −8.9 based on a docking analysis and conformational changes. From this study, this drug will be expected to undergo essential future chemotherapy agent for lung cancer patient.

Structural characterisation of a cyanomesogen by X-ray diffraction and computational methods

The work reports the structural analysis of a cyanomesogen namely, p-cyanobenzylidene p-nonyloxyaniline belonging to the Schiff base class of compounds. It exhibits a partial bilayer/inter-digitated smectic Ad phase and a nematic phase. Polarising Optical Microscopy (POM) and Differential Scanning Calorimetry (DSC) studies have been carried out to study phase sequence and transition temperatures. Micro-structural parameters, like molecular length, spacing, crystallite size, activation energy and micro-strain, have been estimated and pair correlation function (PCF) studies have been done for exhibited mesophases using the X-ray diffraction technique. The study investigated the surface, electronic and optical properties of a compound by computational studies like Density Functional Theory (DFT) and Hirshfeld surface analysis (HSA) methods using the DFT B3LYP/6-311++G(d,p) level basis set. The results reveal that the polarity associated with the polar groups and alkoxy chain lengths, contributes majorly to the mesomorphic characteristics and stability. The Highest Occupied Molecular Orbitals (HOMO) and Lowest Unoccupied Molecular Orbitals (LUMO) analyses of the DFT optimised molecule give electronic and geometrical parameters, atomic charge distribution, polarisability factors which are used to interpret the characteristics of the mesophases and their suitability for electro-optical and non-linear optical applications of the compound in terms of the predicted data.

Solving the Puzzle of the Carbonic Acid Vibrational Spectrum - an Anharmonic Story.

Against the general belief that carbonic acid is too unstable for synthesis, it was possible to synthesize the solid[1,2] as well as gas-phase carbonic acid.[3] It was suggested that solid carbonic acid might exist in Earth's upper troposphere and in the harsh environments of other solar bodies,[4] where it undergoes a cycle of synthesis, decomposition, and dimerization.[5] To provide spectroscopic data for probing the existence of extraterrestrial carbonic acid,[2,6] matrix-isolation infrared (MI-IR) spectroscopy has shown to be essential.[3,4,6-8] However, early assignments within the harmonic approximation using scaling factors impeded a full interpretation of the rather complex MI-IR spectrum of H2CO3. Recently, carbonic acid was detected in the Galactic center molecular cloud,[9] triggering new interest in the anharmonic spectrum.[10] In this regard, we substantially reassign our argon MI-IR spectra based on accurate anharmonic calculations. We calculate a four-mode potential energy surface (PES) at the explicitly correlated coupled-cluster theory using up to triple-zeta basis sets, i. e., CCSD(T)-F12/cc-pVTZ-F12. On this PES, we perform vibrational self-consistent field and configuration interaction (VSCF/VCI) calculations to obtain accurate vibrational transition frequencies and resonance analysis of the fundamentals, first overtones, and combination bands. In total, 12 new bands can be assigned, extending the spectral data for carbonic acid and thus simplifying detection in more complex environments. Furthermore, we clarify disputed assignments between the cc- and ct-conformer.

Chemical space-informed machine learning models for rapid predictions of x-ray photoelectron spectra of organic molecules

Abstract We present machine learning models based on kernel-ridge regression for predicting X-ray photoelectron spec- tra of organic molecules originating from the ionization energies of 1s electrons in carbon (C), nitrogen (N), oxygen (O), and fluorine (F) atoms. We constructed training dataset through high-throughput calculations of core-electron binding energies (CEBEs) for 12,880 small organic molecules in the bigQM7ω dataset, employing the ∆-SCF formalism coupled with meta-GGA-DFT and a variationally converged basis set. The models are cost-effective, as they require the atomic coordinates of a molecule generated using universal force fields while estimating the target-level CEBEs corresponding to DFT-level equilibrium geometry. We explore transfer learning by utilizing the atomic environment feature vectors learned using a graph neural network framework in kernel-ridge regression. Additionally, we enhance accuracy using the ∆-machine learning framework by leveraging inexpensive baseline spectra derived from Kohn–Sham eigenvalues. Upon application to 208 com- binatorially substituted uracil molecules, larger than those in the training set, our analyses reveal that while the models may not yield quantitatively accurate predictions of CEBEs on a molecule-by-molecule basis, they do exhibit a strong linear correlation, which proves valuable for virtual high-throughput screening purposes. We present the dataset and models as the Python module, cebeconf, to facilitate further explorations.

An Approach to Realize Generalized Optimal Motion Primitives Using Physics Informed Neural Networks

Abstract Autonomous manipulation is a challenging problem in field robotics due to uncertainty in object properties, constraints, and coupling phenomenon with robot control systems. Humans learn motion primitives over time to effectively interact with the environment. We postulate that autonomous manipulation can be enabled by basic sets of motion primitives as well, but do not necessitate mimicking human motion primitives. This work presents an approach to generalized optimal motion primitives using physics-informed neural networks. Our simulated and experimental results demonstrate that optimality is notionally maintained where the mean maximum observed final position percent error was 0.564% and the average mean error for all the trajectories was 1.53%. These results indicate that notional generalization is attained using a physics-informed neural network approach that enables near optimal real-time adaptation of primitive motion profiles.

Finite-order method to calculate approximate density matrices in the Fock-space multireference coupled cluster theory

An efficient approach to calculate approximate pure-state and transition reduced density matrices in the framework of the multireference relativistic Fock-space coupled cluster (FS CC) theory is proposed. The method is based on the effective operator formalism and consists of the direct substitution of the FS CC Ansatz for a wave operator into the effective operator expression with the subsequent truncation of expansion at the terms quadratic in cluster amplitudes. The final density matrix is defined by active-space density matrices of different ranks ‘dressed’ with contributions from cluster operators. The method gives a connected expression for pure-state density matrices, provided that the intermediate normalisation condition is fulfilled. Moreover, under some additional assumptions, the connectivity can also be ensured for calculated transition property matrix elements and natural transition spinors. The developed technique allows for fast and accurate calculations of one-particle reduced density matrices for a wide range of electronic states. A pilot application of the new technique to construct averaged atomic natural orbital (ANO) basis sets for fully relativistic electronic structure calculations is presented.

Quantum chemical calculation dataset for representative protein folds by the fragment molecular orbital method

The function of a biomacromolecule is not only determined by its three-dimensional structure but also by its electronic state. Quantum chemical calculations are promising non-empirical methods available for determining the electronic state of a given structure. In this study, we used the fragment molecular orbital (FMO) method, which applies to biopolymers such as proteins, to provide physicochemical property values on representative structures in the SCOP2 database of protein families, a subset of the Protein Data Bank. Our dataset was constructed by over 5,000 protein structures, including over 200 million inter-fragment interaction energies (IFIEs) and their energy components obtained by pair interaction energy decomposition analysis (PIEDA) using FMO-MP2/6-31 G*. Moreover, three basis sets, 6-31 G*, 6-31 G**, and cc-pVDZ, were used for the FMO calculations of each structure, making it possible to compare the energies obtained with different basis functions for the same fragment pair. The total data size is approximately 6.7 GB. Our dataset will be useful for functional analyses and machine learning based on the physicochemical property values of proteins.

Healthcare Costs of Hospitalizations Due to Aspergillosis and 25-Year Trends in Spain, 1997–2021

In the last 40 years, a significant increase in the incidence of lung infections by Aspergillus has been reported. The scarcity of studies that describe the costs of aspergillosis indicates that the economic impact of aspergillosis in the hospital environment is greater than that of other fungal infections. The objective of the study was to evaluate the direct healthcare costs associated with aspergillosis in the Spanish National Health System from 1997 to 2021. A retrospective nationwide longitudinal descriptive study was designed to review hospital records from the Minimum Basic Data Set of patients admitted to hospitals of the National Health System from 1997 to 2021, with a diagnosis of aspergillosis. A total of 44,586 patients were admitted for aspergillosis in the Spanish National Health System. There was a progressive increase in the average annual cost from 1997 to 2012, which reached a maximum peak, EUR 1,395,154.21 (±2,155,192.87). It decreased between 2014 and 2019, but increased again in 2020 and 2021, EUR 28,675.79 (±30,384.12). The Pearson correlation coefficient revealed a weak negative correlation between age and hospital costs and a moderate positive correlation between average length of stay and hospital costs. Our data show that the economic impact of hospitalizations for aspergillosis is significant and increasing at a rate proportionally higher than that of other prevalent diseases. Costs related to Aspergillus infection are associated mainly with respiratory diseases. The results of this economic evaluation may be useful for health authorities to develop a future economic strategy for managing this fungal infection.