Highly efficient implementation of pseudospectral time-dependent density-functional theory for the calculation of excitation energies of large molecules

Computational Multiscale Study of the Interaction Between the PDMS Polymer and Sunscreen-Related Pollutant Molecules

Sunscreen molecules play a critical role in protecting skin from ultraviolet radiation, yet their efficient detection and separation pose challenges in environmental and analytical contexts. In this work, we employ a multilevel modeling approach to investigate the molecular interactions between representative sunscreen molecules and the polydimethylsiloxane (PDMS) polymer, a material widely recognized for its sorbent properties. Our goal is to explore how these interactions can be fine-tuned to facilitate the effective separation of sunscreen molecules in portable membrane inlet mass spectrometry (MIMS) systems, potentially leading to the development of new membrane materials. Using a combination of advanced computational techniques-force field molecular dynamics simulations, semiempirical GFN2-xTB, and density functional theory calculations-we assess the interaction strength and noncovalent interactions of sunscreen molecules, namely oxybenzone, naphthalene, benzo[a]anthracene, avobenzone, and 1,3,5-trichlorobenzene, with PDMS. Additionally, the effect of temperature on the interaction dynamics is evaluated, with the aim of extending the sorbent capacities of PDMS beyond light polar molecules to larger, polar sunscreen compounds. This study provides critical insights into the molecular-level interactions that may guide the design of novel membrane materials for efficient molecular separation.

Quantum chemical package Jaguar: A survey of recent developments and unique features

This paper is dedicated to the quantum chemical package Jaguar, which is commercial software developed and distributed by Schrödinger, Inc. We discuss Jaguar's scientific features that are relevant to chemical research as well as describe those aspects of the program that are pertinent to the user interface, the organization of the computer code, and its maintenance and testing. Among the scientific topics that feature prominently in this paper are the quantum chemical methods grounded in the pseudospectral approach. A number of multistep workflows dependent on Jaguar are covered: prediction of protonation equilibria in aqueous solutions (particularly calculations of tautomeric stability and pKa), reactivity predictions based on automated transition state search, assembly of Boltzmann-averaged spectra such as vibrational and electronic circular dichroism, as well as nuclear magnetic resonance. Discussed also are quantum chemical calculations that are oriented toward materials science applications, in particular, prediction of properties of optoelectronic materials and organic semiconductors, and molecular catalyst design. The topic of treatment of conformations inevitably comes up in real world research projects and is considered as part of all the workflows mentioned above. In addition, we examine the role of machine learning methods in quantum chemical calculations performed by Jaguar, from auxiliary functions that return the approximate calculation runtime in a user interface, to prediction of actual molecular properties. The current work is second in a series of reviews of Jaguar, the first having been published more than ten years ago. Thus, this paper serves as a rare milestone on the path that is being traversed by Jaguar's development in more than thirty years of its existence.

Design of Novel Membranes for the Efficient Separation of Bee Alarm Pheromones in Portable Membrane Inlet Mass Spectrometric Systems

Bee alarm pheromones are essential molecules that are present in beehives when some threats occur in the bee population. In this work, we have applied multilevel modeling techniques to understand molecular interactions between representative bee alarm pheromones and polymers such as polymethyl siloxane (PDMS), polyethylene glycol (PEG), and their blend. This study aimed to check how these interactions can be manipulated to enable efficient separation of bee alarm pheromones in portable membrane inlet mass spectrometric (MIMS) systems using new membranes. The study involved the application of powerful computational atomistic methods based on a combination of modern semiempirical (GFN2-xTB), first principles (DFT), and force-field calculations. As a fundamental work material for the separation of molecules, we considered the PDMS polymer, a well-known sorbent material known to be applicable for light polar molecules. To improve its applicability as a sorbent material for heavier polar molecules, we considered two main factors-temperature and the addition of PEG polymer. Additional insights into molecular interactions were obtained by studying intrinsic reactive properties and noncovalent interactions between bee alarm pheromones and PDMS and PEG polymer chains.

Highly efficient implementation of analytic nonadiabatic derivative couplings within the pseudospectral method

A pseudospectral implementation of nonadiabatic derivative couplings in the Tamm-Dancoff approximation is reported, and the accuracy and efficiency of the pseudospectral nonadiabatic derivative couplings are studied. Our results demonstrate that the pseudospectral method provides mean absolute errors of 0.2%-1.9%, while providing a significant speedup. Benchmark calculations on fullerenes (Cn, n up to 100) using B3LYP achieved 10- to 15-fold, 8- to 17-fold, and 43- to 75-fold speedups for 6-31G**, 6-31++G**, and cc-pVTZ basis sets, respectively, when compared to the conventional spectral method.

Comprehensive Similarity Algorithm and Molecular Dynamics Simulation-Assisted Terahertz Spectroscopy for Intelligent Matching Identification of Quorum Signal Molecules (N-Acyl-Homoserine Lactones)

To investigate the mechanism of aquatic pathogens in quorum sensing (QS) and decode the signal transmission of aquatic Gram-negative pathogens, this paper proposes a novel method for the intelligent matching identification of eight quorum signaling molecules (N-acyl-homoserine lactones, AHLs) with similar molecular structures, using terahertz (THz) spectroscopy combined with molecular dynamics simulation and spectral similarity calculation. The THz fingerprint absorption spectral peaks of the eight AHLs were identified, attributed, and resolved using the density functional theory (DFT) for molecular dynamics simulation. To reduce the computational complexity of matching recognition, spectra with high peak matching values with the target were preliminarily selected, based on the peak position features of AHL samples. A comprehensive similarity calculation (CSC) method using a weighted improved Jaccard similarity algorithm (IJS) and discrete Fréchet distance algorithm (DFD) is proposed to calculate the similarity between the selected spectra and the targets, as well as to return the matching result with the highest accuracy. The results show that all AHL molecular types can be correctly identified, and the average quantization accuracy of CSC is 98.48%. This study provides a theoretical and data-supported foundation for the identification of AHLs, based on THz spectroscopy, and offers a new method for the high-throughput and automatic identification of AHLs.

Large Computational Survey of Intrinsic Reactivity of Aromatic Carbon Atoms with Respect to a Model Aldehyde Oxidase

Aldehyde oxidase (AOX) and other related molybdenum-containing enzymes are known to oxidize the C-H bonds of aromatic rings. This process contributes to the metabolism of pharmaceutical compounds and, therefore, is of vital importance to drug pharmacokinetics. The present work describes an automated computational workflow and its use for the prediction of intrinsic reactivity of small aromatic molecules toward a minimal model of the active site of AOX. The workflow is based on quantum chemical transition state searches for the underlying single-step oxidation reaction, where the automated protocol includes identification of unique aromatic C-H bonds, creation of three-dimensional reactant and product complex geometries via a templating approach, search for a transition state, and validation of reaction end points. Conformational search on the reactants, products, and the transition states is performed. The automated procedure has been validated on previously reported transition state barriers and was used to evaluate the intrinsic reactivity of nearly three hundred heterocycles commonly found in approved drug molecules. The intrinsic reactivity of more than 1000 individual aromatic carbon sites is reported. Stereochemical and conformational aspects of the oxidation reaction, which have not been discussed in previous studies, are shown to play important roles in accurate modeling of the oxidation reaction. Observations on structural trends that determine the reactivity are provided and rationalized.

Theoretical analysis of OLED performances of some aromatic nitrogen-containing ligands

It is well-known that tris(8-hydroxyquinoline) aluminum (Alq3) complex and N,N'diphenyl-N,N'-bis(3-methylphenyl)-1,1'-diphenyl-4,4'-diamine compound (TPD) are widely used as electron transfer material (ETL) and hole transfer material (HTL) in organic light emitting diode (OLED) structure, respectively. Considering the reference materials, in the present work, the OLED performances of some cyclic aromatic structures such as 4,4'azopyridine [AZPY], 4,4'-bipyridine [BIPY], 1,2-bis[4'-(4-methylphenyl)2,2':6'2″-terpyridin6-yl]ethyne (BISTERPY), 5,5'-diamino-2,2'-bipyridine (DABP), dipyrido[3,2-a:2',3'c]phenazine (DPP), 4,7-phenanthroline (PHEN) including nitrogen atom have been theoretically analyzed. It is important to note that B3LYP/6-31G(d) and B3LYP/TZP levels of the theory were taken into account for the calculations about monomeric and dimeric structures, respectively. Additionally, the calculations of the mentioned monomeric form were performed at B3LYP-D3/6-31G, CAM-B3LYP/6-31G and ωB97X-D/6-31G(d) levels. For a detailed theoretical analysis, the reorganization energies (λe and λh), adiabatic and vertical ionization potentials and electron affinities, the effective transfer integrals (Ve and Vh), and the charge transfer rates (We and Wh) of all compounds were computed by means of computational chemistry tools. In the light of calculated parameters, it is determined that these mentioned aromatic cyclic structures will be used in which layers of OLED structure. The results obtained in this study will be helpful in the design and applications of new molecules as OLED materials in the future.

Photocatalytic Activity of the V2O5 Catalyst toward Selected Pharmaceuticals and Their Mixture: Influence of the Molecular Structure on the Efficiency of the Process

Due to the inability of conventional wastewater treatment procedures to remove organic pharmaceutical pollutants, active pharmaceutical components remain in wastewater and even reach tap water. In terms of pharmaceutical pollutants, the scientific community focuses on β-blockers due to their extensive (over)usage and moderately high solubility. In this study, the photocatalytic activity of V2O5 was investigated through the degradation of nadolol (NAD), pindolol (PIN), metoprolol (MET), and their mixture under ultraviolet (UV) irradiation in water. For the preparation of V2O5, facile hydrothermal synthesis was used. The structural, morphological, and surface properties and purity of synthesized V2O5 powder were investigated by scanning electron microscopy (SEM), X-ray, and Raman spectroscopy. SEM micrographs showed hexagonal-shaped platelets with well-defined morphology of materials with diameters in the range of 10−65 µm and thickness of around a few microns. X-ray diffraction identified only one crystalline phase in the sample. The Raman scattering measurements taken on the catalyst confirmed the result of XRPD. Degradation kinetics were monitored by ultra-fast liquid chromatography with diode array detection. The results showed that in individual solutions, photocatalytic degradation of MET and NAD was relatively insignificant (<10%). However, in the PIN case, the degradation was significant (64%). In the mixture, the photodegradation efficiency of MET and NAD slightly increased (15% and 13%). Conversely, it reduced the PIN to the still satisfactory value of 40%. Computational analysis based on molecular and periodic density functional theory calculations was used to complement our experimental findings. Calculations of the average local ionization energy indicate that the PIN is the most reactive of all three considered molecules in terms of removing an electron from it.

Pseudospectral implementations of long-range corrected density functional theory

We have implemented pseudospectral density-functional theory (DFT) with long-range corrected DFT functionals (PS-LRC) in quantum mechanics package Jaguar, and applied it in the calculations of geometry optimizations, dimmer interaction energies, polarizabilities and first-order hyperpolarizabilities, harmonic vibrational frequencies, S₁ and T₁ excitation energies, singlet-triplet gaps, charge transfer numbers, oscillator strengths, reaction barrier heights, electron-transfer couplings, and charge-transfer excitation energies. From our accuracy benchmark analysis, PS grids, PS dealiasing functions, PS atomic corrections, PS multigrid strategy, PS length scales, and PS cutoff scheme perform well in PS DFT with LRC density functionals with very small and ignorable deviations when compared to the conventional spectral (CS) method. The timing benchmark study of S₁ excitation energy calculations of fullerenes (C_n , n up to 540) demonstrates that PS-LRC achieves 1.4-8.4-fold speedups in SCF, 22-32-fold speedups in Tamm-Dancoff approximation, and 6-15-fold speedups in total wall clock time with an average error 0.004 eV of excitation energies compared to the CS method.

Highly efficient implementation of the analytical gradients of pseudospectral time-dependent density functional theory

The accuracy and efficiency of time-dependent density functional theory (TDDFT) excited state gradient calculations using the pseudospectral method are presented. TDDFT excited state geometry optimizations of the G2 test set molecules, the organic fluorophores with large Stokes shifts, and the Pt-complexes show that the pseudospectral method gives average errors of 0.01-0.1 kcal/mol for the TDDFT excited state energy, 0.02-0.06 pm for the bond length, and 0.02-0.12° for the bond angle when compared to the results from conventional TDDFT. TDDFT gradient calculations of fullerenes (C_n, n up to 540) with the B3LYP functional and 6-31G** basis set show that the pseudospectral method provides 8- to 14-fold speedups in the total wall clock time over the conventional methods. The pseudospectral TDDFT gradient calculations with a diffuse basis set give higher speedups than the calculations for the same basis set without diffuse functions included.

Exploring the detailed spectroscopic characteristics, chemical and biological activity of two cyanopyrazine-2-carboxamide derivatives using experimental and theoretical tools

This article represents the spectroscopic and computational studies of two new pyrazine compounds. In order to establish the structure and functional nature of the compounds, we have employed Fourier transformed infrared (FT-IR) and Raman spectra, nuclear magnetic resonance (NMR) spectra, and ultraviolet (UV) absorptions and have compared them with the simulated computational spectra and found that they are in the agreeable range. Simulated hyperpolarisability values are used to obtain the nonlinear optic (NLO) activity of the compound, to be used in organic electronic materials. The charge transfer and related properties was investigated by the simulation of electronic spectrum with time dependent density functional theory (TD-DFT). Natural transition orbitals (NTO) provides information about which region of the molecules are more involved in the electronic transitions and the charge transfer properties for the lowest energy excitation have been analyzed on the basis of electron density variation. Molecular dynamics simulations provide information about the behavior of the molecule in solutions. Frontier orbital analysis and study of various reactivity descriptors like ALIE and Fukui provided deep knowledge on the reactivity side. Molecular docking has been also performed to investigate the interaction between title molecules and exhibits inhibitory activity against Pseudomonas aeruginosa Enoyl-Acyl carrier protein reductase (Fabl).

Modeling MEK4 Kinase Inhibitors through Perturbed Electrostatic Potential Charges

MEK4, mitogen-activated protein kinase kinase 4, is overexpressed and induces metastasis in advanced prostate cancer lesions. However, the value of MEK4 as an oncology target has not been pharmacologically validated because selective chemical probes targeting MEK4 have not been developed. With advances in both computer and biological high-throughput screening, selective chemical entities can be discovered. Structure-based quantitative structure-activity relationship (QSAR) modeling often fails to generate accurate models due to poor alignment of training sets containing highly diverse compounds. Here we describe a highly predictive, nonalignment based robust QSAR model based on a data set of strikingly diverse MEK4 inhibitors. We computed the electrostatic potential (ESP) charges using a density functional theory (DFT) formalism of the donor and acceptor atoms of the ligands and hinge residues. Novel descriptors were then generated from the perturbation of the charge densities of the donor and acceptor atoms and were used to model a diverse set of 84 compounds, from which we built a robust predictive model.

Density Functional Theory Analysis of Deltamethrin and Its Determination in Strawberry by Surface Enhanced Raman Spectroscopy

Deltamethrin is widely used in pest prevention and control such as red spiders, aphids, and grubs in strawberry. It is important to accurately monitor whether the deltamethrin residue in strawberry exceeds the standard. In this paper, density functional theory (DFT) was used to theoretically analyze the molecular structure of deltamethrin, gold nanoparticles (AuNPs) and silver nanoparticles (AgNPs) were used to enhance the surface enhanced Raman spectroscopy (SERS) detection signal. As a result, the theoretical Raman peaks of deltamethrin calculated by DFT were basically similar to the measured results, and the enhancing effects based on AuNPs was better than that of AgNPs. Moreover, 554, 736, 776, 964, 1000, 1166, 1206, 1593, 1613, and 1735 cm⁻¹ could be determined as deltamethrin characteristic peaks, among which only three Raman peaks (736, 1000, and 1166 cm⁻¹) could be used as the deltamethrin characteristic peaks in strawberry when the detection limit reached 0.1 mg/L. In addition, the 500⁻1800 cm⁻¹ SERS of deltamethrin were analyzed by the partial least squares (PLS) and backward interval partial least squares (BIPLS). The prediction accuracy of deltamethrin in strawberry (Rp2 = 0.93, RMSE_p = 4.66 mg/L, RPD = 3.59) was the highest when the original spectra were pretreated by multiplicative scatter correction (MSC) and then modeled by BIPLS. In conclusion, the deltamethrin in strawberry could be qualitatively analyzed and quantitatively determined by SERS based on AuNPs enhancement, which provides a new detection scheme for deltamethrin residue determination in strawberry.

Rapid Determination of Thiabendazole Pesticides in Rape by Surface Enhanced Raman Spectroscopy

Thiabendazole is widely used in sclerotium blight, downy mildew and black rot prevention and treatment in rape. Accurate monitoring of thiabendazole pesticides in plants will prevent potential adverse effects to the Environment and human health. Surface Enhanced Raman Spectroscopy (SERS) is a highly sensitive fingerprint with the advantages of simple operation, convenient portability and high detection efficiency. In this paper, a rapid determination method of thiabendazole pesticides in rape was conducted combining SERS with chemometric methods. The original SERS were pretreated and the partial least squares (PLS) was applied to establish the prediction model between SERS and thiabendazole pesticides in rape. As a result, the SERS enhancing effect based on silver Nano-substrate was better than that of gold Nano-substrate, where the detection limit of thiabendazole pesticides in rape could reach 0.1 mg/L. Moreover, 782, 1007 and 1576 cm⁻¹ could be determined as thiabendazole pesticides Raman characteristic peaks in rape. The prediction effect of thiabendazole pesticides in rape was the best (Rp2 = 0.94, RMSEP = 3.17 mg/L) after the original spectra preprocessed with 1st-Derivative, and the linear relevance between thiabendazole pesticides concentration and Raman peak intensity at 782 cm⁻¹ was the highest (R² = 0.91). Furthermore, five rape samples with unknown thiabendazole pesticides concentration were used to verify the accuracy and reliability of this method. It was showed that prediction relative standard deviation was 0.70–9.85%, recovery rate was 94.71–118.92% and t value was −1.489. In conclusion, the thiabendazole pesticides in rape could be rapidly and accurately detected by SERS, which was beneficial to provide a rapid, accurate and reliable scheme for the detection of pesticides residues in agriculture products.