Basis Set Research Articles

Synthesis, characterization, quantum mechanical calculations and biomedical docking studies on curcumin analogs: 2, 6-(Difurfurylidene) cyclohexanone and 2, 6 – Bis (2,6-Dichloro Benzylidene) Cyclohexanone

The initiation of colorectal cancer is controlled by various factors, including random occurrences and genetic alterations affecting oncogenes and tumor suppressor genes.Curcumin, a significant compound extracted from turmeric, has attracted interest for its robust anticancer properties, particularly regarding its analogs, 2, 6-bisdifurfurylidene cyclohexanone (DFC) and 2, 6-bis (2, 6-dichlorobenzylidene) cyclohexanone (DCC), which were synthesized and assessed for their anticancer efficacy. A combination of spectroscopic techniques and molecular docking methods was utilized to comprehensively evaluate the interaction behaviors of DFC and DCC. The application of density functional theory (DFT) using the B3LYP/6-311G (d, p) basis set facilitated the prediction of spectroscopic properties. The molecular docking investigations conducted using the Glide docking program from Schrodinger Maestro elucidated the interactions of these drugs at the molecular level. In vitro investigations were performed to evaluate the cytotoxic efficacy of the synthesized curcumin analogs. The determined IC50 values revealed that DFC displayed an IC50 of approximately 82 μM, and DCC exhibited a significantly lower IC50 of around 10 μM. This notable disparity highlights the potential of DFC and DCC as a more efficacious cytotoxic agent and further research be conducted on the produced chemicals in the future.

A theoretical study on symmetrical non-fullerene electron acceptors molecules on BDTPT based derivatives to enhance photovoltaic properties of organic solar cells

An innovative and promising approach to developing sustainable energy solutions and promoting an eco-friendly society is the use of organic solar cells. The key component for a solution-processed bulk-heterojunction organic solar cell is the photoactive layer’s embedded donor and acceptor components. This research presents seven modified molecules comprising the A–D–A type structural configuration, entitled A1–A7. All these designed moieties exhibit marvelous outcomes in optoelectronic features, including λmax and band gap, owing to non-fullerene acceptors in the terminal regions. All these compounds are computationally assessed by employing B3LYP at 6-31G (d,p) basis set using chloroform solvent. Compared to the reference molecule, the designed molecules (A1, A2, A4, A5, A6, A7) have reflected breakthrough results. The prerequisite for directing the practical application of designed acceptors is the efficient charge transfer, evidenced by coupling the J61 donor complex with the designed A5 acceptor.

Role of Exact Exchange and Empirical Dispersion in Density Functional Theory-Based Three-Body Noncovalent Interactions.

Total and three-body interaction energies are calculated for a benchmark set of three-body systems using a range of different types of density functional theory (DFT) methods, with the results compared to CCSD(T)/CBS results from the benchmark reference [Phys. Chem. Chem. Phys. 2023, 25, 28621-28637]. Inclusion of Hartree-Fock exchange, via either a global or range-separated hybrid approach or inclusion of empirical dispersion corrections, increases accuracy for total and three-body interactions. Basis set convergence testing shows that the aug-cc-pVTZ basis set is well converged with little to no change seen when using quadruple-ζ basis sets. The accuracy of the DFT methods is similar when calculating interaction energies for both global and local minimum structures. Overall, the CAM-B3LYP-D3BJ, B97D3, and ωB97XD functionals are recommended for calculating three-body interactions.

Improving the Accuracy in the Prediction of Transition-Metal Spin-State Energetics Using a Robust Variation-Based Approach: Density Functional Theory, CASPT2 and MC-PDFT Applied to the Case Study of Tris-Diimine Fe(II) Complexes.

Designing ligands for transition metal complexes with a specified low-spin, high-spin or spin-crossover behavior is challenging. A major advance was recently made by Phan et al. [J. Am. Chem. Soc. 2017, 139, 6437-6447] who showed that the spin state of a homoleptic tris-diimine Fe(II) complex can be predicted from the N-N distance in the free diimine. They could thus predict the change in magnetic behavior on passing from the complexes of 2,2'-bipyridine (bpy), 2,2'-biimidazole (bim) and 2,2'-bis-2-imidazoline (bimz) ligands to those obtained with the modified analogs 4,5-diazafluoren-9-one (dafo), 1,1'-(α,α'-o-xylyl)-2,2'-bisimidazole (xbim) and 2,3,5,6,8,9-hexahydrodiimidazo[1,2-a:2', 1'-c]pyrazine (etbimz), respectively. Theoretically, the challenge lies in the accurate determination of the HS-LS zero-point energy difference ΔEHL°. The issue can be circumvented by using a variation-based approach, wherein ΔEHL° is not directly evaluated but obtained from the estimate of its variation Δ(ΔEHL°) in series of related systems, which include one whose ΔEHL° is accurately known [Phys. Chem. Chem. Phys. 2013, 15, 3752-3763; J. Phys. Chem. A 2022, 126, 6221-6235]. In this study, density functional theory (DFT), second-order multireference perturbation theory in its CASPT2 formulation, multiconfigurational pair DFT (MC-PDFT) and its hybrid formulation (HMC-PDFT) have been applied to the determination of Δ(ΔEHL°) in the pairs of complexes , , , and , . In DFT, we used several semilocal functionals and their global hybrids, as well as their D2, D3, D3BJ and D4 dispersion-corrected forms; and in MC-PDFT, different translated and fully translated functionals. The results are consistent with one another and in very good agreement with experiments. They show small to vanishing dependence on key details of the methods used: namely, the exact-exchange contribution to global hybrids; the ionization potential-electron affinity shift and basis sets used in the CASPT2 calculations; and the active spaces employed for the CASSCF wave functions used in the MC-PDFT and HMC-PDFT calculations. Insights into the change in the spin-state energetics accompanying the ligand exchanges were gained through a complexation energy analysis. Using the accurate CCSD(T) estimate of the HS-LS adiabatic energy difference in [J. Chem. Theory Comput. 2012, 8, 4216-4231], the Δ(ΔEHL°)-approach has been applied to the determination of ΔEHL° in the diimine complexes. The CASPT2 and DFT-D2 methods only give results in agreement with experiments. This suggests for the other methods a limitation in their treatment of dispersion which prevents them from accurately describing the spin-state energetics change accompanying the passing from with the tetragonal arrangement of its nitrile ligands to the tris-diimine complexes with the trigonal packing of their bulkier ligands.

Atmospheric Degradation Mechanism of Isoamyl Acetate Initiated by OH Radicals and Cl Atoms Revealed by Quantum Chemical Calculations and Kinetic Modeling.

Isoamyl acetate is one of the volatile organic compound class molecules relevant to agricultural and industrial applications. With the growing interest in isoamyl acetate applications in industry, the atmospheric fate of isoamyl acetate must be considered. Reaction mechanisms, potential energy profiles, and rate constants of isoamyl acetate reaction with atmospheric relevant oxidant OH radicals and Cl atoms have been obtained from the quantum chemical calculations and kinetic modeling. The geometry optimizations were conducted using M06-2X/6-311++G(3df,3pd) followed by single point-energy calculations at the DLPNO-CCSD(T) method with an extrapolated complete basis set. The rate constants were calculated by solving the master equation. A hydrogen-abstraction reaction dominates the first step of isoamyl acetate degradation, while the addition-substitution reaction plays a small role in the degradation products. The kinetic study was conducted to evaluate the rate constants within a temperature range of 200-400 K. The total rate constants for the isoamyl acetate degradation reactions initiated by the OH radical and Cl atom were determined to be 6.96 × 10-12 and 1.27 × 10-10 cm3 molecule-1 s-1, respectively, under standard temperature and pressure conditions. The product degradation mechanism, ozone formation potential, and atmospheric impacts were discussed.

Synthesis, Structural, DFT studies, Biological Effectiveness, and Molecular Docking Aspects of Multitarget Novel N-(4-Phenylthiazol-2-yl)hydrazine carboxamide Hybrids as COX inhibitors, Anti-TB, and Anti-oxidant activity

We herein report, a series of N-(4-phenylthiazol-2-yl)hydrazinecarboxamide based multitarget compounds 3(a-i) were designed and synthesized. Ten compounds 2, and 3(a-i) were characterized by FTIR, NMR (1H and 13C), and mass analysis and subsequently, tested for their potential COX inhibitors, anti-TB, and anti-oxidant activities. The physicochemical and ADME studies for newly synthesized compounds were disclosed. The DFT calculations were taken for the selected molecules using CAM-B3LYP hybrid functional with a 6–31+g(d) all-electron basis set using the Gaussian 09 package. The compounds 3d, 3e, and 3f recognized outstanding COX-II inhibitions with IC50 values of 0.62, 0.78, and 0.72 μM compared to standard drugs. The compounds 3d, 3e, and 3f showed outstanding anti-TB activity with a MIC value of 0.78μg/mL. The compound 3d, and 3f attested outstanding antioxidant activity at 10 μg/mL with a rate of inhibition of 68.12%, and 68.85% respectively. Finally, the molecular docking studies carried out with cyclooxygenase-2 (PDB ID: 6COX), M. tuberculosis enoyl reductase (INHA) (PDB ID: 4TZK), and cytochrome c peroxidase (PDB ID: 2×08), for all the newly synthesized derivatives. Finally, selected compounds were taken for their molecular dynamic studies.

Jump discontinuities of finite-basis-set exchange–correlation potentials at atomic nuclei

The kinetic energy density of electrons and the gradient of the electron density have pronounced jump discontinuities at the positions of the atomic nuclei in molecules. Certain exact relations then imply that molecular Kohn–Sham exchange–correlation potentials may also be discontinuous at atomic nuclei. Here, we confirm that exchange–correlation potentials derived from Hartree–Fock and correlated wavefunctions within Slater-type basis sets do exhibit such discontinuities. Despite their persistence even in large basis sets, these discontinuities are almost certainly artifacts of basis set finiteness and are expected to disappear in the basis-set limit. The findings imply that imposing electron–nucleus cusp conditions in spherically averaged form may not always be appropriate.

Origin of the vibrational structure of the first absorption band of cis/trans isomeric 1,6-diphenylhexatrienes by (TD)DFT calculations

We present a detailed spectroscopic analysis of the first absorption band of the six possible conformers of 1,6-diphenyl-1,3,5-hexatriene, obtained by changing the configuration of trans double bonds to cis. To this end, we computed the absorption spectra using FCclasses 3.0 code. First, we assessed the performance of PBE0 and CAM-B3LYP density functional theory functionals with different basis sets to reproduce the experimental spectra. Additionally, we considered different models to compute the spectra. PBE0/def2tZVP with an adiabatic hessian model with internal coordinates yields results in very good agreement with experimental data. Subsequently, we analyzed the different contributions of vibronic transition to the spectral structure, correlating ground state conformation with spectral shape, and studied the effect of temperature on the absorption first band.

ExoMol line lists – LIX. High-temperature line list for N2O

Abstract New hot line lists for five isotopologues of N2O are presented, for the parent 14N216O and 4 singly substituted species 14N217O, 14N218O, 14N15N16O and 15N14N16O. The line lists have been computed with the variational program TROVE using a new empirical potential energy surface (PES) and an accurate ab initio dipole moment surface of N2O Ames-1. The PES was obtained by fitting to experimentally derived energies of N2O compiled using the well established MARVEL procedure. Here we also introduce an ‘artificial symmetry group’ Cns(AEM) for an efficient construction of rotation-vibrational basis set of a linear non-symmetric triatomic molecule of the XYZ type. The line lists cover the rotational excitations up to J = 160 and the wavenumber range up to 20000 cm−1. MARVEL energies are also used to improve predicted line positions resulting in excellent agreement with the available experimental spectra, as demonstrated. An extensive comparison with existing line lists for N2O HITRAN, HITEMP, NOSL-296, NOSD-1000 and Ames-296K is provided. The line lists are freely accessible from www.exomol.com.

Synthesis, spectral characterization, anti-cancer activity evaluation, and molecular docking of substituted 1,2-bis(3-butyl-2,6-diphenylpiperidin-4-ylidiene)hydrazone derivatives

A series of 1,2-bis(3‑butyl‑2,6-diphenylpiperidin-4-ylidiene)hydrazone derivatives of the substituted compounds 1–7, were confirmed by the characterization of IR & 1H-, NMR spectrum, together with X-ray crystallography and Mass spectroscopy also incredibly revealed. Docking analysis of the synthesized compounds employed under the protein from PDB code as 5cqh1, 4i40, 5fan, 7nud, 6Ya1 offers more inhibition values from 2D and 3D structures of the compounds, followed by Anti-cancer activity of the carcinoma cell line displayed that, highly potent to be active in nature. The in-vitro performance of -Cl substituted compound 6 (high anticancer activity) was in good correlation with its protein binding behaviour (strong binding). The increasing order of anticancer effect was 5>7>1>2>4>3>6 and the increasing order of binding ability on proteins were 7>3>4>1>2>5>6. The computational investigations were aligned with experimental techniques using Density Functional Theory (DFT) with a 6–31 G (d, p) basis set. Calculated Highest Occupied Molecular Orbital (HOMO) and Lowest Unoccupied Molecular Orbital (LUMO) energies revealed Mulliken charges and charge transfer phenomena occurring within the molecule. Analysis of the Molecular Electrostatic Potential (MEP) unveiled reactive sites, with high electron density regions depicted in red and low electron density regions in blue.

Hyperhardness and hypersoftness of atoms and their ions

ContextThe theory of reactivity based on cDFT has been supplemented with the new method of calculating the atomic and local indices. With the use of previously derived relationship of the electron density gradient to the softness kernel and to the linear response function, we deliver theoretical analysis to obtain significant reactivity indices—the electron density derivatives: local softness and local hypersoftness together with the global hyperhardness index and the derivative of the global softness with respect to the number of electrons. The local derivatives have been applied in the calculation of responses of atoms to perturbation by an external potential by the alchemical approach. The vital role of the local softness has been confirmed; the potential role of the hypersoftness has been indicated.MethodOur original theoretical scheme has been numerically illustrated with the results obtained with electron density calculations with B3LYP method implemented in Gaussian 16 package. The aug-cc-pvqz basis set has been routinely applied, except for the Ca atom (cc-pvqz). Using the pVTZ basis set recommended by Sadlej was necessary for the potassium atom.Graphical

Density functional theory (DFT) study of water autoionization in solvated clusters.

We have implemented a cluster-continuum method using density functional theory to model water clusters and various charged species derived from water. The two aims of this study are to determine the minimal basis required for proper modeling of water autoionization and to determine the minimum number of explicit water molecules required to properly model the energetics of solvation. The thermodynamics of water autoionization converge following a modified power law to deliver chemically accurate values of the Gibbs energy change for autoionization with tractably small clusters. Convergence is slower and not exponential as assumed in previous work. We identify the n = 21 set of clusters as the first effectively bulk water like clusters that can capture the energetic influence of both the first and second solvation shells. In this cluster, a water molecule is encapsulated in the center of a closed shell of other water molecules that hydrogen bond to form five-membered rings. The total energy change for clusters with n ≥ 21 calculated using both the RPBE-D3 and ωB97X-D exchange-correlation functionals with the 6-311+G** basis set is shown to deliver good approximations to the free energy change at 298K. This is true even though neither functional models the individual enthalpy or entropy contributions particularly well.

Anticancer activity, DFT calculation, ADMET study and molecular docking of synthetic Chalcone (2E)-1-(4-Aminophenyl)-3-(3,4,5-trimethoxyphenyl)-2- propen-1-one in breast cancer cell line (MCF-7)

The (2E)-1-(4-Aminophenyl)-3-(3,4,5-trimethoxyphenyl)-2-propen-1-one was prepared by the reaction between 4-aminoacetophenone and 3,4,5, trimethoxy benzaldehyde in ethanolic NaOH. Then in vitro anticancer activity was evaluated against the MCF-7 breast cell line. The molecular properties like bond length, bond angle, molecular electrostatic potential surface (MEPs), Mullikan charge, HOMO-LUMO, and NBO analysis data obtained by DFT/B3LYP method using 6- 31+G(d,p) basis set to get electronic states and molecular parameters of the title molecule. Because of their advantageous pharmacokinetic properties, the synthesized compound was shown to have moderate anticancer activity and could advance to a higher level of the drug development process. Molecular docking studies were utilized to foresee the binding interactions of the title compound with the EGFR tyrosine kinase domain (PDB ID: 1M17). In silico ADME study is also performed to predict the pharmacokinetic profile which expressed good oral drug-like behavior.

Millettia aboensis leaves extract as eco-friendly corrosion inhibitor for mild steel in acidizing solution: From experimental to molecular level prediction

The development of environmentally friendly corrosion inhibitors is becoming more popular since conventional inhibitors are poisonous and non-biodegradable. The anti-corrosive effectiveness of Millettia aboensis leaves extract (MALE) has been assessed in this study using several methods, such as electrochemical measurements, X-ray Photoelectron spectroscopy and scanning electron microscopy on mild steel in acidic solution. Gas chromatography-mass spectrometry (GC-MS) analysis was combined with qualitative and quantitative phytochemical analysis to identify the phytochemicals linked to the activity of the Millettia aboensis extracts and identified key phytoconstituents such as Hexanedioic acid, Phenol derivatives, Octadecanoic acid, and Linolenic acid, which play significant roles in the extract's inhibitory performance. The results showed that the inhibition efficiency (IE%) improved to 88.6 % with the addition of inhibitor concentration from 0.1 g/L to 3.0 g/L and that the corrosion rate drastically reduced from 56.91 mpy to 16.09 mpy. Furthermore, at a greater concentration (3.0 g/L), the Rct values rose from 61.42 Ω cm2 to 176.3 Ω cm2 thus, indicating that the inhibitor molecules were forming a protective film over the metallic surface. Scanning Electron Microscopy (SEM) provided visual evidence of surface morphology, revealing a smoother surface in the presence of the inhibitor, indicative of the protective film formed by the adsorption of organic molecules onto the steel surface. Theoretical calculations using the ωB97XD functional and def2svp basis set further supported the experimental findings, showing that the protonated form of C21H36O4 exhibited the highest interaction energy, correlating with its superior inhibition efficiency on the metal surface.

A theoretical and kinetic study of key reactions between ammonia and fuel molecules, part III: H-atom abstraction from esters by ṄH2 radicals

Hydrogen atom abstraction reactions by ṄH2 radicals play a crucial role in determining the reactivity of ammonia/fuel binary blends. Esters are a typical component of environmentally friendly and economically promising biofuels. The feasibility of the ammonia/biofuel dual-fuel approach has been proven in practical engines. [Energy and Fuels 22 (2008) 2963] and [Int. J. Energy Res. 2023 (2023) 9920670]. ṄH2 radicals play a critical role in the combustion and pyrolysis chemistry of ammonia and N-containing-rich fuels. In ammonia/biofuels hybrid combustion, ṄH2 radicals can react with biofuel molecules in a reaction class that is particularly important especially when sufficient ammonia is blended in order to eliminate NOx emissions. To help unravel the chemistry of ammonia/biofuel blends, a systematic theoretical kinetic study of H-atom abstraction from eleven alky esters of CnH2n+1COOCH3 (n = 1–4), CH3COOCmH2m+1 (m = 1–4), and C2H5COOC2H5, by ṄH2 radicals is performed in this work. The geometry optimization, frequency, and zero-point energy calculations for all related species, as well as the hindrance potential energy surface for low frequency torsional modes in the reactants and transition states, were performed at the M06–2X/6–311++G(d,p) level of theory. Intrinsic reaction coordinate calculations were performed to validate the connections between the transition states and expected minima energy species. The energies of all of the species involved were calculated at the QCISD(T)/cc-pVXZ (X = D, T, Q) and MP2/cc-pVYZ (Y = T, Q) levels of theory and then extrapolated to the complete basis set. Rate constants of 39 reactions were calculated using the Master Equation System Solver (MESS) program in the temperature range of 500 – 2000 K. These rate constants for different H-atom abstraction sites are provided and can be extrapolated to larger esters. The kinetic effects from the functional group are also illustrated by performing detailed comparisons with the previous studies of ṄH2 radical reactions with alkanes, alcohols and ethers.

Experimental and computational insight in molecular interactions of imidazolium based deep eutectic solvent with selected alcohols (C1 and C2)

The experimental thermophysical properties (densities (ρ), speed of sound (u) and refractive (nD) indices) of binary mixtures containing 1:3 mol ratio of [1-butyl-3-dimethylimidazolium chloride to ethylene glycol] [BMIM]Cl: EG], deep eutectic solvent (DES), with selected alcohols (methanol and ethanol) were measured. The measurements of the binary mixtures and the pure solvents were conducted at temperatures ranging from (293.15 to 313.15) K in the interval of 5 K and at ambient pressure using an Anton Paar densitometer. Thermodynamic properties such as excess molar volumes (VmE), intermolecular free length (Lf), isentropic compressibilities (ks), deviation in isentropic compressibilities (Δks) and deviation in refractive indices (ΔnD) were computed from the measured data of physical properties. These excess thermodynamic properties are utilized to investigate the intermolecular interactions occurring between imidazolium DES and methanol or ethanol. Geometry optimization calculations were undertaken using Density Functional Theory (DFT) to model the structure of DES with a 1:1 mol ratio employing the B3LYP-D3 functional and a basis set of 6–31G++**. The DFT results revealed that the interaction energies within the DES in the presence of solvent favoured the formation of the complex. Furthermore, interrogation of the optimized structures at molecular level showed that the interactions were mediated by the electrostatic ionic bonding with the co-solvent playing a significant role in solubility. The excess molar volumes of binary mixtures were correlated using the Lorentz-Lorenz equation, and the densities and refractive indices were predicted using the Lorentz-Lorenz equation. The interaction energies and electrostatic interactions between the different solvents were predicted using quantum mechanics. The experimental data was found to be in good agreement with the theoretical data for the measured properties as well as the computational predictions of the stability of the systems.

Exploiting Non-Abelian Point-Group Symmetry to Estimate the Exact Ground-State Correlation Energy of Benzene in a Polarized Split-Valence Triple-Zeta Basis Set.

Local electronic-structure methods in quantum chemistry operate on the ability to compress electron correlations more efficiently in a basis of spatially localized molecular orbitals than in a parent set of canonical orbitals. However, many typical choices of localized orbitals tend to be related by selected, near-exact symmetry operations whenever a molecule belongs to a point group, a feature which remains largely unexploited in most local correlation methods. The present Letter demonstrates how to leverage a recent unitary protocol for enforcing symmetry properties among localized orbitals to yield a high-accuracy estimate of the exact ground-state correlation energy of benzene (D6h) in correlation-consistent polarized basis sets of both double- and triple-ζ quality. Through an initial application to many-body expanded full configuration interaction (MBE-FCI) theory, we show how molecular point-group symmetry can lead to computational savings that are inversely proportional to the order of the point group in a manner generally applicable to the acceleration of modern local correlation methods.

Quantitative evaluation of noncovalent interactions with negative fragmentation approach including basis set superposition error correction

Abstract We propose a negative fragmentation approach (NFA), including counterpoise (CP) correction to basis set superposition error (BSSE) for quantitatively evaluating intra- and intermolecular noncovalent interactions. Noncovalent interactions are widely found in chemistry and biology and are regarded as essential interactions. However, there are few general methods for evaluating these individual intra- and intermolecular interaction energies because of two issues: (i) difficulty in the evaluation of intramolecular interactions due to the interacting sites connected with covalent bonds and (ii) BSSE affecting the quantitative accuracy of interaction analysis. In our scheme, we overcome the issue (i) using the NFA scheme, which can evaluate intra- and intermolecular interactions as a fragment–fragment interaction of interacting sites, and address the issue (ii) using the CP method. Here, NFA including the CP correction was applied to various molecular systems, providing comparable results for intermolecular interactions to supermolecule calculations with the CP correction and also succeeding in the evaluation of intramolecular interactions and its BSSEs. It is notable that our NFA scheme does not require any particular program or a modification of the program codes. These indicate that many researchers can apply our NFA scheme to various molecular systems.

Molecular Aspects of the Interactions between Selected Benzodiazepines and Common Adulterants/Diluents: Forensic Application of Theoretical Chemistry Methods.

Benzodiazepines are frequently encountered in crime scenes, often mixed with adulterants and diluents, complicating their analysis. This study investigates the interactions between two benzodiazepines, lorazepam (LOR) and alprazolam (ALP), with common adulterants/diluents (paracetamol, caffeine, glucose, and lactose) using infrared (IR) spectroscopy and quantum chemical methods. The crystallographic structures of LOR and ALP were optimized using several functionals (B3LYP, B3LYP-D3BJ, B3PW91, CAM-B3LYP, M05-2X, and M06-2X) combined with the 6-311++G(d,p) basis set. M05-2X was the most accurate when comparing experimental and theoretical bond lengths and angles. Vibrational and 13C NMR spectra were calculated to validate the functional's applicability. The differences between LOR's experimental and theoretical IR spectra were attributed to intramolecular interactions between LOR monomers, examined through density functional theory (DFT) optimization and quantum theory of atoms in molecules (QTAIM) analysis. Molecular dynamics simulations modeled benzodiazepine-adulterant/diluent systems, predicting the most stable structures, which were further analyzed using QTAIM. The strongest interactions and their effects on IR spectra were identified. Comparisons between experimental and theoretical spectra confirmed spectral changes due to interactions. This study demonstrates the potential of quantum chemical methods in analyzing complex mixtures, elucidating spectral changes, and assessing the structural stability of benzodiazepines in forensic samples.

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 \emph{a-priory} only the most dominant configurations. 
This is demonstrated through a set of groundstate calculations of the H$_2$, LiH, BeH$_2$, H$_2$O, NH$_3$ and C$_2$H$_4$ molecules in the sto-3g basis set, performed on IBM quantum devices. We show that the VQ-SCI reaches the full configuration interaction (FCI) 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.